SEARCH · PubMed Health

Results for “Likelihood Functions”

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 505 records · Page 28Linked to original sources

PubMed2005

Rod and cone opsin families differ in spectral tuning domains but not signal transducing domains as judged by saturated evolutionary trace analysis.

The visual receptor of rods and cones is a covalent complex of the apoprotein, opsin, and the light-sensitive chromophore, 11-cis-retinal. This pigment must fulfill many functions including photoactivation, spectral tuning, signal transmission, inactivation, and chromophore regeneration. Rod and cone photoreceptors employ distinct families of opsins. Although it is well known that these opsin families provide unique ranges in spectral sensitivity, it is unclear whether the families have additional functional differences. In this study, we use evolutionary trace (ET) analysis of 188 vertebrate opsin sequences to identify functionally important sites in each opsin family. We demonstrate the following results. (1) The available vertebrate opsin sequences produce a definitive description of all five vertebrate opsin families. This is the first demonstration of sequence saturation prior to ET analysis, which we term saturated ET (SET). (2) The cone opsin classes have class-specific sites compared to the rod opsin class. These sites reside in the transmembrane region and tune the spectral sensitivity of each opsin class to its characteristic wavelength range. (3) The cytoplasmic loops, primarily responsible for signal transmission and inactivation, are essentially invariant in rod versus cone opsins. This indicates that the electrophysiological differences between rod and cone photoreceptors cannot be ascribed to differences in the protein interaction regions of the opsins. SET shows that chromophore binding and regeneration are the only aspects of opsin structure likely to have functionally significant differences between rods and cones, whereas excitatory and adaptational properties of the opsin families appear to be functionally invariant.

Amino Acid Sequence↗

PubMed2004

Molecular adaptation in plant hemoglobin, a duplicated gene involved in plant-bacteria symbiosis.

The evolutionary history of the hemoglobin gene family in angiosperms is unusual in that it involves two mechanisms known for potentially generating molecular adaptation: gene duplication and among-species interaction. In plants able to achieve symbiosis with nitrogen-fixing bacteria, class 2 hemoglobin is expressed at high concentrations in nodules and appears to be a key factor for the achievement and regulation of the symbiotic exchange. In this study, we make use of codon models of DNA sequence evolution with the goal of determining the nature of the selective forces which have driven the evolution of this gene. Our results suggest that adaptive evolution occurred during the period of time following the duplication event (functional divergence) and that a change in the selective pressures arose in class 2 hemoglobin in relation to the acquisition of a symbiotic function.

Adaptation, Biological↗

PubMedDate not supplied

A multiscale optimization approach for the dynamic contour-based boundary detection issue.

We present a new multiscale approach for deformable contour optimization. The method relies on a multigrid minimization method and a coarse-to-fine relaxation algorithm. This approach consists in minimizing a cascade of optimization problems of reduced and increasing complexity instead of considering the minimization problem on the full and original configuration space. Contrary to classical multiresolution algorithms, no reduction of image is applied. The family of defined energy functions are derived from the original (full resolution) objective function, ensuring that the same function is handled at each scale and that the energy decreases at each step of the deformable contour minimization process. The efficiency and the speed of this multiscale optimization strategy is demonstrated in the difficult context of the minimization of a region-based contour energy function ensuring the boundary detection of anatomical structures in ultrasound medical imagery. In this context, the proposed multiscale segmentation method is compared to other classical region-based segmentation approaches such as Maximum Likelihood or Markov Random Field-based segmentation techniques. We also extend this multiscale segmentation strategy to active contour models using a classical edge-based likelihood approach. Finally, time and performance analysis of this approach, compared to the (commonly used) dynamic programming-based optimization procedure, is given and allows to attest the accuracy and the speed of the proposed method.

Algorithms↗

PubMed1999

A generalization of the transmission/disequilibrium test for uncertain-haplotype transmission.

A new transmission/disequilibrium-test statistic is proposed for situations in which transmission is uncertain. Such situations arise when transmission of a multilocus marker haplotype is considered, since haplotype phase is often unknown in a substantial number of instances. Even for single-locus markers, transmission is uncertain if one or both parents are missing. In both these situations, uncertainty may be reduced by the typing of further siblings, whose disease status may be unaffected or unknown. The proposed test is a score test based on a partial score function that omits the terms most influenced by hidden population stratification.

Adult↗

PubMed2003

Evolution of the TCP gene family in Asteridae: cladistic and network approaches to understanding regulatory gene family diversification and its impact on morphological evolution.

