我有一个城市的犯罪数据。通过spatstat
包使用quadratcount()
,我绘制了这个数据集的样方计数。图1显示了每个样方的数据点计数的图像,并通过以下代码创建:
#Quadrat Count Map
q <- quadratcount(xhomicide, 4, 8)
plot(q)
其中xhomicide是先前定义的 Dataframe 。
图1.每个样方的点数图:
图2示出了图1中所示的样方计数的相应强度图。图2是使用以下代码创建的:#添加每个样方图的强度(intensity(q,image=TRUE),main=NULL,las=1)
图2.每个样方的点数的强度图。
手工计算表明,图1和图2都显示了样方计数方法的正确结果。
请注意,图1和图2是在R v4.1.2下创建的。
我将R更新到v4.3.0并安装了spatstat包的当前版本。当在新版本的R下运行上述代码时,我获得了与图1相同的结果(R版本之间每个样方的计数值相同)。但是,当我运行intensity()函数时,结果如图3所示。
图3.使用R v4.3.0绘制的每个样方点数的强度图:
既然两个版本的R的代码是相同的,那么问题是出在新版本的R上,还是spatstat发生了什么变化?我很想知道这个问题是否与新版本的R有关,当前版本的spatstat,或者我的代码不再满足R或spatstat的要求。我找不到任何其他的评论是什么我是造成这个问题。
1条答案
按热度按时间oaxa6hgo1#
这是一个bug。谢谢你提醒我们注意。
此问题已在the github repository上提供的
spatstat.geom
开发版本(版本3.2-1.001
)中得到修复