我想知道如何使用 logmult 包指定 RC(M)_L 模型。例如，我正在尝试使用以下代码复制 Goodman & Hout 1998 中的回归型层效应模型：
library(logmult)
data(Yamaguchi87)
tab<-xtabs(Freq~Son+Father+Country,Yamaguchi87)
rcL.trans(tab,nd=2,layer.effect="none", se="jackknife")
然而，结果与 Goodman & Hout 1998 中表 3 中的 M4* 模型不同。任何人都可以帮助正确指定回归类型的层效应模型。也可以推荐其他 R 包。非常感谢。
参考：Goodman, LA, & Hout, M. (1998)。5. 统计方法和图形显示，用于分析两个定性变量之间的关联在国家、群体之间或随着时间的推移如何不同：一种改进的回归型方法。社会学方法论，28（1），175-230。