r - 如何计算R中列联表的CI

Question

我有从 akxl 列联表计算的 akxl 样本比例表。我想使用 R 计算表格样本比例的 CI。

对于填充了比例的 1 xk 表，有几种方法可以做到这一点。但是 akxl 表呢？根据计算方式，我可以按行或按列的方式对 kxl 比例表执行此操作吗？

例如，我有一个像下面这样的表；

       A     B     C     D
X  0.225 0.637 0.108 0.029
Y  0.455 0.171 0.252 0.122

from a contingency table,

    A  B  C  D
X  23 65 11  3
Y  56 21 31 15

理论上可以按如下方式获得CI吗？

mat <- matrix(c(23, 65, 11, 3, 56, 21, 31, 15),
              ncol = 4, byrow = TRUE)

MultinomCI(mat[1, ],
           conf.level = 0.95,
           method = "sisonglaz")

MultinomCI(mat[2, ],
           conf.level = 0.95,
           method = "sisonglaz")

我想知道可能有更好的方法或功能。

任何意见和答案将不胜感激。

标签： rconfidence-intervalcontingency