Re: [R] graphing plots of plots

2010-08-21 Thread Bernard Leemon
I've now tried Babtiste's code and my reaction is WOW! it shows me how to do just what I need to do. I know enough to follow all the code but it would have taken me a LOOO time to generate it. thank you Babtiste! gary On Sat, Aug 21, 2010 at 2:15 PM, baptiste auguie < baptiste.aug...@g

Re: [R] graphing plots of plots

2010-08-21 Thread Bernard Leemon
many useful suggestions that I'll work on, especially babtiste's detailed code. yes, I want something like Fig 1.7, or 7.18, or 7.22, but where the x,y values are characteristics of the mini-histogram that is plotted. attached (if it makes it through) is what i'm trying to do in R. On Sat, Aug

[R] graphing plots of plots

2010-08-21 Thread Bernard Leemon
I want to make a graph where each element plotted is itself a graph. I can see how to use par(fig=) and viewport to do that, but they require (i think) me to do my own scaling as they are scaled to the graphics window. any advice on which approach I should take (just bite the bullet and do my own

[R] nls() newbie convergence problem

2008-06-04 Thread Bernard Leemon
I'm sure this must be a nls() newbie question, but I'm stumped. I'm trying to do the example from Draper and Yang (1997). They give this snippet of S-Plus code: Specify the weight function: weight < - function(y,x1,x2,b0,b1,b2) { pred <- b0+b1*x1 + b2*x2 parms <- abs(b1*b2)^(1/3) (y-pred)/parms

[R] nls() newbie convergence question

2008-06-04 Thread Bernard Leemon
I'm sure this must be a nls() newbie question, but I'm stumped. I'm trying to do the example from Draper and Yang (1997). They give this snippet of S-Plus code: Specify the weight function: weight < - function(y,x1,x2,b0,b1,b2) { pred <- b0+b1*x1 + b2*x2 parms <- abs(b1*b2)^(1/3) (y-pred)/parms

Re: [R] Using R in a university course: dealing with proposal comments

2008-02-11 Thread Bernard Leemon
Hi Arin, Others have commented wisely an your first issue. As for your 2nd issue, I had my own concerns about using R in undergraduate teaching because I had always used a point-and-click program for that level. I should not have worried. The current generation has been typing on their keyboards

Re: [R] R programming style

2008-02-11 Thread Bernard Leemon
I just got a copy of A First Course in Statistical Programming with R by W. John Braun and Duncan J. Murdoch. Cambridge. at amazon: http://www.amazon.com/First-Course-Statistical-Programming-R/dp/0521694248/ first couple of chapters are base R that most everyone would know before wanting to pro

Re: [R] Do I need to use dropterm()??

2008-02-10 Thread Bernard Leemon
Hi Dani, it would be better to start with a question you are trying to ask of your data rather than trying to figure out what a particular function does. with your variables and model, even if the component terms were not significant, they must in the model or the product of sunlight and aspect wi

Re: [R] R on Mac PRO does anyone have experience with R on such a platform ?

2008-02-09 Thread Bernard Leemon
I have R on all sorts of Macs, including one's a lot whimpier than the one you are describing and it works great on all of them. gary mcclelland colorado On Sat, Feb 9, 2008 at 6:29 PM, Maura E Monville <[EMAIL PROTECTED]> wrote: > I saw there exists an R version for Mac/OS. > I'd like to hear fr

Re: [R] correlation

2008-02-08 Thread Bernard Leemon
It is easy to worry too much about using numbers to represent order when using statistics like the correlation. this little example shows that the correlation is essentially a rank-order correlation itself: > x <- 1:20 > y <- x^2 > cor(x,y) [1] 0.9713482 x and y are definitely not linearly relate

[R] a kinder view of Type III SS

2008-02-07 Thread Bernard Leemon
A young colleague (Matthew Keller) who is an ardent fan of R is teaching me much about R and discussions surrounding its use. He recently showed me some of the sometimes heated discussions about Type I and Type III errors that have taken place over the years on this listserve. I'm presumptive eno

Re: [R] how to calculate chisq value in R

2008-02-07 Thread Bernard Leemon
On Thu, Feb 7, 2008 at 1:06 PM, John Kane <[EMAIL PROTECTED]> wrote: > ?chisq.test > --- jinjin <[EMAIL PROTECTED]> wrote: > > > > > for example, an expression such as chisq(df=1,ncp=0) > > ? > > perhaps pchisq(chisqvalue, df=1, ncp=0) is what you are looking for to evaluate the probability for