我正在使用 R 来分析具有复杂设计的调查数据。因此,我正在使用
survey
包。有没有办法在 R 中获得多向相关矩阵,同时考虑到我的数据的复杂样本设计(即我需要使用采样权重)?
svycor()
包中的jtools
函数仅允许获取皮尔逊相关矩阵。 weightedCorr()
包中的 wCorr
函数仅允许获取由单个权重和变量对(而不是一次整个矩阵)加权的相关性。
提前谢谢您。
您可以使用
survey::svytable
估计人口列联表,然后将其输入到 polycor::polychor
。
如果您需要标准误差,我目前能想到的唯一方法是使用复制权重。也就是说,使用
survey::withReplicates
对每个重复重复上述过程。