cautres {gnm} | R Documentation |
A 4-way contingency table of vote by class by religion in four French elections
data(cautres)
A table of counts, with classifying factors vote
(levels
1:2
), class
(levels 1:6
) and religion
(levels 1:4
) and election
(levels 1:4
).
Bruno Cautres
Cautres, B, Heath, A F and Firth, D (1998). Class, religion and vote in Britain and France. La Lettre de la Maison Francaise 8.
set.seed(1) data(cautres) ## Fit a "double UNIDIFF" model with the religion-vote and class-vote ## interactions both modulated by nonnegative election-specific multipliers doubleUnidiff <- gnm(Freq ~ election:vote + election:class:religion + Mult(Exp(-1 + election), religion:vote) + Mult(Exp(-1 + election), class:vote), family = poisson, data = cautres) ## Deviance should be 133.04 ## Examine the multipliers of the class-vote log odds ratios coefs.of.interest <- grep("Mult2.*election", names(coef(doubleUnidiff))) coef(doubleUnidiff)[coefs.of.interest] ## Mult2.Factor1.election1 Mult2.Factor1.election2 ## -0.5724370 0.1092972 ## Mult2.Factor1.election3 Mult2.Factor1.election4 ## -0.1230682 -0.2105843 getContrasts(doubleUnidiff, coefs.of.interest) ## estimate se ## Mult2.Factor1.election1 -0.3618399 0.2534762 ## Mult2.Factor1.election2 0.3198951 0.1320034 ## Mult2.Factor1.election3 0.0875308 0.1446842 ## Mult2.Factor1.election4 0.0000000 0.0000000 ## Same thing but with election 1 as reference category: getContrasts(doubleUnidiff, rev(coefs.of.interest)) ## estimate se ## Mult2.Factor1.election4 0.3618399 0.2534746 ## Mult2.Factor1.election3 0.4493707 0.2473524 ## Mult2.Factor1.election2 0.6817351 0.2401645 ## Mult2.Factor1.election1 0.0000000 0.0000000