similar to: High breakdown/efficiency statistics -- was RE: Rosner's test [Broadcast]

My second request of the day, sorry to be such a bother. Can you tell me whether Rosner's test for outliers is available in any of the R packages and, if so, which one? I've tried but I can't find it. Thank you very much, Bob Powell

Hello, Could you please tell me how i can get an input from the user in R? C-Ming Mar 29, 2005 --------------------------------- [[alternative HTML version deleted]]

I hope that this job posting does not offend anyone. My sincere apologies if it is inappropriate -- blame me, not my company. -- Bert Gunter Genentech Nonclinical Statistics South San Francisco, CA Genentech has an opening for a fulltime position in its nonclinical biostatistics department, located in South San Francisco, CA. The nonclinical group currently consists of seven MS and PhD

None of Andy's comments) are inconsistent with the point that rlm() and lqs(), if they disagree with lm(), likely offer better places to start than does lm(), in identifying points that should be examined as in some sense outliers. All such methods are to be used, not as crutches, but as sources of insight. Incidentally, the "rlm" class inherits from "lm", and plot.lm()

Thanks Bill. Do you and others think that a link to this guide (or another)should be included in the Posting Guide and/or R FAQ? -- Bert On Tue, May 17, 2011 at 4:07 PM, <Bill.Venables at csiro.au> wrote: > Amen to all of that, Bert. ?Nicely put. ?The google style guide (not perfect, but a thoughtful contribution on these kinds of issues, has avoiding attach() as its very first line.

R 1.7.1 on Windows XP The "normal Q-Q plot" produced by plot.lm() mislabels points when the model is fitted using na.action=na.exclude. Example: x <- 1:50 y <- x + rnorm(50) y[c(5,10,15)] <- NA # insert some NA's y[40] <- 50 # add an outlier plot(lm(y ~ x, na.action=na.omit)) # outlier correctly labeled in all # four plots

Dear R-help, Can some one tell me what could be the problem? I downloaded the R-1.7.1 source and try to compile it on WinXP. "make", "make bitmapdll", "make tcl", "make recommended" all ran w/o problem. However, when I ran "make docs", it failed at creating refman.pdf, with no apparent error that I could see. I have attached the log file

This ought to work: resultdt <- lapply(PGWide[, 240:389], function(x, ...) try(tdt(x, ...))) You can then check the class of each component to see which one failed. Andy From: Farrel Buchinsky > > "Berton Gunter" <gunter.berton at gene.com> wrote in message > news:008601c67097$de1b46e0$5bc4fea9 at gne.windows.gene.com... > > ?try > > > > as in

Hi people, Does anybody know some Density-Based Method for clustering implemented in R? Thanks, Fernando Prass _______________________________________________________

Hi, I am running Samba 3.6.19 and making my way through all the Active Directory documentation. I have joined various Linux boxes to a Windows Server 2008 R2 Active Directory domain. Windows 7 users who log on with their domain accounts can mount the correct shares exported by the Linux boxes without having to supply any additional credentials. All working perfectly! However, I have a

> From: Thomas Lumley [mailto:tlumley at u.washington.edu] > > Please don't send these messages to r-announce. r-announce is for > announcements. Isn't r-announce supposed to be moderated? I posted to r-announce about a minor update in randomForest and was rejected. I was told that it's meant for new or significant upgrades in packages, etc. Andy

Dear R-help, A colleague of mine was running some code on two of our boxes, and noticed a rather large difference in running time. We've so far isolated the problem to the difference between R 1.7.1 and 1.8.0, but not more than that. The exact same code took 933.5 seconds in 1.7.1, and 3594.4 seconds in 1.8.1, on the same box. Basically, the code calls boot() to bootstrap fitting mixture

Dear All, When I run rlm to obtain robust standard errors, my output does not include p-values. Is there any reason p-values should not be used in this case? Is there an argument I could use in rlm so that the output does include p-values? Thanks in advance, Celso [[alternative HTML version deleted]]

This is my first post requesting help to this mailing list. I am new to R. My apologies for any breach in posting etiquette. I am new to this language and just learning my way around. I am attempting to run some sample code and and am confused by the error message: Loading required package: rrcov Error in fitNorm2(fdat[, "FSC-H"], fdat[, "SSC-H"], scalefac = ScaleFactor) :

Does anyone have any ideas how I could do a power calculation for a log rank test. I would like to know what the suggested sample sizes would be to pick a difference when the control to active are in a ratio of 80% to 20%. Thanks Dan -- ************************************************************** Daniel Brewer, Ph.D. Institute of Cancer Research Email: daniel.brewer at icr.ac.uk

I meant to send this to R-help, but it went to Gabor instead... This is my shot at it using diff(), assuming breaks is in non-decreasing order: chop <- function(x, breaks) { s <- sign(x - breaks) if (x <= breaks[1]) { ans <- x } else { numPieces <- sum(s == 1) ans <- diff(c(0, breaks[1:numPieces], x)) } ans } Some tests: >

Hi All Please see the Rprofile file which i have modified as follows and after that when I start R then I see that R says to me "TRUE" for all the packages implying that all loaded at once. But when i try to use commands as simple as help("lm"), it doesnt work nor any of the menu "Packages" is not working. Although the regression using lm ( Y ~ X ) is working

-----Original Message----- From: Ghosh, Sandeep Sent: Thursday, June 30, 2005 5:43 PM To: 'Berton Gunter' Subject: plot legend outside the grid Thanks for the pointers... I managed to get everything to look and feel the way I want except for the legend to plot outside the grid... Thanks for the note on the par, but I'm not able to it to plot outside the plot grid.. dataFrame <-

> > Leaf Sun wrote: > > The histogram is highly screwed to the right, say, the range > > of the vector is [0, 2], but 95% of the value is squeezed in > > the interval (0.01, 0.2). I guess the histogram is as you wrote. See http://web.maths.unsw.edu.au/~tduong/seminars/intro2kde/ for a short explanation. > -----Original Message----- > From: Berton Gunter

Peter You question is difficult to answer without more information about the distribution of your residuals. Different residual patterns call for different transformations to stabilize the variance. One very common form of heterocedasticity is increasing variance with increasing values of an independent predictor, i.e. the variance of the residuals of y=x increase as x increases. In this case a