In the plant subclass Asteridae, bilaterally symmetrical flowers have evolved from a radially symmetrical ancestral phenotype on at least three independent occasions: in the Boraginaceae, Solanaceae, and Lamiales. Development of bilateral flower symmetry has been shown to be determined by the early-acting cycloidea (cyc) and dichotoma (dich) genes in Antirrhinum, a member of the Lamiales. cyc and dich belong to the TCP gene family of putative transcription factors. TCP gene sequences were isolated from 11 Asteridae taxa using an array of degenerate PCR primers. Closely related species exhibiting either ancestral actinomorphic or derived zygomorphic flowers were sampled for each independent origin of bilateral flower symmetry. Cladistic and network-based analyses were performed to establish viable hypotheses regarding the evolution of bilateral symmetry in Asteridae. For the TCP gene family, the use of cladistic phylogenetic analysis to identify orthologous genes is complicated by a paucity of alignable data, frequent gene duplication and extinction, and the possibility of reticulate evolution via intergenic recombination. These complicating factors can be generalized to many regulatory gene families. As an alternative to cladistic analysis, we propose the use of network analysis for the reconstruction of regulatory gene family phylogenetic and functional relationships. Results of analyses support the hypothesis that the origin of bilaterally symmetrical flowers in the Boraginaceae and Solanaceae did not require orthologs or functional analogs of cyc or dich. This suggests that the genetic mechanism that determines bilateral flower symmetry in these taxa is not homologous to that of the Lamiales. Results of analyses are consistent with the hypothesis that the evolution of bilateral floral symmetry in the Lamiales required the origin of a novel gene function subsequent to gene duplication.

Amino Acid Sequence↗

PubMed2006

A simple statistical method for estimating type-II (cluster-specific) functional divergence of protein sequences.

Predicting functional amino acid residues in silico is important for comparative genomics. In this paper, we focus on the issue of how to statistically identify cluster-specific amino acid residues that are related to the functional divergence after gene duplication. We approach this problem using a framework based on site-specific shift of amino acid property (type-II functional divergence), as opposed to site-specific shift of evolutionary rate (type-I functional divergence). An efficient statistical procedure is implemented to facilitate the development of phylogenomic database for cluster-specific residues of large-scale protein families. Our method has the following features: 1) statistical testing of the type-II functional divergence and 2) the site-specific Bayesian profile to measure how amino acid residues contribute to type-II (cluster-specific) functional divergence. Consequently, one may obtain the posterior probability for "functional" cluster-specific residues. Case studies are presented and indicate that radical cluster-specific residues are responsible for most of inferred type-II functional divergence, whereas conserved cluster-specific residues appear less than even those imperfect radical cluster-specific residues to this type of functional divergence.

Amino Acid Sequence↗

PubMed2003

The gamma-frailty Poisson model for the nonparametric estimation of panel count data.

In this article, we study nonparametric estimation of the mean function of a counting process with panel observations. We introduce the gamma frailty variable to account for the intracorrelation between the panel counts of the counting process and construct a maximum pseudo-likelihood estimate with the frailty variable. Three simulated examples are given to show that this estimation procedure, while preserving the robustness and simplicity of the computation, improves the efficiency of the nonparametric maximum pseudo-likelihood estimate studied in Wellner and Zhang (2000, Annals of Statistics 28, 779-814). A real example from a bladder tumor study is used to illustrate the method.

Algorithms↗

PubMed2004

Incorporating individual error rate into association test of unmatched case-control design.

OBJECTIVES: Genotyping error commonly occurs and could reduce the power and bias statistical inference in genetics studies. In addition to genotypes, some automated biotechnologies also provide quality measurement of each individual genotype. We studied the relationship between the quality measurement and genotyping error rate. Furthermore, we propose two association tests incorporating the genotyping quality information with the goal to improve statistical power and inference. METHODS: 50 pairs of DNA sample duplicates were typed for 232 SNPs by BeadArray technology. We used scatter plot, smoothing function and generalized additive models to investigate the relationship between genotype quality score (q) and inconsistency rate (ĩ) among duplicates. We constructed two association tests: (1) weighted contingency table test (WCT) and (2) likelihood ratio test (LRT) to incorporate individual genotype error rate (epsilon(i)), in unmatched case-control setting. RESULTS: In the 50 duplicates, we found q and ĩ were in strong negative association, suggesting the genotypes with low quality score were more likely to be mistyped. The WCT improved the statistical power and partially corrects the bias in point estimation. The LRT offered moderate power gain, but was able to correct the bias in odds ratio estimation. The two new methods also performed favorably in some scenarios when epsilon(i) was mis-specified. CONCLUSIONS: With increasing number of genetic studies and application of automated genotyping technology, there is a growing need to adequately account for individual genotype error rate in statistical analysis. Our study represents an initial step to address this need and points out a promising direction for further research.

