SVCJ程序 求WINBUGS程序

shatian 发表于 2012-12-5 18:48
老师，按你说的方法做了，但是我的这个程序winbugs好像运行不成功，我当时是用OPENBUGS产生的，请问R怎么 ...

那你就要安装package "R2OpenBUGS": Running OpenBUGS from R
语法与 R2WinBUGS大同小异.
你可参考Vignettes: R2OpenBUGS.pdf先自行修改,
解决不了再说.
  http://cran.r-project.org/web/packages/R2OpenBUGS/index.html

epoh 发表于 2012-12-5 19:05
那你就要安装package "R2OpenBUGS": Running OpenBUGS from R
语法与 R2WinBUGS大同小异.
你可参考Vign ...

log中显示DIC monitor not set,是不是我openbugs软件有问题？因为我要的就是DIC的值，在直接运用openbugs时，DIC就运行不了

shatian 发表于 2012-12-5 19:57
log中显示DIC monitor not set,是不是我openbugs软件有问题？因为我要的就是DIC的值，在直接运用openbugs ...

R2OpenBUGS不会有问题的,
等你跑完3000就会有数据了
但就是要时间.

epoh 发表于 2012-11-27 09:12
你下载错了,不是16楼的程序,
16楼的程序是楼主依自己模型自己写的
还在修改中

老師您好.我在24樓重新依您的建議修改我的一些BUGS中不足的部分..但仍是無法成功..是否是寫的太亂或是程度不好太淺..但不論如何再請您幫學生指正一番..感謝..

epoh 发表于 2012-11-16 10:42
1.在c碟,新建文件夹"Bugs",
  作为working.directory
  放进两文件 svcjl.bug,(your code)

老师，我参照你的WINBUGS的程序，写OPENBUGS 的R语言程序，但是总是不对，你能把这个程序再帮我改写成OPENBUGS的吗？

shatian 发表于 2012-12-28 19:47
老师，我参照你的WINBUGS的程序，写OPENBUGS 的R语言程序，但是总是不对，你能把这个程序再帮我改写成OPE ...

哪一个程序?

epoh 发表于 2012-12-28 20:35
哪一个程序?

library(R2WinBUGS)
#Data
n=5038
yy=read.table(file="c:/Bugs/shatian.txt")
r=yy[,1]
summary(r)
data=list("n","r")
#initial values 请自行补上,符号应跟你的不同,请自行更正  
inits = function() {list(mu=,k=,kt=,tauv=,rho=,muy=,lambda=,tauy=,rhoj=,muv=)}
  
parameters <- c("mu","theta","k","sigv","muy","rhoj","sigy","muv","rho","lambda")
#小样先看结果
SVCJ.sim <- bugs(data, inits, parameters.to.save=parameters,"svcj.bug",n.chains=1,
               n.thin=1,n.iter=5000,n.burnin=1500,debug=TRUE,DIC=TRUE,
               bugs.directory="D:/WinBUGS14/",working.directory = "c:/Bugs/")
attach.bugs(SVCJ.sim)
print(SVCJ.sim,digits=4)



这个程序，但是我需要用OPENBUGS，能帮我改写成R2OPENBUGS的吗？

shatian 发表于 2012-12-29 12:57
library(R2WinBUGS)
#Data
n=5038

library(R2OpenBUGS)
#Data
n=5038
yy=read.table(file="c:/Bugs/shatian.txt")
r=yy[,1]
summary(r)
data=list("n","r")
#initial values 请自行补上,符号应跟你的不同,请自行更正  
inits = function() {list(mu=,k=,kt=,tauv=,rho=,muy=,lambda=,tauy=,rhoj=,muv=)}
  
parameters <- c("mu","theta","k","sigv","muy","rhoj","sigy","muv","rho","lambda")
#小样先看结果
SVCJ.sim <- bugs(data, inits, parameters.to.save=parameters,"svcj.txt",n.chains=1,
               n.thin=1,n.iter=5000,n.burnin=1500,debug=TRUE,DIC=TRUE,bugs.seed=9,
               working.directory = "c:/Bugs/")
attach.bugs(SVCJ.sim)
print(SVCJ.sim,digits=4)

epoh 发表于 2012-12-29 13:54
library(R2OpenBUGS)
#Data
n=5038

老师，程序显示错误如下： 参数((bugs.directory = "D:/Program Files/OpenBUGS/OpenBUG322")) 没有用
，这是什么原因呢，我之间自己参考着改的时候，也显示这个问题

shatian 发表于 2012-12-29 14:03
老师，程序显示错误如下： 参数((bugs.directory = "D:/Program Files/OpenBUGS/OpenBUG322")) 没有用
， ...

忘了,不需要bugs.directory
我更正删除了