我在 spatstat 中使用 ppm 函数并寻找关于正交方案中使用的默认虚拟点的文档。 default.dummy 帮助页面显示“如果 nd 缺失，默认值（取决于数据模式 X）由 default.ngrid 计算。”我现在正在寻找有关 default.ngrid 的信息，但找不到任何东西。对于上下文，我将空随机模型与具有协变量的模型进行比较，其中点是类 ppp 的坐标，协变量是窗口中的像素图像。为了生成模型，我使用带有默认参数的 ppm 函数，如下所示：
ppm0 <- ppm(myPoints_ppp ~ 1)
ppm1 <- ppm(myPoints_ppp ~ covariate_A)
我试图了解虚拟点是如何生成的，因为我没有指定 nd 参数。谢谢！
【问题讨论】：
标签： spatstat
该引用已过时。默认虚拟点数由default.n.tiling 确定。这是一个内部的、未记录的函数。
可以参考default.n.tiling的源代码了解具体规则，但这里是一个草图：
目前默认每个维度的最小网格点数是较大的
10 * ceiling(2 * sqrt(npoints(X))/10)
和
spatstat.options('ndummy.min')
这确定了可接受的最小虚拟点数。使用默认值spatstat.options('ndummy.min') = 32 这意味着任何不超过 225 个点的点模式将被赋予一个最小 32 x 32 的虚拟点网格。具有 226 到 400 个点的模式将具有 40 x 40 个虚拟点。具有 401 到 625 个点的模式将具有 50 x 50 个虚拟点，依此类推。
虚拟点的最终数量是通过应用一些其他约束（取决于上下文和用户指定的参数）确定的，并且可能大于上面指定的最小数量。
这个默认值总是可以通过简单地为参数nd指定一个非NULL值来推翻。
我还应该指出，选择默认设置是为了快速产生可接受的结果，而不是产生高度准确的结果。这是必要的，因为 CRAN 对检查包裹的总时间有限制。如果您正在做任何严肃的方法研究，您应该考虑增加spatstat.options('ndummy.min') 或创建自己的正交方案。这也将提高结果的可重复性（因为spatstat 中的默认值可能会更改。）
【讨论】：