# [R-sig-Geo] plotting Ripley fonctions

Alexandre VILLERS alexandre.villers at cebc.cnrs.fr
Thu Oct 9 15:30:30 CEST 2008

```Hey,

You need to transform the K function as follows

Ks<-Kest(object.ppp, border,r=seq(0,50000, by=1000))
plot(Ks\$r, sqrt(Ks\$border/pi)-Ks\$r) if you choose border correction.

Hope that's stuff will make it
Best regards

Alex

P.S.: for a confidence enveloppe, you can use
aleatoire<-runifpoint(length(s\$X), win=aigle.wind, giveup=1000)
Kalea<-Kest(aleatoire,r=seq(0,50000, by=1000))\$border
for (j in 1:99){
aleatoire<-runifpoint(length(s\$X), win=aigle.wind,
giveup=1000)
Kalea<-cbind(Kalea,Kest(aleatoire,r=seq(0,50000,
by=1000))\$border)
}
yKsaleamin<-apply(sqrt(Kalea/pi),1,min, na.rm=T)-Ks\$r
yKsaleamax<-apply(sqrt(Kalea/pi),1,max, na.rm=T)-Ks\$r
lines(Ks\$r,yKsaleamin,lty=2, col="darkgrey", lwd=1.2)
lines(Ks\$r,yKsaleamax,lty=2, col="darkgrey",lwd=1.2)

Lassana Touré a écrit :
> Hi.
> I need some help for plotting L(r) and Lcross so that I can have my
> confidence interval(CI) parallel to X axis.
> The plot I have obtain so far the graphic is not horizontal, then it
> is difficult to interpret strong the aggregation is or the attraction.
> thanks
>
> _______________________________________________
> R-sig-Geo mailing list
> R-sig-Geo at stat.math.ethz.ch
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>
>
> __________ Information from ESET Mail Security, version of virus signature database 3506 (20081009) __________
>
> The message was checked by ESET Mail Security.
> http://www.eset.com
>
>
>     part000.txt - is OK
>
>
>
>

--
Alexandre Villers
PhD Student
Team "Biodiversity"
CEBC-CNRS UPR1934
79360 Beauvoir sur Niort

Phone +33 (0)5 49 09 96 13
Fax   +33 (0)5 49 09 65 26

__________ Information from ESET Mail Security, version of virus signature database 3506 (20081009) __________

The message was checked by ESET Mail Security.
http://www.eset.com

part000.txt - is OK

```