Bias of maximum likelihood estimator of intraclass correlation

Summary A bias correction was derived for the maximum likelihood estimator (MLE) of the intraclass correlation. The bias consisted of two parts: a correction from MLE to the analysis of variance estimator (ANOVA) and the bias of ANOVA. The total possible bias was always negative and depended upon both the degree of correlation and the design size and balance. The first part of the bias was an exact algebraic expression from MLE to ANOVA, and the corrected estimator by this part was ANOVA. It was also shown that the first correction term was equivalent to Fisher's reciprocal bias correction on hisZ scores. The total possible bias of MLE was large for small and moderate samples. Relative biases were larger for small parametric values and vice versa. To ensure a relative bias less than 10% assuming an intraclass correlation of 0.025, which is not unusual in most of the animal genetic studies, the total number of observations (N) should be not less than 500. From a design point of view, minimum bias occurred atn = 2, the minimum family size possible, underN fixed.     

Estimating local interaction from spatiotemporal forest data, and Monte Carlo bias correction

We point out a general problem in fitting continuous time spatially explicit models to a temporal sequence of spatial data observed at discrete times. To illustrate the problem, we examined the continuous time Markov model for forest gap dynamics. A forest is assumed to be apportioned into discrete cells (or sites) arranged in a regular square lattice. Each site is characterized as either a gap or a non-gap site according to the vegetation height of trees. The model incorporates the influence of neighboring sites on transition rate: transition rate from a non-gap to a gap site increases linearly with the number of neighbors that are currently in the gap state, and vice versa. We fitted the model to the spatiotemporal data of canopy height observed at the permanent plot in Barro Colorado Island (BCI). When we used the approximate maximum likelihood method to estimate the parameters of the model, the estimated transition rates included a large bias-in particular, the strength of interaction between nearby sites was underestimated. This bias originated from the assumption that each transition between two observation times is independent. The interaction between sites at local scale creates a long chain of transitions within a single census interval, which violates the independence of each transition. We show that a computer-intensive method, called Monte Carlo bias correction (MCBC), is very effective in removing the bias included in the estimate. The global and local gap densities measuring spatial aggregation of gap sites were computed from simulated and real gap dynamics to assess the model. When the approximate likelihood estimates were applied to the model, the predicted local gap density was clearly lower than the observed one. The use of MCBC estimates, suggesting a strong interaction between sites, improved this discrepancy.     

Estimating rate constants from single ion channel currents when the initial distribution is known

Single ion channel currents can be analysed by hidden or aggregated Markov models. A classical result from Fredkin et al. (Proceedings of the Berkeley conference in honor of Jerzy Neyman and Jack Kiefer, vol I, pp 269–289, 1985) states that the maximum number of identifiable parameters is bounded by 2n_on_c, where n_o and n_c denote the number of open and closed states, respectively. We show that this bound can be overcome when the probabilities of the initial distribution are known and the data consist of several sweeps.     

Fitting parameters of stochastic birth-death models to metapopulation data

Populations that are structured into small local patches are a common feature of ecological and epidemiological systems. Models describing this structure are often referred to as metapopulation models in ecology or household models in epidemiology. Small local populations are subject to demographic stochasticity. Theoretical studies of household disease models without resistant stages (SIS models) have shown that local stochasticity can be ignored for between patch disease transmission if the number of connected patches is large. In that case the distribution of the number of infected individuals per household reaches a stationary distribution described by a birth-death process with a constant immigration term. Here we show how this result, in conjunction with the balancing condition for birth-death processes, provides a framework to estimate demographic parameters from a frequency distribution of local population sizes. The parameter estimation framework is applicable to estimate parameters of disease transmission models as well as metapopulation models.     

Estimating phylogenetic trees from pairwise likelihoods and posterior probabilities of substitution counts

The field of phylogenetic tree estimation has been dominated by three broad classes of methods: distance-based approaches, parsimony and likelihood-based methods (including maximum likelihood (ML) and Bayesian approaches). Here we introduce two new approaches to tree inference: pairwise likelihood estimation and a distance-based method that estimates the number of substitutions along the paths through the tree. Our results include the derivation of the formulae for the probability that two leaves will be identical at a site given a number of substitutions along the path connecting them. We also derive the posterior probability of the number of substitutions along a path between two sequences. The calculations for the posterior probabilities are exact for group-based, symmetric models of character evolution, but are only approximate for more general models.