R rmcorr 不正确的自由度数
R rmcorr incorrect degrees of freedom number
当我使用以下代码分析来自 mrcorr 包的 raz2005 数据的重复测量相关性时：
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 完全相同，我得到 df=400 以上，尽管我的参与者编号是 25。我使用：
mydata.rmc <- rmcorr (participant = participant, measure1 = measure1, measure2 = measure2, dataset = mydata)
根据论文（第 7 页），自由度应该这样计算：
N = n_distinct(raz2005$participant)
k = nrow(raz2005)/N
df = N*(k-1)-1
当我使用以下代码分析来自 mrcorr 包的 raz2005 数据的重复测量相关性时：
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 完全相同，我得到 df=400 以上，尽管我的参与者编号是 25。我使用：
mydata.rmc <- rmcorr (participant = participant, measure1 = measure1, measure2 = measure2, dataset = mydata)
根据论文（第 7 页），自由度应该这样计算：
N = n_distinct(raz2005$participant)
k = nrow(raz2005)/N
df = N*(k-1)-1