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

【问题标题】：Generate r x c contingency tables with bi-variate normal distribution in R在 R 中生成具有双变量正态分布的 r x c 列联表
【发布时间】：2019-12-08 21:05:03
【问题描述】：

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

我尝试从均值为 0 且 sigma 为 diag(2) 的双变量正态分布生成样本。在此之后我不知道如何构建列联表。

library(MASS)

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

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

【问题讨论】：

您能详细说明一下列联表应该是什么样子吗？我想您必须首先将数据放入存储桶中，因为（至少据我所知）列联表是为分类变量设计的。这就是你想要做的吗？
@jludewig ：我只想生成 R x C 列联表，但数据（即）单元格计数应该来自双变量正态分布。
你的表格列名和行名应该是什么（即你的单元格是什么）？
@jludewig ：表格中的单元格应该是正值，单元格的性质是计数及其变量的频率分布。例如，考虑来自此链接的 2 x 2 表：[link] (en.wikipedia.org/wiki/Contingency_table)，该表可以使用 r 中的 rmultinom 函数生成，但我想从二元正态分布生成这样的表。希望这会让你清楚。
使用多项式变量，您正在查看不同类别的计数（将您的数据视为 (A,A,B,A,A,B)），它将定义列联表中的单元格。对于正常数据，您正在查看连续值，因此您的变量可以采用一定范围内的任何值（例如（0.7482、-0.8353、-0.7859、0.9191、0.6021、-0.2223、0.5293、-0.6382））。所以不清楚如何创建列联表，您可能必须先对数据进行分类。你确定这是你想要的吗？

标签： r categorical-data contingency

【解决方案1】：

如果要从持续数据创建竞争表，您将首先要将其转换为分类数据，然后执行类别的计数。您必须定义对您的分析有意义的中断：

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])

【讨论】：

谢谢！你的想法是转变的一种方式。我会根据自己的需要改进它。