search for: berton

What would be nice is to have something like a "robust" task view... Andy From: Berton Gunter > > Many thanks for this Martin. There now are several packages > with what appear to be overlapping functions (or at least > algorithms). Besides those you mentioned, "robust" and > "roblm" are at least two others. Any recommendations about > how o...

Folks: This is a question about the S4 formal class system. Suppose I have a class, 'foo', defined by: setClass('foo',representation(dat='matrix', id='character') ) I wish to define a '[' method for foo that will extract from the 'dat' slot. I would have thought that the following would work, but it doesn't:

How do I manipulate the read.table function to read in only the 2nd column??? [[alternative HTML version deleted]]

Hello R-helpers I want to put a color scale in a plot: I've got an xy plot where the values of the response (z=f(x,y)) is symbolically given by colors (like heat or rainbow color scale) I would like to put such a scale with apprpriates labels in the plot, so as to facilitate the interpretation (like in a finite elements result plot) How is taht possible? Thanks Anne [[alternative

Hi all, I tried to do a variance components using nlme, but I got the following warning mesage ##################################################### not meaningful for factors in: Ops.factor(y[revOrder], Fitted) ###################################################### Can someone point out what is the meaning of this warning message? I tried to look at Ops.factor, but I don't

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?

Dear All, would you please provide your thoughts on the following: let us say I have: a <-c(1,5,8,15,32,69) b <-c(8.5,33) and I would like to extract from "a" the two values that are closest to the values in "b", where the length of this vectors may change but b will allways be shorter than "a". So at the end based on this example I should have the result

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

Press "ESC" is supposed to be able to stop the current computation... but often times it does not really work... I had to kill the R process to stop and restart the session... What is a good way to stop/cancel the current computation? [[alternative HTML version deleted]]

Hi, You could try: A<- matrix(unlist(read.table(text=" 1 2 3 4 5 6 7 8 9 9 8 7 6 5 4 3 2 1 ",sep="",header=FALSE)),ncol=3,byrow=FALSE,dimnames=NULL) library(matrixStats) ?res1<-t(sapply(split(as.data.frame(A),as.numeric(gl(nrow(A),2,6))),colProds)) ?res1 #? [,1] [,2] [,3] #1??? 4?? 10?? 18 #2?? 63?? 64?? 63 #3?? 18?? 10??? 4

...=", homCount, ")"); par(xpd=TRUE); legend(3.8, max(as.vector(dataFrame$marbles_buried)), c(wtLegend, hetLegend, homLegend), col=c('blue', 'red', 'green'), pch=c(1,4,2)); Again thanks so much for all the pointers.. - Sandeep -----Original Message----- From: Berton Gunter [mailto:gunter.berton at gene.com] Sent: Friday, July 01, 2005 7:37 PM To: Ghosh, Sandeep Subject: RE: plot legend outside the grid You are not specifying the x,y position of the legend correctly. You do not need to end each statement with a ";"; you should set par(xpd=TRUE) be...

Hi All I received a raw data set with one record per tennis player (both male and female) and then i cured it by aggregation i.e by 4 age groups, 2 gender levels and 6 income levels. Gender and Income are categorical variables. Please advise me how to use 'R' to model this data set (Actually, i want to know the right regression technique and steps to do that, including removing

Hi all: Newbies (and others!) may find useful the R Reference Card made available by Tom Short and Rpad at http://www.rpad.org/Rpad/Rpad-refcard.pdf or through the "Contributed" link on CRAN (where some other reference cards are also linked). It categorizes and organizes a bunch of R's basic, most used functions so that they can be easily found. For example, paste() is

Can the packages that I've installed in R 2.2.1 be automatically imported to R 2.3.0? Otherwise it is a hassle everytime I update a new version of R... [[alternative HTML version deleted]]

Dear R-Users, how to include error bars within lattice? How should the panel = function(x,y,...){ looks like? Does panel.arrows works here as well? I appreciate any help on this. Regards, Mario AT [[alternative HTML version deleted]]

...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 [mailto:gunter.berton at gene.com] > Sent: Wednesday, November 02, 2005 1:10 PM > To: 'Leaf Sun'; r-help at stat.math.ethz.ch > Subject: [R] Visualizing a Data Distribution -- Was: breaks in hist() > > > Leaf: > > An interesting question concerning graphical p...

I have a Lattice plot in which I want to adjust the number of tick marks used, and I want to have the drawn grid reflect that change. Here is what I'm doing: bwplot(var1 ~ var2, data=df, scales=list(tick.number=10), panel=function(...) { panel.grid(h=0,v=-1,...); panel.stripplot(col="gray40", pch="|", cex=2, ...); panel.bwplot(...); })

Hello, I have data in a frame: 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 I would like them arranged in a single column: 1 2 3 4 . . 8 9 . . 16 etc. I believe this should be possible using reshape, but I can't see how to do it. Thanks for sparing time to help a neophyte, Roger Mason

>On 10/14/2005 1:08 PM, Berton Gunter wrote: >> ?setwd >> e.g. setwd(file.choose()) >> >> BTW, you could have found this on your own via help.search('working >> directory') . Base R has quite good docs -- you should try them first. >> >> -- Bert Gunter > >That won't w...

...ence University of Maryland School of Medicine Division of Gerontology Baltimore VA Medical Center 10 North Greene Street GRECC (BT/18/GR) Baltimore, MD 21201-1524 (Phone) 410-605-7119 (Fax) 410-605-7913 (Please call phone number above prior to faxing) jsorkin at grecc.umaryland.edu >>> Berton Gunter <gunter.berton at gene.com> 8/3/2006 11:56:28 AM >>> I know I'm coming late to this, but ... > > Is someone able to suggest to me a transformation to overcome the > > problem of heterocedasticity? It is not usually useful to worry about this. In my experience...