Soucis avec le package Gstat

[Thomas Galtier]

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]

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]

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]

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]

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]

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]

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]

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]

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]

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]

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]

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]

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]

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]

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)