A simulation approach for power calculation in large cohort studies based on multistate models

Realistic power calculations for large cohort studies and nested case control studies are essential for successfully answering important and complex research questions in epidemiology and clinical medicine. For this, we provide a methodical framework for general realistic power calculations via simulations that we put into practice by means of an R‐based template. We consider staggered recruitment and individual hazard rates, competing risks, interaction effects, and the misclassification of covariates. The study cohort is assembled with respect to given age‐, gender‐, and community distributions. Nested case‐control analyses with a varying number of controls enable comparisons of power with a full cohort analysis. Time‐to‐event generation under competing risks, including delayed study‐entry times, is realized on the basis of a six‐state Markov model. Incidence rates, prevalence of risk factors and prefixed hazard ratios allow for the assignment of age‐dependent transition rates given in the form of Cox models. These provide the basis for a central simulation‐algorithm, which is used for the generation of sample paths of the underlying time‐inhomogeneous Markov processes. With the inclusion of frailty terms into the Cox models the Markov property is specifically biased. An “individual Markov process given frailty” creates some unobserved heterogeneity between individuals. Different left‐truncation‐ and right‐censoring patterns call for the use of Cox models for data analysis. p‐values are recorded over repeated simulation runs to allow for the desired power calculations. For illustration, we consider scenarios with a “testing” character as well as realistic scenarios. This enables the validation of a correct implementation of theoretical concepts and concrete sample size recommendations against an actual epidemiological background, here given with possible substudy designs within the German National Cohort.     

Host-Parasite Interactions and Purifying Selection in a Microsporidian Parasite of Honey Bees

To clarify the mechanisms of Nosema ceranae parasitism, we deep-sequenced both honey bee host and parasite mRNAs throughout a complete 6-day infection cycle. By time-series analysis, 1122 parasite genes were significantly differently expressed during the reproduction cycle, clustering into 4 expression patterns. We found reactive mitochondrial oxygen species modulator 1 of the host to be significantly down regulated during the entire infection period. Our data support the hypothesis that apoptosis of honey bee cells was suppressed during infection. We further analyzed genome-wide genetic diversity of this parasite by comparing samples collected from the same site in 2007 and 2013. The number of SNP positions per gene and the proportion of non-synonymous substitutions per gene were significantly reduced over this time period, suggesting purifying selection on the parasite genome and supporting the hypothesis that a subset of N. ceranae strains might be dominating infection.     

Effect of feeding boxes on the behavior of stereotyping amur tigers (Panthera tigris altaica) in the Zurich Zoo,Zurich, Switzerland

The stereotyped pacing shown by the two Amur tigers in the Zurich Zoo was hypothesized as being caused by permanently frustrated appetitive foraging behavior. Several electrically controlled feeding boxes were installed and access to each box was possible only twice a day for 15 min at semi‐random times. The boxes had to be opened actively by the tigers. Two trials were carried out: one with solitary confinement, and one with paired confinement. During box feeding, the female's stereotyped pacing was significantly reduced from 16% (solitary confinement, conventional feeding) and 7% (paired confinement, conventional feeding) to 1% (solitary confinement) and less than 0.01% (paired confinement) of the daily observed time. The female's sleeping increased significantly in both solitary and paired confinement. The male only showed a significant reduction in stereotyped pacing behavior when kept with the female (conventional feeding: 10%; box feeding: <0.01% of the daily observed time). On days with a box‐feeding regime in paired confinement, the male spent 25% (83 min) of the observed time with active behavior at the feeding boxes. The results support the hypothesis that permanently frustrated appetitive foraging behavior causes stereotyped pacing in adult tigers. Zoo Biol 21:573–584, 2002. © 2002 Wiley‐Liss, Inc.