Displaying 1 result from an estimated 1 matches for "fev20s".
Did you mean:
fev20ns
2012 Jul 09
1
R to winbugs interface
...[1]
FEVdataBUGS=list("n","r","beta0","C0inv","a","b","FEV","Age","Smoke")
#creates list of imputs for winbugs code
parameters=c("beta","tau","meanFEVs","meanFEVns","FEV20s","FEV20ns")
#parameters to be monitored by winbugs
inits=list(list(tau=1,beta=c(0,0,0,0),FEV20s=2,FEV20ns=2)) # giving initial
values
FEV.fit=bugs(FEVdataBUGS,inits,parameters,"FEVWBModel.txt",n.chains=1,n.iter=60000,n.thin=1,n.burnin=10000)
print(FEV.fit,digits=3)
attach.b...