# [R] Help with

Andrew Criswell arc at arcriswell.com
Thu Oct 10 01:42:56 CEST 2002

```Hello All:

I hope I can get someone interested in this problem:

Agresti in "Analysis of Categorical Data," p. 289, applies a "row and column
effects model" to analyze a two-dimensional cross-classification of ordinal
data.

He got his results in either SAS or GLIM. Is there a way to replicate his
results with R?

He claims the RC model fits well with G^2(RC) = 3.57 with df = 8.  The ML
estimates for the row scores are (-1.11, -1.12, -0.37, 0.03, 1.01, 1.82) and
the column estimates are (-1.68, -0.14, 0.14, 1.41) and beta = 0.17

Below, is the data:

tbl.0 <- expand.grid(SOCIO = letters[1:6],
MENTAL = c('well', 'mild', 'moderate',
'impaired))

COUNTS <- c( 64, 57, 57, 72, 36, 21,
94, 94, 105, 141, 97, 71,
58, 54, 65, 77, 54, 54,
46, 40, 60, 94, 78, 71)

data.0 <- data.frame(tbl.0, COUNTS)

Thanks,
ANDREW

-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !)  To: r-help-request at stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._

```