similar to: support for POSIXct classes in image.default axes

Questions on image: 1) How can I put labels on the x and y axes? 2) How to tell it to use e.g. 32 grey levels (not some colour map)? 3) How to know what the legend is (i.e. each grey level = what z value)? Thanks very much for any help. BTW I was wondering if persp was on the To Do list. That would be great! Bill Simpson

I noticed some undocumented and inconsistent behavior in device_Raster when a plot is produced with reflected axes such as: image(volcano, xlim = c(1,0), useRaster = TRUE) image(volcano, ylim = c(1,0), useRaster = TRUE) The `pdf` device will perform horizontal and vertical reflections, while `quartz` will ignore the transformations when plotting to the screen, but when plotting to a

On Thu, 4 Aug 2016 at 17:53 Michael Sumner <mdsumner at gmail.com> wrote: > On Thu, 4 Aug 2016 at 11:17 Paul Murrell <paul at stat.auckland.ac.nz> wrote: > >> Hi >> >> Just to clarify, I think this IS a problem with grid.path() as well as >> polypath(). >> >> > Hi, oh dear - sorry about that > > I appreciate the deeper explanation, I

Hello, I'm having trouble with using the alpha channel for transparency with lines with lattice levelplots. If I use transparency via the alpha argument to rgb to overplot lines on levelplot the transparent colour affects all of the region colours in the plot. Can anyone explain why the difference in region colours? #### Warning: this code attempts to create PDF files in working directory

hi, i am using par(mrow=c(6,6)) function to get 6x6 plots on one screen. the problem that i am having is that the axis tick labels are far away from the ticks and going into previous plots (see attached figure). i need to know how can i reduce the distance between the ticks and their values (y axis values).the part of the code that i am using (after reading in the data) to create the top row is:

Hi everyone I want a square scatterplot with abline(0,1) going exactly through the SW and NE corners. By "square" I mean that the plotting region is exactly square, and that the axis limits are identical. x <- 1:20 y <- x+rep(c(-1,1),10) lims <- range(c(x,y)) None of the following do this: plot(x,y) ; abline(0,1) #not square plot(x,y,asp=1);abline(0,1) #diagonal

Dear All, I was looking for the color index in?image function,?such as from topo.colors(n) and etc. but still never found it. For instance, from the help menu. ########################################### # Volcano data visualized as matrix. Need to transpose and flip # matrix horizontally. image(t(volcano)[ncol(volcano):1,]) # A prettier display of the volcano x <- 10*(1:nrow(volcano)) y

Full_Name: Fran?ois Birgand Version: 2.9.0 OS: windows xp Submission from: (NULL) (152.1.16.161) When I type this sequence: x11() Ylim<-c(0,14) Xlim<-c(0,14) plot(0,0,xlim=Xlim,ylim=Ylim,col="white",main="",xlab="",ylab="",bty="n",xaxt="n",yaxt="n",xaxs="i",yaxs="i",asp=1)

I ran into a weird, to me at least, problem, and hoping someone can shed some light into it. In a nutshell, there seems to be some problem when one calls plot with a formula, from within another function, using ... to pass arguments, and one of those arguments being xlim (and only xlim shows this problem). Here is an example: > plotw <- function(obj,...) { + plot(k~j,

Full_Name: Mark Wall Version: 1.6.0 OS: linux Submission from: (NULL) (63.251.119.254) I want to plot a histogram of a *subset* of some data: >t = c(0:9) >hist(t,right=FALSE,breaks=10,xlim=c(0,5),xaxs="i") This means I should plot a histogram from 0 to 5 with breaks at 1,2,3,4. This should produce exactly 5 bars of frequency=1. Instead I get 5 and 1/4 bars. I do not want the

Hi Just to clarify, I think this IS a problem with grid.path() as well as polypath(). For the example you give, grid.path() diverts to drawing a polygon (because there is no 'id' specified), and the NAs in 'x' generate two separate polygons, which get drawn one on top of the other. The correct analogy to the polypath() example is ... x2 <- matrix(x[!is.na(x)], ncol=2)

Hi,   I would like to distinguish my plots using colors but I got error message.  How do I correct that?     plot(ecdf(z), main ="CDF for observed and simulated weighted sum",type="l",lwd=2,col="blue", xlab="Weighted sum (mm)", ylab="Cumulative Percent", xlim=c(0,15), xaxs ='i', yaxs ='i',ylim=c(0,1)) par(new=TRUE)

Hello, if I interrupt Sweave while it's processing a file it seemingly leaves an open sink connection that hides printed output. Can this be changed to reset the sink on exit? I've been baffled by this for years. This is seen in Windows (R Under development (unstable) (2013-10-20 r64082)) and an older Linux (R version 3.0.0 (2013-04-03)). Run the code below in two parts with a manual

Hi! Whenever I save a graphic in ps/eps format generated by filled.contour or image, an undesired grid is added to it (not visible on the X11 screen). For example: postscript("volcano.eps") filled.contour(volcano,col=gray(seq(0,1,,50)),levels=seq(min(volcano),max(volcano),,50)) dev.off() Any ideia how to eliminate this grid? Thanks, Rachel [[alternative HTML version deleted]]

Dear R-helpers, I am having difficulty making R plot only integer labels on the x-axis of a simple graph. I want to plot the median values of a score on each of three occasions. Non-integer occasions are impossible. But, R keeps labelling the x-axis with half-occasions, despite my attempts to stop this using the "xaxs" and "xaxp" parameters of 'plot'. p1=c(1,2,3);

Dear R users, The parameters xaxs and yaxs (function par, package graphics) seem not to work with the function hist (package graphics), even when the parameters xlim and ylim are defined. Is there any way to make yaxs="i" and xaxs="i" work properly with the function hist, mainly to produce histograms that "touch" the horizontal axis? The R documentation and the R

Hi I am using the ggplot2 package for the volcano plot and I am using the following code for the same: g = ggplot(data=data, aes(x=data[11], y=-log10(data[12]), colour=threshold)) + + geom_point(alpha=0.4, size=1.75) + + opts(legend.position = "none") + + xlim(c(-10, 10)) + ylim(c(0, 15)) + + xlab("log2 fold change") + ylab("-log10 p-value") data[11] is a

When plotting a reversed time axis with the base graphics system, R fails to label the axes. Example: times <- c(Sys.time() - 100, Sys.time()) plot(times, 1:2, xlim = times) # correctly labels the x-axis plot(times, 1:2, xlim = rev(times)) # fails to label the x-axis axis.POSIXct(1, times) # still fails to add labels dates <- c(Sys.Date() - 10, Sys.Date()) plot(dates, 1:2, xlim = dates)

I have some code as shown below. Basically, I would like three boxplots to be set next to each other with no ylabels on the two "inner" plots, and I want the same x axis range on all three. However, it seems like boxplot does not respect the xlim setting. I've tried the various ways I thought would work (par, boxplot(...xlim=)) but none of them seem to work. I then tried plot.window,

Hi Folks, I've been digging for the solution to this for several hours now. If there is a solution, it must be one of the worst "needle-in-a-haystack" examples in R documentation! Essentially, I want to make an x-y plot in which the X-axis really is the X-axis (i.e. its vertical position is at y=0), and the Y-axis really is the Y-axis (i.e. its horizontal position is at x=0).