similar to: summary: Call SAS from R

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")

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

Hi, According to help(sub), the ^ should match the zero-length string at the beginning of a string: sub('^','var',1:3) # "1" "2" "3" sub('$','var',1:3) # "1var" "2var" "3var" # This generates what I expected from the first case: sub('^.','var',11:13) # "var1" "var2"

Hi All, I am new to post question on this list. I apologize if this question is too easy or irrelevant. I am doing a simulation study and I need to read a data file that can be easily read by SAS. So, what I try to do is to execute SAS in R and then read the output of SAS to R. I try the following code but it didn't work. system("c:\\program files\\sas institute\\v8\\sas.exe

Hi All, I'd like to plot something like http://www.nawcwpns.navy.mil/~weather/mugu/mesodata/analysis.html Looking through the galleries at http://addictedtor.free.fr/graphiques/allgraph.php http://r-spatial.sourceforge.net/gallery/ http://fawn.unibw-hamburg.de/cgi-bin/Rwiki.pl?GraphGallery demo(graphics) I did not find a function to plot a 2d field on a matrix. I did find mention of a

Hi all, I'd like to quantize a variable to map it into a limited set of integers for use with a colormap. "image" and filled.contour" do this mapping inside somewhere, but I'd like to choose the colors for plotting a set of polygons. Is there a pre-existing function that does something like this well? i.e., is capable of using 'breaks'?

Hi, I'd like to produce a marginless or zero margin plot so that the pixel coordinates represent the mathematics. xy<-data.frame(x=c(0,1,1,0,0),y=c(0,1,0,0,1)) png('junk.png',width=300,height=300) par(mar=c(0,0,0,0)) plot(xy$x,xy$y,xlim=c(0,1),ylim=c(,1)) dev.off() The resultant file has about a 10 pixel margin around these lines, and I'm not sure what parameter or function

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")

Hi All In trying to correlate some tide gauge data I need to deal with varying timezones. From the documentation on strptime, it seemed that the tz variable might have some effect on the conversion, but I'm not seeing an effect. > strptime("20061201 1:02 PST",format="%Y%m%d %H:%M",tz="PST")+0 [1] "2006-12-01 01:02:00 EST" >

I have some disconnected boundary data from a finite element ocean model and I'd like to make a plot. Maptools looks promising, but since my data is not in a shapefile or a map, I'm unclear on what the best way to approach the problem. >geom[1:10,] lon lat depth 1 -75.42481 35.58192 16.172 2 -75.40726 35.58567 18.045 3 -75.41351 35.60312 17.333 4 -75.38888 35.58959

I don't want to gripe here, but I would like to know the best method for submitting suggestions and patches to the documentation. It looks like for many help items, we need to track from the html file which help.search("*")/help("*") presents, back to the authoritiative documentation file from which they are generated. It seems the installation process discards much of

Suppose I have x<-data.frame(v1=1:4, v2=c(2,4,NA,7), v3=rep(1,4), v4=LETTERS[1:4],v5=rep('Z',4)) or a much larger frame, and I wish to test for and remove the constant numeric columns. I made: is.constant<-function(x){identical(min(x),max(x))} and apply(x,2,is.constant) # Works for numerics x[,-which(apply(x,2,is.constant))] I'd really like to be able to

Is there a tool for reading a graphics file into an object? I might not be looking with the correct vocabulary, but I'm finding lots of references to producing graphics from R, but not any for inputting graphics into R. I'd like to use a jpeg image and have the data available in R for analysis. I could convert the image to another format, PPM perhaps, and parse that in, GRASS seems like

>>>>> "David" == David Forrest <drf5n@maplepark.com> >>>>> on Tue, 22 Mar 2005 15:02:20 -0600 (CST) writes: David> According to help(sub), the ^ should match the David> zero-length string at the beginning of a string: yes, indeed. David> sub('^','var',1:3) # "1" "2" "3"

One quirk to be watched is the use of connections to the DBMS via the Unix sockets vs ports. The PostgreSQL driver bundled with unixODBC will use Unix sockets to `localhost', but this driver seems unreliable (see the ChangeLog). The current driver will only use a TCP/IP port, and to use that needs postmaster started with the -i flag (which is not the default) and with tcp/ip

Hello, I'd like to do something like: n<-100 zz<-cbind(rnorm(n),rnorm(n),floor(runif(n)*3+1)) colors<-c("red","green","blue") plot(zz,col=colors(zz[3])) and have a matrix of scatterplots colored by class. The above does not work, of course, but I'm not sure exactly what function I'm looking for. Thank you for your time, Dave, -- Dave

Is there a package for stacks with pushing and popping? I'd really like to do something like push(par(no.readonly=TRUE)) some stuff.... par(pop()) It seems like it wouldn't be difficult, and that someone may have already implemented a set of stack functions, and I wouldn't like to duplicate it. Thanks for your time, Dave -- Dave Forrest (434)924-3954w(111B) (804)642-0662h

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 -

I know browsing is an old topic, but I've not seemed to find out what is causing my problem. Browsing seems to work from the command line, but not graphically. Complaint: inside of windows, mug does not show up on the browse list, but allows connections. Environment: se-ntserver is a NT Server on domain SE-NTNET 128.143.6.* mug is linux on 128.143.6.* (running samba 2.0.6) vmware

Okay, been working with tripack, seems the most mature package for this. Got it to work well with their test data set - data(tritest). When i tried random numbers to explore further, i am getting some results that don't reconcile. example run this: library(tripack) y <- runif(100) x <- runif(100) vm <- voronoi.mosaic(x,y) plot(vm) par(new=T) plot(x,y,col='blue') when