Re: [R] Least Median Square Regression

2016-10-08 Thread Enrico Schumann
On Sat, 08 Oct 2016, Bryan Mac writes: > I am confused reading the document. > > I have installed and added the package (MASS). > > What is the function for LMS Regression? > In MASS, it is 'lqs'. But the vignette provides a code example for how to compute 'manually' an LMS-regression, i.e. ho

Re: [R] Least Median Square Regression

2016-10-08 Thread Bert Gunter
Well, first of all, note that there is no "lms" method for the stats package's lm() function. You can't just make stuff up, you know! And second, ?lmsreg -- after loading MASS via library(MASS), if you haven't already done this after your install -- is what you want. Other than ?lmsreg and what E

Re: [R] Least Median Square Regression

2016-10-08 Thread Bryan Mac
I am confused reading the document. I have installed and added the package (MASS). What is the function for LMS Regression? Bryan Mac bryanmac...@gmail.com > On Oct 8, 2016, at 6:17 AM, Enrico Schumann wrote: > > On Sat, 08 Oct 2016, Bryan Mac writes: > >> Hi R-help, >> >> How do you p

Re: [R] Least Median Square Regression

2016-10-08 Thread Enrico Schumann
On Sat, 08 Oct 2016, Bryan Mac writes: > Hi R-help, > > How do you perform least median square regression in R? Here is what I have > but received no output. > > LMSRegression <- function(df, indices){ > sample <- df[indices, ] > LMS_NAR_NIC_relation <- lm(sample$NAR~sample$NIC, data = samp

Re: [R] Least Median Square Regression

2016-10-08 Thread ruipbarradas
Hello, Use package quantreg, function rq(). install.packages("quantreg") ?rq Hope this helps, Rui Barradas Citando Bryan Mac : Hi R-help, How do you perform least median square regression in R? Here is what I have but received no output. LMSRegression <- function(df, indices){ sampl

[R] Least Median Square Regression

2016-10-08 Thread Bryan Mac
Hi R-help, How do you perform least median square regression in R? Here is what I have but received no output. LMSRegression <- function(df, indices){ sample <- df[indices, ] LMS_NAR_NIC_relation <- lm(sample$NAR~sample$NIC, data = sample, method = "lms") rsquared_lms_nar_nic <- summary(