[程式] [R]計算配對後相關係數的對數

看板Statistics作者 (critical thinking )時間13年前 (2010/10/10 15:42), 編輯推噓1(1016)
留言17則, 4人參與, 最新討論串1/1
------------------------------------------------------------------------ [軟體程式類別]:R [程式問題]:假設Raw data在不同變項尚有缺漏值,今計算相關係數矩陣,如何得知 配對成幾對? [軟體熟悉度]:高(1年以上) [問題敘述]:如前所述 [程式範例]: mydata=matrix( c(1,2,NA,4,5,NA,7,8,9,10), 5,2) cor(mydata,use="pairwise") [,1] [,2] [1,] 1 NA [2,] 2 7 [3,] NA 8 [4,] 4 9 [5,] 5 10 由範例可知,應可配對成三對,但R的cor指令僅會輸出相關矩陣結果,並不會像SPSS 一樣可輸出N或SD,有無相關指令可以查詢? 謝謝R前輩協助解惑! ----------------------------------------------------------------------------- -- ※ 發信站: 批踢踢實業坊(ptt.cc) ◆ From: 220.136.216.217

10/10 17:19, , 1F
u can use the funtion named cor.test()
10/10 17:19, 1F

10/10 18:32, , 2F
感謝樓上回答,但研究一下似乎還不太知道如何得知!
10/10 18:32, 2F

10/10 18:55, , 3F
手動生產一個 mydata.2 , 把含有 NA 的 rows 消除不行嗎?
10/10 18:55, 3F

10/10 19:09, , 4F
回覆樓上先進,剛lin5有回覆可以使用na.omit作用類似你的
10/10 19:09, 4F

10/10 19:10, , 5F
的作法,不過因為我的資料是3000*12的矩陣,如此處理還是不
10/10 19:10, 5F

10/10 19:10, , 6F
方便,SPSS雖可以直接產生幾對相關係數但因係採R做模擬,
10/10 19:10, 6F

10/10 19:11, , 7F
10次,轉成SPSS處理不方便,所以才就教各位!
10/10 19:11, 7F

10/10 19:59, , 8F
這樣的話就抓df去算也可以
10/10 19:59, 8F

10/10 21:29, , 9F
好的 謝謝大家回覆喔!
10/10 21:29, 9F

10/10 22:50, , 10F
不懂你的邏輯ㄟ... 所以你的意思是 , 消除含有 NA 的 row
10/10 22:50, 10F

10/10 22:51, , 11F
是可行的辦法 , 但是你不會 , 是這樣嗎?
10/10 22:51, 11F

10/10 23:03, , 12F
mydata[!rowSums(is.na(mydata)),] 試試看這行(for winXP)
10/10 23:03, 12F

10/10 23:05, , 13F
3000*12 實在是一個超級小的矩陣~~~ (煙)
10/10 23:05, 13F

10/11 08:43, , 14F
回覆樓上先進,上述做法我大概會做,只是稍嫌麻煩!
10/11 08:43, 14F

10/11 08:44, , 15F
除了g大所提供的方法外,我也提供一個解決方式
10/11 08:44, 15F

10/11 08:44, , 16F
N=length(na.omit(missing001)[,1])
10/11 08:44, 16F

10/11 12:50, , 17F
na.omit() 還滿好用的 , good
10/11 12:50, 17F
文章代碼(AID): #1CiMvSfh (Statistics)