search for: hubers

Hello: I found the "hubers" function in MASS library is NOT working on the following data: > a <- >c(7.19,7.19,7.19,9.41,6.79,9.41,7.19,9.41,1.64,7.19,7.19,7.19,7.19,1.422,7.19,6.79,7.19,6.79,7.19,7.19,4.44,6.55,6.79,7.19,9.41,9.41,7.19,7.19,7.19,7.19,1.64,1.597,1.64,7.19,1.422,7.19,6.79,9.38,7.19,1.64,7.19,...

We have a choice when calculating the Huber location estimate: > set.seed(221205) > y <- 7 + 3*rt(30,1) > library(MASS) > huber(y)$mu [1] 5.9117 > coefficients(rlm(y~1)) (Intercept) 5.9204 I was surprised to get two different results. The function huber() works directly with the definition whereas rlm() uses iteratively reweighted least squares. My surprise is

How to get robust M-estimators of multivariate scatter using Huber's psi? Which package/function should I look into? Ideally, I hope I can self-define thresholds of Huber's psi function. Thanks a lot!!! -- View this message in context: http://www.nabble.com/How-to-get-robust-M-estimator-of-multivariate-scatter-using-Huber%27s-psi--tp20585755p20585755.html Sent from the R help mailing

Hi Everyone, Given the scenario I have, I was wondering if anyone would be able to give me a hind on how to get the results from hubers() in a more efficient way. I have an outcome on an array [N x S x D]. I also have a factor (levels 1,2,3) stored on a matrix N x S. My objective is to get "mu" and "sigma" for each of the N rows (outcome) stratified by the factor (levels 1, 2 and 3) for each of the D &q...

R-experts, I have fitted many different lines. The fast-tau estimator (yellow line) seems strange to me?because this yellow line is not at all in agreement with the other lines (reverse slope, I mean the yellow line has a positive slope and the other ones have negative slope). Is there something wrong in my R code ? Is it because the Y variable is 1 vector and should be a matrix ? Here is the

You need to pay attention to the documentation more closely. If you don't know what something means, that is usually a signal that you need to study more... in this case about the difference between an input variable and a design (model) matrix. This is a concept from the standard linear algebra formulation for regression equations. (Note that I have never used RobPer, nor do I regularly

Good day, The use of "rounding" also doesn't make sense. If The number is halfway between two integers, it is rounded to the nearest even integer. > round(2.5) [1] 2 -------------------------------------- Dario Strbenac University of Sydney Camperdown NSW 2050 Australia

Hi, I am finding the recently added [1] functionality of embedding raster images into plots on R devices very useful! Thanks to Paul Murrell and others for providing that. I noted that in https://svn.r-project.org/R/trunk/src/library/grDevices/src/devPS.c a macro is defined: #define MAX_RASTERS 64, and consequently, I get Error in grid.Call.graphics("L_raster", x$raster, x$x, x$y,

FWIW, I suspect this is related to the function R_unif_index that was introduced in src/main/RNG.c around revision 72356, or the way this function is used in do_sample in src/main/random.c. 20.9.18 08:19, Wolfgang Huber scripsit: > Besides wording of the documentation re truncating vs rounding, there is > something peculiar going on with the fractional part of n: > > >

Hello, first of all: I'm new to R and have only used SPSS befor this (which can't do this at all...). I'm trying to output some regression results to latex. The regressions are normal OLS and I'm trying to output the results with standardized regression coefficients and t-statistics based on "Huber-White sandwich estimator for variance". The final result should be

I've used optim to analyze some data I have with good results, but need to correct the var-cov matrix for possible effects of clustering of observations (respondents) in small groups (non-independence). Is there any function to adjust the matrix? I heard some time ago that the vcovHC function would have a cluster capability added to it, but I don't see that in my fairly recent version.

Since now many browsers support (ECMA/Java-)scripted SVG, I am wondering whether there are already any examples of inserting R code into SVG documents (or a Javascript canvas?) either directly, or perhaps more likely through a JavaScript layer, to dynamically generate graphics or make them interactive? I am aware of the excellent packages gridSVG and SVGAnnotation, which facilitate making

...de is run z.score<-function (group) { for (t in levels(group$Country)) { # this will give t countries y<-subset(group,factor(Country)==t) #particular analyte over all countries equal to y #calculate overall huber mean and sigma for a particular analyte over all countries ov.mu<-hubers(group$X)$mu ov.sigma<-hubers(group$X)$s #define arrays #p.mu<-array() #p.sigma<-array() #z.value<-array() #calculate huber mean and sigma for given analyte (defined by group) by selected country (y) p.mu<-hubers(y$X)$mu #p.sigma<-hubers(y$X)$s #calcul...

...de is run z.score<-function (group) { for (t in levels(group$Country)) { # this will give t countries y<-subset(group,factor(Country)==t) #particular analyte over all countries equal to y #calculate overall huber mean and sigma for a particular analyte over all countries ov.mu<-hubers(group$X)$mu ov.sigma<-hubers(group$X)$s #define arrays #p.mu<-array() #p.sigma<-array() #z.value<-array() #calculate huber mean and sigma for given analyte (defined by group) by selected country (y) p.mu<-hubers(y$X)$mu #p.sigma<-hubers(y$X)$s #calcul...

Hi all, I'm trying to fit a nonlinear regression by "nlrob": model3=nlrob(y~a1*x^a2,data=transient,psi=psi.bisquare, start=list(a1=0.02,a2=0.7),maxit=1000) However an error message keeps popping up saying that the function psi.bisquare doesn't exist. I also tried psi.huber, which is supposed to be the default for nlrob: model3=nlrob(y~a1*x^a2,data=transient,psi=psi.huber,

Dear all, is there a good way to create SVG plots with R whose elements have titles (tooltips) or act as hyperlinks? I am using the RSvgDevice package, which works great - but it doesn't seem to support the notion that plot objects have titles or are act as hyperlinks, so I am helping myself by giving the objects funny unique colors and then postprocessing the .svg file. I wonder

hi guys and gals ... How are you all ... i have to do something in robust regression by R programm , and i have some problems as following: *the first :* suppose w(r) =1/(1 r^2) and r <- c(7.01,2.07,7.061,5.607,8.502,54.909,12.222) and i want to exclude some values from r so that (abs(r)>4.9 )... after ,i want to used (w) to get on coefficients beta0 and beta1 (B1 <-

Hi Is it possible to subset an n-dimensional array by a vector of n dimensions? E.g. assume I have > x <- array(1:24, dim=2:4) > x[1,1,2] [1] 7 > dims <- c(1,1,2) I would like a function that I can supply x and dims as parameters to, and have it return 7. Also, I would like to do something like: > x[1,1,] [1] 1 7 13 19 > dims2<- c(1,1,NA) And have a function of x and

On 31/03/2018 11:57 AM, varin sacha via R-help wrote: > Dear R-experts, > > Here below my reproducible R code. I want to add many straight lines to a plot using "abline" > The last fit (fast Tau-estimator, color yellow) will not appear on the plot. What is going wrong ? > Many thanks for your reply. > It's not quite reproducible: you forgot the line to create

Dear Everybody! I want to export data to LaTeX. As I want to employ the data as freely as possible I want to avoid the xtable-command and instead generate some List like \MyOwnPrettyCommand{Adam}{Auer}{17} \MyOwnPrettyCommand{Bertram}{Bauer}{14} \MyOwnPrettyCommand{Christoph}{Huber}{75} \MyOwnPrettyCommand{Damian}{Dorfer}{69} \MyOwnPrettyCommand{Emanuel}{Eder}{43} with \MyOwnPrittyCommand