Soucis avec le package Gstat

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

Thomas Galtier
Messages : 93
Enregistré le : 14 Mai 2007, 12:08

Soucis avec le package Gstat

Messagepar Thomas Galtier » 21 Mai 2007, 09:18

Bonjour tout le monde,

est-ce que quelqu'un a déjà utilisé Gstat avec R ?

Là le souci que j'ai c'est que sur la doc, en exemple ils disent de taper par exemple :

data(zinc): 'zinc.eas', x=1,y=2,v=3;

Le truc c'est que dans data() j'ai pas de variables zinc par exemple.
COmment le creer ? Qu'est-ce qu'un fichier .eas ?

Merci par avance,

Thomas

Logez Maxime
Messages : 2897
Enregistré le : 26 Sep 2006, 11:35

Messagepar Logez Maxime » 21 Mai 2007, 09:39

Bonjour,

Il semblerait que le jeu de données utilisé ne soit pas "zinc" mais "meuse", qui contient bien une variable zinc.

Maxime

Thomas Galtier
Messages : 93
Enregistré le : 14 Mai 2007, 12:08

Messagepar Thomas Galtier » 21 Mai 2007, 10:19

Merci Maxime,

Mais comment faire avec un autre jeux de données ?
Enfin un jeu de données perso... creer quelque chose pour ecrire par exemple : data (temp): ... ?

encore merci ;)

Nicolas Péru
Messages : 1408
Enregistré le : 07 Aoû 2006, 08:13

Messagepar Nicolas Péru » 21 Mai 2007, 11:21

il ne sert à rien de vouloir créer ce genre d'objet. Data() est fait pour appeler les données contenues dans R. Ce n'est pas ça qui rend le jeu de données compatible avec gstat.
Je crois (de mémoire) qu'il faut que tu crées un objet de type gstat avec la fonction gstat() et après tu pourras utiliser la plupart des fonctions contenues dans le package gstat.

Thomas Galtier
Messages : 93
Enregistré le : 14 Mai 2007, 12:08

Messagepar Thomas Galtier » 21 Mai 2007, 11:27

Merci Nicolas,

En fait je ne sais pas creer un objet de type gstat ... enfin d'un point de vu code. J'ai un fichier .txt avec trois colonnes qui sont trois variables differentes et j'aimerai implanter ça dans gstat car en gros j'ai ma longitude, ma latitude et mon jeu de donnée ... mais d'un point de vu code j'y comprend rien.
Pourtant j'ai la doc R et une doc Gstat, mais les deux que j'ai sont conradictoires.
J'ai déjà reussi à lire le fichier (de façon pas très simple d'ailleurs) ou meme creer une data.frame mais creer un objet gstat ... pas reussi

Merci encore à vous ;)

Nicolas Péru
Messages : 1408
Enregistré le : 07 Aoû 2006, 08:13

Messagepar Nicolas Péru » 21 Mai 2007, 11:35

Donne nous le code que tu as employé pour essayer de créer ton objet gstat et le message d'erreur qui s'y rapporte. Je ne suis pas certains de pouvoir t'aider mais ça ne coute rien d'essayer et ensuite peut être que quelqu'un de plus compétent pourra te donner une réponse mais il faut être plus précis dans ta demande en donnant des éléments sur tes données, ce sur quoi tu es bloqué exactement etc... ;)

Thomas Galtier
Messages : 93
Enregistré le : 14 Mai 2007, 12:08

Messagepar Thomas Galtier » 21 Mai 2007, 11:43

Bon voilà mes données :

Long Lat Temp
1.6 50.73 28.063
1.83 50.13 32.024
2.33 51.05 26.849
3.1 50.56 31.984
3.15 50.21 32.651
-1.93 49.71 23.19
-1.48 49.65 25.2
-0.45 49.18 33.008
0.06 49.5 29.389
-0.76 49.35 23.491
0.16 49.36 29.991
1.18 49.38 33.77
1.23 49.01 35.079
-1.5 49.51 27.13
1.1 49.93 25.579
0.08 49.53 29.694
2.11 49.46 35.286
..............................

que j'ai reussi à stocker grace à cette commande :

Code : Tout sélectionner

mydata_temp<-read.table("mydata_temp.txt");


Donc c'est bienstocker dans R.

Après j'ai testé ça :

