similar to: Hoaglin Outlier Method

>>>>> Serguei Sokol <sokol at insa-toulouse.fr> >>>>> on Tue, 30 May 2017 16:01:17 +0200 writes: > Le 30/05/2017 ? 09:33, Martin Maechler a ?crit : ... >> However, even after the patch, The example from the SO >> post differs from the result of Richie Cotton's >> function... > The explanation is quite simple.

>>>>> 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

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?

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

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]]

I have a function of two arguments, (m, n) that returns a list with the same structure as the function makefoo below. When m=1, it returns a list of vectors, each of length n. When m>1, it returns a list of square matrices, each of size m x m. makefoo <- function(m, n) { if (m==1) { A <- sample(1:100, n) B <- sample(1:100, n) C <- sample(1:100, n)

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

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

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

On Sat, Nov 16, 2013 at 11:14 PM, EGO.II-1 <eoconnor25 at gmail.com> wrote: > On 11/16/2013 10:04 PM, Michael B Allen wrote: >> Broke affinity for irq 27 > Found this online....don't know if it pertains to your issue, but check > it out. > > https://www.centos.org/forums/viewtopic.php?t=3941 Unfortunately updating the BIOS did not fix the problem (even though the

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

As it turns out, with the proprietary firmware, the VP3 and VP4 interfaces are identical. Furthermore, this is all already implemented for nvc0. So these patches (a) move the easily sharable bits of the nvc0 implementation into the nouveau directory, and then (b) implement the other parts in nv50. The non-shared parts are still largely copies, but there are some differences, not the least of which

Hello, I have samba 3.20 running in test on Solaris 8 and 9 beautifully it is curently configured as a Domain Client it authenticates using winbind and nsswitch libraries using ADS with windows 2003 server. I have krb5 and ldap working just fine but I now find myself with a mess of bugs and questions and I could really use some of your expertise in the matters so here they go, Thank you soooo

Hi all, I'm using LLVM 3.0, for which I've filed following bug http://llvm.org/bugs/show_bug.cgi?id=12130. I'm trying to solve this problem by myself digging into LLVM sources. It seems that problem that I'm experiencing is related to presence or absence of NSW attribute on a 'mul'. Considering following code: define void @t2(double* %x) { L.entry: %a = alloca [2 x

Hi all, I'm using LLVM 3.0, for which I've filed following bug http://llvm.org/bugs/show_bug.cgi?id=12130. I'm trying to solve this problem by myself digging into LLVM sources. It seems that problem that I'm experiencing is related to presence or absence of NSW attribute on a 'mul'. Considering following code: define void @t2(double* %x) { L.entry: %a = alloca [2 x

Hi all, I'm using LLVM 3.0, for which I've filed following bug http://llvm.org/bugs/show_bug.cgi?id=12130<https://webmail-eu.st.com/owa/redir.aspx?C=f2882890c0064e71ba4c8238dc7fd399&URL=http%3a%2f%2fllvm.org%2fbugs%2fshow_bug.cgi%3fid%3d12130> . I'm trying to solve this problem by myself digging into LLVM sources. It seems that problem that I'm experiencing is related to

Hello all! I have a problem to extract values greater that for example 1820. I try this code: x[x[,1]>1820,]->x1 Please help me! Thank you! The data structure is: structure(c(2.576, 1.728, 3.434, 2.187, 1.928, 1.886, 1.2425, 1.23, 1.075, 1.1785, 1.186, 1.165, 1.732, 1.517, 1.4095, 1.074, 1.618, 1.677, 1.845, 1.594, 1.6655, 1.1605, 1.425, 1.099, 1.007, 1.1795, 1.3855, 1.4065, 1.138, 1.514,

Background: Our subnet:128.2.224.* samba machine 128.2.224.1 (tmc1) AIX 4.2 using precompiled binaries NT Server 128.2.224.6 (tmc6) nt 4.0 workstation 128.2.224.37 (tmc37) smbd/nmbd version 2.0.5a smb.conf ----------------------------------------------- workgroup = TMC guest ok = yes hosts allow = 128.2.224. dns proxy = no #interfaces = 128.2 security = SERVER password server = 128.2.224.6 [tmp]