全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 R语言论坛
2215 2
2010-05-12
程序代码如下所示,为什么赋初值时出错,显示this chain contains uninitialized variables
model {
   for (i in 1:n){
     Y[i]~dnorm(s[i],tauC)
     s[i]<-theta[1]/(1+theta[2]*exp(-theta[3]*t[i]))}
    theta[1]~dnorm(mu[1],tau[1])
    theta[2]~dnorm(mu[2],tau[2])
    theta[3]~dnorm(mu[3],tau[3])
    tauC~dgamma(0.001,0.001)
    sigmaC<-1/sqrt(tauC)
  for (k in 1:3) {
    mu[k]~dnorm(0,0.000001)
    tau[k]~dgamma(0.001,0.001)}
}
list(t=7)
t[] Y[]
1 2
2 5
3 10
4 20
5 30
6 40
7 44
END
INITS
list(theta1=1,theta2=1,theta3=1,tauc=1)
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

全部回复
2010-5-12 19:57:05
model {
   for (i in 1:n){
     Y[i]~dnorm(s[i],tauC)
     s[i]<-theta[1]/(1+theta[2]*exp(-theta[3]*t[i]))}
    theta[1]~dnorm(mu[1],tau[1])
    theta[2]~dnorm(mu[2],tau[2])
    theta[3]~dnorm(mu[3],tau[3])
    tauC~dgamma(0.01,0.01)
    sigmaC<-1/sqrt(tauC)
  for (k in 1:3) {
    mu[k]~dnorm(0,0.000001)
    tau[k]~dgamma(0.01,0.01)}
}
#Data
list(n=7,t=c(1,2,3,4,5,6,7),Y=c(2,5,10,20,30,40,44))

#Init
list(theta=c(1,1,1),mu=c(1,1,1),tau=c(1,1,1),tauC=1)
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2010-5-12 20:16:36
非常感谢你!
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群