Salut, Sans exemples et codes reproductibles c'est difficile de te proposer une solution. Normalement la fonction xtable du package éponyme permet juste d'afficher et exporter des tables ( data.frame ) en HTML ou LaTeX, donc je ne vois trop comment tu construis tes tables avec xtable . De manière gé...
Comment se présentent vos résultats ? Quels packages ont été utilisés pour l'analyse (gstat ?) ? Avec un bout de code vous aurez rapidement une solution parce que dans l’état actuel, votre question est trop vague.
Salut, Désolé pour le retard, ci-joint une approche qui utilise le package raster et rgdal . Par contre, il faut avoir le driver HDF4 pour ta version de gdal installé sur ta machine. Si jamais tu ne l'as pas et que tu n'arrives pas à compiler une version de gdal avec HDF4, tu peux utiliser le Modis ...
Je vois, serait il possible d'avoir un exemple de fichier hdf en lien (dropbox ou autre) ? Le serveur qui permet de les télécharger est en maintenance (http://nsidc.org/data/MOD10A1).
C'est quoi comme données hdf ? Il ne s'agit pas forcement de données de types hdf5 parfois et très souvent c'est hdf4 (e.g certains produits MODIS). As tu un exemple de fichier pour rendre tout cela plus concret ? Tu es sur Windows ou Unix (Linux, MAC OS) ?
Je ne comprends pas le problème car l'objet xx1 est bien un SpatialPolygon donc tu peux faire: ## Faire la carte plot(xx1) A mon humble avis, il faut d'abord revoir les bases du logiciel R et comprendre comment on manipule les objets. Sinon ce code marche plutôt bien, tu peux t'en inspirer pour fair...
Tu peux aussi utiliser la fonction gsub et sa capacité à faire de la "capture de groupes". x <- c("5 rue de la Source 85150 LA MOTHE ACHARD", "Les Bazinieres 85220 SAINT REVEREND", "RESIDENCE LES BOSQUETS 42 BD MOURAIN DU PATIS 85300 CHALLANS") gsub("(.+\...
Salut, Je sais pas exactement ce que tu veux faire comme traitement, mais si tu veux faire un moyenne de tes 12 fichiers et en faire une carte le plus simple c'est d'utiliser le package raster. Ci-joint un exemple sur 2 fichiers ( à adapter sur tes 12 fichiers): library(raster) library(ncdf4) l <- l...
Salut, Tu peux utiliser la fonction shQuote shQuote(c("a", "b", "c")) [1] "'a'" "'b'" "'c'" Donc tu peux par exemple faire une serie de commande du style: x1 <- shQuote(levels(VarA)[1]) fm1 <- paste0("I(",VarAA,"==",x1,&...
Salut, Avec un si grand nombre de données (relativement aux observations), il peut être plus intéressant d'utiliser les méthodes de pénalisation de style L1 (lasso), L2 (ridge) ou encore une combinaison des deux. Tu peux essayer le package penalizedLDA [1] qui implémente la méthode présenté dans l'e...
Ah ok, je vois. Tu as deux options, tu peux utiliser sprintf ou paste0 pour construire ta formule puis l'utiliser dans la fonction. Tu peux par exemple essayer un truc du genre: fm <- paste0("I(", VarA, " == 1)") svyciprop(reformulate(fm),dw) Ou fm <- printf("I(%s == 1)"...
Tu compliques un peu les choses, j'ai bien dit que reformulate prends en paramètre une chaine de caractère. Donc la commande ci-dessous devrait marcher: svyciprop(reformulate("I(VarA == 1)") ,dw) Si tu as des doutes, cherches a comprendre la différence entre mon code et le tien. Amicalement,
Salut, Tu peux utiliser la fonction reformulate pour construire une formule, elle transforme les chaines de caractères en formule Il suffira juste de changer le second paramètre de ta fonction et mettre la nom de la variable sous forme de caractère. J'ai pas tes données pour tester mais je pense que...
Sans les données, c'est dur de faire un diagnostic. Quand tu utilise par exemple, la base de données Mars (data(Mars)), les commandes marchent bien. Essayes de voir la différence entre ce jeu de données et le tien pour commencer.