search for: abunn

Hi, How can I plot two scatterplots on the same, one panel? I have two times series (price data for two goods) for the same period. Many thanks, -R [[alternative HTML version deleted]]

Hi, I've been doing a lot of CA modeling lately and am now wanting to make some checkerboard plots in R. Let's say I have a matrix: > is.matrix(junk) [1] TRUE > junk [,1] [,2] [,3] [,4] [,5] [1,] 0 0 1 0 0 [2,] 0 1 1 1 0 [3,] 0 1 0 0 1 [4,] 0 1 1 1 1 [5,] 0 1 0 0 0 > and I want to make a

Martin Maechler wrote: > More generally, please, please, everyone : > > Replace > if (something == TRUE) > with if (something) > and > if (something.or.other == FALSE) > with if (!something.or.other) Amen! And right on! And you tell 'em! And a few other exclamations, several of which are not suitable for a family program. That sort of syntax,

Hi, I have been trying for a while to use cbind and rbind to add a row and column to the same table but seem to be able only to add one OR the other??? Any help would be most welcome, Tony Evans Australia [[alternative HTML version deleted]]

Thank you all for your help :-) Much Much Aprreciated........ I tried with the "jitter" command and it worked :-) Cheers, Latha > -----Original Message----- > From: Latha Raja > Sent: Monday, February 07, 2005 3:21 PM > To: 'r-help at stat.math.ethz.ch' > Subject: Need your help with my R plot > Importance: High > > Hi, > > I am using R to

Just want to offer my congratulations to the JGR developers as the recepient of the 2005 Chambers Award. Great job, guys!! http://stats.math.uni-augsburg.de/JGR/ [Now, could JGR be updated to work with 2.1.0 (or be made R version independent, please... 8-)] Best, Andy

Helpers. An instrument sends me data that is mostly nonlinear. I have a group of functions to manipulate this data so that it is useful to the user. One of them involves a nls model that called to fit a line to the data and returns the fits. This works well 99 out of 100 times. Occasionally, the sensor sends some bad data to which a nls model cannot be fit. When that happens, the nls function

Although this should be trivial, I'm having a spot of trouble. I want to make a lattice plot of the format y1+y2 ~ x | grp but then fit a lm to each y variable and add an abline of those models in different colors. If the xyplot followed y~x|grp I would write a panel function as below, but I'm unsure of how to do that with y1 and y2 without reshaping the data before hand. Thoughts

OK. Another amateur question. I have a list with attributes on pine trees, like the stem's location, a logical value set to T if it's alive, some parameters for growth, diameter, etc. The tree list has another list in it which is a new data type for me. I want to make a new list that retains all the live trees. That is where Living == T. Here's the summary of the list: >

With 2.1 on Windows XP SP2 (32 bit) I also get no title in a png plot, so I can reproduce this bug. R is installed from pre-built binary. No title appears when I run this: png("foo.png") plot(1:10, main = "foo") dev.off() The jpeg device behaves as expected: jpeg("foo.jpg") plot(1:10, main = "foo") dev.off() Thoughts on this?

Suppose I have a list where I want to extract only the elements that occur in every component. For instance in the list foo I want to know that the numbers 2 and 3 occur in every component. The solution I have seems unnecessarily clunky. TIA, Andy foo <- list(x = 1:10, y=2:11, z=1:3) bar <-unlist(foo) bartab <- table(bar) as.numeric(names(bartab)[bartab==length(foo)])

Hi, I have a matrix not unlike this: foo <- matrix(,5,5) foo[5,1] <- 1 foo[1:3,2] <- 1 foo[3:4,3] <- 1 foo[4:5,4] <- 1 foo[2:4,5] <- 1 foo [,1] [,2] [,3] [,4] [,5] [1,] NA 1 NA NA NA [2,] NA 1 NA NA 1 [3,] NA 1 1 NA 1 [4,] NA NA 1 1 1 [5,] 1 NA NA 1 NA I want to get a vector that is the column numbers as sorted

I've ported somebody else's rather cumbersome Matlab model to R for colleagues that want a free build of the model with the same type of I/O. The Matlab model reads a text file with the initial parameters specified as: C:\Data\Carluc\Rport>more Params.R # Number of years to simulate nYears = 50; # Initial year for graphing purposes year0 = 1970; # NPP/GPP ratio (cpp0 unitless) fnr =

Hi all, I'm porting a library I wrote in S-plus over to R. All but one of the functions work fine. The only one that doesn't is really giving me a headache. I'm admitting defeat. Can anybody help? Object example3 has been read in using read.table. The statement works fine in S-Plus. When I run the same statement in R the column Prefix, which is read in using the as.character

Hi there - I have been trying to generate some simple stats and save the results to a file. My data looks like this: x y z exp 0 3 5 1 2 11 10 1 4 4 5 1 7 6 4 1 11 1 2 2 5 7 1 2 3 3

All Is there an R command to produce error bars around a plotted point. Crawley's book uses the command error.bar() but my version of R rejects it. Must be an S+ command(?). Thanks REX [[alternative HTML version deleted]]

Dear R list: I modelled changes in a variable (mconc) over time (d) for individuals (replicate) given one of three treatments (treatment) using: mconc.lme <- lme(mconc~treatment*poly(d,2), random=~poly(d,2)|replicate, data=my.data) summary(mconc.lme) shows that the linear coefficient of one of the treatments is significantly different to zero, viz. Value Std.Error

Hi all (really probably just Deepayan): In the plot below I want to add text on either side of each violin plot that indicates the number of observations that are either positive or negative. I'm trying to do this with ltext() and I've also monkeyed about with panel.text(). The code below is generally what I want but my calls to ltext() are wrong and I'm not sure how to fix them.

Hello, Is there anyone who could give me an example of how to import an ascii grid (i.e. ArcGIS exported raster) into R. I want to use it with gstat but don't know the appropriate import routine. Thanks very much for your help. Regards, femke Femke Reitsma Graduate Student (ABD) Geography Department 2181 LeFrak Hall University of Maryland College Park, MD 20742 Phone: 301-405-4121

...help at stat.math.ethz.ch mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide! > http://www.R-project.org/posting-guide.html > ------------------------------ Message: 14 Date: Wed, 10 Nov 2004 09:39:09 -0500 From: "Andy Bunn" <abunn at whrc.org> Subject: RE: [R] Does something like partition.rpart() exist? To: "Simone Vantini" <vantini at mate.polimi.it>, <R-help at stat.math.ethz.ch> Message-ID: <NEBBIPHDAMMOKDKPOFFIKEOGCLAA.abunn at whrc.org> Content-Type: text/plain; charset="iso-8859-1&...