similar to: color of opposite sign values in filled.contour

Hi all I'm plotting colored contour maps using filled.contour. My levels are very unevenly spaced, with, say, high resolution in the small numbers but ranges that can be an order of magnitude or two larger in absolute value compared to where the action takes place. Aside from transforming the data, is there a way to control the color spacing in the key to the right of the map? Right now I

Here is the answer to an old question of mine, regarding maps in R Hope this help. claudia ---------- Forwarded message ---------- Date: Mon, 25 Jun 2001 13:48:05 -0400 From: Kieran Healy <kjhealy at Princeton.EDU> To: Claudia Tebaldi <tebaldi at ucar.edu> Subject: Re: [R] maps in R Hi Claudia -- you can get the maps() library from here:

Hi all Just this morning I upgraded to R 2.12.0 (for Mac OS X 10.6.4). All went well until I needed to run a help() or help.search() in my session, which I'm running within Emacs (ESS 5.3.7). Say I need help with the command 'density'. When I type help(density) or ?density the ESS help buffer opens, it is titled *help[R](density)* but it contains only a couple of lines saying,

I'm running R-1.4.1 on a linux debian (version 2.2) box. I've tried to upgrade to R-1.6.0, and the execution of make halts with the following error messages: make[4]: Entering directory `/home/tebaldi/R/R-1.6.0/src/library/methods' dumping R code in package 'methods' Error in testRversion(descfile) : This package has not been installed properly See the Note in ?library

Hi all I have a dataset whose response is a categorical variable of the ordinal scale type (6 levels). I'm interested in building classification models, and I'm wondering if there is something implemented in R (or its packages) that I'm not aware of, to treat the ordinal scale measurements straightforwardly. I can think of the alternative of building a conditional (hierarchical)

Hello I have several replicates (30) of a -- short -- time series (19 obs) and I'd like to estimate the first few lags autocorrelation function. Is there an immediate way of handling multiple realizations when computing the ACF in R? I'm using R-1.4.1 on a Linux Debian box. Thanks in advance claudia ------------------------------------------------------------------------- claudia

Dear R-helpers, I have lat/lon coordinates of regularly spaced grid points, about 4Km apart, covering the entire US continental region. I would like to mask this rectangular grid in order to extract all and only the grid points within a specific region. Today I want to extract Montana, say, from this grid, and I am hoping to somehow use the returned value of the function

There is probably an elegant and efficient way to do this -- perhaps already implemented, and some of you knows it. Myself, I'm getting stuck in long loops ... Here is the question: I have a long (80,000 + obs) binary time series. For each of its elements I'd like to compute and save in a vector "the number of steps 'til the first switch of state happens". For example,

Hi I'm trying to print several plots to a multi-page postscript file. Some plots are output of basic commands, some are produced by "trellis" commands (from the package lattice). Right now I'm not able to get a straightforward black and white color scheme for the latter kind: when I open the postscript file the "trellis" plots are invisible, or parts of them are

Hi all I'm running lm()/leaps() kind of computations within a loop, and I submitted the entire thing by sourcing a file of R code in interactive mode. I'd like to interrupt the evaluation of this (it's taking way too long) but I am not able to. I'm running R-1.3.0 from an emacs window (ESS-5.0) and I'm on a linux box, (Debian). I tried Control-c Control-c (nothing happens),

Hello I'm trying to render in 3D what I usually plot by image(), or image.plot() from the library fields, followed by a map("world",add=TRUE) type of command. More concretely, I have a field of temperature values, for a given geographic area, and I would like to plot a 3D surface whose x and y axes consist of longitude and latitude values and whose height and color-coding correspond

Hello I'm upgrading from R-1.4.1 to R-1.5.1, locally on my Linux PC running Debian. after running 'configure' this is what I get (just in case it helps figuring out what follows later...) R is now configured for i686-pc-linux-gnu Source directory: . Installation directory: /usr/local C compiler: gcc -D__NO_MATH_INLINES -mieee-fp -g -O2 C++

Dear R colleagues, I am trying to plot some geophysical data as a filled contour on a continent map and so far the guidance from the R-help archives has been invaluable. The only bit that still eludes me is the colour key (legend) coming with filled.contour: I prefer to generate my own colour palette, mainly based on the quantiles of tenths of the data in order to capture the whole range (of

Dear R-users, I haven't found a way in the searchable archive to overplot a contour (lines) over a surface. I have a (n,m) matrix that represents sea surface temperature that I have plotted using image.plot(), filled.contour() or image(). I would like to overplot this image with some contour lines of mixed layer depth values(same size matrix). How can I do this? Any help is appreciated,

Hello, I am trying to plot SADIE red-blue plots of cluster indicies using filled.contour. I want a plot which only has three bins for the data: <-1.5, -1.5 - 1.5, >1.5, but I am having trouble getting there. example X1 X2 X3 X4 X5 1 -5 -4.5 1.0 4.5 6 2 -3 -2.0 1.2 -1.0 3 3 0 0.0 0.0 -0.5 -1 4 -2 -3.0 1.0 1.5 3 5 -6 -2.0 0.5 3.0 2 example<-as.matrix(example)

> From: Martin Maechler <maechler at stat.math.ethz.ch> > To: royle at bearmtn-e0.cgd.ucar.edu > CC: lgygax at access.unizh.ch, r-help at stat.math.ethz.ch > Subject: contour labelling [was "Re: [R] image ()"] > Mime-Version: 1.0 (generated by tm-edit 7.106) > > >>>>> "Jeffrey" == Jeffrey A Royle <royle at bearmtn-e0.cgd.ucar.edu>

Dear all, As a newbie in R I would like to do the following (simple?) thing: to plot a filled.contour plot over a map showing country boundaries (e.g. for Europe) What i do is: map('worldHires',xlim=c(-10,40),ylim=c(35,70),boundary = TRUE,border=0.1) map.axes() filled.contour(mslp, zlim=c(1000,1020),color.palette = colorRampPalette(c("blue", "white",

Hello I'm having trouble figuring out how to use the output of "segmented()" with a new set of predictor values. Using the example of the help file: ??set.seed(12) xx<-1:100 zz<-runif(100) yy<-2+1.5*pmax(xx-35,0)-1.5*pmax(xx-70,0)+15*pmax(zz-.5,0)+rnorm(100,0,2) dati<-data.frame(x=xx,y=yy,z=zz) out.lm<-lm(y~x,data=dati) o<-## S3

Hi, I am plotting values of log(hazards ratio) as a function of two predictors, using the plotting function filled.contour(). Here is a simple simulated example of this: x <- seq(0,1, length=20) y <- seq(0,1, length=20) z <- outer(x,y, function(x,y) x^2 + y^2 ) zmat <- matrix( rexp(n=400, rate = z+0.001), 20, 20) filled.contour(x, y, log(zmat),

Hi I have integer values representing unix time, stored into a simple numeric vector, no POSIX* class inheritance of any kind. Something like crazylookingtimes_c(993340800.00, 993254400.00, 994809600.00) I'd like to convert them into a format like YYDDMMHHMMSS but I cannot find out how, at least not in the documentation on DateTimeClasses that I've skimmed. Can anybody help, please?