在R中计算Spearman系数时，如何考虑聚类？

Question

我有一组人用血液和头发来测量他们的药物浓度(也就是说，每个人都有三个用血样测量的值，还有三个用头发样本测量的值)。我想计算两个测量值之间的Spearman系数，但我不知道如何解释个人内部的重复测量。在R中有办法做到这一点吗？

id<-rep(c(1:100),times=3) ##id variable
df1<-data.frame(id)
df1$var1 <- sample(500:1000, length(df1$id))  ##measurement1
df1$var2 <- sample(500:1000, length(df1$id))  ##measurement2
cor.test(x=df1$var1, y=df1$var2, method = 'spearman') ## this doesn't account for clustering within individuals

谢谢!

Answer 1

也许R包'rmcorr‘提供了您正在寻找的功能。该包帮助计算重复的度量相关性：

install.packages("rmcorr")

rmcorr::rmcorr(participant = id, measure1 = var1, measure2 = var2, dataset = df1)