作者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