Création d'une boucle avec R

[invitec4d2e6d6]

Bonjour,

Je voulais savoir comment créer une boucle avec le Logiciel R.

Je dispose d'un table de données ce présentant sous cette forme

Année Effectif
1999 12
1999 13
1999 10
2000 9
2000 8
2000 16
... ...

Mon tableau comprenant 30 ans de mesures, je voulais savoir comment créer une boucle qui me permettrais d'avoir la somme des effectifs pour chaque année, afin de gagner du temps.

Merci d'avance
-----

[invitec5eb4b89]

Juste sur le petit exemple que tu as donné :

tab <-  data.frame(Annee=rep(1999:2000,c(3,3)),Effectif=c(12,13,10,9,8,16))
tab 
un <- unique(tab$Annee)
somme <- c()

for (i in seq_along(un) ) somme[i] <- sum(tab$Effectif[tab$Annee==un[i]])
somme

[invitec4d2e6d6]

Merci pour l'info.

Cordialement.

Falindor

[invitee3e4ae46]

Étant nouvelle sur le forum, bonjour à tous,

J'ai le même problème que Falindor : je ne sais pas faire une boucle sous R !!

J'ai suivi le dégagement d'oxygène pendant 50 jour sur 100 échantillons. Tous les échantillons n'ont pas été mesurés le même nombre de fois (1 mesure, 1 jour).

Mon tableau s'organise de la sorte :

Ech . O2 . Jour
1 . 0,1 . 1
1 . 0,4 . 5
1 . NA . 13
... ...
100 . 0,3 . 50

J'aimerais calculer le dégagement gazeux cumulé total pour chaque échantillon et récupérer les résultats dans un tableau.

La fonction pour calculer l'aire total (soit le dégagement gazeux cumulé total) est :
sum(lintegrate(Jour, O2, xint=Jour)) package tis

Comment faire une boucle pour m'éviter de faire ce calcul 100 fois.

Je remercie d'avance mon/mes sauveur(s).

Aneh

[A voir en vidéo sur Futura]

[invite986312212]

tu commences par créer une fonction qui fait le calcul pour un seul échantillon.
disons
aire <- function(x)sum(lintegrate(x$Jo ur,x$O2, xint=x$Jour))

et ensuite tu fais (en supposant que ton tableau s'appelle "donnees")

resultat <- tapply(donnees, donnees$Ech,aire)

[invitee3e4ae46]

tu commences par créer une fonction qui fait le calcul pour un seul échantillon.
disons
aire <- function(x)sum(lintegrate(x$Jo ur,x$O2, xint=x$Jour))

et ensuite tu fais (en supposant que ton tableau s'appelle "donnees")

resultat <- tapply(donnees, donnees$Ech,aire)

Quand je tape ces lignes de codes, le message suivant s'affiche :

Erreur dans tapply(donnees, donnees$Ech, aire) : 
  les arguments doivent avoir même longueur

Que dois-je faire, une boucle pour remplacer le tapply ?

[invite986312212]

ah c'est vrai...

tu peux toujours faire une boucle

si la variable Ech contient les nombres 1,2,... (dans n'importe quel ordre)

resultat<-rep(NA,length(unique(donnees$E ch)))

for (e in unique(donnees$Ech))
resultat[e]<-aire(donnees[donnees$Ech==e,])

si Ech contient un code unique qui repère chaque échantillon, mais pas nécessairement les premiers entiers, il faut remplacer Ech par as.numeric(as.factor(Ech)) par exemple (il y a d'autres façons de faire)
où aire est définie comme précédemment, mais il y a moyen de vectoriser ça, je vais y penser.