On this Picostat.com statistics page, you will find information about the occupationalStatus data set which pertains to Occupational Status of Fathers and their Sons. The occupationalStatus data set is found in the datasets R package. You can load the occupationalStatus data set in R by issuing the following command at the console data("occupationalStatus"). This will load the data into a variable called occupationalStatus. If R says the occupationalStatus 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 occupationalStatus R data set. The size of this file is about 223 bytes.
Occupational Status of Fathers and their Sons
Cross-classification of a sample of British males according to each
subject's occupational status and his father's occupational status.
table of counts, with classifying factors
origin (father's occupational status; levels
destination (son's occupational status; levels
Goodman, L. A. (1979)
Simple Models for the Analysis of Association in Cross-Classifications
having Ordered Categories.
J. Am. Stat. Assoc., 74 (367), 537–552.
The data set has been in package gnm and been provided by the
require(stats); require(graphics)plot(occupationalStatus)## Fit a uniform association model separating diagonal effects
Diag <- as.factor(diag(1:8))
Rscore <- scale(as.numeric(row(occupationalStatus)), scale = FALSE)
Cscore <- scale(as.numeric(col(occupationalStatus)), scale = FALSE)
modUnif <- glm(Freq ~ origin + destination + Diag + Rscore:Cscore,
family = poisson, data = occupationalStatus)summary(modUnif)
plot(modUnif) # 4 plots, with warning about h_ii ~= 1
Dataset imported from https://www.r-project.org.