使用 R 检验比例的统计差异

Question

我正在尝试解决以下问题：

一个人可以分为 GroupA、GroupB 或 GroupC。

我想知道属性 1（或属性 2）如何影响这些组中的观察比例。请注意，attribute1:attribute2 具有 1:N 的关系。 Attribute1 有五个可能的值，A,B,C,D,E，而 attribute2 有两个可能的值：A,B。

模拟数据：

obsGroupA <- round(runif(40, 240, 63535))
obsGroupB <- round(runif(40, 2478, 95063))
obsGroupC <- round(runif(40, 3102, 104799))
propGroupA <- obsGroupA/(obsGroupA + obsGroupB + obsGroupC)
propGroupB <- obsGroupB/(obsGroupA + obsGroupB + obsGroupC)
propGroupC <- obsGroupC/(obsGroupA + obsGroupB + obsGroupC)
#propGroupA + propGroupB + propGroupC
attributeA <- c("A", "B", "C", "D", "E")[runif(40, 1, 5)]
attributeB <- ifelse(attributeA %in% c("A", "B", "E"), "A", "B")  

模型尝试：

#y <- cbind(obsGroupA, obsGroupB, obsGroupC)
y <- cbind(propGroupA, propGroupB, propGroupC)
model <- glm(y ~ attributeA)

我收到以下错误：

x[good, , drop = FALSE] 中的错误：（下标）也是逻辑下标 长

知道如何在 R 中进行统计测试吗？
任何对正确统计测试的参考也将不胜感激。

谢谢。

Answer 1

嗯，你应该先看看已经评论过的回归分析。你在理解那里有一些问题。但是，这就是你想要的：

obsGroupA <- round(runif(40, 240, 63535))
obsGroupB <- round(runif(40, 2478, 95063))
obsGroupC <- round(runif(40, 3102, 104799))
propGroupA <- obsGroupA/(obsGroupA + obsGroupB + obsGroupC)
propGroupB <- obsGroupB/(obsGroupA + obsGroupB + obsGroupC)
propGroupC <- obsGroupC/(obsGroupA + obsGroupB + obsGroupC)
#propGroupA + propGroupB + propGroupC
attributeA <- c("A", "B", "C", "D", "E")[runif(40, 1, 5)]
attributeB <- ifelse(attributeA %in% c("A", "B", "E"), "A", "B")

y <- data.frame(propGroupA, propGroupB, propGroupC,attributeA,attributeB)
model <- glm(propGroupA ~ attributeA ,data=y )
summary(model)