R：兩個隨機變量的列聯表

Question

我有兩個Beta變量B1（8,5）和B2（4,7）我已經爲它們中的每一個產生了1000個樣本。

如果我不得不使用一個2路的應急表來測試對兩組樣品的獨立性，

1）會是適合使用應急表？

2）任何指針正確地處理？

（I創建一個頻率表，但我看到的是這樣的 - X的所有樣品排列成行，和Y中的所有樣品佈置成列，並且每個單元寫入0）

Answer 1

貝塔隨機變量可以從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