标记变异函数计算距离规范

Question

在spatstat包中，为什么markvario函数不允许使用大于窗口长度1/4的距离来进行方差计算？

(X，更正=c(“各向同性”，“里普利”，“翻译”)，r=空，方法=“密度”，.，normalise=FALSE)

参数r是一个数字向量。应该评估标记变异函数γ(R)的参数r的值。

云杉数据集的窗口长度为200米，但方差图显示的距离仅为50米，甚至指定r到0,200。

r=seq(0,200，by=0.5)) 云杉标记变异函数

Answer 1

函数markvario返回类fv (函数值)的对象。它实际上是一个具有附加属性的data.frame。一个属性告诉plot命令(分派给plot.fv)推荐的绘图范围。如果您查看打印输出的底部，您将看到r运行在0到200之间，但建议的范围是0到50。使用xlim参数覆盖默认值：

library(spatstat)
rslt <- markvario(longleaf,r=seq(0,200, by=0.5))
rslt
#> Function value object (class 'fv')
#> for the function r -> gamma(r)
#> ...........................................................................
#>       Math.label               
#> r     r                        
#> theo  {gamma[]^{iid}}(r)       
#> trans {hat(gamma)[]^{trans}}(r)
#> iso   {hat(gamma)[]^{iso}}(r)  
#>       Description                                       
#> r     distance argument r                               
#> theo  theoretical value (independent marks) for gamma(r)
#> trans translation-corrected estimate of gamma(r)        
#> iso   Ripley isotropic correction estimate of gamma(r)  
#> ...........................................................................
#> Default plot formula:  .~r
#> where "." stands for 'iso', 'trans', 'theo'
#> Recommended range of argument r: [0, 50]
#> Available range of argument r: [0, 200]
#> Unit of length: 1 metre
plot(rslt, xlim = c(0, 200))

​

​