[R-sig-Geo] plotting Ripley fonctions

[Alexandre VILLERS]

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