A program for maximum likelihood estimation and likelihood ratio tests in one-locus ABO-like systems based upon population phenotype data.

Phenotypes in an ABO-like system of a number of genetically-independent persons from a number of populations are supposed to be observed. The program which is written in FORTRAN calculates maximum likelihood estimates of gene frequencies and their standard errors in each population and in the populations taken together. Furthermore the program calculates expected values and likelihood ratio and goodness of fit chi-square tests of Hardy-Weinberg equilibrium. If several subpopulations are pooled together a likelihood ratio test of homogeneity is performed.     

Estimating odds ratios in genome scans: an approximate conditional likelihood approach

In modern whole-genome scans, the use of stringent thresholds to control the genome-wide testing error distorts the estimation process, producing estimated effect sizes that may be on average far greater in magnitude than the true effect sizes. We introduce a method, based on the estimate of genetic effect and its standard error as reported by standard statistical software, to correct for this bias in case-control association studies. Our approach is widely applicable, is far easier to implement than competing approaches, and may often be applied to published studies without access to the original data. We evaluate the performance of our approach via extensive simulations for a range of genetic models, minor allele frequencies, and genetic effect sizes. Compared to the naive estimation procedure, our approach reduces the bias and the mean squared error, especially for modest effect sizes. We also develop a principled method to construct confidence intervals for the genetic effect that acknowledges the conditioning on statistical significance. Our approach is described in the specific context of odds ratios and logistic modeling but is more widely applicable. Application to recently published data sets demonstrates the relevance of our approach to modern genome scans.     

Neurofibromatosis-1: a maximum likelihood estimation of mutation rate

Summary Methods of classical segregation analysis were applied to a sample of 129 sibships with one or more individuals affected by neurofibromatosis-1 (NF-1). The sample consists only of subjects with NF-1; all the probands had been referred for genetic counselling because of café-au-lait spots, and a diagnostic protocol was invariably applied. No deviation from the segregation ratio expected for a fully penetrant Mendelian dominant gene was observed. A maximum likelihood estimate of the proportion of sporadic cases was obtained, and the mutation rate was estimated to be 6.5×10^-5 gametes per generation (95% CI 5.0–8.1).     

A maximum likelihood approach to correlation dimension and entropy estimation

To obtain the correlation dimension and entropy from an experimental time series we derive estimators for these quantities together with expressions for their variances using a maximum likelihood approach. The validity of these expressions is supported by Monte Carlo simulations. We illustrate the use of the estimators with a local recording of atrial fibrillation obtained from a conscious dog.     

Bayesian and maximum likelihood estimation of genetic maps

There has recently been increased interest in the use of Markov Chain Monte Carlo (MCMC)-based Bayesian methods for estimating genetic maps. The advantage of these methods is that they can deal accurately with missing data and genotyping errors. Here we present an extension of the previous methods that makes the Bayesian method applicable to large data sets. We present an extensive simulation study examining the statistical properties of the method and comparing it with the likelihood method implemented in Mapmaker. We show that the Maximum A Posteriori (MAP) estimator of the genetic distances, corresponding to the maximum likelihood estimator, performs better than estimators based on the posterior expectation. We also show that while the performance is similar between Mapmaker and the MCMC-based method in the absence of genotyping errors, the MCMC-based method has a distinct advantage in the presence of genotyping errors. A similar advantage of the Bayesian method was not observed for missing data. We also re-analyse a recently published set of data from the eggplant and show that the use of the MCMC-based method leads to smaller estimates of genetic distances.     

Methods of estimation in log odds ratio regression models

McCullagh's (1984, Journal of the Royal Statistical Society, Series B 46, 250-256) approximation to the conditional maximum likelihood estimator in log odds ratio regression models is shown to have negligible asymptotic bias unless the odds ratios are large and the sample sizes in individual 2 X 2 tables are very small. In application to two sets of case-control data, it yields results virtually indistinguishable from those of the conditional analysis. A generalization of the Mantel-Haenszel estimator proposed by Davis (1985, Biometrics 41, 487-495) does not approximate the conditional results nearly as well.     

Multiple imputation for missing values through conditional Semiparametric odds ratio models

Multiple imputation is a practically useful approach to handling incompletely observed data in statistical analysis. Parameter estimation and inference based on imputed full data have been made easy by Rubin's rule for result combination. However, creating proper imputation that accommodates flexible models for statistical analysis in practice can be very challenging. We propose an imputation framework that uses conditional semiparametric odds ratio models to impute the missing values. The proposed imputation framework is more flexible and robust than the imputation approach based on the normal model. It is a compatible framework in comparison to the approach based on fully conditionally specified models. The proposed algorithms for multiple imputation through the Markov chain Monte Carlo sampling approach can be straightforwardly carried out. Simulation studies demonstrate that the proposed approach performs better than existing, commonly used imputation approaches. The proposed approach is applied to imputing missing values in bone fracture data.     

Estimation of a common odds ratio in case-control studies of familial aggregation

A new estimator of a common odds ratio is proposed for case-control studies of familial aggregation. The proposed estimator is a modification to the usual Mantel-Haenszel estimator that relies on an empirical adjustment for the within-family clustering which is typical of such designs. A simulation study shows that the estimator tends to have smaller mean squared error than the unmodified Mantel-Haenszel estimator under conditions likely to arise in practice. The construction of confidence intervals is also discussed.