R中使用卡方检验(GOF)的正确方法

Question

我正在尝试确定今天的数据是否在四个类别上与昨天的数据真的不同。

我统计的数据是：

data <-data.frame(yesterday=c(10741, 1575, 174, 2),
              today = c(11987, 1705, 211, 2), 
              row.names = c("a", "b", "c", "unknown"))

> data
        yesterday today
a           10741 11987
b            1575  1705
c             174   211
unknown         2     2

因此，我使用stats包中的卡方进行测试，如下所示：

stats::chisq.test(x = data$yesterday, y = data$today)

结果是：

Pearson's Chi-squared test

data:  data$yesterday and data$today
X-squared = 12, df = 9, p-value = 0.2133

我的问题是，我假设这应该是相同的：

stats::chisq.test(data)

但是你可以看到结果是完全不同的。

Pearson's Chi-squared test

data:  data
X-squared = 1.3846, df = 3, p-value = 0.7092

so....which是使用此测试比较来自同一数据集的两个样本的正确方法吗？

Answer 1

我假设问题在于您在第一种情况下对列联表的列应用了chisq.test，而函数期望x和y是因数。因此，提供列联表的版本应该是正确的，至少它对应于文档中的示例