SEARCH · PubMed Health

Results for “EM algorithm”

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 397 records · Page 22Linked to original sources

PubMed2003

Genetic variability of the marine mussel Mytilus galloprovincialis assessed using two-dimensional electrophoresis.

Two-dimensional electrophoresis (2-DE) has been used to measure the degree of genetic variability of the marine mussel Mytilus galloprovincialis. Genetic polymorphisms were detected in 33 of a total of 86 polypeptides scored among the most abundant proteins from foot samples in 38 individuals. Estimates of average heterozygosity were 0.101+/-0.018 and 0.114+/-0.021 in a natural and a cultured population, respectively, from the NW of the Iberian Peninsula. These are the highest estimates of average heterozygosity reported by 2-DE in an animal species to date. We consider that these data throw open the question of the level of genetic variability detectable by two-dimensional electrophoresis. Multilocus genotype data were used to infer haplotypic frequencies by means of the EM algorithm in order to detect linkage disequilibrium between loci coding abundant proteins. Significant associations were found in 22.7% of the 406 two-locus pairs analysed. Also, clusters of loci in which all pairwise combinations exhibit statistically significant associations were detected and physical linkage between some of these loci is postulated from the linkage disequilibrium data.

Animals↗

PubMed2005

Sequencing drug response with HapMap.

The information about how DNA sequence varies across the human genome is crucial for unravelling the genetic basis of drug response. A haplotype map, or HapMap, intended to reveal such a variation pattern, has been recently developed by the International HapMap Consortium. Here, we present a conceptual model for directly characterizing specific DNA sequence variants that are responsible for drug response based on the haplotype structure provided by HapMap. Our model is developed in the maximum likelihood context, incorporated by clinically meaningful mathematical functions that model drug response and implemented with the EM algorithm. Our model is employed to a pharmacogenetic study of cardiovascular disease with 107 patients. We found that the haplotype constituted by allele Gly16 (G) at codon 16 and allele Glu27 (G) at codon 27 genotyped within the beta2AR candidate gene exhibits a different effect on heart rate curve from the rest haplotypes. Parents with the diplotype consisting of two copies of haplotype GG are more sensitive in heart rate to increasing dosages of dobutamine than those with other haplotypes. This model provides a powerful tool for elucidating the genetic variants of drug response and ultimately designing personalized medications based on each patient's genetic constitution.

Algorithms↗

PubMed1999

The additive genetic gamma frailty model for linkage analysis of age-of-onset variation.

Age of onset is a key factor in the linkage analysis of many complex diseases. Current methods in nonparametric linkage analysis are mainly concentrated on the affected relative pairs or affected family members with age of onset information either ignored or taken into account by specifying age-dependent penetrances for liability classes. On the other hand, gamma frailty models were developed in the biostatistics literature to model familial aggregation of age of onset. However, these frailty models cannot be used directly for linkage analysis. This paper extends the gamma frailty model by incorporating inheritance vector information and provides a semiparametric approach for linkage testing. For a given inheritance vector at the putative disease locus, we construct an additive genetic gamma frailty for each individual within a nuclear family and use the Cox proportional hazard model to model age of onset. We derive the conditional hazard ratio parameter for sib pairs and define a likelihood ratio based LOD score statistic under our model. The EM algorithm is used for estimating the parameters and the maximum likelihood functions. Simulated data sets are used to illustrate these new statistical methods.

Age Factors↗

PubMed2003

Accelerated gene counting for haplotype frequency estimation.

Current implementations of the EM algorithm for estimating haplotype frequencies from genotypes on proximal loci require computational resources that grow as nh2k, where n is the number of individuals genotyped and h is the number of haplotypes possible on k loci. For diallelic loci hk=2k. We present an approach whose computational requirement grows as n2t where t is the largest number of loci at which an individual in the sample is heterozygous. The method is illustrated by haplotype frequency estimation from a sample of 45 individuals genotyped at 26 single nucleotide polymorphisms in the PIK3R1 gene.

1-Phosphatidylinositol 4-Kinase↗

PubMed2003

Association of single nucleotide polymorphisms of the bile salt export pump gene with intrahepatic cholestasis of pregnancy.

