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@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
_______________________________________________
R-sig-Geo mailing list
R-sig-Geo@stat.math.ethz.ch
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
