PubMed Health⌕ Search

SEARCH · PubMed Health

Results for “Monte Carlo Method”

Explore indexed PubMed citations for clinical trials, systematic reviews and public health research. Read source abstracts and follow each citation to its original PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 829 records · Page 46Linked to original sources

Air-ground interface correction factors for gamma emitters in air.

For isotropic point photon sources above an air-ground interface, the spectral energy fluences at the ground have been calculated by the Monte Carlo method for source energies in the range 50 keV-10 MeV. The source heights and the horizontal source-detector distances varied from 0.05 to 5 mean free paths in air. Air-ground interface correction factors have been calculated and then analytically approximated. The greatest deviations from unity were found for 100-keV sources located on the ground. Near the source, kerma free in air is enhanced by 30%, and far from the source it is reduced by a factor of 3.

Air Pollution, Radioactive↗

Determination of HPGe detector response using MCNP5 for 20-150 keV X-rays.

The Monte Carlo method was used to calculate the efficiency, escape and Compton curves of a planar high-purity germanium detector (HPGe) in the 20-150 keV energy. These curves were used for the determination of photons spectra produced by an X-ray machine in order to allow a precise characterization of photon beams applied to medical diagnosis. The detector was modeled with the MCNP5 computer code and validated by comparison with experimental data. The air kerma calculated after the spectra stripping was compared with ionization chamber measurements.

Germanium↗

The accuracy of single-seed dose superposition for I-125 implants.

The Monte Carlo method was used to study perturbations of single I-125 seed dose distributions created by the presence of one or three neighboring seeds for the case of seeds immersed in a water phantom. Perturbation factors were determined within the geometric shadow of neighboring seeds for two-seed designs, four-seed spacings, and several choices of dose point. The results were compared to dose estimates obtained by the simple superposition of single-seed data for one- and two-plane implants. Some significant differences were found.

Brachytherapy↗

Correction factors for water-proofing sleeves in kilovoltage x-ray beams.

This paper investigates the effect of the waterproofing sleeve on the calibration of kilovoltage photon beams (50-300 kV). The sleeve effect correction factor, ps has been calculated using the Monte Carlo method as the ratios of the air kerma in an air cavity of a cylindrical chamber without the waterproofing sleeve to that with a sleeve. Three sleeve materials have been studied, PMMA, nylon and polystyrene. The calculations were carried out using the EGS4 (Electron Gamma Shower version 4) code system with the application of a correlated-sampling variance-reduction technique. The results show that the sleeve correction factor for 1-mm thick nylon and polystyrene sleeves, ps varies from 0.992 to 1.000 and from 0.981 to 1.000, respectively, for the same beam quality range. The ps factor varies with sleeve thickness, beam quality and phantom depth. No significant dependence of the ps factor on field size and source-surface distance has been found. Measurements for PMMA, nylon and polystyrene sleeves of various thicknesses have also been carried out and show excellent agreement with Monte Carlo calculations.

Biophysical Phenomena↗

Detecting pseudoknots and other local base-pairing structures in RNA sequences.

The current version of RNAFOLD now has the capability for conducting comparative studies relative to nonbifurcating hairpins contained in the global structures produced by the dynamic programming and Monte Carlo methods. It also has the capability of both finding and comparing hairpins and pseudoknots independently of global structures. The efficacy of these increased capabilities has been tested for select families of sequences, and the results thus far indicate favorable utility. Under consideration is a further extension designed to incorporate pseudoknots within global structures. Written in the C language, RNAFOLD and its companion program, GENALIGN, for doing multiple alignments in the comparison of pseudoknots and hairpins, is available for UNIX systems on standard 1/2-inch, 9-track tape or on SUN tape cartridges.

Base Sequence↗

Transfer of 137Cs from soil to grass--analysis of possible sources of uncertainty.

On the basis of experimental data on the fallout of 137Cs from the Chernobyl accident, a statistical analysis was made of possible values of the coefficients of transfer via the soil-grass pathway. Model calculations of the possible radionuclide concentration in milk are performed by the Monte Carlo method using the probability distribution function obtained. It is shown that not only our limited knowledge and the possible diversity of the external conditions, but also the very nature of the phenomenon under consideration, make it desirable to switch from a deterministic to a stochastic formulation of the model. The fact that the dependence of the soil-to-grass radionuclide transfer coefficients on the initial concentration in soil might be non-linear in character is an indication of the limitation of the conventional linear model.

