search for: hoaglin

...user of R so please bear with me. I have reviewed some R books, FAQs and such but the volume of material is great. I am in the process of porting my current SAS and SVS Script code to Lotus Approach, R and WordPerfect. My question is, can you help me determine the best R method to implement the Hoaglin Outlier Method? It is used in the Appendix A and B of the fo llowing link. http://trb.org/publications/nchrp/nchrp_w71.pdf The sample data from Appendix A for determining outliers in R: T314Data <- structure(list(Lab = as.integer(c(1:60)), X = c(4.89, 3.82, 2.57, 2.3, 2.034, 2, 1.97, 1.85, 1.8...

...sen for this thread, "... does not produce the correct Tukey line", so I think we should get better. Apart from Richie / my MMline() function, I've also noticed that ACSWR :: resistant_line() exists. However "the literature" (see references below), notably the two with Hoaglin, strongly recommends smarter iterations, and -- lo and behold! -- when this topic came up last (for me) in Dec. 2014, I did spend about 2 days work (or more?) to get the FORTRAN code from the 1981 - book (which is abbreviated the "ABC of EDA") from a somewhat useful OCR scan into compila...

>>>>> Serguei Sokol <sokol at insa-toulouse.fr> >>>>> on Mon, 29 May 2017 15:28:12 +0200 writes: > Sorry, I have seen it too late that we had different tab > width in the original file and my editor. Here is the > patch with all white spaces instead of mixing tabs and > white spaces. thank you - it still gives quite a few

Hello R-helpers.. I want to ask about how I can sample data sets without having the infinite numbers coming out. For example, set.seed(1234) a<-rnorm(15,0,1) b<-rnorm(15,0,1) c<-rnorm(15,0,1) d<-rnorm(15,0,36) After come out with the sample, I need to do a transformation (by Hoaglin, 1985) for each data set. Actually I need to measure the skewness and kurtosis, that's why I need to do the transformation. After transformation, there will be 'Inf' value in my data sets and I cannot proceed with the next step where I need to compute the trimmed mean and sum square of...

Hi List, Just wondering where the documentation exists for the statistics which makeup the bwplot. I'm guessing that if R is like similar products that the graph is constructed as The median is the filled circle. The box surrounding the filled circle depicts the 25th and 75th quartile. The range of values is given by the dotted lines (?whiskers?) outside of each box, and possible

Can anyone comment or point me to a discussion of the pros and cons of robust regressions, vs. a more "manual" approach to trimming outliers and/or "normalizing" data used in regression analysis?

Le 31/05/2017 ? 17:30, Serguei Sokol a ?crit : > > More thorough reading revealed that I have overlooked this phrase in the > line's doc: "left and right /thirds/ of the data" (emphasis is mine). Oops. I have read the first ref returned by google and it happened to be tibco's doc, not the R's one. The layout is very similar hence my mistake. The latter does not

...art from the bug which typically shows for non-equidistant x[]}. We could also consider to eventually add a new 'method = <string>' argument to line() one version of which would continue to compute the current solution, another would compute the one corresponding to Velleman & Hoaglin (1981)'s FORTRAN implementation (which had to be corrected for some infinite-loop cases!)... not in the close future though Given all this discussions here, I think I should commit what I currently have ASAP. Martin > x[DELETED ATTACHMENT line.c.patch2, plain text]

Hello R-Help, Could anybody tell me how the boxplot-function calculates the upper and lower extremes of the notch contained in $conf which I assume is the confidence interval? Is it reliable for data which is not normally distributed? If not, how can I calculate and boxplot a specific confidence interval for not normally distributed data in R (increasing the sample size does not normalize the

Thomas, First of all, thanks for the help, but it isn't exactly what I'm looking for. smooth() doesn't perform the smooth the way I want it to do. I want, precisely, the 4253H method. R doesn't give the option to do that. Thanks, Henrique.

Hi, trying to do a robudt regression of a two-way linear model, I keep getting the following error: > lqs(obs ~ y + s -1,method="lms", contrasts=list(s=("contr.sum"))) Error: lqs failed: all the samples were singular Robust regression with M-estimators works (also regular least square fits, of course): rlm.formula(formula = obs ~ y + s - 1, method = "M",

Hello, I'm learning about robust M-estimators right now and had settled on the "Huber Proposal 2" as implemented in MASS, but further reading made clear, that at least 2 further weighting functions (Hampel, Tukey bisquare) exist. In a post from B.D. Ripley going back to 1999 I found the following quote: >> 2) Would huber() give me results that are similar (i.e., close

I would like to do some fairly basic stem-and-leaf displays in R. I am aware (I might even say painfully aware) of stem(base) and have tried it. That's why I'm hoping someone has a usable stem- and-leaf display for R so that I don't have to write my own. r-project.org > Search > R Site Search > "stem and leaf display" finds nothing. I also tried the mail archive

For boxplot(), is it possible to pass in a parameter to change the default way that the 1st and 3rd quartiles are computed? (specifically, I'd like to use type 6 described in the quantile function). Also, what are the options for how outliers are computed, and how can one change them? Thank you [[alternative HTML version deleted]]

...Smooth of a time serie To: " 'Henrique Patr?cio Sant'Anna Branco' " <hpsbranco at superig.com.br>, r-help at stat.math.ethz.ch Message-ID: <3A822319EB35174CA3714066D590DCD50205C7CB at usrymx25.merck.com> Content-Type: text/plain; charset=iso-8859-1 Vellman & Hoaglin's "ABC of EDA" book has listing of Fortran program for that (and other) smoother. You can try to load that into R. Another thing you can try is to port things in the "smoothers" collection on StatLib's S section. That also seems to contain the 3253H smoother. Andy &...