class kde() --> conversion --> data.frame X,Y,Z

Postez ici vos questions, réponses, commentaires ou suggestions - Les sujets seront ultérieurement répartis dans les archives par les modérateurs

Modérateur : Groupe des modérateurs

Nicolas Bergeron
Messages : 17
Enregistré le : 13 Déc 2007, 16:18

class kde() --> conversion --> data.frame X,Y,Z

Messagepar Nicolas Bergeron » 28 Fév 2008, 21:32

Bonjour aux utilisateurs R,

J'ai un problème pour convertir un objet de class "ker" (library(ks)) vers un data frame qui donne, pour chaque ligne, la localisation X,Y et la valeur Z (fonction de densité) estimé par la fonction ker().

Voici le début:

library(ks)

HpiScale <- Hpi(x=xy, pilot="samse", pre="scale")

J'ai pour résultat ceci:

> str(fhat.HpiScale)
List of 5
$ x : num [1:23, 1:2] -304031 -305444 -304087 -305393 -305270 ...
..- attr(*, "dimnames")=List of 2
.. ..$ : chr [1:23] "426" "427" "428" "429" ...
.. ..$ : chr [1:2] "x" "y"
$ eval.points:List of 2
..$ : num [1:151] -307405 -307343 -307282 -307220 -307158 ...
..$ : num [1:151] 743923 743977 744031 744085 744139 ...
$ estimate : num [1:151, 1:151] 5.08e-13 8.23e-13 1.30e-12 2.00e-12 2.99e-12 ...
$ H : num [1:2, 1:2] 149287 27563 27563 147516
$ binned : logi FALSE
- attr(*, "class")= chr "kde"

La fonction a calculer la matrix "$estimate" pour chacune des combinaisons $eval.points, soit x et y où "$estimate" a été évalué. Je voudrais produire un data.frame qui contiendrait pour chaque ligne, la combinaison X,Y,Z, soit:

xyz <- data.frame(X= fhat.HpiScale$eval.points[1],
Y= fhat.HpiScale$eval.points[2],
Z= ??) # Matrix ??

Puis-je utiliser lapply ???? Je suis un peu confus

Je pourrais ainsi importer le tout dans GRASS avec le module r.in.xyz.

javascript:emoticon(':roll:')

POuvez-vous me donner une piste de solution s-v-p?


Merci!
NicB

Retourner vers « Questions en cours »

Qui est en ligne

Utilisateurs parcourant ce forum : Aucun utilisateur enregistré et 1 invité