remplacement des NA par la moyenne de la ligne.

[antoine rivierre]

Bonjour.
Afin d'effectuer une ACP avec une matrice remplie de trous je cherche à remplacer les NA par la valeur moyenne de la ligne.
Comme les données sont centrées et réduites avant l'analyse j'imagine que ceci ne doit pas influencer le résultat de l'ACP??
J'en profite donc également pour vous demander si mathématiquement et statistiquement parlant ceci se justifie. Dans le cas contraire que me conseillez vous??

Merci pour votre aide.

[Renaud Lancelot]

Code : Tout sélectionner

> m <- matrix(rnorm(12), ncol = 3)
> m[c(1, 6, 11)] <- NA
> m
           [,1]       [,2]       [,3]
[1,]         NA  0.3547303  0.5882354
[2,]  0.3930100         NA  0.6428642
[3,]  0.7740637 -0.2386702         NA
[4,] -0.8935756  0.7266868 -0.1749101
>
> fn <- function(x){
+           m. <- mean(x, na.rm = TRUE)
+           x[is.na(x)] <- m.
+           x
+           }
>
> t(apply(m, 1, fn))
           [,1]       [,2]       [,3]
[1,]  0.4714829  0.3547303  0.5882354
[2,]  0.3930100  0.5179371  0.6428642
[3,]  0.7740637 -0.2386702  0.2676967
[4,] -0.8935756  0.7266868 -0.1749101

Cette procédure n'est pas idéale. Elle ne biaise pas la moyenne, mais elle biaise la variance (en la réduisant). Il est préférable d'utiliser des procédures d'imputations multiples. Il y a des packages R pour ça, mais je ne les ai pas utilisés.

[antoine rivierre]

merci beaucoup pour la fonction. Avez vous une idée du nom des packages??

merci Antoine

[Renaud Lancelot]

Code : Tout sélectionner

RSiteSearch("imputation")