similar to: plot.formula

Could I suggest the following extension to plot.formula: plot(cbind(y1,y2) ~ x, ...) should plot (y1 against x) and (y2 against x) on the same plot. The default y axis limits would be determined by the range of c(y1,y2). This would be pretty handy sometimes, replacing 4 lines of code. The current plot.formula evaluates cbind(y1,y2), which is a matrix, so plot.formula looks for

Peter Dalgaard <p.dalgaard at biostat.ku.dk> writes: > > Barry Rowlingson <B.Rowlingson at lancaster.ac.uk> writes: > > > Liaw, Andy wrote: > > > Stupid me: fell into this trap: > > > > > >>0 == 0 == 0 > > > [1] FALSE > > > > > > > Ouch! > > > > Python's comparison operators don't have this

From: Martin Maechler <maechler at stat.math.ethz.ch> MIME-Version: 1.0 Content-Type: text/plain; charset=us-ascii Content-Transfer-Encoding: 7bit Message-ID: <17739.46076.735981.117358 at stat.math.ethz.ch> Date: Fri, 3 Nov 2006 22:26:20 +0100 To: Barry Rowlingson <B.Rowlingson at lancaster.ac.uk> Cc: Sarah Goslee <sarah.goslee at gmail.com>, r-help at

Full_Name: Albrecht Gebhardt Version: 0.90.0 OS: Linux + Tru64 Submission from: (NULL) (143.205.61.72) This is where example(krige) from library(sgeostat) stops with R 0.90.0: krige> image(grid$x, grid$y, grid$krige$zhat, add = T) Error in length(x) == nrow(z) : comparison (1) is possible only for vector types in previous versions of R the above image() call worked, because image() did the

Full_Name: Barry Rowlingson Version: 1.6.0 OS: RH8 i386 Submission from: (NULL) (148.88.136.205) Strange things happen if I premultiply a difftime() object with a number. Example: > d1 <- difftime(Sys.time(),Sys.time()) > d2 <- 1 * difftime(Sys.time(),Sys.time()) > d3 <- difftime(Sys.time(),Sys.time()) * 1 > d1 Time difference of 0 secs - thats fine > d2 [1] 0

gReetings, I want to manipulate a formula object, containing the name "." so that "." is replaced by a desired (arbitrary) expression. What is a safe way to do this? Adrian Baddeley

Barry, Would you mind providing the necessary patch for this behavior? Thanks, Greg > -----Original Message----- > From: Barry Rowlingson [mailto:B.Rowlingson@lancaster.ac.uk] > Sent: Monday, November 04, 2002 8:04 AM > To: r-devel@stat.math.ethz.ch > Subject: Re: longjmp - was: seemingly random "nesting of > readline input" > warnings > > > Luke

I see some value in Duncan?s proposal to implement this as an extra package instead of a change to base R, if only to see if the idea has legs. I?m minded to do so myself using your suggestion, but is there a particular reason why you recommend using the remotes package instead of devtools? The latter seems to have the same functions I would need, and I believe it is more widely installed that

Hi, Sorry if this is already known... contour() sometimes draws contour lines that intersect. Is there a temporary fix? A dataset which causes problems is at http://www.maths.uwa.edu.au/~adrian/dumpdata.R If you try just source("dumpdata.R") image(huh) contour(huh) the 100 x 100 matrix 'huh' contains an hourglass-shaped region of values around 0.8. The contour plot

Hi, I am fitting a weighted least square regression and trying to compute SSE,SST and SSReg but I am not getting SST = SSReg + SSE and I dont know what I am coding wrong. Can you help please? xnam <-colnames(X) # colnames Design Matrix fmla1 <- as.formula(paste("Y ~",paste(xnam, collapse=

Wise and merciful R-helpers: I want to equip a data frame with an attribute which specifies how to plot some of the columns. Up to now we have been doing this by giving the data frame a `formula' attribute, that can be passed to plot.formula. For example dat <- data.frame(x=1:100,y=runif(100),z=100:1) attr(dat, "plotme") <- (z ~ x) ...... ......

Hello, I'd like to take all possible sub-summands of a vector in the quickest and most efficient way possible. By "sub-summands" I mean for each sub-vector, take its sum. Which is to say: if I had the vector x<-1:4 I'd want the "sum" of x[1], x[2], etc. And then the sum of x[1:2], x[2:3], etc. And then...so on. The result would be: 1 2 3 4 2 5 7 6 9 10 I can

I am working on a two stage threshold quantile regression model in R, and my aim is to estimate the threshold of the reduced-form equation (call it rhohat), and the threshold of the structural equation (call it qhat), in two stages. On the first stage, i estimate rhohat by quantile regression and obtain the fitted values. I use these fitted values to estimate qhat on the second stage. The code is

Hi - the list of mirror sites given at www.r-project.org/ also mentions http://mirror.au.r-project.org/pub/CRAN The correct URL seems to be http://cran.planetmirror.com/ regards Adrian Baddeley -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-devel mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info",

On Mon, Mar 8, 2010 at 8:46 PM, Barry Rowlingson <b.rowlingson at lancaster.ac.uk> wrote: > On Mon, Mar 8, 2010 at 6:44 PM, Ma Ismail - NewYork-MEAG-NY > <ima at meag-ny.com> wrote: >> Hi, >> >> A few of the developers on our Quant team are using R for data calculation and to generate a [snip] > ?I've noticed a lot of financial corporates getting into R

Hello, All: The vignette with the sos package used "upquote.sty", required for R Journal when it was published in 2009. Current CRAN policy disallows "upquote.sty", and I've so far not found a way to pass "R CMD check" with sos without upquote.sty. I changed sos.Rnw per an email exchange with Prof. Ripley without solving the problem; see below. The

Dear all, Is there a reason, why there is no +-operator for aes (i.e. uneval) objects (as there is for themes and gg objects)? I had a couple of cases where such an operator would be useful, for instance to combine the result of aes and aes_string in functions. Any flaws with the following proposition: `+.uneval` <- function(e1, e2) { dup <- names(e1) %in% names(e2) if (any(dup)) {

Hi Peter, The ssconvert tool (part of gnumeric) is very good at converting spreadsheets to csv-files. There is a wrapper in the "gnumeric" package on cran. Cheers, Thomas > Date: Fri, 4 Oct 2013 09:08:50 +0100 > From: Barry Rowlingson <b.rowlingson at lancaster.ac.uk> > To: Peter Maclean <pmaclean2011 at yahoo.com> > Cc: "r-help at r-project.org"

Thank you for interesting examples. I would find useful to document this behavior also in `?mean`, while `+` operator is also affected, the `sum` function is not. For mean, NA / NaN could be handled in loop in summary.c. I assume that performance penalty of fix is the reason why this inconsistency still exists. Jan On Mon, Jul 2, 2018 at 8:28 PM, Barry Rowlingson < b.rowlingson at

On Thu, Jul 26, 2018 at 1:43 PM <luke-tierney at uiowa.edu> wrote: > > On Thu, 26 Jul 2018, G?bor Cs?rdi wrote: > > > On Thu, Jul 26, 2018 at 12:25 PM Barry Rowlingson > > <b.rowlingson at lancaster.ac.uk> wrote: > >> > >> On Thu, Jul 26, 2018 at 12:22 AM, G?bor Cs?rdi <csardi.gabor at gmail.com> wrote: > >>> I am trying to