Erreur dans dataframe *tmp* remplacement de données x a y

Je suis un débutant dans R. Voici un code très simple où je suis en train d'essayer de sauver la durée de vie résiduelle:

# Create variables for child's EA:

dat$cldeacdi <- rowMeans(dat[,c('cdcresp', 'cdcinv')],na.rm=T)
dat$cldeacu <- rowMeans(dat[,c('cucresp', 'cucinv')],na.rm=T)

# Create a residual score for child EA:

dat$cldearesid <- resid(lm(cldeacu ~ cldeacdi, data = dat))

J'obtiens le message suivant:

Error in `$<-.data.frame`(`*tmp*`, cldearesid, value = c(-0.18608488908881,  : 
  replacement has 366 rows, data has 367

J'ai cherché pour cette erreur, mais ne pouvais pas trouver quelque chose qui pourrait résoudre ce problème. En outre, j'ai créé exactement le même code pour la maman de l'EA, et il a sauvé le résiduel très bien, avec pas d'erreurs. Je serais reconnaissant si quelqu'un pouvait m'aider à résoudre ce problème.

OriginalL'auteur Marishka Usacheva | 2017-11-10

  1. 8

    J'ai le sentiment que vous avez NAs dans vos données. Regardez cet exemple:

    #mtcars data set
test <- mtcars
#adding just one NA in the cyl column
test[2, 2] <- NA

#running linear model and adding the residuals to the data.frame
test$residuals <- resid(lm(mpg ~ cyl, test))
Error in `$<-.data.frame`(`*tmp*`, "residuals", value = c(0.382245430809409,  : 
  replacement has 31 rows, data has 32

    Comme vous pouvez le voir il en résulte une erreur similaire à la vôtre.

    Comme une validation:

    length(resid(lm(mpg ~ cyl, test)))
#31
nrow(test)
#32

    Cela se produit parce que lm sera exécuté na.omit sur l'ensemble de données avant l'exécution de la régression, donc si vous avez des lignes avec NA ces seront supprimées résultant en moins de résultats.

    Si vous exécutez na.omit sur votre dat ensemble de données (c'est à dire dat <- na.omit(dat) au début de votre code, puis votre code devrait fonctionner.

    Merci beaucoup! Vous étiez sur le point droit!!
    Heureux de vous aider. N'hésitez pas à accepter si cela a aidé.
    En fait, je viens de remarquer que quand je lance na.omettre efface mon ensemble de données, comme dans - 0 observations de mes variables... Pourquoi est-ce?
    na.omit va supprimer la rangée complète de données.cadre si il y a au moins une colonne de la rangée avec un NA. lm le fait aussi que depuis qu'il appelle na.omit en interne. Assurez-vous d'exclure les colonnes que vous n'avez pas besoin à la régression avant d'exécuter na.omit .

    OriginalL'auteur LyzandeR