On this Picostat.com statistics page, you will find information about the VADeaths data set which pertains to Death Rates in Virginia (1940). The VADeaths data set is found in the datasets R package. You can load the VADeaths data set in R by issuing the following command at the console data("VADeaths"). This will load the data into a variable called VADeaths. If R says the VADeaths data set is not found, you can try installing the package by issuing this command install.packages("datasets") and then attempt to reload the data. If you need to download R, you can go to the R project website. You can download a CSV (comma separated values) version of the VADeaths R data set. The size of this file is about 176 bytes.
Death Rates in Virginia (1940)
Death rates per 1000 in Virginia in 1940.
A matrix with 5 rows and 4 columns.
The death rates are measured per 1000 population per year. They
are cross-classified by age group (rows) and
population group (columns). The age groups are: 50–54, 55–59,
60–64, 65–69, 70–74 and the population groups are Rural/Male,
Rural/Female, Urban/Male and Urban/Female.
This provides a rather nice 3-way analysis of variance example.
Molyneaux, L., Gilliam, S. K., and Florant, L. C.(1947)
Differences in Virginia death rates by color, sex, age,
and rural or urban residence.
American Sociological Review, 12, 525–535.
McNeil, D. R. (1977)
Interactive Data Analysis.
n <- length(dr <- c(VADeaths))
nam <- names(VADeaths)
d.VAD <- data.frame(
Drate = dr,
age = rep(ordered(rownames(VADeaths)), length.out = n),
gender = gl(2, 5, n, labels = c("M", "F")),
site = gl(2, 10, labels = c("rural", "urban")))
coplot(Drate ~ as.numeric(age) | gender * site, data = d.VAD,
panel = panel.smooth, xlab = "VADeaths data - Given: gender")
summary(aov.VAD <- aov(Drate ~ .^2, data = d.VAD))
opar <- par(mfrow = c(2, 2), oma = c(0, 0, 1.1, 0))
Dataset imported from https://www.r-project.org.