BACKGROUND: We determined whether genetic variability in the gene encoding the bile salt export pump (BSEP) contributes to individual differences in susceptibility to the development of intrahepatic cholestasis of pregnancy (ICP). METHODS: The study involved 57 affected and 115 healthy control pregnant women who were genotyped for two single nucleotide polymorphisms (SNPs) in the BSEP gene. Chi-square analysis was used to assess genotype and allele frequency differences between the cholestatic and control groups. In addition, single locus analysis was expanded to pair of loci haplotype analysis to examine the estimated haplotype frequencies of the two SNPs, of unknown phase, among the cholestatic and control groups. Estimated haplotype frequencies were assessed using the maximum-likelihood method, employing an expectation-maximization (EM) algorithm. RESULTS: The genotype and allele frequency distribution of the two intragenic SNPs in the ICP and control groups revealed significant evidence of association with the exon 28 SNP (P=0.04 and P=0.02, respectively). In addition, a borderline allele association was noted with the intron 19 SNP (P=0.08). Although the overall distribution of estimated haplotypes of intron 19 and exon 28 SNPs did not differ between the ICP and control groups, the most common haplotype, A-G, was significantly overrepresented in the ICP group (P=0.02), at an odds ratio of 1.73 (95% CI: 1.08-2.74). CONCLUSIONS: The use of two intragenic SNPs in both single locus and haplotype analyses of association suggests that the BSEP gene is a susceptibility gene in intrahepatic cholestasis of pregnancy.

ATP Binding Cassette Transporter, Subfamily B, Mem↗

PubMed1995

Estimation of the shelf-life of drugs with mixed effects models.

This paper proposes a normal mixed effects model for stability analysis. An EM algorithm is developed to compute the maximum likelihood estimates of regression coefficients of the fixed effects and random effects, and variance components. The likelihood ratio test is used for the preliminary testing of batch-to-batch variation. An example from a marketing stability study is given to illustrate the proposed procedure.

Algorithms↗

PubMed2003

A likelihood-based, counterfactual approach to accounting for treatment failures in clinical trials.

Consider a two-armed, placebo-controlled trial in which subjects may experience treatment failure. For ethical reasons, it is necessary to administer emergency or rescue medications for such subjects. However, the emergency medications may bias the set of response measurements. When analyzing the data from clinical trials, the standard approach is to perform an intent-to-treat (ITT) analysis, wherein the data are analyzed according to treatment assignment. Secondary statistical analyses that supplement the ITT analysis can be performed to account for the impact of treatment failures and emergency medications. A likelihood-based, counterfactual approach to supplemental analyses uses the expectation maximization (EM) algorithm for parameter estimation and a likelihood ratio test to test the equality of the placebo and experimental treatment means for subjects who would not fail under either treatment assignment. A simulation study is performed to assess the operating characteristics of the likelihood ratio test. An example from the Asthma Clinical Research Network (ACRN) is used to draw comparisons between the standard ITT procedure and the developed supplemental analysis.

Clinical Trials as Topic↗

PubMed2000

Segregation analysis in Shwachman-Diamond syndrome: evidence for recessive inheritance.

Shwachman-Diamond syndrome is a rare disorder of unknown cause. Reports have indicated the occurrence of affected siblings, but formal segregation analysis has not been performed. In families collected for genetic studies, the mean paternal age and mean difference in parental ages were found to be consistent with the general population. We determined estimates of segregation proportion in a cohort of 84 patients with complete sibship data under the assumption of complete ascertainment, using the Li and Mantel estimator, and of single ascertainment with the Davie modification. A third estimate was also computed with the expectation-maximization (EM) algorithm. All three estimates supported an autosomal recessive mode of inheritance, but complete ascertainment was found to be unlikely. Although there are no overt signs of disease in adult carriers (parents), the use of serum trypsinogen levels to indicate exocrine pancreatic dysfunction was evaluated as a potential measure for heterozygote expression. No consistent differences were found in levels between parents and a normal control population. Although genetic heterogeneity cannot be excluded, our results indicate that simulation and genetic analyses of Shwachman-Diamond syndrome should consider a recessive model of inheritance.

Abnormalities, Multiple↗

PubMed1991

Parallel simulated annealing for emission tomography.

