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.
How to Save the residuals of an LM object greater or less than a certin value to an R object?
3 messages · Faiz Rasool, Jeff Newmiller, Peter Dalgaard
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.
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