1

我正在尝试遵循包 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)
4

2 回答 2

3

该文件位于R_LIBS/fanplot/model/您R_LIBS的包库的位置。

在示例代码中，该文件是通过以下方式找到的system.file()：

system.file("model", "ew.txt", package = "fanplot")

然后下一行将模型写入工作目录ew.txt中的文件。目前尚不清楚您想要哪个。如果随包装一起找到它，如上所示。或者从 CRAN 下载源代码压缩包（来自cran 上的包页面。ew.txtR_LIBStar.gz

于 2012-08-02T21:09:23.523 回答
0

当您在交互式会话中运行此代码时，使用getwd来查看您的工作目录是什么。这应该是包含已保存文本文件的目录。

于 2012-08-02T20:55:00.037 回答