faiz rasool
2018-Feb-23 04:09 UTC
[R] How to Save the residuals of an LM object greater or less than a certin value to an R object?
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.
Jeff Newmiller
2018-Feb-23 04:56 UTC
[R] How to Save the residuals of an LM object greater or less than a certin value to an R object?
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.
Peter Dalgaard
2018-Feb-23 11:40 UTC
[R] How to Save the residuals of an LM object greater or less than a certin value to an R object?
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
Seemingly Similar Threads
- How to Save the residuals of an LM object greater or less than a certin value to an R object?
- Two questions, first about contingency tables, and second about table () and data.frame (), from a visually impaired user.
- Problem in performing goodness of fit test in R.
- How to Save the residuals of an LM object greater or less than a certin value to an R object?
- couple of how-to-do it in R questions regarding corelations and mean and SD of likert items