作者cecepipi (pepe)
看板Statistics
標題[問題] R資料擷取問題Q^Q....
時間Sun Mar 24 04:22:06 2013
如果是跟統計軟體有關請重發文章
如果跟論文有關也煩請您重發文章
文章類別是為了幫助大家搜尋資料與解答，造成不便之處請見諒
小弟大約接觸R一～三個月！目前欲分析一組資料的時候遇到一個奇怪的問題＠＠..
我有一大組（五萬多筆）資料想要做t test之後去找 p-value < 0.001的數
我首先用了一個loop
aa <- rep(0,54674)
I <- 1:54674
for(i in I)
{
aa[i] <- t.test(data[ ,i+5]~data$SCQ_grp,alternative="two.sided")$p.value
}
去得到我想要的每一組資料p value
然後在寫一個
bb <- rep(0.54674)
for(i in I)
{
bb[i] <- t.test(data[ ,i+5]~data$SCQ_grp,alternative="two.sided")$statistic
}
去得到我的t值
我的問題是＠＠當我得到這個數據想做篩選的時候
cc <- which(aa<0.001) (我想找出所有比0.001小的數值）
我得到的結果是我aa的組別號碼,不是p-value
我該怎麼得到組別號碼＋p-value呢qq？還有如果有辦法能把我的t值一併合起來嘛qq？
我目前想到的方法只有aa[i],bb[i]去一個一個找＠＠可是這樣我要找4000多次qq
想問版上各位大大有沒有一個方法可以快速得到呢？
謝謝您耐心的看完～如有大大知道該如何處理的話
煩請您不吝提點小弟一些方法～感激不盡！！！
--
▲ ∕ ￣￣＼
＝◥ /
▅▅▅ ﹨
查理~ 給我腎臟~
皿 ● i \﹨
◢
▼ i \\﹨
咦 不可以嗎?
ξ
CC
|￣￣i ζ
ξ| | | |
◥ ψropkcat
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 208.102.218.186
1F：推 andrew43:b <- t.test(...) 則看看 b$statistic 和 b$p.value 就有 03/24 05:42
2F：→ andrew43:你可以學一下 ?t.test 裡的 "Value" 一節是什麼. 03/24 05:43
3F：推 blackalert:cc是你的index，利用它放回去找就好 => aa[cc], bb[cc] 03/24 08:44
4F：→ blackalert:要合併可以用cbind 03/24 08:45