The R Project SVN R

Rev

Rev 61153 | Rev 88548 | Go to most recent revision | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed

Rev 61153 Rev 61168
Line 37... Line 37...
37 
n <- length(dr <- c(VADeaths))
 37 
n <- length(dr <- c(VADeaths))
38 
nam <- names(VADeaths)
 38 
nam <- names(VADeaths)
39 
d.VAD <- data.frame(
 39 
d.VAD <- data.frame(
40 
 Drate = dr,
 40 
 Drate = dr,
41 
 age = rep(ordered(rownames(VADeaths)), length.out = n),
 41 
 age = rep(ordered(rownames(VADeaths)), length.out = n),
42 
 gender= gl(2,5,n, labels= c("M", "F")),
 42 
 gender = gl(2, 5, n, labels = c("M", "F")),
43 
 site =  gl(2,10,  labels= c("rural", "urban")))
 43 
 site =  gl(2, 10, labels = c("rural", "urban")))
44 
coplot(Drate ~ as.numeric(age) | gender * site, data = d.VAD,
 44 
coplot(Drate ~ as.numeric(age) | gender * site, data = d.VAD,
45 
       panel = panel.smooth, xlab = "VADeaths data - Given: gender")
 45 
       panel = panel.smooth, xlab = "VADeaths data - Given: gender")
46 
summary(aov.VAD <- aov(Drate ~ .^2, data = d.VAD))
 46 
summary(aov.VAD <- aov(Drate ~ .^2, data = d.VAD))
47 
opar <- par(mfrow = c(2,2), oma = c(0, 0, 1.1, 0))
 47 
opar <- par(mfrow = c(2, 2), oma = c(0, 0, 1.1, 0))
48 
plot(aov.VAD)
 48 
plot(aov.VAD)
49 
par(opar)
 49 
par(opar)
50 
}
 50 
}
51 
\keyword{datasets}
 51 
\keyword{datasets}