A method for implementing simulated annealing in parallel to speed up the execution of emission tomography (ET) image reconstruction is presented. A high degree of parallelism can be attained by using a parallel-acceptance partitioning strategy, in which perturbations to subsets of the estimate are evaluated in parallel. However because the point spread function in ET imaging systems is globally dependent, processors cannot update the current estimate independently. Consequently, processors must be synchronized each time a perturbation is accepted to avoid introducing error. This can produce excessive communications overhead, especially when the acceptance rate is high. In this paper an energy function is constructed to reduce the synchronization requirements by using a reformulation of the log-likelihood function from the expectation maximization (EM) algorithm. The approach is to change the global dependence in the energy function from the current estimate to the estimate generated during the last iteration. The synchronization requirements for guaranteed convergence are then significantly reduced from once per acceptance to once per iteration. This parallel implementation on 54 Inmos T800 transputers connected in a ring topology resulted in execution times that were almost 50 times faster than on a VAX 8600.

Algorithms↗

PubMed1992

Truncation artifact suppression in cone-beam radionuclide transmission CT using maximum likelihood techniques: evaluation with human subjects.

Transverse image truncation can be a serious problem for human imaging using cone-beam transmission CT (CB-CT) implemented on a conventional rotating gamma camera. If this problem can be solved, CB-CT will be useful for attenuation compensation of SPECT images. This paper presents a reconstruction method to reduce or eliminate the artifacts resulting from the truncation. The method uses a previously published transmission maximum likelihood EM algorithm, adapted to the cone-beam geometry. The reconstruction method is evaluated qualitatively using three human subjects of various dimensions and various degrees of truncation. For the two smaller subjects, with moderate truncation, the maximum likelihood method is very successful, nearly eliminating the artifacts seen with conventional filtered backprojection of truncated geometries. The use of an expanded reconstructed space, which contains the entire transverse slice of the subject, is necessary for optimal truncation removal. For the largest subject investigated, the truncation was substantial, and the artifacts were only partially removed by the maximum likelihood reconstruction. Nonetheless, the images were qualitatively superior to those obtained with filtered backprojection. An added elliptical support prior moderately increased the rate of convergence, and helped to force a reasonable body contour.

Adult↗

PubMed1993

The convergence of object dependent resolution in maximum likelihood based tomographic image reconstruction.

Study of the maximum likelihood by EM algorithm (ML) with a reconstruction kernel equal to the intrinsic detector resolution and sieve regularization has demonstrated that any image improvements over filtered backprojection (FBP) are a function of image resolution. Comparing different reconstruction algorithms potentially requires measuring and matching the image resolution. Since there are no standard methods for describing the resolution of images from a nonlinear algorithm such as ML, we have defined measures of effective local Gaussian resolution (ELGR) and effective global Gaussian resolution (EGGR) and examined their behaviour in FBP images and in ML images using two different measurement techniques. For FBP these two resolution measures are equal and exhibit the standard convolution behaviour of linear systems. For ML, the FWHM of the ELGR monotonically increased with decreasing Gaussian object size due to slower convergence rates for smaller objects. For the simple simulated phantom used, this resolution dependence is independent of object position. With increasing object size, number of iterations and sieve size the object size dependence of the ELGR decreased. The FWHM of the EGGR converged after approximately 200 iterations, masking the fact that the ELGR for small objects was far from convergence. When FBP is compared to a nonlinear algorithm such as ML, it is recommended that at least the EGGR be matched; for ML this requires more than the number of iterations (e.g., < 100) that are typically run to minimize the mean square error or to satisfy a feasibility or similar stopping criterion. For many tasks, matching the EGGR of ML to FBP images may be insufficient and >> 200 iterations may be needed, particularly for small objects in the ML image because their ELGR has not yet converged.

Algorithms↗

PubMed1995

A Monte Carlo investigation of dual-energy-window scatter correction for volume-of-interest quantification in 99Tcm SPECT.

Using Monte Carlo simulation of 99Tcm single-photon-emission computed tomography (SPECT), we investigate the effects of tissue-background activity, tumour location, patient size, uncertainty of energy windows, and definition of tumour region on the accuracy of quantification. The dual-energy-window method of correction for Compton scattering is employed and the multiplier which yields correct activity for the VI as a whole calculated. The model is usually a sphere containing radioactive water located within a cylinder filled with a more dilute solution of radioactivity. Two simulation codes are employed. Reconstruction is by ML-EM algorithm with attenuation compensation. The scatter multiplier depends only slightly on the sphere location or the cylinder diameter. It also depends little on whether correction is before or after reconstruction. At low background level, it changes with VOI size, but not at higher background. For a geometrical VOI, it is 1.25 at zero background, decreases sharply to 0.56 for equal concentrations, and is 0.44 when the background concentration is very large. Quantification is accurate (less than 9% error) if the test background is reasonably close to that used in setting the universal scatter-multiplier value, or if the test backgrounds are always large and so is the universal-value background, but not if the test backgrounds cover a large range of values including zero. Results largely agree with those from experiment after the experimental data with background is re-evaluated with prejudice.

