首页 > 解决方案 > 有没有办法在创建线段的随机模式时指定线段交叉的数量?

问题描述

我一直在使用该spatstat程序包来确定点模式是聚类的、随机的还是规则的,方法是将其与在完全空间随机性 (CSR) 下生成的最近邻距离的相对频率分布进行比较。代码如下(Ypppx,y坐标的对象):

Y<-rpoint(60, 2, fmax=NULL, win=unit.square(),giveup=1000, verbose=FALSE,
       nsim=1, drop=TRUE)
envelope(Y, Gest, nsim = 999, nrank = 25, global=FALSE, fix.n=TRUE)

但是,我刚刚意识到随机点必须分布在窗口内多边形网络的交点处;它们不能与完全的空间随机性进行比较,因为它们仅限于多边形网络段的交叉点。有什么方法可以模拟网络模式的空间随机性?

我设法创建随机线段并在每个交叉点添加点:

l<-rpoisline(4, win=owin())
i<-selfcrossing.psp(l)
par(mfrow=c(1,2))
plot(l)
plot(i)

但是,我不确定如何将其与envelope()功能集成。我还需要交叉次数保持不变。那么无论如何要指定rpoisline()函数中的交叉点数吗?

标签: rpoissonspatstat

解决方案


推荐阅读