R中Chisquared离差检验的误差

Question

我对R非常陌生，正在尝试确定包含~500k点的形状文件是否是随机分布的。不管我输入什么参数，我都会得到相同的错误信息。在

以下是我从一开始就有的：

> library(spatstat)

spatstat 1.40-0       (nickname: ‘Do The Maths’) 
For an introduction to spatstat, type ‘beginner’
> as.ppp(area)
marked planar point pattern: 500000 points
Mark variables: 
[1] OBJECTID   Encoded_Ti Time_      Filter     Category   Severity      Action_    Hit_Count  Profile    Encoded_So Source_IP  Source_Por Encoded_De
[14] Dest_IP    Dest_Port  VLAN_Tag   Source_Cou Source_Reg Source_Cit     Source_Lat Source_Lon Dest_Count Dest_Regio Dest_City  Dest_Latit Dest_Longi
window: rectangle = [-159.964, 178.417] x [-46.4, 70.6349] units
Warning message:
some mark values are NA in the point pattern x 
> quadrat.test(area)
Error in UseMethod("quadrat.test") : 
no applicable method for 'quadrat.test' applied to an object of class          "c('SpatialPointsDataFrame', 'SpatialPoints', 'Spatial')"


<bytecode: 0x0000000024e7a660>
<environment: namespace:spatstat>

所以我的下一个尝试是：

^{pr2}$

我最后的尝试是：

> quadrat.test(X, nx = 20, ny = 20)
Error in rectquadrat.countEngine(X$x, X$y, tess$xgrid, tess$ygrid) : 
xbreaks do not span the actual range of x coordinates in data
> quadrat.test(X, nx = 20, ny = 20, xbreaks= NULL, ybreaks = NULL)
Error in rectquadrat.countEngine(X$x, X$y, tess$xgrid, tess$ygrid) : 
xbreaks do not span the actual range of x coordinates in data
> data(X)
Warning message:
In data(X) : data set ‘X’ not found
> quadrat.test(X)
Error in rectquadrat.countEngine(X$x, X$y, tess$xgrid, tess$ygrid) : 
xbreaks do not span the actual range of x coordinates in data
> quadrat.test(X, 10)
Error in rectquadrat.countEngine(X$x, X$y, tess$xgrid, tess$ygrid) : 
xbreaks do not span the actual range of x coordinates in data

就像我说的，我在这方面很新，而且只有业余的python经验，但是我是一个研究生，需要在项目中使用这个函数。任何帮助都将不胜感激。在

干杯

Answer 1

您应该知道R（几乎）永远不会改变您的输入变量，所以您需要为您的命令分配一个输出并从那里开始工作。特别是如果area是一个变量，它可以通过as.ppp转换为ppp对象，则应为结果命名：

X <- as.ppp(area)

然后您可以将quadratcount或quadrat.test函数应用于新创建的ppp对象：

^{pr2}$

对于~500k点，20×20的方格可能看起来是一个非常大的空间尺度，用于样方计数，但这当然取决于您的具体设置。在

由于您还没有提供一个可复制的示例，我不能说这些命令是否适用于您的设置，但它们确实可以处理适当的数据。在