Re: [R] help change auto.key colors in xyplot

2009-02-13 Thread Deepayan Sarkar
On Fri, Feb 13, 2009 at 11:22 AM, Dieter Menne dieter.me...@menne-biomed.de wrote: Dimitri Liakhovitski ld7631 at gmail.com writes: the code below works just fine to produce a dotplot. However, I am not successful changing the color of the lines in the legend (auto.key). If I add col=..., it

Re: [R] help change auto.key colors in xyplot

2009-02-13 Thread hadley wickham
On Fri, Feb 13, 2009 at 1:22 PM, Dieter Menne dieter.me...@menne-biomed.de wrote: Dimitri Liakhovitski ld7631 at gmail.com writes: the code below works just fine to produce a dotplot. However, I am not successful changing the color of the lines in the legend (auto.key). If I add col=..., it

Re: [R] Optimizing Multiple Models...any suggestions?

2009-02-13 Thread Paul Heinrich Dietrich
Hi Max, Thanks for the suggestion, that's exactly what I was looking for. Thanks again. Paul -- View this message in context: http://www.nabble.com/Optimizing-Multiple-Models...any-suggestions--tp21979556p22006494.html Sent from the R help mailing list archive at Nabble.com.

[R] second axis title orientation

2009-02-13 Thread Jörg Groß
Hi, I have added a second y-axis via axis() to a plot. Then I tried to add an y axis title. But the axis() function seems to have no argument like ylab. and if I add a title with title() the text is centered at the first axis, not the second defined one. Is there a way to add an

Re: [R] second axis title orientation

2009-02-13 Thread Gabor Grothendieck
There is an example in: library(zoo) example(plot.zoo) See ?plot.zoo On Fri, Feb 13, 2009 at 6:19 PM, Jörg Groß jo...@licht-malerei.de wrote: Hi, I have added a second y-axis via axis() to a plot. Then I tried to add an y axis title. But the axis() function seems to have no argument like

[R] PCA functions

2009-02-13 Thread glenn
Hi All, would appreciate an answer on this if you have a moment; Is there a function (before I try and write it !) that allows the input of a covariance or correlation matrix to calculate PCA, rather than the actual data as in princomp() Regards Glenn [[alternative HTML

[R] question about tool collisions

2009-02-13 Thread Carl Witthoft
Recently I got introduced to two packages: {seewave} and {audio} . Turns out they both have a tool to call a system audio tool, and in both cases the name of the tool is play(). Naturally these two tools do slightly different things with different arguments. So, what should a user do, and

Re: [R] second axis title orientation

2009-02-13 Thread Alexandre Swarowsky
?mtext On Sat, 2009-02-14 at 00:19 +0100, Jörg Groß wrote: Hi, I have added a second y-axis via axis() to a plot. Then I tried to add an y axis title. But the axis() function seems to have no argument like ylab. and if I add a title with title() the text is centered at the first

Re: [R] question about tool collisions

2009-02-13 Thread Bert Gunter
See the warn.conflicts arguments of require() and library(). People thought about this a long time ago. -- Bert Gunter -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Carl Witthoft Sent: Friday, February 13, 2009 4:05 PM To:

Re: [R] Bootstrap or Wilcoxons' test?

2009-02-13 Thread Murray Cooper
Charlotta, I'm not sure what you mean when you say simple linear regression. From your description you have two groups of people, for which you recorded contaminant concentration. Thus, I would think you would do something like a t-test to compare the mean concentration level. Where does the

Re: [R] question about tool collisions

2009-02-13 Thread Marc Schwartz
on 02/13/2009 06:05 PM Carl Witthoft wrote: Recently I got introduced to two packages: {seewave} and {audio} . Turns out they both have a tool to call a system audio tool, and in both cases the name of the tool is play(). Naturally these two tools do slightly different things with

[R] Outlier Detection for timeseries

2009-02-13 Thread Pele
Hello R users, Can someone tell if there is a package in R that can do outlier detection that give outputs simiilar to what I got from SAS below. Many thanks in advance for any help! Outlier Details

Re: [R] Bootstrap or Wilcoxons' test?

2009-02-13 Thread David Winsemius
I must disagree with both this general characterization of the Wilcoxon test and with the specific example offered. First, we ought to spell the author's correctly and then clarify that it is the Wilcoxon rank-sum test that is being considered. Next, the WRS test is a test for differences

Re: [R] Bootstrap or Wilcoxons' test?

2009-02-13 Thread Murray Cooper
First of all, sorry for my typing mistakes. Second, the WRS test is most certainly not a test for unequal medians. Although under specified models it would be. Just as under specified models it can be a test for other measures of location. Perhaps I did not word my explanation correctly, but I

[R] several ifelse problems...

2009-02-13 Thread kathie
Dear R users, From the code below, I try to compute y value. (In fact, y looks like a trapezoid) -- x - seq(0,1,.01) y - ifelse(abs(x-.5)=0.3,0, ifelse(abs(w-.5)=0.4,-1, ifelse((0.1w

Re: [R] Bootstrap or Wilcoxons' test?

2009-02-13 Thread Daniel Malter
Hi Charlotta, to be more constructive toward your goal. If you bootstrap the regression when the regression is ill-specified, the bootstrap may not help you. Further, a test as difficult as a regression does not seem to be necessary in your case. A t-test if your dependent variable is

Re: [R] Outlier Detection for timeseries

2009-02-13 Thread stephen sefick
what are you looking for? I am not familiar with SAS as I am a poor scientist. I am not promising anything, but if you were to tell me what you wanted - method etc. I may know of a package or something that would work. Stephen On Fri, Feb 13, 2009 at 9:16 PM, Pele drdi...@yahoo.com wrote:

[R] superscript

2009-02-13 Thread David Douterlungne
Dear R-users. I'm struggeling to fix the superscript of a label of a figure axis. For some reason R doesn't recognize the hat symbol. plot(1,1,xlab=ligth intensity (PAR),ylab=expression(mass Pteridium rhizomes (gr/0.25m^2))) A very similiar scriptline does not give any problem

Re: [R] superscript

2009-02-13 Thread Daniel Moreira
Try: plot(1,1,xlab=ligth intensity (PAR),ylab=expression(mass Pteridium rhizomes (gr/0.25*m^2))) Daniel Moreira, MD Research Associate Duke University Medical Center DUMC 2626, MSRB-I Room 455 571 Research Drive Durham, North Carolina 27710 Telephone:

Re: [R] several ifelse problems...

2009-02-13 Thread Berwin A Turlach
Dear Kathie, On Fri, 13 Feb 2009 21:08:25 -0800 (PST) kathie kathryn.lord2...@gmail.com wrote: Dear R users, [snip] However, even though the results show that y=8 for x=0.11, when x=0.11, actual y value is -0.9. And, y=-0.8 for x=0.88. I cannot understand the above results. It may help

Re: [R] Bootstrap or Wilcoxons' test?

2009-02-13 Thread David Winsemius
The Wilcoxon rank sum test is not plain and simple a test equality of distributions. If it were such, it would be able to test for differences in variance when locations were similar. For that purpose it would, in point of fact, be useless. Compare these simple situations w.r.t. the WRS:

<    1   2