> #load library > library(gnm) > > # data > clogg<-structure( + c(44,59,23,27,258, + 11,41,11,8,57, + 38,147,13,16,105, + 62,293,27,27,110), + dim = c(5, 4), + dimnames = structure(list( + att_abort = c("Error", "Yes, Yes, Yes", "Yes, Yes, No", "Yes, No, No", "No, No, No"), + att_sex = c("Always wrong", "Almost always wrong", "Sometimes wrong", "Not wrong"))), + class = "table") > > # model > nullModel <- gnm(Freq ~ att_abort + att_sex, family = poisson, data = clogg, verbose = FALSE) > nullModel Call: gnm(formula = Freq ~ att_abort + att_sex, family = poisson, data = clogg, verbose = FALSE) Coefficients: (Intercept) att_abortYes, Yes, Yes 3.8344 1.2481 att_abortYes, Yes, No att_abortYes, No, No -0.7394 -0.6867 att_abortNo, No, No att_sexAlmost always wrong 1.2295 -1.1666 att_sexSometimes wrong att_sexNot wrong -0.2534 0.2333 Deviance: 236.3379 Pearson chi-squared: 222.6175 Residual df: 12 > > rc2 <- gnm(Freq ~ att_abort + att_sex + instances(Mult(att_abort, att_sex), 2), family = poisson, data = clogg) Initialising Running start-up iterations.. Running main iterations....................................................... Done Warning message: In strsplit(specialTerms, ", inst = |,? ?\\)$", perl = TRUE) : perl = TRUE is only fully implemented in UTF-8 locales > > rc2 <- gnm(Freq ~ att_abort + att_sex + Mult(att_abort, att_sex), family = poisson, data = clogg) Initialising Running start-up iterations.. Running main iterations...... Done Warning message: In strsplit(specialTerms, ", inst = |,? ?\\)$", perl = TRUE) : perl = TRUE is only fully implemented in UTF-8 locales > rc2 Call: gnm(formula = Freq ~ att_abort + att_sex + Mult(att_abort, att_sex), family = poisson, data = clogg) Coefficients: (Intercept) att_abortYes, Yes, Yes att_abortYes, Yes, No 3.76936 1.13345 -0.74598 att_abortYes, No, No att_abortNo, No, No att_sexAlmost always wrong -0.69899 1.11114 -1.01442 att_sexSometimes wrong att_sexNot wrong Mult(., att_sex).att_abortError -0.11838 0.28147 -0.21653 Mult(., att_sex).att_abortYes, Yes, Yes Mult(., att_sex).att_abortYes, Yes, No Mult(., att_sex).att_abortYes, No, No -3.62635 0.62784 0.90373 Mult(., att_sex).att_abortNo, No, No Mult(att_abort, .).att_sexAlways wrong Mult(att_abort, .).att_sexAlmost always wrong 3.05984 0.22156 0.05569 Mult(att_abort, .).att_sexSometimes wrong Mult(att_abort, .).att_sexNot wrong -0.05131 -0.13927 Deviance: 5.557435 Pearson chi-squared: 5.445697 Residual df: 6 > > getContrasts(rc2, pickCoef(rc2, ").att_abort"), ref="first", scaleRef="first", scaleWeights="unit") $covmat Mult(., att_sex).att_abortError Mult(., att_sex).att_abortYes, Yes, Yes Mult(., att_sex).att_abortYes, Yes, No Mult(., att_sex).att_abortError 0.243016437 -0.05218705 -0.10259587 Mult(., att_sex).att_abortYes, Yes, Yes -0.052187051 0.04163889 -0.04130263 Mult(., att_sex).att_abortYes, Yes, No -0.102595865 -0.04130263 0.41915092 Mult(., att_sex).att_abortYes, No, No -0.092035029 -0.01169584 -0.16910856 Mult(., att_sex).att_abortNo, No, No 0.003477522 0.05748423 -0.09209758 Mult(., att_sex).att_abortYes, No, No Mult(., att_sex).att_abortNo, No, No Mult(., att_sex).att_abortError -0.09203503 0.003477522 Mult(., att_sex).att_abortYes, Yes, Yes -0.01169584 0.057484230 Mult(., att_sex).att_abortYes, Yes, No -0.16910856 -0.092097577 Mult(., att_sex).att_abortYes, No, No 0.38306375 -0.098551727 Mult(., att_sex).att_abortNo, No, No -0.09855173 0.116143402 $qvframe Estimate Std. Error Mult(., att_sex).att_abortError 0.0000000 0.0000000 Mult(., att_sex).att_abortYes, Yes, Yes -0.6912988 0.1512059 Mult(., att_sex).att_abortYes, Yes, No 0.1711848 0.1840496 Mult(., att_sex).att_abortYes, No, No 0.2271176 0.1761647 Mult(., att_sex).att_abortNo, No, No 0.6642435 0.1013186 $relerrs NULL $modelcall gnm(formula = Freq ~ att_abort + att_sex + Mult(att_abort, att_sex), family = poisson, data = clogg) attr(,"class") [1] "qv" > > getContrasts(rc2, pickCoef(rc2, ").att_sex"), ref="first", scaleRef="first", scaleWeights="unit") $covmat Mult(att_abort, .).att_sexAlways wrong Mult(att_abort, .).att_sexAlmost always wrong Mult(att_abort, .).att_sexAlways wrong 5.218724e-04 1.239207e-05 Mult(att_abort, .).att_sexAlmost always wrong 1.239207e-05 9.382739e-04 Mult(att_abort, .).att_sexSometimes wrong 1.526107e-05 1.072107e-05 Mult(att_abort, .).att_sexNot wrong 1.555223e-05 1.504050e-05 Mult(att_abort, .).att_sexSometimes wrong Mult(att_abort, .).att_sexNot wrong Mult(att_abort, .).att_sexAlways wrong 1.526107e-05 1.555223e-05 Mult(att_abort, .).att_sexAlmost always wrong 1.072107e-05 1.504050e-05 Mult(att_abort, .).att_sexSometimes wrong 4.053890e-04 4.897128e-05 Mult(att_abort, .).att_sexNot wrong 4.897128e-05 3.114257e-04 $qvframe Estimate Std. Error Mult(att_abort, .).att_sexAlways wrong 0.0000000 0.00000000 Mult(att_abort, .).att_sexAlmost always wrong -0.3442358 0.06058244 Mult(att_abort, .).att_sexSometimes wrong -0.5663083 0.03293807 Mult(att_abort, .).att_sexNot wrong -0.7488635 0.02974229 $relerrs NULL $modelcall gnm(formula = Freq ~ att_abort + att_sex + Mult(att_abort, att_sex), family = poisson, data = clogg) attr(,"class") [1] "qv" >