簡體   English   中英

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=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.

 
粵ICP備18138465號  © 2020-2024 STACKOOM.COM