R：两个随机变量的列联表答案

【问题标题】：R : contingency table for two random variablesR：两个随机变量的列联表
【发布时间】：2014-02-11 07:21:11
【问题描述】：

我有两个 Beta 变量 B1(8, 5) 和 B2(4, 7)，我为每个变量生成了 1000 个样本。

如果我必须使用 2-Way 列联表来检验两组样本的独立性，

1) 使用列联表是否合适？

2) 有正确方法的指针吗？

（我创建了一个频率表，但我看到的只是这个 - X 的所有样本按行排列，Y 的所有样本按列排列，每个单元格都写 0）

【问题讨论】：

标签： r statistics probability

【解决方案1】：

Beta 随机变量可以连续取 0 到 1 之间的任何值。所以简单的列联表没有多大意义。

您可以查看协方差或绘图，或对数据进行分箱，然后查看列联表。也许像

> set.seed(1)
> B1 <- rbeta(1000, shape1=8, shape2=5)
> B2 <- rbeta(1000, shape1=4, shape2=7)
> cov(B1,B2)
[1] 0.0003400774
> plot (B1,B2)
> CT <- table(cut(B1,4), cut(B2,4))
> print(CT)

                (0.0518,0.246] (0.246,0.44] (0.44,0.635] (0.635,0.829]
  (0.214,0.401]             15           30           11             3
  (0.401,0.587]             77          173           83            12
  (0.587,0.774]            106          231          126            20
  (0.774,0.96]              25           54           30             4
> chisq.test(CT)

        Pearson's Chi-squared test

data:  CT 
X-squared = 2.4747, df = 9, p-value = 0.9816

Warning message:
In chisq.test(CT) : Chi-squared approximation may be incorrect

【讨论】：

这对理解为什么它不起作用以及协方差更好的原因非常有帮助（整个代码不是我要求的，但感谢它，以及努力！）。谢谢！