作者discount (critical thinking )
看板Statistics
標題[程式] [R]計算配對後相關係數的對數
時間Sun Oct 10 15:42:18 2010
------------------------------------------------------------------------
[軟體程式類別]: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
1F：推 Prozac:u can use the funtion named cor.test() 10/10 17:19
2F：→ discount:感謝樓上回答,但研究一下似乎還不太知道如何得知! 10/10 18:32
3F：→ gsuper:手動生產一個 mydata.2 , 把含有 NA 的 rows 消除不行嗎? 10/10 18:55
4F：→ discount:回覆樓上先進，剛lin5有回覆可以使用na.omit作用類似你的 10/10 19:09
5F：→ discount:的作法,不過因為我的資料是3000*12的矩陣,如此處理還是不 10/10 19:10
6F：→ discount:方便，SPSS雖可以直接產生幾對相關係數但因係採R做模擬, 10/10 19:10
7F：→ discount:10次，轉成SPSS處理不方便，所以才就教各位! 10/10 19:11
8F：→ lin15:這樣的話就抓df去算也可以 10/10 19:59
9F：→ discount:好的 謝謝大家回覆喔! 10/10 21:29
10F：→ gsuper:不懂你的邏輯ㄟ... 所以你的意思是 , 消除含有 NA 的 row 10/10 22:50
11F：→ gsuper:是可行的辦法 , 但是你不會 , 是這樣嗎? 10/10 22:51
12F：→ gsuper:mydata[!rowSums(is.na(mydata)),] 試試看這行(for winXP) 10/10 23:03
13F：→ gsuper:3000*12 實在是一個超級小的矩陣~~~ (煙) 10/10 23:05
14F：→ discount:回覆樓上先進,上述做法我大概會做,只是稍嫌麻煩! 10/11 08:43
15F：→ discount:除了g大所提供的方法外,我也提供一個解決方式 10/11 08:44
16F：→ discount:N=length(na.omit(missing001)[,1]) 10/11 08:44
17F：→ gsuper:na.omit() 還滿好用的 , good 10/11 12:50