Accidents↗

Thermodynamical properties of simple models of protein-like heteropolymers.

The lattice approximation of a heteropolymer chain as a model of a single polypeptide was used in the computer simulation. The residues of a model polypeptide were represented by the chain of alpha-carbons located on a very flexible [310] lattice. The force field that mimic the intramolecular interactions contained the long-range contact potential between the residues and the local preferences in forming helical structures. The chain consisted of two types of residues that had different hydrophobicity. The simulations were performed by means of the Replica Exchange Monte Carlo method combined with the Histogram method. The series of simulations were carried out to investigate the influence of both components of the force field on the transition temperature and the characteristics of the coil-to-globule transition. The properties of low-temperature ordered structures were determined. The thermodynamical description of the model chain was also given. The phase transition was found to be sharp and cooperative for longer chains and strong helical potential. The collapsed globule contained the strongly hydrophobic residues inside the globule while the remaining residues were mainly located close to the globule surface.

Asparagine↗

Temperature distribution in selective laser-tissue interaction.

Selective photothermal interaction using dye enhancement has proven to be effective in minimizing surrounding tissue damage and delivering energy to target tissue. During laser irradiation, the process of photon absorption and thermal energy diffusion in the target tissue and its surrounding tissue are crucial. Such information allows the selection of proper operating parameters such as dye concentrations, laser power, and exposure time for optimal therapeutic effect. Combining the Monte Carlo method for energy absorption and the finite difference method for heat diffusion, the temperature distributions in target tissue and surrounding tissue in dye enhanced laser photothermal interaction are obtained. Different tissue configurations and dye enhancement are used in the simulation, and different incident beam sizes are also used to determine optimum beam sizes for various tissue configurations. Our results show that the algorithm developed in this study could predict the thermal outcome of laser irradiation. Our simulation indicates that with appropriate absorption enhancement of the target tissue, the temperature in the target tissue and in the surrounding tissue can be effectively controlled. This method can be used for optimization of lesion treatment using laser photothermal interactions. It may also provide guidance for laser immunotherapy in cancer treatment, since the immunological responses are believed to be related to tissue temperature changes.

Animals↗

Quasi Monte Carlo-based isotropic distribution of gradient directions for improved reconstruction quality of 3D EPR imaging.

In continuous wave (CW) electron paramagnetic resonance imaging (EPRI), high quality of reconstructed image along with fast and reliable data acquisition is highly desirable for many biological applications. An accurate representation of uniform distribution of projection data is necessary to ensure high reconstruction quality. The current techniques for data acquisition suffer from nonuniformities or local anisotropies in the distribution of projection data and present a poor approximation of a true uniform and isotropic distribution. In this work, we have implemented a technique based on Quasi-Monte Carlo method to acquire projections with more uniform and isotropic distribution of data over a 3D acquisition space. The proposed technique exhibits improvements in the reconstruction quality in terms of both mean-square-error and visual judgment. The effectiveness of the suggested technique is demonstrated using computer simulations and 3D EPRI experiments. The technique is robust and exhibits consistent performance for different object configurations and orientations.

Algorithms↗

Room scatter studies in the air kerma strength standardization of the Amersham CDCS-J-type 137Cs source: a Monte Carlo study.

Corrections for room scatter, [kSC (d, H)]RM, as a function of source-to-detector distance, d, and source-to-floor-height, H, and also departure from constant room scatter. ks1, have been computed for rooms of various sizes using Monte Carlo methods for air kerma strength standardization of the Amersham CDCS-J-type 137Cs brachytherapy source. These corrections will also be applicable to any type of 137Cs source that may be considered for standardization. It was found that, depending upon the relative position of the source with respect to the surrounding concrete scattering surfaces (side walls, floor and ceiling) and different set of d values, the assumption of constant room scatter overestimated the air kerma strength, Sk, by between 0.2% and 0.6%.

Air↗

Comparison of a finite-element multigroup discrete-ordinates code with Monte Carlo for radiotherapy calculations.

Radiotherapy calculations often involve complex geometries such as interfaces between materials of vastly differing atomic number, such as lung, bone and/or air interfaces. Monte Carlo methods have been used to calculate accurately the perturbation effects of the interfaces. However, these methods can be computationally expensive for routine clinical calculations. An alternative approach is to solve the Boltzmann equation deterministically. We present one such deterministic code, Attila. Further, we computed a brachytherapy example and an external beam benchmark to compare the results with data previously calculated by MCNPX and EGS4. Our data suggest that the presented deterministic code is as accurate as EGS4 and MCNPX for the transport geometries examined in this study.

