Question

我正在尝试按照包的粉丝图中的示例，但我找不到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)

Answer 1

该文件位于R_LIBS/fanplot/model/，其中R_LIBS是您的软件包库的位置。

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

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

然后下一行将模型写入工作目录中的文件ew.txt。目前尚不清楚您想要哪个ew.txt。如果包中附带的那个如上所述在R_LIBS中找到它。或者从CRAN下载源tarball（package page on cran中的tar.gz。

Answer 2

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

无法在R中的包fanplot中找到运行示例的ew.txt文件

2 个答案: