我在R中使用spatstat在线性网络上生成点。我已成功在道路网络上生成泊松分布式泊松过程。
现在当我想使用'nsim'时,为了实现随机过程的多个实现,我得到一个错误(使用这个数据)http://biogeo.ucdavis.edu/data/diva/rds/NLD_rds.zip
roads<-readShapeSpatial("NL_rd.shp")
spatstat.roads<-as.psp(roads)
final_roads<-as.linnet(spatstat.roads)
abc<-rpoislpp(0.02/5000,final_roads, nsim=2)
plot(abc)
当我使用该nsim时,我收到如下错误
Error in (function (..., check = TRUE, promote = TRUE, demote = FALSE) :
Some arguments of solist() are not 2D spatial objects
虽然它适用于一个实现,但我想绘制两个实现来直观地比较随机性。所以我不明白一个实现是如何到达而不是两个和非二维空间对象的错误对话。这对我来说有点神秘。
请您澄清一下这个错误以及解决方法。
感谢。
答案 0 :(得分:1)
这是一个错误。
很快会在the development version of spatstat,版本1.47-0.018
及更高版本中修复。
与此同时,您可以使用replicate
获得所需的结果,例如
abc <- replicate(2, rpoislpp(0.02/5000,final_roads), simplify=FALSE)
感谢您引起我们的注意。