Body Burden↗

Stopping power data for high-energy photon beams.

Spencer-Attix stopping power ratios for the dosimetry of high-energy photon beams used in radiation therapy have been calculated using the Monte Carlo method. The stopping power ratios are calculated in a more consistent way than previously and are given as a function of the attenuation properties of the beam. The dependence of the stopping power ratio on the electron contamination of the beam as well as on depth and field size has also been investigated. Results are compared with stopping power ratios recommended in different dosimetry protocols and to experimental results. The agreement with most dosimetry protocols is within about one per cent and with recent experimental data is better than half a per cent.

Evaluation Studies as Topic↗

Monte Carlo tests for associations between disease and alleles at highly polymorphic loci.

In an association analysis comparing cases and controls with respect to allele frequencies at a highly polymorphic locus, a potential problem is that the conventional chi-squared test may not be valid for a large, sparse contingency table. However, reliance on statistics with known asymptotic distribution is now unnecessary, as Monte Carlo simulations can be performed to estimate the significance level of any test statistic. We have implemented a Monte Carlo method for four 'chi-squared' test statistics, three of which involved combination of alleles, and evaluated their performance on a real data set. Combining rare alleles to avoid small expected cell counts, and considering each allele in turn against the rest, reduced the power to detect a genuine association when the number of alleles was very large. We should either not combine alleles at all, or combine them in such a way that preserves the evidence for an association.

Alleles↗

Functional form and risk adjustment of hospital costs: Bayesian analysis of a Box-Cox random coefficients model.

While risk-adjusted outcomes are often used to compare the performance of hospitals and physicians, the most appropriate functional form for the risk adjustment process is not always obvious for continuous outcomes such as costs. Semi-log models are used most often to correct skewness in cost data, but there has been limited research to determine whether the log transformation is sufficient or whether another transformation is more appropriate. This study explores the most appropriate functional form for risk-adjusting the cost of coronary artery bypass graft (CABG) surgery. Data included patients undergoing CABG surgery at four hospitals in the midwest and were fit to a Box-Cox model with random coefficients (BCRC) using Markov chain Monte Carlo methods. Marginal likelihoods and Bayes factors were computed to perform model comparison of alternative model specifications. Rankings of hospital performance were created from the simulation output and the rankings produced by Bayesian estimates were compared to rankings produced by standard models fit using classical methods. Results suggest that, for these data, the most appropriate functional form is not logarithmic, but corresponds to a Box-Cox transformation of -1. Furthermore, Bayes factors overwhelmingly rejected the natural log transformation. However, the hospital ranking induced by the BCRC model was not different from the ranking produced by maximum likelihood estimates of either the linear or semi-log model.

Bayes Theorem↗

Comparison of measured and calculated dose rates in water near I-125 and Ir-192 seeds.

Recent theoretical and experimental work indicates that currently accepted 125I dosimetry data may overestimate dose in water at 1 cm by 10%-24%. Among the most comprehensive measurements are those of the NCI-sponsored brachytherapy contract participants. Absolute dose rates in water calculated by the Monte Carlo method have been compared with the NCI dose measurements about 125I and 192Ir seeds embedded in solid-water phantoms. The photon transport code allows realistic geometric simulation of the complex internal seed structure, the National Institute of Standards and Technology air-kerma strength standardization geometry, and the dose measurement setup. When the appropriate measurement medium and geometry are assumed, agreement between theory and measurement is excellent, within 3% at 1 cm and averaging 3% at larger distances. However, the data do not support the water equivalence of solid water at 125I energies indicating that solid-water measurements underestimate 125I specific dose-rate constants in water by 4.3%. Because of its higher ratio of absorption to scatter, 125I dose distributions measured in solid water are less penetrating (by 35% at 10 cm) than those measured in liquid water. For model 6711, model 6702, and steel-clad 192Ir seeds, Monte Carlo calculations yielded specific dose-rate constants (assuming liquid water medium) of 0.877, 0.932, and 1.122 cGy cm2 h-1 per unit air-kerma strength, respectively. For 125I, currently accepted values are 18% and 11% larger for the two seed models.

Humans↗

A spectral unfolding method to determine source depth distribution.

