Bonjour,
PASSAGE m'est inconnu.
Les modules de statistiques spatiales d'ADE-4 classique ne sont
(toujours) pas passés sous R. Il existe sous R le package spatstat plus
particulièrement spécialisé dans la simulation de processus ponctuels.
Il permet également des analyses exploratoires dans un esprit voisin de
celui d'ADS. Inconvénient : spatstat ne permet pas le calcul
d'intervales de confiance et ne propose la correction des effets de bord
par la méthode de Ripley que pour des domaines de forme rectangulaire.
Il y aurait donc un certain intérêt à envisager le passage des routines
ADS sous R. Toute bonne volonté sera la bienvenue ...
Pour ce qui est de la transformation de K(r) calculée par spatstat en
L(r), c'est très simple (en reprenant l'exemple BPoirier de la pile de
données d'ADE-4 classique) :
library(spatstat)
# lire le fichier des coordonnées des hêtres :
BP_Beech<-read.table("BP_Beech.txt")
# créer une fenêtre d'échantillonnage rectangulaire de 110 x 90 m :
BP_Rect<-owin(c(0,110),c(0,90))
# créer un objet de classe "ppp" associant le semis et la fenêtre
d'échantillonnage :
BP_BeechRect<-ppp(BP_Beech[,1],BP_Beech[,2],window=BP_Rect)
# visualiser le semis :
plot(BP_BeechRect)
Contrairement à ADS, spatstat n'élimine pas automatiquement les points
situés à l'extérieur du domaine d'étude. Lancer une analyse de
BP_BeechRect provoquerait une erreur d'éxecussion.
# réduire le semis aux points contenus dans la fenêtre d'échantillonnage
rectangulaire :
BP_BeechRect2<-subset.ppp(BP_BeechRect,window=BP_Rect)
# visualiser le nouveau semis :
plot(BP_BeechRect2)
# calcul de K(r) pour des distances de 0, 1, 2, …, 25 m avec correction
locale des effets de bord par la méthode de Ripley :
BP_BeechRect2.K<-Kest(BP_BeechRect2,seq(0,25,1),correction="Ripley")
# visualiser les résultats :
plot(BP_BeechRect2.K)
# pour une meilleure lisibilité, on peut calculer L(r) :
BP_BeechRect2.L<-BP_BeechRect2.K
BP_BeechRect2.L$theo<-sqrt(BP_BeechRect2.L$theo/pi)- BP_BeechRect2.L$r
BP_BeechRect2.L$iso<-sqrt(BP_BeechRect2.L$iso/pi)- BP_BeechRect2.L$r
plot(BP_BeechRect2.L)
On peut faire le même genre de manip pour calculer la fonction intertype
L12(r) à partir de Kcross{spatstat}
Cordialement,
-- *********************************************** Raphael Pelissier UMR AMAP TA40/PS2 34398 Montpellier cedex 5 France Tel. +33 (0)4 67 61 75 23 Fax. +33 (0)4 67 61 56 68 RP Home Page http://pelissier.free.fr/ ************************************************Alain Bellido a écrit :
> Bonjour > > J'ai utilise le module Ripley de ADE "classique", ca marche bien. > Le pb c'est que je ne trouve pas les memes resultats avec le logiciel > PASSAGE (c'est surement passage qui se trompe...) > Qqun a t il deja eu ce pb ? > > Histoire de me mettre au gout du jour je cherche a calculer K et L de > Ripley sous ADE-R, mais je ne trouve rien. > Il faut dire que je debute sous R ! > Il y a bien aussi le package spatstat mais apparemment il ne calcule > que K. > > Merci de votre aide > Bien cordialement > Alain Bellido >
This archive was generated by hypermail 2b30 : Tue Mar 08 2005 - 17:00:51 MET