[程式] R迴圈跟結果
[軟體程式類別]:
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
→
05/14 21:00, , 1F
05/14 21:00, 1F
→
05/14 21:17, , 2F
05/14 21:17, 2F
→
05/14 21:21, , 3F
05/14 21:21, 3F
→
05/14 21:31, , 4F
05/14 21:31, 4F
→
05/15 16:08, , 5F
05/15 16:08, 5F
→
05/15 17:25, , 6F
05/15 17:25, 6F
→
05/15 18:16, , 7F
05/15 18:16, 7F