similar to: US 2004 Elections map

A few questions about maps... (1) How can I find a listing of the internal data sets that map() from the maps library contains? For example, "usa", "county", "state", "nz" all work. Are there any others? (2) Is there an easier, more generalized way to produce this (http://www.ai.rug.nl/~hedderik/R/US2004/ ) type of plot than this

I'm trying to get Sweave running for automatic report generation, and it seems to run fine when just using verbatim output. However, I've ran into a problem with xtable. I would like to print the following matrix using xtable: > dim(counts) [1] 19 15 All columns are filled with real/integer numbers > 0 and < 1000. Just typing xtable(counts) gives correct LaTeX output, but

A couple of months ago, probably using an older version of R, R used to run the following code just fine: library("foreign") data.exp1 <- as.data.frame(read.spss("dataDef.sav")) Issuing the same commands now (after starting R using --vanilla), gives me the following behavior: > library("foreign") > x <- read.spss("dataDef.sav") Error

Hedderik van Rijn <hedderik@cmu.edu> writes: > If the second argument to unsplit is not a simple vector (but a "list > containing multiple lists"), the function seems to have some problems. > > Given a slight modification of the examples in help(split): > > > xg <- split(x,list(g1=g,g2=g)) > > unsplit(xg,list(g1=g,g2=g)) > [1] -0.7877109

Hi, [some background] I have a dataset which describes a number of subjects doing a "scientific discovery". That is, they have to discover the rules underlying a particular domain. To do so, they have to set the levels of 5 variables which leads to a certain outcome. To identify what kind of "experiments" the subjects did, I want to combine these levels into one variable

Dear all, I've got a quite large dataframe (stor) with rows subject and rt (reaction time). I would like to split the reaction times per subject into 6 bins of equal size. Right now, I'm using the following code: bindata <- function(rt) { bindata <- rep(-1,length(rt)) binwidth <- length(rt)/6 bindata[order(rt)[(0*binwidth)+1:(1*binwidth)]] <- 1

Hi: How do I find average across three columns such as weight1, weight2, weight3 (20 people, each person has three different measures)? They have NA's also. Thanks Tom, Univ of Pittsburgh, ttrut at yahoo.com

I use AfterStep on Linux. I made an R.xpm file (attached) and put this icon in the "wharf". To invoke R I click on the R icon. I did not find an offical R.xpm, and so I am submitting this one for your collective consideration.... The idea behind the icon is that it looks like a scatterplot, which is "statistical". My main beefs with my R.xpm are that it doesn't have a

On Wed, Jan 7, 2009 at 11:05 AM, ajay ohri <ajayohri at yahoo.com> wrote: > > FYI..not a R -Help Topic, buy I dont know which list to post discussions like this. > Regards, > Ajay > > ---------- Forwarded message ---------- > From: ajay ohri <ajayohri at yahoo.com> > Date: Wed, Jan 7, 2009 at 10:46 AM > Subject: Re: R "Threatens" SAS, According to

I would be very thankfull for such a solution, too .... "R.E. Darnell" wrote: > > Would introducing a "Are you sure?" question after kill the *R* ( or > whatever) buffer with a hasty C-x k be a popular default > "feature". How can I include this for my own use? > > Thanks > > Ross Darnell --

After some fiddling around with the tapply command, I discovered that the factors (the INDEX argument) given to tapply must be specified in fastest-cycling first order. The following code shows how I discovered my error: (R version 1.2.1) -o-o-o-o-o- x <- as.data.frame(list(data=c(-9,0,3,1,-9,1,0,-9,0,3,1,-9,1,0), subj=c(rep(1,7),rep(2,7)),

Hi, I have an array of 2d node locations and an array triangles, and would like to plot something like a image or persp. An example of doing it with rgl is: library(ncdf) library(rgl) # wget http://www.maplepark.com/~drf5n/extras/teapot.nc teapot<-open.ncdf("teapot.nc") z<-get.var.ncdf(teapot,"tris") l<-get.var.ncdf(teapot,"locations")

Hedderik van Rijn <hedderik@cmu.edu> writes: > > Thanks for pointing this out. I think your code should work, although > > I'm slightly uneasy about actually modifying f, so how about > > Just curious, why are you uneasy about that? Does it have side effects? It is largely due superstition, but if you change f before calling split<-() you have split() and

Hi, I'd like to partition a 2d probability density function into regions of equal probability. It is straightforward in the 1d case, like qnorm(seq(0,1,length=5)) but for 2d I'd need more constraints. Any suggestions for how to approach this? Is seems like a spatial sampling problem but I'm not sure where to look. Thanks for your time, Dave -- Dr. David Forrest drf at

Hello all, It seems that the next improvement to the R Graph Gallery is categorization of the graphics, that way each graph will be easier to find. That step should be done *carefully* if we want to avoid the opposite side-effect : graph not reachable through the categories. That's why the wisdom of the R community is required. Graphics will be classified in : - categories -

Using netscape 4.03 (or 2.x or 3.x)/ R-0.62.1 with SPARC Solaris 2.4: when netscape is running for help, graphics in R quit with "Color allocation error". Closing netscape removes the problem. netscape does not interfere with other X11 imaging programs including ghostscript and ImageMagick except that it increases the contrast in my mailer windows slightly. (netscape does block

The following code produces an eps file with the tops of each of the ylabs clipped off. par(mfrow=c(2,2)) plot(runif(10), ylab="Function(Lengthy Expression)",xlab="Prediction") plot(runif(10), ylab=expression(Delta * Beta^2),xlab="Prediction") plot(runif(10), ylab="Function(Lengthy Expression)",xlab="Prediction")

[This email is either empty or too large to be displayed at this time]

Dear all, I'm trying to plot some data over a geographic map. I'm reading the documentation (?map)... The problem is here: *map(database = "world", regions = ".", exact = FALSE, boundary = TRUE,* * interior = TRUE, projection = "", parameters = NULL, orientation = NULL,* * fill = FALSE, col = 1, plot = TRUE, add = FALSE, namesonly = FALSE,* * xlim =

My friend Veslot proposed me the Raven ! http://www.teteamodeler.com/allopass/images/corbeau.jpg As a matter of fact the intelligence of this bird is comparable to one of a monkey, although its brain is close to a reptile brain. R is quite the same : small, compact, but so clever. This is a canadian website explaining (in french) the R-aven intelligence