[問題] R資料擷取問題Q^Q....
如果是跟統計軟體有關請重發文章
如果跟論文有關也煩請您重發文章
文章類別是為了幫助大家搜尋資料與解答，造成不便之處請見諒
小弟大約接觸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
推
03/24 05:42, , 1F
03/24 05:42, 1F
→
03/24 05:43, , 2F
03/24 05:43, 2F
推
03/24 08:44, , 3F
03/24 08:44, 3F
→
03/24 08:45, , 4F
03/24 08:45, 4F