similar to: Multiple plots per window

Warranty on Accuracy, Precision, Legality, ... of R in Research (These questions may well have been raised.) What is the implied warranty of using R for research & publications, consulting, etc.? Alternately, how does one obtain such a warranty? Your answers will be much appreciated. Perhaps you can point me to some websites which discussed this subject in the past. Thanks & regards

I'm trying to draw an arrow with a curved shaft on the graph as a straight one looks messy on a detailed graph. I've looked in arrows but it doesn't seem to give an option. larrows doesn't look much more promising. I had a look in the archive and couldn't find anything. Any thoughts? Thanks Paul -- View this message in context:

This was my initial attempt at creating a title on a graph of the R squared value: x<-rnorm(10) y<-rnorm(10) plot(x,y, main=paste(expression(R^2)," = ",round(summary(lm(y~ x))$r.squared, digits=3), sep="")) I've read various other posts that say expression needs to be taken outside the paste, but I can't seem to get it work as the following fails plot(x,y,

> An orphaned package? anyone in Switzerland know if there's an alternative? Note the email. I guess CRAN-R should be notified. regards, Bob C > The original message was received at Fri, 6 Apr 2007 16:32:57 -0700 > from vayu.arc.nasa.gov [143.232.122.22] > > ----- The following addresses had permanent fatal errors ----- > <juerg.schmidli at env.ethz.ch> >

A simple question - I have a small dataset to read in and want to copy and paste part from Excel and paste it into an R script file without creating more files saving it as a .txt/.csv and then reading that in. I want to read in 3 columns e.g. 1 2.5 3.4 1 2.3 3.1 1 2.6 3.9 2 2.9 2.8 2 2.6 2.9 2 2.7 2.9 3 2.3 3.3 3 2.4 3.0 3 2.7 3.2 I thought I could use scan() but don't know how to extend it

Hi - I'm trying to avoid using a 'for' loop due to inefficiency and instead use a function (and ultimately tapply as I'm working on a matrix) but I can't figure out how to get 'function' to take the variables as anything other than vectors for example aa<-0 x<-1:4 test.fun<-function(x) {aa<-(x*x +aa) return(aa)} test.fun(1:4) This code

Hi, The default value for lag.max in ACF implementation is 10*log10(N) There several publications recommending setting lag.max to: - N/4 (Box and Jenkins, 1970; Chatfield, 1975; Anderson, 1976; Pankratz, 1983; Davis, 1986; etc.) - sqrt(N)+10 (Cryer, 1986) - 20<=N<=40 (Brockwell and Davis) Why R uses 10*log10(N) as a default? Please, give me a reference to a book or article where the

I'm trying to assign NAs to values that satisfy certain conditions (more complex than shown below) and it gives the right result, but breaks the loop having done the first one viz: new<-c(rep(5,4),6) for (i in 1:6) {new[new[i]>5.5][i]<-NA} gives the correct result, though an error message appears which causes a break if it's in a loop. If I can get rid of the error message and

Simple problem - I want the ylab to automatically pick up x1 rather than having to define x1 in the plot statement. x1<-c(1.2,2,3);x2<-c(1,2.1,2.6) y<-x1 plot(1:3,y, ylab="x1") There must be a way of accessing the name x1 somehow, but unfortunately I don't know how to search for it. Any help would be great, -- View this message in context:

Am trying to produce a graph which prints out well in black and white using ggplot2. I have the following example set up nicely, but want to shade the red bars in one pattern and the blue in another so they print out clearly. I tried changing colours to 1 light, 1 dark, but then the overlapping colour looks virtually identical to the darker one. I noted the option density in barplot, but

If you redefine your NAs as below to be detected as some arbitrary large number, then the code should work through. Any 5's left in your dataset can be replaced just as easily by NAs again. Not elegant, but effective. site <- c("s1", "s1", "s1", "s2","s2", "s2") pref <- c(1, 2, 3, 1, 2, 3) R1 <- c(NA, NA, 1, NA,NA,NA) R2

The pursuit projection packages ppr is an excellent contribution to R. It is great for one-to-three ridge fits, often somewhat intuitive, and for multi-ridge fits, where it at least describes a lot of variance. Like many folk, I need to report the fits obtained from ppr to the greater, outside, non-R world. It is fairly obvious how to use the terms alpha and beta to report on

Hi - I can access R from Tinn-R by going to Options->Main->Application/R and setting the search path, but each time I exit Tinn-R I have to redefine the search path. Is there no way of fixing that directory as default? I have installed R under its default directory C:/Program Files/R/R-2.7.1 and Tinn under a variety of different places to try to rectify the problem though currently under

Dear all, given two numeric vectors x and y, the ACF(x) at lag k is cor(x(t),x(t+k)) while the CCF(x,y) at lag k is cor(x(t),y(t-k)). See below for a simple example. > set.seed(1) > x <- rnorm(10) > y <- rnorm(10) > x [1] -0.6264538 0.1836433 -0.8356286 1.5952808 0.3295078 -0.8204684 0.4874291 0.7383247 0.5757814 -0.3053884 > y [1] 1.51178117 0.38984324

I'm interested in plotting a y with an x factor as the combination of 2 factors and colour with respect to a third, which the code below does with interaction.plot(). However, this is because I redefine the x to be 1 factor. Is there a way of getting it to plot without redefining it, and ideally to not join up the lines BETWEEN levels a and b, but just join those between after and before for

Hello, what are they meaning? It could be something that would show a threshhold above which the result is indicating different meanings then just random noise. But there is no description on the definition of those lines, so it means nothing, if it is not clearly defined. Where can I find a detailed definiton? Ciao, Oliver

I'm trying to plot a graph where the axes go through 0,0, rather than around it combined with a box round the graph, so x<-0:10;y<-0:10 plot(x,y) gives me a box but doesn't go through the point 0,0, but stays at a distance. In trying to circumvent this problem, I wrote plot(x,y) axis(1, pos=c(0,0));axis(2,pos=c(0,0)) which gives me axes that go through the origin, but then I

Hello, I've come up with some challenges with my process that are a bit too complicated for the mailing list. Is there anyone out there, preferably a real "statistician", who is willing to consult with me via phone/email for a few hours. I'm happy to pay you for your time. Thanks, -Noah

Hi - I am trying to substitute for "the_other_y" in the code below. I want y2 and y3 to be there when i is 1, y1 and y3 to be there when i is 2 and y1 and y2 to be there when i is 3. I'm sure it's to do with what format the data should be in and I've tried alldata[,-i], but it fits all the columns of alldata except i rather than each column one at a time. I've tried

Hi, Are there some texts about time series that teatches it with R? And site 'bout time series "for dummies" , can anyone indicate them for me? thanks a lot Rafael Bertola Undergraduate in Statistics at UNICAMP - Brasil -------------------------------------------------- --------------------------------------------------