search for: whrc

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

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

...qrt(tau.sp) tau ~ dunif(0,100) sigma.sppl<-1/sqrt(tau) #sd of species x plot random effect } ----------------------------------------- John Poulsen, PhD Assistant Scientist Woods Hole Research Center 149 Woods Hole Road Falmouth, MA 02540-1644 USA Tel: 508.540.9900 x168 www.whrc.org [[alternative HTML version deleted]]

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?

Hello, I am trying to make a xyplot plot with points that are different symbols. I want to call the symbol type (pch) from a column in my dataframe. Here is a simplified example. In my real example I also have groups, which I have not included here. This example doesn't change the symbols or colors. Any help you can provide would be appreciated. Thanks, John x<-c(1:12)

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)])

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 =

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.

...tat.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" I...

This seems like a FAQ, but I can't figure it out. I have a mv ts object: R > tsp(pg) [1] 1982 2003 1 R > dim(pg) [1] 22 12 and a univariate ts: R > tsp(rw) [1] 1690 1996 1 Yet, when I try to intersect them: R > tsp(ts.intersect(rw, pg)) [1] 1982 2176 1 the process goes awry. How to I get rw and pg to be one ts that runs from 1982 to 1996 and has 13 univariate time

How can I specify the repositories for upgrade()? > x <- packageStatus(repositories = "http://cran.us.r-project.org//bin/windows/contrib/2.0") > upgrade(x, ask = FALSE) Error in update[, 3] : incorrect number of dimensions x, the object of class "packageStatus", prints and summarizes fine. I also ran install.packages() without incident. >

I have a very simple query with regard to summarizing the number of factors present in a certain snippet of a data frame. Given the following data frame: foo <- data.frame(yr = c(rep(1998,4), rep(1999,4), rep(2000,2)), div = factor(c(rep(NA,4),"A","B","C","D","A","C")), org = factor(c(1:4,1:4,1,2))) I want to get

Hi, I would like to skip a preset number of values when reading in a binary file (readBin). Is there any way to do this? kind regards, and thanks in advance for the help. Pieter Beck [[alternative HTML version deleted]]

Hello, I want to build a function to call up a column of a data.frame by the names of the columns. I have column names that are sequentially named (col1, col2, etc.). How do I change a character expression into something that will be understood as a data.frame column. For example: example<-data.frame(cbind(col1=1:10, col2=21:30, col3=41:50)) call.fun<-function(t){

Preamble: Sorry for being dense. Now that that's done, here are my questions. I want to put a polygon on a plot of a time series. I'm going to add lines from a smooth.spline interpolation and other annotation to it. But here's the general idea: # Start code n <- 121 dat <- rnorm(n) plot(dat, type="n") polygon(c(1:n,n:1), c(dat,c(rep(0,n))), col =

Hi, Somebody asked me to make sure that all the machines running the in our lab (XP and Linux, both running 2.0) have R installed and that A) All the packages are installed and B) kept up-to-date. Obediently, I began to modify a shared Rprofile so that once a week it checks for new packages and updates to the current version of the installed packages on CRAN. Sounds simple enough. Plus some