前段时间,有个朋友,问我如何用R解决生态学的部分问题,主要是空间生态学的问题。这对我来说,已经跨行了,后来咨询一位生态学钟博士,他给了一点代码,在此基础上,我做了一些分析。关于spatstat包的中文帖子好像没有,于是,写点用法,供有需要的网友参考。
1 数据读取
library(spatstat)
setwd("C:/Users/yzhlin/Desktop/LS")
df
df1
2 径级结构ggplot(df1s,aes(x=Dbh))+
geom_histogram(binwidth=3,color='black',fill='white')
3 空间分布格局
3.1 空间分布点图sp
range(df1s$Xcoor.60),range(df1s$Ycoor.60))
spxy
ggplot(spxy,aes(x,y))+geom_point()
3.2 点格局分析
(1)Ripley’s K(r)函数k
plot(k)
(2)L函数L
plot(L)
(3)双变量的Ripley’s K(r)函数
df1sab
df1sab$S