# 点過程の計算例
# 
# インストールされていなければ、以下のパッケージをインストールする必要がある。
# install.packages("spatstat")					#Spatial Point Pattern analysis, model-fitting and simulation
# install.packages("sm")							#kernel smoothing methods
# install.packages("akima")							#Interpolation of irregularly spaced data
# install.packages("gam")							#Generalized Additive Models
# インストールされているパッケージの一覧は次の命令で得ることができる。
# library()

library(spatstat)	# ライブラリの呼び出し

# マークなし点過程
x <- runif(50)	# 一様分布の乱数を50点発生
y <- runif(50)	# 一様分布の乱数を50点発生
pp <- ppp(x, y, c(0,1), c(0,1))	# 点過程オブジェクトの作成、範囲はxyともに0から1まで
plot(pp)	# 散布図の作成
k <- Kest(pp)	# K関数の計算
pc <- pcf(pp)	# pair correlationの計算
plot(k)	# K関数のグラフ
plot(k, sqrt(iso/pi)-r ~ r)	# L関数のグラフ
plot(pc)	# pair correlationのグラフ

# マーク付き点過程
pp <- ppp(x, y, c(0,1), c(0,1), marks=rep(factor(c("a","b")),25))	# 点過程オブジェクトの作成、aとbが交互にある。
plot(pp)	# aとbを含めた散布図の作成
plot(split(pp))	# aとbを分けた散布図の作成
k <- Kmulti(pp, pp$marks=="a", pp$marks=="b")	# K関数の計算
pc <- pcf(pp)	# pair correlationの計算
plot(k)	# K関数のグラフ(aとbの関係)
plot(k, sqrt(iso/pi)-r ~ r)	# L関数のグラフ(aとbの関係)
plot(pc)	# pair correlationのグラフ(aとbの関係)