[R] How to Save the residuals of an LM object greater or less than a certin value to an R object?

Peter Dalgaard pdalgd at gmail.com
Fri Feb 23 12:40:54 CET 2018


Also,

which( abs( stdresiduals ) > 2.5 )

will tell you which of the standardized residuals are bigger than 2.5 in absolute value. It returns a vector of indices, as in

> set.seed(1234)
> x <- rnorm(100)
> which (abs(x) > 2.5)
[1] 62
> x[62]
[1] 2.548991


-pd

> On 23 Feb 2018, at 05:56 , Jeff Newmiller <jdnewmil at dcn.davis.CA.us> wrote:
> 
> Residuals are stored as a numeric vector. The R software comes with a document "Introduction to R" that discusses basic math functions and logical operators that can create logical vectors:
> 
> abs( stdresiduals ) > 2.5
> 
> It also discusses indexing using logical vectors:
> 
> stdresiduals[ abs( stdresiduals ) > 2.5 ]
> 
> Note that in most cases it is worth going  the extra step of making your example reproducible [1][2][3] because many problems arise from issues in the data or in code that you don't think is broken. 
> 
> [1] http://stackoverflow.com/questions/5963269/how-to-make-a-great-r-reproducible-example
> 
> [2] http://adv-r.had.co.nz/Reproducibility.html
> 
> [3] https://cran.r-project.org/web/packages/reprex/index.html (read the vignette)
> -- 
> Sent from my phone. Please excuse my brevity.
> 
> On February 22, 2018 8:09:49 PM PST, faiz rasool <faiz7r at gmail.com> wrote:
>> Dear list members,
>> 
>> I want to  save residuals above or less than a certain value to an R
>> object. I have performed a multiple linear regression, and now I want
>> to find out which cases have a residual of above + 2.5 and – 2.5.
>> 
>> Below I provide the R  commands I have used.
>> 
>> Reg<-lm(a~b+c+d+e+f) # perform multiple regression with a as the
>> dependent variable.
>> 
>> Residuals<-residuals(reg) # store residuals to an R object.
>> stdresiduals<-rstandard(reg) #save the standardized residuals.
>> #now I want to type a  command that will save the residuals  of
>> certain range to an object.
>> 
>> I realize that by plotting  the data I can quickly  see the  residuals
>> outside a certain boundtry, however, I am totally blind, so visually
>> inspecting the plot is not possible for me.
>> 
>> Thank you for any help.
>> 
>> Regards,
>> Faiz.
>> 
>> ______________________________________________
>> R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see
>> 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.
> 
> ______________________________________________
> R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see
> 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.

-- 
Peter Dalgaard, Professor,
Center for Statistics, Copenhagen Business School
Solbjerg Plads 3, 2000 Frederiksberg, Denmark
Phone: (+45)38153501
Office: A 4.23
Email: pd.mes at cbs.dk  Priv: PDalgd at gmail.com



More information about the R-help mailing list