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