Cette page appartient aux archives web de l'EPFL et n'est plus tenue à jour.
This page belongs to EPFL's web archive and is no longer updated.

Exercice 1, série 11

Bonjour,

il me semble que lorsque j'utilise la fonction bisquare, ma regression robuste fonctionne. Cependant quand je veux plotter cette regression, je reçois le message d'erreur suivant:

rfit2 <- rlm(calls~year, phones, method="M", psi=psi.bisquare)

plot(rfit2)

 

Erreur dans lm.influence(x, do.coef = FALSE) : 
  non-NA residual length does not match cases used in fitting
 
je ne comprends pas...
 
 
Posted by Mahé Fellay on Sunday 15 December 2013 at 11:05
Comments
Your regression is probably working fine, but you cannot use the plot command with the output of the rlm function. To plot the regression line, you should use the abline command instead as instructed in the exercise sheet.
Posted by Mikael Kuusela on Sunday 15 December 2013 at 15:25
Thank's!
Posted by Mahé Fellay on Sunday 15 December 2013 at 15:31