R Gibbs区域相互作用模型中的误差

我正在使用spatstat v.1.59-0构建一些点过程模型，但我在一些验证工具上遇到了问题，特别是

effectfun

和剩余图

parres

该模型的组成部分是由52个考古遗址的点位置组成的

rings\u pp

，以及转换为

spatstat

像素图像的DEM

elev

。我试图评估吉布斯点过程模型与区域交互参数（

AreaInter

）和

elev

我使用以下代码将模型与数据相匹配：

---

rr1再次感谢您提出的明确问题

这是一个bug，可以通过简单的示例演示

fit <- ppm(cells ~ x, AreaInter(0.07))
plot(effectfun(fit, se.fit=TRUE))
plot(parres(fit))

fit再次感谢您提出的明确问题

这是一个bug，可以通过简单的示例演示

fit <- ppm(cells ~ x, AreaInter(0.07))
plot(effectfun(fit, se.fit=TRUE))
plot(parres(fit))

fit感谢您提出了一个包含许多细节的好问题。您是否可以提供数据，或者尝试使用模拟数据集或使用内置数据集重新创建问题？调试可以在本地重现的问题要容易得多。你好，鲁巴克博士。谢谢你的回复。我无法用其中一个示例数据集（bei）重现问题，我认为不可能将数据上载到stackoverflow。DEM协变量相当大（380MB）。但是，如果您同意的话，我很乐意将数据和R脚本通过电子邮件发送给您。如果可能的话，请在Dropbox上发布一个公共链接，供任何人下载和测试该示例。在接下来的一周或更长的时间里，我将没有时间运行代码并调试它，但我可能会在某个时候开始。我无法接收380 MB的电子邮件。我想我的限制是10MB。谢谢你提出一个很好的问题，有很多细节。您是否可以提供数据，或者尝试使用模拟数据集或使用内置数据集重新创建问题？调试可以在本地重现的问题要容易得多。你好，鲁巴克博士。谢谢你的回复。我无法用其中一个示例数据集（bei）重现问题，我认为不可能将数据上载到stackoverflow。DEM协变量相当大（380MB）。但是，如果您同意的话，我很乐意将数据和R脚本通过电子邮件发送给您。如果可能的话，请在Dropbox上发布一个公共链接，供任何人下载和测试该示例。在接下来的一周或更长的时间里，我将没有时间运行代码并调试它，但我可能会在某个时候开始。我无法接收380 MB的电子邮件。我想我的限制是10MB。嗨，巴德利博士。好的，非常感谢！你和鲁巴克和特纳的R书太棒了。有第二版的计划吗？嗨，巴德利博士。好的，非常感谢！你和鲁巴克和特纳的R书太棒了。有第二版的计划吗？
fit <- ppm(cells ~ x, AreaInter(0.07))
plot(effectfun(fit, se.fit=TRUE))
plot(parres(fit))