Bonjour, vous avez moyen de faire la même chose en beaucoup plus simple avec la fonction dredge() du package MuMIn : modele.complet <- glm(nat~dbh+decay+site,family=binomial(link="logit")) # Ce qui a l'air d'être le modèle "complet" dans votre cas, mais vous pouvez inclure des in...
En restant sur votre example et en utilisant le package diagram, on peut faire un truc comme ça : par(mar=par("mar")+c(3,0,-3,1)) plot(1,xlim=c(0,10),ylim=c(0,1)) mtext(text=c(1149,855,628,459,325,214,132,80,29,2,0),side=1,at=c(0:10),line=4.5,cex=0.8) mtext(text=c(1149,855,628,459,325,214,...
Bonjour, une solution fourbe et sans doute pas optimale en utilisant notamment textrect() du package diagram : plot(1,xlab="") # Toutes les commandes additionnelles du graphe ici par(new=TRUE,mar=c(0,0,0,0)) plot(1,xlim=c(0,1),ylim=c(0,1),axes=FALSE,type="n") textrect(c(0.53,0.02...
Vu comment le code a été écrit, je pense que la solution la plus simple est de créer un sous-jeu de données avec seulement les 436 lignes qui vous intéressent, et d'ajouter maxpos comme colonne supplémentaire. Puis de travailler là-dessus : data2 <- data[indmax,] data2$maxpos <- maxpos continue <- T...
Y a un truc que je pige pas, votre jeu de données réel n'a pas l'air fichu comme celui que vous avez donné en exemple. Pourquoi les variables indmax et maxpos ne sont pas dans le tableau 'data' ? Est-ce qu'indmax fait référence seulement à certaines lignes du tableau, et que la boucle ne doit tourne...
Si le code est bien adapté à votre jeu de données, il devrait fonctionner même si c'est long. Mais en effet le jeu de données est assez grand donc je ne suis pas surpris que ça mette un certain temps, voire un temps certain... Personnellement je ne vois pas de moyen facile d'accélérer tout ça, mais ...
Si j'ai bien compris l'idée du code : - On regarde si l'écart est d'un jour ou non (donc 24h hein ?) - si écart < 1 jour on supprime la ligne du jeu de données - sinon on ne touche pas à la ligne et on l'a va comparer la prochaine valeur avec la ligne que l'on a pas touché Oui c'est globalement ça....
Bonjour, je tente un truc : jeu <- data.frame(date_heure=as.POSIXct(c("2016-01-01 01:00","2016-01-01 10:00", "2016-01-02 03:00","2016-01-02 08:00","2016-01-03 01:00")),Max=c(9,7,5,6,8)) jeu date_heure Max 1 2016-01-01 01:00:00 9 2 2016-01-01 10:00:00...
Bonjour Jérémie, par curiosité, y a-t-il un problème avec cette fonction ? Je vais mettre à jour le package dans quelques jours donc autant corriger s'il y a un bug. Maxime
C'est un simple problème d'importation du tableau. Vous avez des virgules qui séparent vos unités des décimales, or R est anglo-saxon donc le séparateur décimal est le point. Ajoutez dec="," à read.table().
la fonction dudi.pca() ne le permet pas en tout cas. Regardez la fonction pca() du package mixOmics pour une ACP avec des données manquantes (pas trop quand même, hein).