林學菜鳥—R語言點格局分析

本人是R語言菜鳥一枚，有不足的地方希望大佬們指教。這段時間忙著寫畢業論文，因為論文的一部分內容涉及到樹種的分佈格局，所以自己鼓搗了一段時間。下面是我用R語言（Rstudio）做的一些工作（下列4圖分別是：原始資料；輸出結果；K分佈格局；G分佈格局）：

help("Kest")
# 幫助函式,對搞清楚包、函式的作用及使用方法賊有用
install（spatstat）
安裝spatstat包
require(spatstat)
#載入spatstat包
setwd("E:/cycl")
#建立工作環境


cycl<-read.csv("cycl.csv",header=T)
head(cycl)
#資料讀取
cycl.ppp<-ppp(cycl$X,cycl$Y,window=owin(xrange=c(0,24.5),yrange=c(0,17.2)))
#該函式建立點模式資料集
cycl.envelope<-envelope(cycl.ppp,fun=Kest,nsim=199)
cycl.envelope<-envelope(cycl.ppp,fun=Gest,nsim=199)
#Monte-Carlo檢驗(Envelope)；k(r)org(r);nism(迭代次數)

write.csv(cycl.envelope, "cyclenvelope.csv")
# 將envelope()函式計算的結果輸出到一個.csv檔案裡，可以根據輸出結果自己作圖


plot(cycl.envelope,xlab="尺度Scale (m)",ylab="K值K(r)",main="b",xlim=c(0,6))
# xalb:x軸標題；ylab:y軸標題；main:圖示題；xlim(調整x軸範圍，即尺度範圍)

plot(cycl.envelope,xlab="尺度Scale (m)",ylab="G值G(r)",main="b",xlim=c(0,15))
# xlim(調整x軸範圍，即尺度範圍)


plot(cycl.ppp)
# # 樹種空間分佈圖
plot(envelope(cycl.ppp,Kest))
# # K分佈格局
plot(envelope(cycl.ppp,Gest))
# # G分佈格局