r - 错误:在“模拟”中:缺少参数“参数”,没有默认值
问题描述
嘿,我最近正在学习包“盛况”。而且我不知道为什么我的代码会出现这样的错误,来源来自论文https://kingaa.github.io/pomp/vignettes/pompjss.pdf
library(pomp)
# state function
gompertz.proc.sim <- function(x, t, params, delta.t, ...) {
eps <- exp(rnorm(n = 1, mean = 0, sd = params["sigma"]))
S <- exp(-params["r"] * delta.t)
setNames(params["K"]^(1 - S) * x["X"]^S * eps, "X")
}
# obs
gompertz.meas.sim <- function(x, t, params, ...){
setNames(rlnorm(n = 1, meanlog = log(x["X"]), sdlog = params["tau"]), "Y")
}
# density
gompertz.meas.dens <- function(y, x, t, params, log, ...) {
dlnorm(x=y["Y"], meanlog = log(x["X"]), sdlog = params["tau"], log = log )
}
# build our 'pomp' object
gompertz <- pomp(data = data.frame(time = 1:100, Y = NA), times = "time",
rprocess = discrete_time(step.fun = gompertz.proc.sim, delta.t = 1),
rmeasure = gompertz.meas.sim, t0 = 0
)
# assign values to params
theta <- c(r = 0.1, K = 1, sigma = 0.1, tau = 0.1, X.0 = 1)
# simulate values for our obs
gompertz <- simulate(gompertz, params = theta )
```
解决方案
您的代码将适用于旧版本的 pomp。但是,pomp 目前有一个更新的版本。请检查此链接https://kingaa.github.io/pomp/vignettes/upgrade_guide.html
请注意,JSS 中的论文是针对旧盛况的。
推荐阅读
- python - cx_oracle + 线程池先插入调用函数慢
- maven - 带有 google_checks 和 4 个空格 indentSize 的 Maven Checkstyle 插件
- c# - 我应该如何填充我的 201 Created 响应?
- bash - 如何提取简单文本并存储到文件中?
- javascript - Typescript,通用可变参数工厂函数返回元组
- python - 在 Pandas 中加入两列字典
- python-3.x - 将鼠标悬停在 Tkinter 中的按钮上时如何停止颜色变化?
- css - 如何跟踪和查找开发工具中所做更改的源代码?
- react-bootstrap - React.Children.only 期望接收单个 React 元素子元素。但文件中只有独生子女?
- c++ - 如何将 QDockWidgets 添加为选项卡?