ecespa (version 1.1-1)

figuras: Artificial point data.

Description

The three different point patterns in the figure 3.1 of De la Cruz (2007)

Usage

data(fig1)
data(fig2)
data(fig3)

encoding

latin1

References

De la Cruz, M. 2007. M�todos para analizar datos puntuales. En: Introducci�n al An�lisis Espacial de Datos en Ecolog�a y Ciencias Ambientales: M�todos y Aplicaciones (eds. Maestre, F. T., Escudero, A. y Bonet, A.), pp 000-000. Asociaci�n Espa�ola de Ecolog�a Terrestre, Universidad Rey Juan Carlos y Caja de Ahorros del Mediterr�neo, Madrid.{Metodos para analizar datos puntuales. En: Introduccion al Analisis Espacial de Datos en Ecologia y Ciencias Ambientales: Metodos y Aplicaciones (eds. Maestre, F. T., Escudero, A. y Bonet, A.), pp 000-000. Asociacion Espa�ola de Ecologia Terrestre, Universidad Rey Juan Carlos y Caja de Ahorros del Mediterraneo, Madrid.}

Examples

data(fig1)

data(fig2)

data(fig3)

# transform to ppp format of spatstat with function haz.ppp:

fig1.ppp <- haz.ppp(fig1)

fig2.ppp <- haz.ppp(fig2)

fig3.ppp <- haz.ppp(fig3)

#Analyses as in Fig.3.2 of De la Cruz (2007). First, compute function K:

cosa1 <- Kest(fig1.ppp)

# Plot different estimators. 
# Fig. 3.2a: 

par("mar"=par("mar")+c(0,1,0,0))

plot(cosa1, col=c(1,0,0,1), lwd=c(2,2,2,2), lty=c(1,1,1,2), 
	 main="")

# Fig. 3.2b:

plot(cosa1, sqrt(./pi)-r~r, col=c(1,0,0,1), lwd=c(2,2,2,2), 
         lty=c(1,1,1,2), main="", ylab="L(r)")

# Fig. 3.2c:

plot(cosa1, .-(pi*r^2)~r, col=c(1,0,0,1), lwd=c(2,2,2,2), 
         lty=c(1,1,1,2), main="", ylab=expression(K(r)-pi*r^2))

# Fig. 3.2d:

plot(cosa1,(./(pi*r^2))-1~r, col=c(1,0,0,1), lwd=c(2,2,2,2), 
	 lty=c(1,1,1,2), main="",
         ylab=expression((K(r)/pi*r^2)-1))

## Analyses as in fig. 3.7 of De la Cruz (2007).
## First, compute function K and pointwise envelopes:

cosa1.env <- envelope(fig1.ppp, Kest)

cosa2.env <- envelope(fig2.ppp, Kest)

cosa3.env <- envelope(fig3.ppp, Kest)

## Plot function L with pointwise envelopes:

plot(cosa1.env,sqrt(./pi)-r~r, lwd=c(1,1,2,2), 
	 lty=c(1,1,3,3), col=c(1,1,1,1), xlab="r", 
         ylab="L(r)", main="", ylim=c(-2,2))

## Add simultaneous envelopes of Ripley (+-1.68 *sqrt(A)/N):

abline(h=1.68*sqrt(fig1.ppp$w$area)/fig1.ppp$n, 
         lty=2, lwd=2)

abline(h=-1.68*sqrt(fig1.ppp$w$area)/fig1.ppp$n, 
         lty=2, lwd=2)

## Plot function L with pointwise envelopes:

plot(cosa2.env,sqrt(./pi)-r~r, lwd=c(1,1,2,2), 
         lty=c(1,1,3,3), col=c(1,1,1,1), xlab="r", 
         ylab="L(r)", main="")

## Add simultaneous envelopes of Ripley:

abline(h=1.68*sqrt(fig2.ppp$w$area)/fig2.ppp$n, 
         lty=2, lwd=2)

abline(h=-1.68*sqrt(fig2.ppp$w$area)/fig2.ppp$n, 
         lty=2, lwd=2)

## Plot function L with pointwise envelopes:

plot(cosa3.env,sqrt(./pi)-r~r, lwd=c(1,1,2,2), 
         lty=c(1,1,3,3), col=c(1,1,1,1), xlab="r", 
         ylab="L(r)", main="")

## Add simultaneous envelopes of Ripley:

abline(h=1.68*sqrt(fig3.ppp$w$area)/fig3.ppp$n, 
         lty=2, lwd=2)

abline(h=-1.68*sqrt(fig3.ppp$w$area)/fig3.ppp$n, 
         lty=2, lwd=2)