尝试通过 simulate.mppm 模拟点模式会导致错误

Question

我正在尝试基于适合多种多类型点模式数据集的模型来模拟点模式。但是，我不能运行 simulate.mppm 而不引发错误。我附上了下面的代码（使用 optimbase_1.0-9 和 spatstat_1.62-2）：

library(spatstat)
library(optimbase)
flusubset <- flu[1:4]
typelist <- lapply(lapply(flusubset$pattern, marks), levels)
num_marks <- length(typelist[[1]])
iradii <- 50*ones(num_marks,num_marks)
hradii <- 3*ones(num_marks,num_marks)
Int <- anylist()
for (i in 1:dim(flusubset)[1]) {
  Int[[i]] <- MultiStraussHard(iradii=iradii, hradii=hradii)
}
Int <- as.hyperframe(Int)
multmodel <- mppm(pattern ~ 1, data=flusubset, interaction=Int)
simulate(multmodel)

每当我运行这段代码时，我都会遇到错误：

Error in check.nvector(w, nrow(x), things = "rows of x") : The length of ‘w’ (=16744) should equal the number of rows of x (=15424)
9. stop(whinge)
8. check.nvector(w, nrow(x), things = "rows of x")
7. sumouter(mom, lam * wt * glmsub)
6. vcmGibbs(object, ..., what = what, err = err)
5. vcov.mppm(object, what = "fisher", err = "null")
4. vcov(object, what = "fisher", err = "null")
3. subfits(object)
2. simulate.mppm(multmodel)
1. simulate(multmodel)

感谢您的帮助！

Answer 1

这是 spatstat 当前版本中的错误。

我已在 GitHub repository 的开发版本 spatstat 1.62-2.011 中修复了它。包含错误修复的 spatstat 的下一个 public 版本将于 2020 年 2 月发布。

对于 spatstat 中的错误，请 post 在 GitHub 存储库的问题页面上提出问题。

尝试通过 simulate.mppm 模拟点模式会导致错误

Trying to simulate point patterns via simulate.mppm causes error

r

spatstat