Code : Tout sélectionner

library(gstat)
Le chargement a nécessité le package : sp
> gstat(mydata_temp):'mydata_temp.txt';
Erreur dans inherits(formula, "formula") : l'argument "formula" est manquant, avec aucune valeur par défaut
> gstat(mydata_temp):'mydata_temp.txt',x=1,y=2,v=3;
Erreur : syntax error, unexpected ',', expecting '\n' or ';' dans "gstat(mydata_temp):'mydata_temp.txt',"


En fait le code pour data() mais avec gstat à a place.

Merci encore à toi ;)

t merci par avance à chaque personne qui pourra m'aider ou essayera en tout cas c'est super sympa :)

Thomas[/list]

Nicolas Péru
Messages : 1408
Enregistré le : 07 Aoû 2006, 08:13

Messagepar Nicolas Péru » 21 Mai 2007, 12:05

tape

Code : Tout sélectionner

?gstat
et tu as l'exemple avzec les données meuse pour la création d'un objet gstat...ce n'est pas le code que tu emplois

Thomas Galtier
Messages : 93
Enregistré le : 14 Mai 2007, 12:08

Messagepar Thomas Galtier » 21 Mai 2007, 12:09

Merci nicolas,

Mais le souci dans leur exemple c'est que ça fait appel quand meem à la data(meuse) or moi je prends mes données dans un fichier texte en gros.

Logez Maxime
Messages : 2897
Enregistré le : 26 Sep 2006, 11:35

Messagepar Logez Maxime » 21 Mai 2007, 12:13

Re,


Peut importe que dans leur exemple il fasse appel a data(meuse) la synthaxe data = "data.frame", fonctionne aussi normalement, donc tu dois pouvoir te servir de mydata_temp.

Maxime

Nicolas Péru
Messages : 1408
Enregistré le : 07 Aoû 2006, 08:13

Messagepar Nicolas Péru » 21 Mai 2007, 12:19

En fait pour être plus précis, le fait de taper data(meuse) appelle l'objet meuse contenu dans la librairie Gstat.
Dans cet objet meuse tu as tout un ensemble de chose dont un tableau de données qui est un objet de class data.frame. Si tu fais class(mydata_temp) tu verras qu'il s'agit aussi d'un objet de class data.frame Donc tu peux faire les mêmes choses avec ton tableau que ce qui est fait dans l'exemple sur meuse en adaptant avec ton nom de tableau comme te le dis maxime.

Thomas Galtier
Messages : 93
Enregistré le : 14 Mai 2007, 12:08

Messagepar Thomas Galtier » 21 Mai 2007, 12:20

Merci maxime,

Mais en fait je ne vois pas du tout comment faire.
Ici l" id" je ne sais pas comment il s'apelle dans mon fichier texte, ou meme dans ma data.frame si je la recharge.

Le reste, c'est surtout des options donc ça avec la commande ? gstat j'ai ce qu'il me faut mais le coup de charger les données j'y arrive pas ...

Thomas Galtier
Messages : 93
Enregistré le : 14 Mai 2007, 12:08

Messagepar Thomas Galtier » 21 Mai 2007, 12:22

Merci nicolas j'avais pas vu ta reponse :)

Donc en fait au lieu de taper data= meuse, je tape : data=class(mydata_temp) ça devrait passer ?

J'essaye au passage de mieux comprendre l'histoire de classe, objet, etc ...

Logez Maxime
Messages : 2897
Enregistré le : 26 Sep 2006, 11:35

Messagepar Logez Maxime » 21 Mai 2007, 12:25

re,

Pourquoi te compliquer la vie, data dois faire référence a des données si tu tapes data = class(mydata_temp), ça ne marchera pas car class(mydata_temp) te donne un vecteur de taille 1 : "data.frame" ! Il te suffit de taper data = mydata_temp.

Maxime

Thomas Galtier
Messages : 93
Enregistré le : 14 Mai 2007, 12:08

Messagepar Thomas Galtier » 21 Mai 2007, 12:33

Et pour l' "ID" je fais comment pour lui donner un nom ou l'apeller à travers la fonction gstat( id =, formula =, locations = ... ,data=mydata_temp)


Retourner vers « Questions en cours »

Qui est en ligne

Utilisateurs parcourant ce forum : goel fenech et 1 invité