作者d2389758 (花!)
看板Statistics
标题[程式] R回圈跟结果
时间Mon May 14 19:23:49 2012
[软体程式类别]:
R
[程式问题]:
应该是资料处理
[软体熟悉度]:
新手(不到1个月)
[问题叙述]:
我有data(phenotype)要分别对所有genotype做回归
跑出来的结果我想要deviance跟0.025quant/0.5quant/0.975quant
可是我不知道怎麽把他们写进loop
查了很久也找不到QQ 可能是没用到关键字或是问题太简单了
单一次跑是这样
dataf1<-data.frame(cbind(ID,y,genodat[,1]))
names(dataf1) <- c("ID","y","SNP1")
model1 = y ~ 1 + SNP1 + f(ID, model="generic0", Cmatrix=AInv)
fit1 = inla(model1, data=dataf1,
control.data = list(hyper = param.data),
control.compute=list(dic=T))
summary(fit1)
假设我要跑25个 然後把跑出来的结果(deviance跟0.025quant/0.5quant/0.975quant)
用成table或matrix都可
所以我想loop应该是
for(i in 1:25){
datafi<-data.frame(cbind(ID,y,genodat[,i]))
names(datafi) <- c("ID","y","SNP")
modeli = y ~ 1 + SNP + f(ID, model="generic0", Cmatrix=AInv)
fiti = inla(modeli, data=datafi,
control.data = list(hyper = param.data),
control.compute=list(dic=T))
}
可是我不知道怎麽把那些值叫出来(写进loop)
麻烦大家帮我解答了
谢谢
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 193.170.62.226
1F：→ obarisk:做空的list放进去！ 05/14 21:00
2F：→ gsuper:新手应该不会用list吧 05/14 21:17
3F：→ d2389758:a没错...感谢楼上...刚刚在爬文研究怎麽弄list 05/14 21:21
4F：→ gsuper:就学学呗~~没多难 05/14 21:31
5F：→ moonbat:不好意思，想请问一下原PO"inla"这个是什麽? 05/15 16:08
6F：→ d2389758:inla是我用的package! 05/15 17:25
7F：→ gsuper:inlainlainlainla 05/15 18:16