use the function resid instead of lm.object$residuals to extract the
residuals:
> resid(lm(a~b, na.action=na.exclude))
1 2 3 4
5 6
-2.533445e-17 4.222409e-17 -8.444818e-18 -8.444818e-18
NA NA
> lm(a~b, na.action=na.exclude)$residuals
1 2 3 4
-2.533445e-17 4.222409e-17 -8.444818e-18 -8.444818e-18
On Tue, Jul 1, 2008 at 12:07 AM, Georg Ehret <[EMAIL PROTECTED]> wrote:
> Dear R community, Can I please get some advice on the following: I wish
> to obtain a list of residuals, padded by NAs for NAs in my source data. I
> tried several options of "na.action", but did not succeed...
>
> Example: I would like to get "0,0,0,0,NA,NA":
> > a
> [1] 1 2 3 4 NA 9
> > b
> [1] 3 4 5 6 7 NA
> > lm(a~b)$residuals
> 1 2 3 4
> 0 0 0 0
> > lm(a~b,na.action=na.exclude)$residuals
> 1 2 3 4
> 0 0 0 0
>
> Thank you and best regards!
> Georg.
> **************************
> Georg Ehret
> Johns Hopkins
> Baltimore, USA
>
> [[alternative HTML version deleted]]
>
> ______________________________________________
> R-help@r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide
> http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
>
[[alternative HTML version deleted]]
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