Case-Control Studies↗

PubMed1992

Longitudinal studies with continuous responses.

The analysis of serial measurements obtained in longitudinal studies plays an increasingly prominent role in applied research. The last few years have seen the development of many new techniques for carrying out analyses, including computer software. These methods can be used in a variety of standard problems, including repeated measures and cross-over designs, as well as growth curve analyses. We review these new methods, their application, and available computer packages. Data from a longitudinal study of lung function is used to illustrate the methods.

Adolescent↗

PubMed2006

Adaptive evolution of Hox-gene homeodomains after cluster duplications.

BACKGROUND: Hox genes code for homeodomain-containing transcription factors that function in cell fate determination and embryonic development. Hox genes are arranged in clusters with up to 14 genes. This archetypical chordate cluster has duplicated several times in vertebrates, once at the origin of vertebrates and once at the origin of gnathostoms, an additional duplication event is associated with the origin of teleosts and the agnanths, suggesting that duplicated Hox cluster genes are involved in the genetic mechanisms behind the diversification of vertebrate body plans, and the origin of morphological novelties. Preservation of duplicate genes is promoted by functional divergence of paralogs, either by subfunction partitioning among paralogs or the acquisition of a novel function by one paralog. But for Hox genes the mechanisms of paralog divergence is unknown, leaving open the role of Hox gene duplication in morphological evolution. RESULTS: Here, we use several complementary methods, including branch-specific dN/dS ratio tests, branch-site dN/dS ratio tests, clade level amino acid conservation/variation patterns, and relative rate ratio tests, to show that the homeodomain of Hox genes was under positive Darwinian selection after cluster duplications. CONCLUSION: Our results suggest that positive selection acted on the homeodomain immediately after Hox clusters duplications. The location of sites under positive selection in the homeodomain suggests that they are involved in protein-protein interactions. These results further suggest that adaptive evolution actively contributed to Hox-gene homeodomain functions.

Adaptation, Biological↗

PubMed2006

Application of maximum likelihood to direct methods: the probability density function of the triple-phase sums. XI.

The maximum-likelihood method is applied to direct methods to derive a more general probability density function of the triple-phase sums which is capable of predicting negative values. This study also proves that maximization of the origin-free modulus sum function S yields, within the limitations imposed by the assumed approximations, the maximum-likelihood estimates of the phases. It thus represents the formal theoretical justification of the S function that was initially derived from Patterson-function arguments [Rius (1993). Acta Cryst. A49, 406-409].

Journal Article↗

PubMed2005

Fisher information matrix of the Dirichlet-multinomial distribution.

In this paper we derive explicit expressions for the elements of the exact Fisher information matrix of the Dirichlet-multinomial distribution. We show that exact calculation is based on the beta-binomial probability function rather than that of the Dirichlet-multinomial and this makes the exact calculation quite easy. The exact results are expected to be useful for the calculation of standard errors of the maximum likelihood estimates of the beta-binomial parameters and those of the Dirichlet-multinomial parameters for data that arise in practice in toxicology and other similar fields. Standard errors of the maximum likelihood estimates of the beta-binomial parameters and those of the Dirichlet-multinomial parameters, based on the exact and the asymptotic Fisher information matrix based on the Dirichlet distribution, are obtained for a set of data from Haseman and Soares (1976), a dataset from Mosimann (1962) and a more recent dataset from Chen, Kodell, Howe and Gaylor (1991). There is substantial difference between the standard errors of the estimates based on the exact Fisher information matrix and those based on the asymptotic Fisher information matrix.

Abnormalities, Drug-Induced↗

PubMed2002

Estimation of kinetic parameters without input functions: analysis of three methods for multichannel blind identification.

Compartment modeling of dynamic medical image data implies that the concentration of the tracer over time in a particular region of the organ of interest is well modeled as a convolution of the tissue response with the tracer concentration in the blood stream. The tissue response is different for different tissues while the blood input is assumed to be the same for different tissues. The kinetic parameters characterizing the tissue responses can be estimated by multichannel blind identification methods. These algorithms use the simultaneous measurements of concentration in separate regions of the organ; if the regions have different responses, the measurement of the blood input function may not be required. Three blind identification algorithms are analyzed here to assess their utility in medical imaging: eigenvector-based algorithm for multichannel blind deconvolution; cross relations; and iterative quadratic maximum-likelihood (IQML). Comparisons of accuracy with conventional (not blind) identification techniques where the blood input is known are made as well. Tissue responses corresponding to a physiological two-compartment model are primarily considered. The statistical accuracies of estimation for the three methods are evaluated and compared for multiple parameter sets. The results show that IQML gives more accurate estimates than the other two blind identification methods.

