楼主: lipj
1800 0

WinBUGS编译时无响应求助! [推广有奖]

教授

已卖:1105份资源

学术权威

3%

还不是VIP/贵宾

-

威望
1
论坛币
22372 个
通用积分
2226.6850
学术水平
61 点
热心指数
104 点
信用等级
55 点
经验
13322 点
帖子
2262
精华
0
在线时间
6620 小时
注册时间
2008-10-5
最后登录
2026-1-5

楼主
lipj 在职认证  发表于 2014-8-1 22:20:13 |AI写论文
1论坛币
附件是EGARCH的winbugs模型,求帮助!!!
model
{
  # Measurement Equation
  ysigma20 <- 1/(logsigma20)
  
  y0 ~ dnorm(alpha, ysigma20)
  ymean[1] <- alpha + beta * y0
  y[1] ~ dnorm(ymean[1], ysigma2[1])
  for(t in 2:T){
  
    ymean[t] <- alpha + beta * y[t-1]
ysigma2[t] <- 1/(sigma2[t])
    y[t] ~ dnorm(ymean[t], ysigma2[t])
  }
  
  # State Equation
  logsigma20 <- mu
  
  logsigma2[1] <- mu/(1-phi)
  
  sigma[1] <- exp(logsigma2[1]/2)
  
  sigma2[1] <- pow(sigma[1], 2)
  
  ysigma2[1] <- 1/(sigma2[1])
  
  logsigma2[2] <- mu + phi * logsigma2[1] + psi1 * ((y[1]-alpha-beta * y0)/sigma[1]) + psi2 * abs((y[1]-alpha-beta * y0)/sigma[1])
  
  sigma[2] <- exp(logsigma2[2]/2)

  sigma2[2] <- pow(sigma[2], 2)
  for(t in 3:T){
  
    logsigma2[t] <- mu + phi * logsigma2[t-1] + psi1 * ((y[t-1]-alpha-beta * y[t-2])/sigma[t-1]) + psi2 * abs((y[t-1]-alpha-beta * y[t-2])/sigma[t-1])
  
    sigma[t] <- exp(logsigma2[t]/2)

sigma2[t] <- pow(sigma[t], 2)
  
  }
  # Piror Distribution
  alpha ~ dnorm(0, 1)
  beta ~ dnorm(0, 1)
  mu ~ dunif(0, 1)
  phi ~ dunif(0, 1)
  psi1 ~ dunif(0, 1)
  psi2 ~ dunif(0, 1)
}


关键词:winbugs WINBUG BUGS bug Win 模型
为了美好的未来生活而奋斗!

您需要登录后才可以回帖 登录 | 我要注册

本版微信群
加好友,备注jltj
拉您入交流群
GMT+8, 2026-1-6 04:14