[英]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=more than 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
聲明:本站的技術帖子網頁,遵循CC BY-SA 4.0協議,如果您需要轉載,請注明本站網址或者原文地址。任何問題請咨詢:yoyou2525@163.com.