Algorithms↗

PubMed2005

The estimation of distributions and the minimum relative entropy principle.

Estimation of Distribution Algorithms (EDA) have been proposed as an extension of genetic algorithms. In this paper we explain the relationship of EDA to algorithms developed in statistics, artificial intelligence, and statistical physics. The major design issues are discussed within a general interdisciplinary framework. It is shown that maximum entropy approximations play a crucial role. All proposed algorithms try to minimize the Kullback-Leibler divergence KLD between the unknown distribution p(x) and a class q(x) of approximations. However, the Kullback-Leibler divergence is not symmetric. Approximations which suppose that the function to be optimized is additively decomposed (ADF) minimize KLD(q||p), the methods which learn the approximate model from data minimize KLD(p||q). This minimization is identical to maximizing the log-likelihood. In the paper three classes of algorithms are discussed. FDA uses the ADF to compute an approximate factorization of the unknown distribution. The factors are marginal distributions, whose values are computed from samples. The second class is represented by the Bethe-Kikuchi approach which has recently been rediscovered in statistical physics. Here the values of the marginals are computed from a difficult constrained minimization problem. The third class learns the factorization from the data. We analyze our learning algorithm LFDA in detail. It is shown that learning is faced with two problems: first, to detect the important dependencies between the variables, and second, to create an acyclic Bayesian network of bounded clique size.

Algorithms↗

PubMed2004

A fast algorithm for functional mapping of complex traits.

By integrating the underlying developmental mechanisms for the phenotypic formation of traits into a mapping framework, functional mapping has emerged as an important statistical approach for mapping complex traits. In this note, we explore the feasibility of using the simplex algorithm as an alternative to solve the mixture-based likelihood for functional mapping of complex traits. The results from the simplex algorithm are consistent with those from the traditional EM algorithm, but the simplex algorithm has considerably reduced computational times. Moreover, because of its nonderivative nature and easy implementation with current software, the simplex algorithm enjoys an advantage over the EM algorithm in the dynamic modeling and analysis of complex traits.

Algorithms↗

PubMed1996

A Bayesian approach to Weibull survival models--application to a cancer clinical trial.

In this paper we outline a class of fully parametric proportional hazards models, in which the baseline hazard is assumed to be a power transform of the time scale, corresponding to assuming that survival times follow a Weibull distribution. Such a class of models allows for the possibility of time varying hazard rates, but assumes a constant hazard ratio. We outline how Bayesian inference proceeds for such a class of models using asymptotic approximations which require only the ability to maximize the joint log posterior density. We apply these models to a clinical trial to assess the efficacy of neutron therapy compared to conventional treatment for patients with tumours of the pelvic region. In this trial there was prior information about the log hazard ratio both in terms of elicited clinical beliefs and the results of previous studies. Finally, we consider a number of extensions to this class of models, in particular the use of alternative baseline functions, and the extension to multi-state data.

Bayes Theorem↗

PubMed2001

The estimation of statistical parameters for local alignment score distributions.

The distribution of optimal local alignment scores of random sequences plays a vital role in evaluating the statistical significance of sequence alignments. These scores can be well described by an extreme-value distribution. The distribution's parameters depend upon the scoring system employed and the random letter frequencies; in general they cannot be derived analytically, but must be estimated by curve fitting. For obtaining accurate parameter estimates, a form of the recently described 'island' method has several advantages. We describe this method in detail, and use it to investigate the functional dependence of these parameters on finite-length edge effects.

Algorithms↗

PubMed1999

Estimation in independent observer line transect surveys for clustered populations.

This paper introduces a framework for animal abundance estimation in independent observer line transect surveys of clustered populations. The framework generalizes an approach given in Chen (1999, Environmental and Ecological Statistics 6, in press) to accommodate heterogeneity in detection caused by cluster size and other covariates. Both parametric and nonparametric estimators for the local effective search widths, given the covariates, can be derived from the framework. A nonparametric estimator based on conditional kernel density estimation is proposed and studied owing to its flexibility in modeling the detection functions. A real data set on harbor porpoise in the North Sea is analyzed.

Animals↗