创建一个只有两个变量的交叉表矩阵答案

【问题标题】：Creating a Crosstabulation matrix with just two variables创建一个只有两个变量的交叉表矩阵
【发布时间】：2014-01-19 20:54:36
【问题描述】：

如果问题过于简单，我深表歉意，但我现在才开始研究 R 和统计相关的问题。

我正在使用来自内置 infert 数据集的数据，对于分配，我需要构建一个 交叉表矩阵，其中列将是 induced 变量，而行将成为spontaneous 变量。

这是一个矩阵示例，供不熟悉infert的人参考：

> infert
    education age parity induced case spontaneous stratum pooled.stratum
1      0-5yrs  26      6       1    1           2       1              3
2      0-5yrs  42      1       1    1           0       2              1
3      0-5yrs  39      6       2    1           0       3              4
4      0-5yrs  34      4       2    1           0       4              2
5     6-11yrs  35      3       1    1           1       5             32
6     6-11yrs  36      4       2    1           1       6             36
7     6-11yrs  23      1       0    1           0       7              6
8     6-11yrs  32      2       0    1           0       8             22
9     6-11yrs  21      1       0    1           1       9              5
10    6-11yrs  28      2       0    1           0      10             19
11    6-11yrs  29      2       1    1           0      11             20
...
247   12+ yrs  29      1       0    0           1      82             43
248   12+ yrs  23      1       0    0           1      83             40

之后，我还需要对结果进行卡方列联表检验，看看这两个变量是否相关。

问题是，由于我缺乏经验，我不知道从哪里开始。谷歌没有产生相关结果，我找不到任何关于构建交叉表矩阵的方法（）或教程。

但是，一个非常相关的页面还包含我发现的卡方测试方法的信息是 this，尽管它使用外部包 gmodels。

有什么方法可以通过内置的 R 函数做到这一点，还是我必须得到gmodels？

【问题讨论】：

@MarkHeckmann 我已经在我的问题中链接到同一个页面。下次发表评论时请仔细阅读。
对您的刺激感到抱歉。也许我应该添加：“在您链接的页面上都有描述。”：attach(infert); mytable <- table(spontaneous, induced); chisq.test(mytable)

标签： r matrix

【解决方案1】：

您可以使用table 创建列联表：

freq <- table(infert[c("induced", "spontaneous")])
#        spontaneous
# induced  0  1  2
#       0 67 47 29
#       1 45 16  7
#       2 29  8  0

现在，您可以使用chisq.test 来测试两个变量是否独立：

chisq.test(freq)

#  Pearson's Chi-squared test
#
# data:  freq
# X-squared = 18.198, df = 4, p-value = 0.001129

【讨论】：