R中绘制样方计数强度图:版本更新后结果异常排查
问题:spatstat包intensity()函数生成的样方计数强度图在R新版本下异常
我使用某城市犯罪点数据集,通过spatstat包的quadratcount()函数绘制样方计数图:
旧环境(结果验证正常)
在R v4.1.2及对应版本spatstat环境中,执行以下代码生成的样方点数量分布图和样方计数强度图,经手动计算验证结果均正确:
# 绘制样方计数图 q <- quadratcount(xhomicide, 4, 8) plot(q)
# 绘制样方计数强度图 plot(intensity(q, image=TRUE), main=NULL, las=1)
注:xhomicide为预先定义的点模式对象(原描述标注为数据框,应为笔误)。
新环境(强度图结果异常)
将环境升级至R v4.3.0并安装最新版spatstat后,执行相同代码出现以下情况:
quadratcount()输出的样方计数结果与旧版本完全一致;intensity(q, image=TRUE)生成的样方计数强度图结果异常。
我需要排查该问题的原因:是R版本升级导致、spatstat包版本更新导致,还是代码写法不符合新版本要求?目前未找到相关问题说明。
内容的提问来源于stack exchange,提问作者Fritz
相关产品推荐
相关产品推荐

