Bonjour, En utilisant des expressions régulières ? Un peu de lecture : http://www.duclert.org/Aide-memoire-R/Le-langage/Chaines-de-caracteres.php Un exemple rapide: vecteur=c("1222133123123") positions_sequences <- gregexpr(c("12"), vecteur) nb_sequences <...
Bonjour, Je souhaiterais convertir une date Matlab en format date R. Par exemple, sous MATLAB la valeur de ma date est: 13735196586055 soit le 11-Jul-2013 05:14:18 Ci-dessous le code que j'utilise: data <- readMat(filename_matlab) date.time <- as.Date(((data$time[1,]/86400)+7...
Finalement, je passe par la fonction Rgshhs du package maptools. Il faut télécharger le fichier binaire ici http://www.soest.hawaii.edu/pwessel/gshhg/gshhg-bin-2.3.0.zip puis le dézipper. Et ensuite utiliser le code suivant : cote<-Rgshhs("C:/Users/.../Desktop/GSHHS/gshhg-bin-2.3.0/gshhs_f....
Bonjour, A partir de mon jeu de données (points de mesures effectuées au large des côtes marocaines) je souhaiterais : 1- Cartographier les points avec des valeurs associées en ajoutant le polygone du continent, en distinguant les mesures de jour et de nuit. 2- Calculer la distance entre chaque poin...
Un petit up sur le sujet. Avec la fonction predict du package adehabitat, je pouvais remettre les 34310 observations sur une grille de dimension 360*180 sst_asc <- asc.from.raster(raster(sst_climato)) chl_asc <- asc.from.raster(raster(chloro_climato)) predictors_enf...
Bonjour, En effet, il y a des NA en entrée de la fonction dudi.pca: > str(predictors_enfa) Formal class 'SpatialPixelsDataFrame' [package "sp"] with 7 slots ..@ data :'data.frame': 39642 obs. of 2 variables: .. ..$ SST: num [1:39642] NA NA NA NA NA NA NA NA NA NA ... .. ..$ CHL: nu...
Bonjour, Je souhaite utiliser les fonctions d'ENFA du package 'adehabitatHS'. Une erreur apparaît avec la fonction dudi.pca concernant la présence de NA. Cf. le code ci-dessous: predictors_enfa <- stack(raster(sst_climato),raster(chloro_climato)) predictors_enfa <- as(pre...
Bonjour, Petit update de ce post. La solution trouvée pour l'interprétation du modèle a été de passer par une visualisation en 2D par l'utilisation de l'option 'plot.type="contour"' de la fonction vis.gam du package mgcv : vis.gam(model,view=c("CHLA","SST"),...
Finalement j'ai opté pour la formule suivante : gam(W ~ te(X, log(Y)),data=dataset,family="binomial") En revanche, l'interprétation du graphique résultat du modèle est assez complexe... http://www.casimages.com/img.php?i=14011406325417122511899120.png Si quelqu'un a...
Bonjour, Je teste actuellement un GAM avec la formule suivante : gam(W~s(X,by=log(Y)),data=dataset,family="binomial") Est-ce que le fait d'introduire le paramètre "by" correspond d'une manière à une interaction X*log(Y) ? Il y a peu de documentation sur ce...