The depth distribution in tissue of a gamma-emitting isotope can be estimated by observing the energy spectrum of photons emerging from the body. Monte Carlo methods were used to compute surface energy spectra for 99Tcm point sources at various depths. From these noise-free results a discrete response matrix was constructed to relate the input source depth to the corresponding output spectrum. Observed spectra with noise were then simulated independently by Monte Carlo calculations assuming various measurement times and source strengths. The source distribution was determined by solving the discrete inverse unfolding problem by the Method of Regularisation. The effective depth resolution is seen to depend strongly upon measurement time. The results show that limited depth resolution can be obtained from a single view of the object region.

Monte Carlo Method↗

EUD-based radiotherapy treatment plan evaluation: incorporating physical and Monte Carlo statistical dose uncertainties.

The purpose of this work is to quantify the impact of dose uncertainty on radiobiologically based treatment plan evaluation. Dose uncertainties are divided into two categories: physical and statistical. Physical dose uncertainty is associated with the systematic and/or random errors incurred during treatment planning and/or delivery. The dose uncertainty associated with Monte Carlo calculated dose distributions is deemed statistical and noted as artificial with respect to the actual delivered dose. We will refer to all dose uncertainties that arise from either calculation or delivery as stochastic. Both physical and statistical dose uncertainties are considered at the intra- and inter-voxel levels. To account for voxel dose uncertainty, we calculate the mean survival fraction (SF) for the random dose deposition. Mathematically, the expression for the mean survival fraction is identical to that used by Niemierko (1997 Med. Phys. 24 103-10) in defining equivalent uniform dose (EUD). To distinguish between spatial and probabilistic dose variations, we define equivalent stochastic dose (ESD) as a voxel dose that gives the mean expected survival fraction for the randomly deposited dose. For a probability density function f(D), that represents the probabilistic voxel dose, SF(ESD) can be calculated by convolving SF(D) with f(D). In the case where the probability density function follows a Gaussian distribution, an analytic expression is derived for SF(ESD). The derived expression is verified using the Monte Carlo method and ESD values calculated with varied radiosensitivities for cases of 60 and 70 Gy at 2 Gy per fraction. The analytic expression is also extended to account for a multi-voxel dose distribution that incorporates a spatial dose heterogeneity. The results show that survival fraction increases with an increased dose uncertainty. This reduction depends on radiobiological parameters attributed to tissue and tumour. For tissue, ESD drops to 55% of the mean physical dose when the dose has a 10% intra- and inter-voxel dose uncertainty and inhomogeneity.

Body Burden↗

Polymorphisms in the +252(A/G) lymphotoxin-alpha and the -308(A/G) tumor necrosis factor-alpha genes and susceptibility to chronic periodontitis in a Czech population.

BACKGROUND: Periodontitis is an inflammatory disease that leads to irreversible attachment loss, bone destruction and eventually tooth loss. Tumor necrosis factor (TNF), a pluripotent proinflammatory cytokine that is able to induce tissue destruction and bone resorption, has been implicated in the pathogenesis of periodontal disease. METHODS: In this study, we investigated an association between chronic periodontitis and two previously described bi-allelic polymorphisms in the TNF locus: a G to A transition at position -308 in the 5'promoter region of the TNF-alpha gene and an NcoI restriction fragment length polymorphism (RFLP) in the first intron (position +252A/G) of the lymphotoxin alpha (LT-alpha) gene. Genomic DNA was obtained from 132 patients with chronic periodontitis together with 114 age- and gender-matched unrelated control subjects. RESULTS: The TNF-alpha (-308G/A) polymorphism itself showed no association with chronic periodontitis, whereas the frequency distribution of the LT-alpha (+252A/G) genotypes showed statistically significant differences between patients and the reference group. The proportion of individuals carrying the LT-alpha 1/1 genotype was significantly lower in the group of patients with chronic periodontitis (0.8%) than in the control group (8.8%) (P < 0.0094, Pcorr < 0.05). However, the significant differences in the frequencies of the combined genotypes (TNF-alpha and LT-alpha) between the control and the patient groups were found using a simulation by applying the Monte-Carlo method (P < 0.01). CONCLUSION: Our data suggest that combined genotypes composed of the TNF-alpha and LT-alpha gene polymorphisms may influence the susceptibility to chronic periodontitis. We also showed that, comparing the two genes, the 1/1 genotype of the NcoI polymorphism in the first intron of the LT-alpha gene is a more informative marker and it may be one of the protective genetic factors against chronic periodontitis in our population.

Adenine↗