请问: 在R里调用WinBUGS软件,通过函数bugs,例如:
bugs(data,inits,model.file="univt.bug",
parameters=c("mu","tau","df","lambda")
其中的参数model.file是要先在WinBUGS软件下写关于模型的代码:
model{
for(i in 1:N){
y ~ dt(mu,tau,df)
}
mu ~ dnorm(0.0,1.0E-6)
tau ~ dgamma(0.1,0.01)
df ~ dunif(2,50)
lambda <- sqrt(1/tau)
}
然后存储为univt.bug文件吗?