R:两个随机变量的列联表

Beta随机变量可以连续取值0到1。因此,一个简单的应急表没有多大意义。

你可以查看协方差,或者一个图表,或者数据的bin,然后查看一个列联表。也许是这样

> 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