[R] lm and NA

[Rolf Turner]

On 1/07/2008, at 9:30 AM, Peter Alspach wrote:

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

	But he (Georg) did indeed use na.action=na.exclude !!!

	So that's not the problem/solution.

	The difference between what he did and what you did is
	essentially that you (correctly) used the ***resid()*** function to
	extract the residuals.  Whereas Georg used (whatever)$residuals
	--- as one is repeatedly advised ***NOT*** to do!

		cheers,

			Rolf Turner
>
> 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.
>
> ______________________________________________
> 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.


######################################################################
Attention:\ This e-mail message is privileged and confid...{{dropped:9}}