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

Question

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

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

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

我收到以下错误：

Error in x[good, , drop = FALSE] : (subscript) logical subscript too long

知道如何在 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)

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

Test for statistical difference in proportions using R

r

statistics

glm