在 R 中生成具有双变量正态分布的 r x c 列联表

Question

我想使用 R 从双变量正态分布生成列联表。使用 rmultinom 和其他的多名义分布生成表的一种方法是 r2dtable，但我想生成交叉分类数据使用具有不同相关结构的二元正态分布。

我尝试从双变量正态分布生成样本，均值为 0，西格玛为 diag(2)。之后就不知道构造列联表了

library(MASS)

x<-mvrnorm(n=100, mu=c(0,0), Sigma=diag(2))

请帮我解决这个问题。谢谢

Answer 1

如果您想从连续数据中创建意外事件 table，您首先必须将其转换为分类数据，然后对类别进行计数。您必须定义对您的分析有意义的中断：

library(MASS)

x<-mvrnorm(n=100, mu=c(0,0), Sigma=diag(2))

x <- apply(x,FUN = function(y) {cut(y,breaks = c(-Inf,-0.5,0.5,Inf))},MARGIN = 2)

table(x[,1],x[,2])

在 R 中生成具有双变量正态分布的 r x c 列联表

Generate r x c contingency tables with bi-variate normal distribution in R

r

categorical-data

contingency