Discussion :

[sassélesas]

Bonjour,
Me voila avec un nouveau problème, je souhaiterais faire une moyenne pondérée par groupe, seulement j'ai un sérieux souci.
Exemple:

NOMVAR	FLUX	pond
FC	220782663	1
FC	344521891	1
FC	222356	1.07
TG	400814730	1.32
TG	229027	1.20
TE	58879013	1
TE	632149870	1.031
TE	456210	1.19
TE	87953	1

avant de faire la pondération, j'ai supprimé toutes les observations de TE et TG, et je n'ai laissé que les FC, comme suit:

NOMVAR	FLUX	pond
FC	220782663	1
FC	344521891	1
FC	222356	1.07

après ça, j'ai essayé de faire une moyenne pondérée avec une fonction que j'ai trouvée sur ce forum:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
tab_etu= tapply(seq_along(tab_fin$FLUX),tab_fin$Nomvar
function(xx){return(weighted.mean(x=tab_fin$FLUX[xx],w=tab_fin$pond))})

Le problème:
1. j'ai dû transposer la table que j'ai obtenue parce qu'il a considéré les colonnes comme des lignes ;
2. j'ai eu les variables TE et TG dans les résultats mais qui ont des valeurs manquantes.

NOMVAR	FLUX_pond
FC	123654
TG	NA
TE	NA

Une solution?
Merci d'avance

[ced]

Le problème pourrait venir du fait que NOMVAR est un facteur, et qu'il possède dans ses levels TG et TE.
Pour le vérifier, que donne :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

levels(tab_fin$Nomvar)

Si ça renvoie FC, TG et TE, il faut que vous limitiez les niveaux du facteur à FC seulement, comme suit :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

levels(tab_fin$Nomvar) <- c("FC")

Sinon, donnez-nous le résultat du code suivant :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

str(tab_fin)

[sassélesas]

Voilà le code d'erreurs de la fonction levels(tab_fin$Nomvar) <- c("FC"):

Error in `levels<-.factor`(`*tmp*`, value = "FC") :
number of levels differs

C'est peut-être dû au fait que le nomvar n'est pas exactement FC, mais il y a plusieurs FC du genre FC_1236547 (le FC ou TG OU TE sont suivis d'une série de chiffres) ( j'avais simplifié dans l'exemple).
Pour le résultat de ce code str(tab_fin) est le suivant:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
'data.frame':   2298 obs. of  4 variables:
 $ Codbanq           : int  6 6 6 6 6 6 6 6 6 6 ...
 $ Nomvar            : Factor w/ 446 levels "MFCDTA_0010E6NREM1ETTETTES5ATAT",..: 124 138 151 165 179 89 100 113 94 105 ...
 $ MTOBS             : num  2.21e+08 3.45e+08 3.45e+08 4.01e+08 4.01e+08 ...
 $ pond_corrigee_8097: num  1 1 1 1 1 1 1 1 1 1 ...