Humans↗

PubMed1997

Validation of the central-ray approximation for attenuated depth-dependent convolution in quantitative SPECT reconstruction.

In order to model photon attenuation and detector resolution variation as a depth-dependent convolution for efficient reconstruction of quantitative SPECT, a central-ray approximation is necessary. This work investigates the impact of the approximation upon reconstruction accuracy and computational efficiency. A patient chest CT image was acquired and converted into an object-specific attenuation map. From a segmentation of the map, an emission thorax phantom was constructed with a cardiac insert. To generate a system-specific resolution-variant kernal, a point source was measured at several depths from the surface of a low-energy, high-resolution, parallel-hole collimator of a SPECT system. Projections of parallel-beam geometry were simulated from the phantom, the map, and the kernel on an elliptical orbit. Reconstruction was performed by the ML-EM algorithm with and without the central-ray approximation. The approximation cuts down dramatically (more than 100 fold) the computing time with a negligible loss (less than 1%) of reconstruction accuracy.

Algorithms↗

PubMed1997

Scatter compensation methods in 3D iterative SPECT reconstruction: a simulation study.

Effects of different scatter compensation methods incorporated in fully 3D iterative reconstruction are investigated. The methods are: (i) the inclusion of an 'ideal scatter estimate' (ISE); (ii) like (i) but with a noiseless scatter estimate (ISE-NF); (iii) incorporation of scatter in the point spread function during iterative reconstruction ('ideal scatter model', ISM); (iv) no scatter compensation (NSC); (v) ideal scatter rejection (ISR), as can be approximated by using a camera with a perfect energy resolution. The iterative method used was an ordered subset expectation maximization (OS-EM) algorithm. A cylinder containing small cold spheres was used to calculate contrast-to-noise curves. For a brain study, global errors between reconstruction and 'true' distributions were calculated. Results show that ISR is superior to all other methods. In all cases considered, ISM is superior to ISE and performs approximately as well as (brain study) or better than (cylinder data) ISE-NF. Both ISM and ISE improve contrast-to-noise curves and reduce global errors, compared with NSC. In the case of ISE, blurring of the scatter estimate with a Gaussian kernel results in slightly reduced errors in brain studies, especially at low count levels. The optimal Gaussian kernel size is strongly dependent on the noise level.

Algorithms↗

PubMed1999

Megavoltage CT on a tomotherapy system.

A megavoltage computed tomography (MVCT) system was developed on the University of Wisconsin tomotherapy benchtop. This system can operate either axially or helically, and collect transmission data without any bounds on delivered dose. Scan times as low as 12 s per slice are possible, and scans were run with linac output rates of 100 MU min(-1), although the system can be tuned to deliver arbitrarily low dose rates. Images were reconstructed with clinically reasonable doses ranging from 8 to 12 cGy. These images delineate contrasts below 2% and resolutions of 3.0 mm. Thus, the MVCT image quality of this system should be sufficient for verifying the patient's position and anatomy prior to radiotherapy. Additionally, synthetic data were used to test the potential for improved MVCT contrast using maximum-likelihood (ML) reconstruction. Specifically, the maximum-likelihood expectation-maximization (ML-EM) algorithm and a transmission ML algorithm were compared with filtered backprojection (FBP). It was found that for expected clinical MVCT doses enough imaging photons are used such that little benefit is conferred by the improved noise model of ML algorithms. For significantly lower doses, some quantitative improvement is achieved through ML reconstruction. Nonetheless, the image quality at those lower doses is not satisfactory for radiotherapy verification.

Algorithms↗

PubMed2005

A Monte Carlo and physical phantom evaluation of quantitative In-111 SPECT.

