时间戳

首页 > 解决方案 > R rmcorr 不正确的自由度数

r - R rmcorr 不正确的自由度数

问题描述

当我分析包中的raz2005数据mrcorr以进行重复测量时，与以下代码相关：

brainvolage.rmc <- rmcorr (participant = Participant, measure1 = Age, measure2 = Volume, dataset = raz2005)

正如论文https://www.ncbi.nlm.nih.gov/pmc/articles/PMC5383908/pdf/fpsyg-08-00456.pdf中给出的那样，我得到 df=71，这是正确的，因为参与者的数量是72 在这个数据中。

然而，当我对我的数据使用相同的代码时，我的数据结构raz2005与

mydata.rmc <- rmcorr (participant = participant, measure1 = measure1, measure2 = measure2, dataset = mydata)

标签： rcorrelation

解决方案

根据论文（第 7 页），自由度应按如下方式计算：

N = n_distinct(raz2005$participant)
k = nrow(raz2005)/N
df = N*(k-1)-1

推荐阅读