На самом деле список, который возвращает chisq.test
, Xsq
, также содержит residuals
и stdres
объекты, и это также tables
. Только expected
не является таблицей. Но у меня нет проблем с использованием rbind
или cbind
для их объединения. Мы можем даже сделать массив следующим образом:
Results <- array(cbind(Xsq$observed, Xsq$expected), dim=c(2, 3, 2),
dimnames=list(gender=c("F", "M"), party=c("Democrat", "Independent",
"Republican"), c("Observed", "Expected")))
options(digits=5)
ftable(Results, row.vars=1, col.vars=2:3)
# party Democrat Independent Republican
# Observed Expected Observed Expected Observed Expected
# gender
# F 762.00 703.71 327.20 319.78 468.00 533.71
# M 484.00 542.29 239.00 246.42 477.00 411.29
ftable(Results, row.vars=1:2, col.vars=3)
# Observed Expected
# gender party
# F Democrat 762.00 703.71
# Independent 327.20 319.78
# Republican 468.00 533.71
# M Democrat 484.00 542.29
# Independent 239.00 246.42
# Republican 477.00 411.29