我正在尝试按照fanplot包中的一个示例进行操作,但我找不到ew.txt文件。如果你能帮助我,我将不胜感激。
data(ew)
plot(ew)
## Not run:
plot(ew.mcmc)
##
##Create your own (longer) MCMC sample
##
# data
r <- diff(log(ew))
r <- ts(exp(r) - 1, start(ew))
library("R2OpenBUGS")
# write model file:
ew.bug <- dget(system.file("model", "ew.txt", package = "fanplot"))
write.model(ew.bug, "ew.txt")
# take a look:
file.show("ew.txt")
# run openbugs
ew.mcmc <- bugs(data = list(N = length(r), H = 25, r = c(r), p = c(ew)),
inits = list(list(a = exp(12), psi = 0.5, itau2 = 0.5)),
param = c("alpha", "psi", "tau", "r.new", "p.new", "y.sim"),
model = "ew.txt",
n.iter = 11000, n.burnin = 1000, n.chains = 1, n.thin = 1)
## End(Not run)发布于 2012-08-03 05:09:24
该文件位于R_LIBS/fanplot/model/中,其中R_LIBS是包库的位置。
在示例代码中,可以通过system.file()找到该文件
system.file("model", "ew.txt", package = "fanplot")然后,下一行将模型写入working目录中的文件ew.txt。目前还不清楚您想要哪个ew.txt。如果包附带的包可以在上面的R_LIBS中找到它。或者从CRAN下载源代码tarball (从package page on cran下载tar.gz。
发布于 2012-08-03 04:55:00
当您在交互式会话中运行此代码时,使用getwd查看您的工作目录。这应该是包含保存的文本文件的目录。
https://stackoverflow.com/questions/11785520
复制相似问题