Size and Power of Test Statistics for Gene Correlation in 2 × 2 Contingency Tables

Monte Carlo simulation of size and power of two proposed tests for linkage disequilibrium between two genes each with two alleles were investigated. Results were compared with two commonly used statistics, the correlation coefficient r and the log-odds ratio tests. Depending on the sign of the linkage disequilibrium, the new tests were found to be more powerful than either of the correlation or log-odds ratio tests. However, on average (positive and negative linkage disequilibrium) the Chi-square test using the correlation coefficient was to a small extent more powerful than the other tests.     

Improved exact confidence intervals for the odds ratio in two independent binomial samples

For two independent binomial samples, the usual exact confidence interval for the odds ratio based on the conditional approach can be very conservative. Recently, Agresti and Min (2002) showed that the unconditional intervals are preferable to conditional intervals with small sample sizes. We use the unconditional approach to obtain a modified interval, which has shorter length, and its coverage probability is closer to and at least the nominal confidence coefficient.     

Estimating the Weighting Components of a Health Quality Index

This paper introduces a maximum likelihood estimation model that uses patients' pretherapy characteristics to determine the weighting components of multipe outcomes of quality health. The resulting health quality index' can then be included in cost-effectiveness analyses to compare alternative treatment plans. A Monte Carlo simulation is performed to test the model's robustness and power.     

Comparison of statistical methods for the analysis of limiting dilution assays

Summary This study reports the results of a critical comparison of five statistical methods for estimating the density of viable cells in a limiting dilution assay (LDA). Artificial data were generated using Monte Carlo simulation. The performance of each statistical method was examined with respect to the accuracy of its estimator and, most importantly, the accuracy of its associated estimated standard error (SE). The regression method was found to perform at a level that is unacceptable for scientific research, due primarily to gross underestimation of the SE. The maximum likelihood method exhibited the best overall performance. A corrected version of Taswell's weighted-mean method, which provides the best performance among all noniterative methods examined, is also presented.     

Testing the Equality of Dependent Variances

The likelihood ratio test for testing equality of vgE;2 correlated variables is developed. In general, evaluation of the test statistic involves the iterative optimization of a likelihood function with 1 + v(v – 1)/2 parameters. The explicit form of the test statistic is derived in the bivariate case, and an iterative algorithm for determining the maximum likelihood estimates is suggested. A limited Monte Carlo study determines the behavior of the proposed procedure under the null hypothesis and variety of parameter values.     

Test Statistics for Three-Factor Interaction in 2 × 2 × 2 Contingency Tables

A Monte Carlo simulation was conducted in order to determine the size and power of two proposed tests (the covariance and correlation tests) for three-factor interaction in 2 × 2 × 2 contingency tables. Results were compared to the log-odds ratio test statistic. Simulation showed the correlation test to be more conservative than the covariance test, but less so than the log-odds ratio test. However, the correlation test was the most powerful among the three tests.     

Importance of anisotropy in detachment rates for force production and cargo transport by a team of motor proteins

Many cellular processes are driven by collective forces generated by a team consisting of multiple molecular motor proteins. One aspect that has received less attention is the detachment rate of molecular motors under mechanical force/load. While detachment rate of kinesin motors measured under backward force increases rapidly for forces beyond stall‐force; this scenario is just reversed for non‐yeast dynein motors where detachment rate from microtubule decreases, exhibiting a catch‐bond type behavior. It has been shown recently that yeast dynein responds anisotropically to applied load, i.e. detachment rates are different under forward and backward pulling. Here, we use computational modeling to show that these anisotropic detachment rates might help yeast dynein motors to improve their collective force generation in the absence of catch‐bond behavior. We further show that the travel distance of cargos would be longer if detachment rates are anisotropic. Our results suggest that anisotropic detachment rates could be an alternative strategy for motors to improve the transport properties and force production by the team.     

Pedigree-based assignment tests for reversing coyote (Canis latrans) introgression into the wild red wolf (Canis rufus) population

The principal threat to the persistence of the endangered red wolf (Canis rufus) in the wild is hybridization with the coyote (Canis latrans). To facilitate idengification and removal of hybrids, assignment tests are developed which use genotype data to estimate identity as coyote, 1/4, 1/2, 3/4 or full red wolf. The tests use genotypes from the red wolves that founded the surviving population and the resulting pedigree, rather than a contemporary red wolf sample. The tests are evaluated by analysing both captive red wolves at 18 microsatellite loci, and data simulated under a highly parameterized, biologically reasonable model. The accuracy of assignment rates are generally high, with over 95% of known red wolves idengified correctly. There are, however, tradeoffs between ambiguous assignments and misassignments, and between misidengifying red wolves as hybrids and hybrids as red wolves. These result in a compromise between limiting introgression and avoiding demographic losses. The management priorities and level of introgression determine the combination of test and removal strategy that best balances these tradeoffs. Ultimately, we conclude that the use of the assignment tests has the capacity to arrest and reverse introgression. To our knowledge, the presented approach is novel in that it accounts for genetic drift when the genotypes under analysis are temporally separated from the reference populations to which they are being assigned. These methods may be valuable in cases where reference databases for small populations have aged substantially, pedigree information is available or data are generated from historical samples.     

最大似然法和贝叶斯推论与似然比检验探讨泽泻目系统发育关系

应用最大似然法(ML)、贝叶斯推论(BI)、邻接法(NJ)和似然比检验(hLRTs)进行泽泻目分子系统学研究。所用的rbcL基因序列代表了泽泻目14科46属以及作为外类群的6相关属。研究结果表明,*等级制似然比检验表明泽泻目rbcL序列最适合的DNA进化模型为GTR+I+G,最大似然法、贝叶斯法和邻接法构建的系统发育树拓扑结构相似,没有显著的差异,但贝叶斯树支持率较高;泽泻目为一单系类群,由两个主要谱系分支构成,深层分布格局由5个主要分支构成。基于分子系统发育树,文中对泽泻目科间、水鳖科+茨藻科、泽泻科+花蔺科+黄花蔺科、和Cymodoeaceae complex的系统发育关系进行了讨论。研究结果还表明,泽泻目系统发育关系可能还需要更多的证据进一步的澄清。     

Test of the Inverse Monte Carlo Method for the Calculation of Interatomic Potential Energies in Atomic Liquids

Abstract

The principle purpose of this paper is to demonstrate the use of the Inverse Monte Carlo technique for calculating pair interaction energies in monoatomic liquids from a given equilibrium property. This method is based on the mathematical relation between transition probability and pair potential given by the fundamental equation of the “importance sampling” Monte Carlo method. In order to have well defined conditions for the test of the Inverse Monte Carlo method a Metropolis Monte Carlo simulation of a Lennard Jones liquid is carried out to give the equilibrium pair correlation function determined by the assumed potential. Because an equilibrium configuration is prerequisite for an Inverse Monte Carlo simulation a model system is generated reproducing the pair correlation function, which has been calculated by the Metropolis Monte Carlo simulation and therefore representing the system in thermal equilibrium. This configuration is used to simulate virtual atom displacements. The resulting changes in atom distribution for each single simulation step are inserted in a set of non-linear equations defining the transition probability for the virtual change of configuration. The solution of the set of equations for pair interaction energies yields the Lennard Jones potential by which the equilibrium configuration has been determined.