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)$residuals1 2 3 4 0 0 0 0> lm(a~b,na.action=na.exclude)$residuals1 2 3 4 0 0 0 0 Thank you and best regards! Georg. ************************** Georg Ehret Johns Hopkins Baltimore, USA [[alternative HTML version deleted]]
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.
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)$residuals1 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 <georgehret@gmail.com> 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]]