如何在spatstat中使用Lest()函数绘制Ripley‘s L
如何在spatstat中使用Lest()函数绘制Ripley‘s L
提问于 2019-06-12 14:13:28
我有一个包含多个点的csv文件。数据包括经度、纬度和日期。我根据季节(春、夏、秋、冬)对我的数据进行了分类。我需要计算每个赛季的雷普利L。我使用了R中的spatstat包中的Lest()。我可以使用这个函数绘制Ripley‘s L。但是每个小区(季节)的x轴(距离)的范围是不同的。如何在同一范围(距离)内制作所有四个地块(季节)？
sa <- readOGR(".", "boundary") # Don't add the .shp extension
W <- as(sa, "owin")
Summer.ppp <- as(Sum_UTM, "ppp")
Autumn.ppp <- as(Autumn_UTM, "ppp")
Winter.ppp <- as(Winter_UTM, "ppp")
Spring.ppp <- as(Spring_UTM, "ppp")
Summer <- Lest(Summer.ppp, correction="Ripley")
plot(Summer , xlab="distance (m)", ylab="K(d)", cex.lab= 2, cex.axis=2)
plot(Summer, . -r ~ r, ylab=expression(hat("L")), xlab="distance (m)",cex.lab= 2, cex.axis=2)
Autumn <- Lest(Autumn.ppp, correction="Ripley")
plot(Autumn, xlab="distance (m)", ylab="K(d)", cex.lab= 2, cex.axis=2)
plot(Autumn, . -r ~ r, ylab=expression(hat("L")), xlab="distance (m)",cex.lab= 2, cex.axis=2)
Winter <- Lest(Winter.ppp, correction="Ripley")
plot(Winter, xlab="distance (m)", ylab="K(d)", cex.lab= 2, cex.axis=2)
plot(Winter, . -r ~ r, ylab=expression(hat("L")), xlab="distance (m)",cex.lab= 2, cex.axis=2)
Spring <- Lest(Spring.ppp, correction="Ripley")
plot(Spring , xlab="distance (m)", ylab="K(d)", cex.lab= 2, cex.axis=2)
plot(Spring, . -r ~ r, ylab=expression(hat("L")), xlab="distance (m)",cex.lab= 2, cex.axis=2)回答 1
Stack Overflow用户
发布于 2019-06-13 04:36:21
你写你有纬度和经度，但同时你的数据集被命名为"UTM"，这将表明你有平面x，y坐标。到底是哪一个？在spatstat中，您只能使用平面坐标。您不需要声明您正在使用哪些包。上面转换为ppp的as函数可能来自maptools？在这种情况下，我认为ppp的窗口只是设置为坐标的边界框。您有一个名为W的对象，您可能希望将其用作ppp的窗口。你可以这样做：
Window(Spring.ppp) <- W正如注释中所提到的，您可以简单地在plot调用中指定xlim = c(lower,upper)，其中lower和upper是您想要的x轴上的限制的数值。
页面原文内容由Stack Overflow提供。腾讯云小微IT领域专用引擎提供翻译支持
原文链接：
https://stackoverflow.com/questions/56555692
复制相关文章
相似问题