Fitting complex random effect models with standard software using data augmentation: application to a study of male and female fecundability

Abstract

We discuss fitting of a complex random effect model using a standard statistical package (Stata) to carry out block-wise Gibbs sampling within a multiprocessor computing environment. The application involves a dataset concerning artificial insemination by donor (AID). Success or failure at each of 12 100 menstrual cycles is modelled with a mixed model with random effects due to woman, conception attempt within woman, semen donor, donation within donor, and physician who carries out the treatment. Given the availability of software within Stata to fit a model with single random effect, the full model can be fitted by an alternating imputation algorithm (Clayton and Rasbash, 1999) implemented with five copies of Stata running on separate processors and communicating via disk files.

Keywords

fecundability frailty models Gibbs sampling parallel algorithms random effects statistical computing

Get full access to this article

View all access options for this article.

References

Breslow N , Clayton D (1993) Approximate inference in generalized linear mixed models . Journal of the American Statistical Association, 88, 9-25 .

Clayton DG (1994) Some approaches to the analysis of recurrent event data . Statistical Methods in Medical Research, 3, 244-262 .

Clayton DG , Ecochard R (1997) Artificial insemination by donor: discrete time survival data with crossed and nested random effects . In Lin D , Fleming T eds. Proceedings of the First Seattle Symposium in Biostatistics: Survival Analysis. New York: Springer Verlag, 99-122 .

Clayton DG , Rasbash J (1999) Estimation in large crossed random effect models by data augmentation . Journal of the Royal Statistical Society (A), 162, 425-436 .

Conoway MR (1990) A random effects model for binary data . Biometrics, 46, 317-328 .

Crouchley R (1995) A random effects model for ordered categorical data . Journal of the American Statistical Association, 90, 489-498 .

Ecochard R , Clayton DG (1998) Multilevel modelling of conception in artificial insemination by donor . Statistics in Medicine, 17, 1137-1156 .

Ecochard R , Clayton DG (2000) Multivariate parametric random effect regression models for fecundability studies . Biometrics, 56, 1023-1029 .

Gelfand AE , Smith AFM (1990) Sampling based approaches to calculating marginal densities . Journal of the American Statistical Association, 85, 398-409 .

10.

Gini C (1924) Premières recherches sur la fécondabilité de la femme . Proceedings of the International Mathematical Congress, 2, 889-892 .

11.

Goldstein H , Healy MR (1995) The graphical presentation of a collection of means . Journal of the Royal Statistical Society A, 158, 175-177 .

12.

Goldstein H , Spiegelhalter DJ (1996) League tables and their limitations: Statistical issues in comparisons of institutional performance . Journal of the Royal Statistical Society A, 159, 385-443 .

13.

Marshall EC , Spiegelhalter DJ (1988) Reliability of league tables of in vitro fertilisation clinics: retrospective analysis of live birth rates . British Medical Journal, 316, 1701-1704 .

14.

Rasbash J , Goldstein H (1994) Efficient analysis of mixed hierarchical and crossed-classified random structures using a multilevel model . Journal of Educational and Behavioral Statistics, 19, 337-350 .

15.

Raudenbush SW (1994) A crossed random effects model for unbalanced data with applications in cross-sectional and longitudinal research . Journal of Educational Statistics, 18, 321-349 .

16.

Scheike TH , Jensen TK (1997) A discrete survival model with random effects: an application to time to pregnancy . Biometrics, 53, 318-329 .

17.

Schafer J (1997) Analysis of incomplete multivariate data. London: Chapman and Hall .

18.

Spiegelhalter DJ , Thomas A , Best NG (1996) Computation on Bayesian graphical models . Bayesian Statistics, 5, 407-425 .

19.

Tanner MA , Wong WH (1987) The calculation of posterior distributions by data augmention . Journal of the American Statistical Association, 82, 528-550 .