在R中是否有一种方法可以在一个应急表中指定所有不同的三方组合？
在R中是否有一种方法可以在一个应急表中指定所有不同的三方组合？
提问于 2022-09-06 06:18:01
我有许多变量(5)，我想检查所有的交互项(响应上的a_b，响应上的a_c，等等)关于多项反应。如下所示，响应有三个级别：
multinom(response~(a+b+c+d+e)^2)我希望创建应急表，以查看我的预测器和响应之间是否存在某种级别的分离。我可以手动写一些东西来获取每个组合的指示数，如下所示：
xtabs(~response+a+b)
xtabs(~response+a+c)
xtabs(~response+a+d)
...但这似乎很麻烦。是否有一种方法可以输出所有可能的三种方式组合，我可以在某种类型的表函数中使用类似于上述回归公式的公式输入？
谢谢,
回答 1
Stack Overflow用户
回答已采纳
发布于 2022-09-06 18:00:54
下面是一个使用内置UCBasmissions数据集的示例
(somedata <- as.data.frame(UCBAdmissions))
(collection <- setdiff(names(somedata), "Freq"))
(combinations <- combn(x = collection,
m=2,
simplify = FALSE))
(all_results <- lapply(combinations, \(x){xtabs(data=somedata,
formula = as.formula(paste("~Freq",x[1],x[2],sep = " + ")))}))页面原文内容由Stack Overflow提供。腾讯云小微IT领域专用引擎提供翻译支持
原文链接：
https://stackoverflow.com/questions/73617490
复制相关文章
相似问题
腾讯云开发者