Accurate estimation of the 3D in vivo activity distribution is important for dose estimation in targeted radionuclide therapy (TRT). Although SPECT can potentially provide such estimates, SPECT without compensation for image degrading factors is not quantitatively accurate. In this work, we evaluated quantitative SPECT (QSPECT) reconstruction methods that include compensation for various physical effects. Experimental projection data were obtained using a GE VH/Hawkeye system and an RSD torso phantom. Known activities of In-111 chloride were placed in the lungs, liver, heart, background and two spherical compartments with inner diameters of 22 mm and 34 mm. The 3D NCAT phantom with organ activities based on clinically derived In-111 ibritumomab tiuxetan data was used for the Monte Carlo (MC) simulation studies. Low-noise projection data were simulated using previously validated MC simulation methods. Fifty sets of noisy projections with realistic count levels were generated. Reconstructions were performed using the OS-EM algorithm with various combinations of attenuation (A), scatter (S), geometric response (G), collimator-detector response (D) and partial volume compensation (PVC). The QSPECT images from the various combinations of compensations were evaluated in terms of the accuracy and precision of the estimates of the total activity in each organ. For experimental data, the errors in organ activities for ADS and PVC compensation were less than 6.5% except the smaller sphere (-11.9%). For the noisy simulated data, the errors in organ activity for ADS compensation were less than 5.5% except the lungs (20.9%) and blood vessels (15.2%). Errors for other combinations of compensations were significantly (A, AS) or somewhat (AGS) larger. With added PVC, the error in the organ activities improved slightly except for the lungs (11.5%) and blood vessels (3.6%) where the improvement was more substantial. The standard deviation/mean ratios were all less than 1.5%. We conclude that QSPECT methods with appropriate compensations provided accurate In-111 organ activity estimates. For the collimator used, AGS was almost as good as ADS and may be preferable due to the reduced reconstruction time. PVC was important for small structures such as tumours or for organs in close proximity to regions with high activity. The improved quantitative accuracy from QSPECT methods has the potential for improving organ dose estimations in TRT.

Algorithms↗

PubMed2001

Detection of recombination in DNA multiple alignments with hidden Markov models.

Conventional phylogenetic tree estimation methods assume that all sites in a DNA multiple alignment have the same evolutionary history. This assumption is violated in data sets from certain bacteria and viruses due to recombination, a process that leads to the creation of mosaic sequences from different strains and, if undetected, causes systematic errors in phylogenetic tree estimation. In the current work, a hidden Markov model (HMM) is employed to detect recombination events in multiple alignments of DNA sequences. The emission probabilities in a given state are determined by the branching order (topology) and the branch lengths of the respective phylogenetic tree, while the transition probabilities depend on the global recombination probability. The present study improves on an earlier heuristic parameter optimization scheme and shows how the branch lengths and the recombination probability can be optimized in a maximum likelihood sense by applying the expectation maximization (EM) algorithm. The novel algorithm is tested on a synthetic benchmark problem and is found to clearly outperform the earlier heuristic approach. The paper concludes with an application of this scheme to a DNA sequence alignment of the argF gene from four Neisseria strains, where a likely recombination event is clearly detected.

Algorithms↗

PubMed2003

RMDP: a dedicated package for 131I SPECT quantification, registration and patient-specific dosimetry.

The limitations of traditional targeted radionuclide therapy (TRT) dosimetry can be overcome by using voxel-based techniques. All dosimetry techniques are reliant on a sequence of quantitative emission and transmission data. The use of (131)I, for example, with NaI or mIBG, presents additional quantification challenges beyond those encountered in low-energy NM diagnostic imaging, including dead-time correction and additional photon scatter and penetration in the camera head. The Royal Marsden Dosimetry Package (RMDP) offers a complete package for the accurate processing and analysis of raw emission and transmission patient data. Quantitative SPECT reconstruction is possible using either FBP or OS-EM algorithms. Manual, marker- or voxel-based registration can be used to register images from different modalities and the sequence of SPECT studies required for 3-D dosimetry calculations. The 3-D patient-specific dosimetry routines, using either a beta-kernel or voxel S-factor, are included. Phase-fitting each voxel's activity series enables more robust maps to be generated in the presence of imaging noise, such as is encountered during late, low-count scans or when there is significant redistribution within the VOI between scans. Error analysis can be applied to each generated dose-map. Patients receiving (131)I-mIBG, (131)I-NaI, and (186)Re-HEDP therapies have been analyzed using RMDP. A Monte-Carlo package, developed specifically to address the problems of (131)I quantification by including full photon interactions in a hexagonal-hole collimator and the gamma camera crystal, has been included in the dosimetry package. It is hoped that the addition of this code will lead to improved (131)I image quantification and will contribute towards more accurate 3-D dosimetry.

Humans↗