similar to: Symbolic substitution of expressions into expressions

I'd like to be able to modify axlab in (C) below so that 'Inf' is replaced by the infinity symbol. y <- rnorm(40) breaks <- c(-Inf, -1, 1, Inf) x <- cut(y, breaks=breaks) plot(unclass(x), y, xaxt="n", xlab="") ## A: The following gives the axis labels "(-Inf, 1]", etc. axis(1, at=1:3, labels=expression("(-Inf,-1]", "(-1,1]",

expression() accepts multiple expressions as arguments, thus: plot(1:2, 1:2) legend("topleft", expression(y == a * x^b, "where "* paste(y=="wood; ", x=="dbh"))) Is there a way to do this when values are to be substituted for a and b? i.e., the first element of the legend argument to legend() becomes,

The inst/doc directory of the DAAG package has 6 files coral551.spot, ... that are around 0.85 MB each. It would be useful to be able to zip then, but that as matters stand interferes with the use of the Sweave file that uses them to demonstrate input of expression array data that is in the "spot" format. They do not automatically get unzipped when required. I have checked that

In this connection, note the following > a4 <- 4 > plotThis <- bquote(alpha==.(a), list(a=a4)) > do.call(plot, list(1:10, main=do.call(expression, c(plotThis)))) > do.call(plot, list(1:10, main=do.call(expression, plotThis))) Error in do.call(expression, plotThis) : second argument must be a list > ## Whereas plotThis has class "call", c(plotThis) has class

I propose the following enhancements and changes to plot.lm(), the most important of which is the addition of a Residuals vs Leverage plot. (1) A residual versus leverage plot has been added, available by specifying which = 5, and not included as one of the default plots. Contours of Cook's distance are included, by default at values of 0.5 and 1.0. The labeled points, if any, are those

I am looking for information on experimentation with the use of R in the teaching of statistics and science in schools. Any leads would be very welcome. I am certain that there is such experimentation. I've made this inquiry on r-sig-teaching, with no response. John. John Maindonald email: john.maindonald at anu.edu.au phone : +61 2 (6125)3473 fax : +61 2(6125)5549 Centre

Can anyone suggest an equivalent, for lattice graphs, of the base graphics argument bty="l"? NB that I am leaving off the box around the strip, with a strip function: stripfun <- function(which.given,which.panel, factor.levels=as.expression(levlist), ...){ panel.text(x=0, y=0.5, lab = as.expression(levlist[which.panel[which.given]]),

In Section 1.4.2 of "Writing R Extensions" %\VignetteEngine{knitr::knitr} should be %\VignetteEngine{knitr::knit} > sessionInfo() R version 3.0.2 (2013-09-25) Platform: x86_64-apple-darwin10.8.0 (64-bit) Is this sort of thing best reported here, or is a huge report in order? John Maindonald email: john.maindonald at anu.edu.au phone : +61 2 (6125)3473 fax

<<< Thus factor analysis is in essence a model for the correlation matrix of x, ? = ?'? + ? >>> This should surely be ? = ??' + ? Also line 3 under ?Details? says <<< for a p?element row-vector x, ? >>> x is here surely a column vector, albeit the transpose of a row vector from the data matrix. cf page 322 of ?Modern Applied Statistics with S?, 4th

I have written a function to plot data which will be used for various different chemistries. A simplified version is: plot_data <- function(risk,levels,chem,sd2,measure){ plot(risk, levels,main=paste ("per", sd2, measure, "\n in usual", chem)) } The problem is with the title. This works fine if the variable "chem" is just text, but if it is an expression then

Does anyone know any more than is in the following press release about REvolution Computing and their commercialization of R? http://www.intel.com/capital/news/releases/080122.htm "Intel Capital, the global investment arm of Intel Corporation, today announced that it has invested in the Series A financing of REvolution Computing, creator of parallel computing software for computational

The following is 'semicolon.Rnw' > \SweaveOpts{engine=R, keep.source=TRUE} > > <<xycig-A, eval=f, echo=f>>= > library(SMIR); data(bronchit); library(KernSmooth) > @ % > > Code for panel A is > <<code-xycig-A, eval=f, echo=t>>= > <<xycig-A>> > @ % Sweave("semicolon") yields the following 'semicolon.tex'

rug() may add integer axis labels when called with a POSIXlt object as argument. dtimes <- c("09/29/2007 12:54", "09/30/2007 00:14", "10/01/2007 00:14", "10/02/2007 00:14", "10/03/2007 00:14", "10/04/2007 00:14", "10/05/2007 00:14", "10/06/2007 00:14", "10/07/2007

You need to give the model formula that gave your output. There are two sources of variation (at least), within and between locations; though it looks as though your analysis may have tried to account for this (but if so, the terms are not laid out in a way that makes for ready interpretation. The design is such (two locations) that you do not have much of a check that effects are consistent over

In the following, 'probability=T' is ignored, i.e. the vertical axis is labelled in frequencies. However hist(islands, freq=F) achieves the desired effect. > data(islands) > hist(islands, probability=T) # Vertical axis gives frequencies > hist(islands, freq=F) # Expected result --please do not edit the information below-- Version: platform = i386-pc-mingw32

I have a file 4lmetc.Rnw, intended for inclusion in a LaTeX document, that starts: \SweaveOpts{engine=R, keep.source=TRUE} \SweaveOpts{eps=FALSE, prefix.string=snArt/4lmetc} The attempt to process the file through Sweave generates the error: > Sweave("4lmetc") Writing to file 4lmetc.tex Processing code chunks ... 1 : keep.source term verbatim Error in file(srcfile$filename, open =

I am doing [R version 3.1.2 (2014-10-31) -- "Pumpkin Helmet?; Platform: x86_64-apple-darwin10.8.0 (64-bit)] > R CMD build DAAGviz > R CMD check DAAGviz_1.0.3.tar.gz Without a .Rinstignore file, I get: <<< The following files should probably not be installed: ?figs10.pdf?, ?figs11.pdf?, ?figs12.pdf?, ?figs13.pdf?, ?figs14.pdf?, ?figs5.pdf?, ?figs6.pdf?, ?figs9.pdf?

I get the same style of path as Hadley. This is on Windows 7 Home Premium with SP1. I start R by clicking on the R-2.31.0 icon. I'd assumed that it was a change that came with R-2.13.0! (On 32-bit Windows XP, which I have just checked, I do indeed get the 8.3 paths.) > R.home() [1] "C:/Programs/R/R-2.13.0" > sessionInfo() R version 2.13.0 (2011-04-13) Platform:

I am puzzled by the warning message in the output below. It appears whether or not I fit the seasonal term (but the precise point of doing this was to fit what is effectively a second seasonal term). Is there some deep reason why AR parameters ("Warning message: some AR parameters were fixed: ...") should somehow intrude into the fitting of a model that has only MA terms? >

The following code demonstrates an annoyance with plot.lm(): library(DAAGxtras) x11(width=3.75, height=4) nihills.lm <- lm(log(time) ~ log(dist) + log(climb), data = nihills) plot(nihills.lm, which=5) OR try the following xy <- data.frame(x=c(3,1:5), y=c(-2, 1:5)) plot(lm(y ~ x, data=xy), which=5) The "Cook's distance" text overplots the label for the point with the