data(M3IR5)
sim <- expand.grid(dimnames(M3IR5))
sim$value <- as.vector(M3IR5)
require(lattice)
# Figure 16
xyplot(log(value) ~ dev, groups=origin, data=sim, t="l", auto.key=list(space="right"))
sim$origin <- as.numeric(sim$origin)
sim$dev <- as.numeric(sim$dev) - 1
sim$cal <- sim$origin + sim$dev - 1
# Table 3.1
summary(lm(log(value) ~ cal + dev, data=sim))
crm2 <- t(matrix(c(11073, 6427, 1839, 766, 14799, 9357, 2344, NA, 15636, 10523, NA,NA,16913, NA,NA,NA), nrow=4))
dimnames(crm2) <- list(origin=0:3, dev=0:3)
crm <- expand.grid(dimnames(crm2))
crm$value <- as.vector(crm2)
crm <- na.omit(crm)
Run the code above in your browser using DataLab