[R] lm and NA

Peter Alspach PAlspach at hortresearch.co.nz
Mon Jun 30 23:30:46 CEST 2008


Georg

?lm suggests "'na.exclude' can be useful"

Thus:

> round(resid(lm(c(1:4,NA,9)~c(3:7,NA), na.action=na.exclude)), 5)
 1  2  3  4  5  6 
 0  0  0  0 NA NA 

HTH ....

Peter Alspach


> -----Original Message-----
> From: r-help-bounces at r-project.org 
> [mailto:r-help-bounces at r-project.org] On Behalf Of Georg Ehret
> Sent: Tuesday, 1 July 2008 9:07 a.m.
> To: r-help
> Subject: [R] lm and NA
> 
> 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 at 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.
> 

The contents of this e-mail are privileged and/or confidential to the named
 recipient and are not to be used by any other person and/or organisation.
 If you have received this e-mail in error, please notify the sender and delete
 all material pertaining to this e-mail.



More information about the R-help mailing list