为什么跳跃那么多？！----mcmc 方法估计跳跃扩散模型的问题-经管之家官网！

坛友互助群

扫码加入各岗位、行业、专业交流群

MCMC to estimate a jump-diffusion model

We use Winbugs to do MCMC through R, the code is below:

# winbugs program#
modle {
for (i in 1:N) {
y ~dnorm(r,tau) # tau corresponding to
r<-mu+Xi*J # Xi corresponding to #
Xi[t]~dnorm(Jmu,Jtau)
J[t]~dbern(lamda)
}
Sigma<-1/sqrt(tau)
JSigma<-1/sqrt(Jtau)
mu~dnorm(0,1) # prior distributions for parameters
tau~dgamma(2, 2)
Jmu~dnorm(0,1)
Jtau~dgamma(2,2)
lamda~dbeta(2,2)

}

# R code with interface to Winbugs

library(R2WinBUGS)
y=read.table("c:/data.txt") #import data#
N=length(N)
Xi0=rep(0,N) # Give initial values to state variables#
J0=rep(0,N)
data=list("N","y")
# Give initial values to the parameters for winbugs #
inits = function() {list ( tau=2,
Jtau=2,
mu=0,
Jmu=0,
lamda = 0.05, Xi=Xi0,J=J0)}
parameters<- c("Sigma","JSigma","mu","Jmu","lamda","Xi","J")
# quote the winbugs#
jump.sim<-bugs(data,inits,parameters,"c:/mcmc/model.txt",n.chains=1,n.iter=1000)

无论怎么调整初始值和先验参数，跳跃的强度 lamda 始终在0.2以上，请问为什么 是不是我模型有问题？data 已经上传：）期待指教：）

扫码或添加微信号：坛友素质互助

「经管之家」APP：经管人学习、答疑、交友，就上经管之家！
免流量费下载资料----在经管之家app可以下载论坛上的所有资源，并且不额外收取下载高峰期的论坛币。
涵盖所有经管领域的优秀内容----覆盖经济、管理、金融投资、计量统计、数据分析、国贸、财会等专业的学习宝库，各类资料应有尽有。
来自五湖四海的经管达人----已经有上千万的经管人来到这里，你可以找到任何学科方向、有共同话题的朋友。
经管之家（原人大经济论坛），跨越高校的围墙，带你走进经管知识的新世界。
扫描下方二维码下载并注册APP