A two-stage model for multiple time series data of counts

We propose a two-stage model for time series data of counts from multiple locations. This method fits first-stage model(s) using the technique of iteratively weighted filtered least squares (IWFLS) to obtain location-specific intercepts and slopes, with possible lagged effects via polynomial distributed lag modeling. These slopes and/or intercepts are then taken to a second-stage mixed-effects meta-regression model in order to stabilize results from various locations. The representation of the models from the stages into a combined mixed-effects model, issues of inference and choices of the parameters in modeling the lag structure are discussed. We illustrate this proposed model via detailed analysis on the effect of air pollution on school absenteeism based on data from the Southern California Children's Health Study.     

A semiparametric model for the analysis of recurrent-event panel data

In many longitudinal studies, interest focuses on the occurrence rate of some phenomenon for the subjects in the study. When the phenomenon is nonterminating and possibly recurring, the result is a recurrent-event data set. Examples include epileptic seizures and recurrent cancers. When the recurring event is detectable only by an expensive or invasive examination, only the number of events occurring between follow-up times may be available. This article presents a semiparametric model for such data, based on a multiplicative intensity model paired with a fully flexible nonparametric baseline intensity function. A random subject-specific effect is included in the intensity model to account for the overdispersion frequently displayed in count data. Estimators are determined from quasi-likelihood estimating functions. Because only first- and second-moment assumptions are required for quasi-likelihood, the method is more robust than those based on the specification of a full parametric likelihood. Consistency of the estimators depends only on the assumption of the proportional intensity model. The semiparametric estimators are shown to be highly efficient compared with the usual parametric estimators. As with semiparametric methods in survival analysis, the method provides useful diagnostics for specific parametric models, including a quasi-score statistic for testing specific baseline intensity functions. The techniques are used to analyze cancer recurrences and a pheromone-based mating disruption experiment in moths. A simulation study confirms that, for many practical situations, the estimators possess appropriate small-sample characteristics.     

On power and sample size calculations for likelihood ratio tests in generalized linear models

A direct extension of the approach described in Self, Mauritsen, and Ohara (1992, Biometrics 48, 31-39) for power and sample size calculations in generalized linear models is presented. The major feature of the proposed approach is that the modification accommodates both a finite and an infinite number of covariate configurations. Furthermore, for the approximation of the noncentrality of the noncentral chi-square distribution for the likelihood ratio statistic, a simplification is provided that not only reduces substantial computation but also maintains the accuracy. Simulation studies are conducted to assess the accuracy for various model configurations and covariate distributions.     

Homogeneity/Heterogeneity Hypotheses for Standardized Mortality Ratios Based on Minimum Power‐divergence Estimators

This paper analyzes the power divergence estimators when homogeneity/heterogeneity hypotheses among standardized mortality ratios (SMRs) are taken into account. A Monte Carlo study shows that when the standard mortality rate is not external, that is it is estimated from the sample data, these estimators have a good performance even for small sample sets and in particular the minimum chi‐square estimators have a better behavior compared to the classical maximum likelihood estimators. In order to make decisions under homogeneity/heterogeneity hypotheses of SMRs we propose some test‐statistics which consider the minimum power divergence estimators. Through a numerical example focused on SMRs of melanoma mortality ratios in different regions of the US, a homogeneity/heterogeneity study is illustrated.     

Testing Homogeneity of Two Zero‐inflated Poisson Populations

The problem of testing treatment difference in the occurrence of a safety parameter in a randomized parallel‐group comparative clinical trial under the assumption that the number of occurrence follows a zero‐inflated Poisson (ZIP) distribution is considered. Likelihood ratio tests (LRT) for homogeneity of two ZIP populations are derived under the hypotheses that (i) there is no difference in inflation parameters, (ii) there is no difference in non‐zero means; and (iii) there is no difference in both inflation parameters and non‐zero means. Approximate formulas for sample size calculation are also obtained for achieving a desired power for detecting a clinically meaningful difference under the corresponding alternative hypotheses. An example concerning the assessment of the gastrointestinal (GI) safety in terms of the number of erosion counts of a newly developed compound for the treatment of osteoarthritis and rheumatoid arthritis is given for illustration purpose (© 2009 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Sufficiency of the number of segregating sites in the limit under finite-sites mutation

We show that the number of segregating sites is a sufficient statistic for the scaled mutation parameter (θ) in the limit as the number of sites tends to infinity and there is free recombination between sites. We assume that the mutation parameter at each site tends to zero such than the total mutation parameter (θ) is constant in the limit. Our results show that Watterson’s estimator is the maximum likelihood estimator in this case, but that it estimates a composite parameter which is different for different mutation models. Some of our results hold when recombination is limited, because Watterson’s estimator is an unbiased, method-of-moments estimator regardless of the recombination rate. The quantity it estimates depends on the details of how mutations occur at each site.     

Bayesian Predictive Probability Functions for Count Data that are Subject to Misclassification

We develop three Bayesian predictive probability functions based on data in the form of a double sample. One Bayesian predictive probability function is for predicting the true unobservable count of interest in a future sample for a Poisson model with data subject to misclassification and two Bayesian predictive probability functions for predicting the number of misclassified counts in a current observable fallible count for an event of interest. We formulate a Gibbs sampler to calculate prediction intervals for these three unobservable random variables and apply our new predictive models to calculate prediction intervals for a real‐data example. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

On a probabilistic model for the numerical estimation of nocturnal migration of birds

The study of nocturnal bird migration by cone methods of observation has a century-long history but has continued to be used up to the present. To describe the flux and estimate the number of passing birds a probabilistic model is proposed. This model is based on the concept of dynamic Poisson ensemble of points in appropriate phase space and has two parameters. One is scalar and the other one is functional. We constructed consistent estimations of these parameters and discuss their use for the numerical estimation of the flux of birds observed in a narrow light cone generated by the bright lunar disk and formed by an open angle of telescope. Selection on the same type of birds was suggested as the necessary condition for the model application. Ground speed of each bird was introduced into the model as a new but obligatory value determining the quantification of the flux of bird.