| EVS {cmm} | R Documentation |
European Values Study 1999/2000, Views on Women's Roles http://www.europeanvaluesstudy.eu/
The data are tabulated in Bergsma, Croon, and Hagenaars (2009, Table 5.1a). Section 5.1.2 in Bergsma, Croon and Hagenaars (2009).
data(EVS)
A data frame with 960 observations on the following variables.
SAERWEuropean Values Study 1999/2000
Bergsma, W. P., Croon, M. A., & Hagenaars, J. A. P. (2009). Marginal models for dependent, clustered, and longitudinal categorical data. New York: Springer.
data(EVS)
# Table SAERW
var = c("S", "A", "E", "R", "W");
dim = c(2, 3, 3, 3, 2);
# matrices for table SA
at1 = MarginalMatrix(var, c("S", "A"), dim);
bt1 = ConstraintMatrix(c("S", "A"), list(c("S"), c("A")), c(2, 3));
# matrices for table SAER
at2 = MarginalMatrix(var, c("S", "A", "E", "R"), dim);
bt2 = ConstraintMatrix(c("S", "A", "E", "R"), list(c("S", "A", "E"), c("S", "R"), c("A", "R")), c(2, 3, 3, 3));
# matrices for table SAERW: constraints
at3 = MarginalMatrix(var, c("S", "A", "E", "R", "W"), dim);
bt3 = ConstraintMatrix(c("S", "A", "E", "R", "W"), list(c("S", "A", "E", "R"), c("S", "W"), c("A", "W"), c("E", "W"), c("R", "W")), c(2, 3, 3, 3, 2));
# matrix for table SAERW: design matrix
x = DesignMatrix(c("S", "A", "E", "R", "W"), list(c("S", "A", "E", "R"), c("S", "W"), c("A", "W"), c("E", "W"), c("R", "W")), c(2, 3, 3, 3, 2));
# model1: specification using only constraints
at = rbind(at1, at2, at3);
bt = DirectSum(bt1, bt2);
bt = DirectSum(bt, bt3);
model1 = list(bt, "log", at);
# model2: same as model1 but using both constraints and a design matrix to specify a loglinear model for the joint distribution
at = rbind(at1, at2);
bt = DirectSum(bt1, bt2);
model2 = list(list(bt, "log", at), x);
nkps3 = MarginalModelFit(EVS, model2, MaxError = 10.^-25,
MaxSteps = 1000,
ShowProgress = 10,
MaxStepSize = .7 );