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