aids {gamlss.data} | R Documentation |
The quarterly reported AIDS cases in the U.K. from January 1983 to March 1994 obtained from the Public Health Laboratory Service, Communicable Disease Surveillance Centre, London.
data(aids)
A data frame with 45 observations on the following 3 variables.
The counts y can be modelled using a (smooth) Poisson regression model in time x with the quarterly effects i.e. cs(x,df=7)+qrt. Overdispersion persists, so use a Negative Binomial distribution of type I or II. The data also can be used to find a break point in time, see Rigby and Stasinopoulos (1992).
Public Health Laboratory Service, Communicable Disease Surveillance Centre, London.
Stasinopoulos, D.M. and Rigby, R. A. (1992). Detecting break points in generalized linear models. Computational Statistics and Data Analysis, 13, 461–471.
data(aids) attach(aids) plot(x,y,pch=21,bg=c("red","green3","blue","yellow")[unclass(qrt)]) detach(aids)