Script R to prepare ahs.RData.

data-raw/ahs.R

+##----------------------------------------------------------------------
+## Prepare de the data set.
+
+setwd("/home/walmes/GitLab/mcglm/data-raw")
+
+ahs <- read.table("ahs.txt", header=TRUE, sep="\t")
+str(ahs)
+
+ahs$sex <- factor(ahs$sex, labels=c("male", "female"))
+str(ahs)
+
+xtabs(~chcond1+chcond2, data=ahs)
+
+ahs$chcond <- factor(with(ahs, chcond1+2*chcond2),
+                     labels=c("otherwise", "not limited", "limited"))
+xtabs(~chcond, data=ahs)
+ahs$chcond1 <- NULL
+ahs$chcond2 <- NULL
+ahs$id <- NULL
+
+names(ahs)[c(1:9, 15, 10:14)]
+ahs <- ahs[, c(1:9, 15, 10:14)]
+str(ahs)
+
+save(ahs, file="../data/ahs.RData")
+
+##----------------------------------------------------------------------
+## Include in the @examples
+
+library(lattice)
+library(latticeExtra)
+
+data(ahs, package="mcglm")
+str(ahs)
+
+xt <- xtabs(~age+sex, data=ahs)
+mosaicplot(xt)
+
+xt <- xtabs(~age+chcond, data=ahs)
+mosaicplot(xt)
+
+useOuterStrips(
+    combineLimits(
+        xyplot(Ndoc+Nndoc+Nadm+Nhosp+Nmed~age|sex,
+               outer=TRUE, data=ahs,
+               jitter.x=TRUE, amount=0.01,
+               type=c("p", "a"),
+               scales=list(y=list(relation="free")),
+               ylab="Number or occurences",
+               xlab="Age (years/100)")
+    )
+)
+
+useOuterStrips(
+    combineLimits(
+        xyplot(Ndoc+Nndoc+Nadm+Nhosp+Nmed~income|sex,
+               outer=TRUE, data=ahs,
+               jitter.x=TRUE, amount=0.01,
+               type=c("p", "a"),
+               scales=list(y=list(relation="free")),
+               ylab="Number or occurences",
+               xlab="Age (years/100)")
+    )
+)
+
+##----------------------------------------------------------------------