\name{muscatine} \alias{muscatine} \docType{data} \title{Muscatine} \description{This example is a subset of data from the Muscatine Coronary Risk Factor Study, a longitudinal study of coronary risk factors in school children from Muscatine (Iowa, USA). } \usage{data(muscatine)} \format{ A data frame with 156 observations on the following 7 variables. \describe{ \item{\code{id}}{identifies de number of the individual profile. This vector contains observation of 52 individuals.} \item{\code{obese}}{a numeric vector that identify the obesity status (1="yes", 0="no") of a child at each occasion. } \item{\code{sex}}{a factor with levels \code{1} for "female" and \code{0} for "male". } \item{\code{time}}{a numeric vector (1,2,3) indicating the observed time points.} \item{\code{counts}}{a numeric vector indicating the number of times that each profile is replicated.} } } \details{The data set presented by Fitzmaurice, Laird and Lipsitz (1994) contains records on 1014 children who were 7-9 years old in 1977 and were examined in 1977, 1979 and 1981. Height and weight were measured in each survey year, and those with relative weight greater than 110% of the median weight in their respective stratum were classified as obese. The binary response of interest is whether the child is obese (1) or not (0). However, many data records are incomplete, since not all children participate in all the surveys. This data set was also analyzed by Azzalini (1994).} \source{ Fitzmaurice, G. M., Laird, N. M. and Lipsitz, S. R. (1994). Analyzing incomplete longitudinal binary responses: a likelihood based approach. \emph{Biometrics}, 38, 602-612.} \references{Azzalini, A. (1994). Logistic regression for autocorrelated data with application to repeated measures. \emph{Biometrika}, 81, 767-775.} \examples{ str(muscatine) # we decompose the time effect in orthogonal components muscatine$time1 <- c(-1, 0, 1) muscatine$time2 <- c(1, -2, 1) # second order Markov Chain without random effects musc2 <- bild(obese~(time1+time2)*sex, data=muscatine, time="time1", aggregate=sex, trace=TRUE, dependence="MC2") summary(musc2) getAIC(musc2) getLogLik(musc2) } \keyword{datasets}