R的调查包中的Pearson相关系数

很抱歉,如果这是非常明显的,但我看不出如何在调查包中的两个变量之间进行简单的Pearson关联.我的数据有层次,所以它相当于在apistrat中为api00和api99找到r.

library(survey)
data(api)

dstrat <- svydesign(id=~1,strata=~stype, weights=~pw, data=apistrat, fpc=~fpc)

我敢肯定必须有一个简单的方法,使用svyvar或svyglm或其他东西,但我看不到它？

您可以使用svyvar估计方差 - 协方差矩阵,然后将其缩放到相关性:

library(survey)
data(api)

dstrat <- svydesign(id=~1,strata=~stype, weights=~pw, data=apistrat, fpc=~fpc)
v <- svyvar(~api00+api99, dstrat)

as.matrix(v)
cov2cor(as.matrix(v))

这适用于任何数量的相关性和任何设计.