similar to: persp + postscript/pdf (PR#1419)

Hi, I am using the persp function to plot 3D surfaces, but the plots have "little white lines" when I print them to a pdf file (visible in Acrobat, Foxit, Evince, Xpdf and Gimp). This does not happen when I create png or tiff images. Here is some sample code: pdf("test.pdf") x <- seq(0,1,length=101) f <- dnorm(x, 0, 0.25) z <- c() for(i in 1:100) z <-

I am having some trouble getting the persp() package to change the x and y axis on a 3d plot. It defaults to the [0,1] interval and when I try to change it I get errors. Example: This works: ------------ D <- c(1,2,3,4,5,6,7,8,9,10) M <- c(11,12,13,14,15,16,17,18,19,20) DM <- cbind(D,M) persp(DM, theta = 40, phi = 30, expand = 0.5, col = "lightblue", ltheta = 120,

Hello all, I am trying to use computer modern fonts in my r grahics. I tried to do, as described here: http://www.stat.auckland.ac.nz/~paul/R/CM/CMR.html but unfortunately, it does not work. First of all I downloaded the cm-lgc package and the AFM and PFB-files from the page and put them in my R working directory, so far, so good. Then I tried to run the following code: > sn <-

Dear All, I have been trying to to do this for a few days now. I can generate persp plot OK and can generate it in colour using the command persp(x,y,z,col=terrain.colors(20)) However, I was wondering if it is possible to shade the 3d surface like a contour plot. i.e. black for large z, white for small z, say Thanks in advance Colin

Hello folks, I am a beginner R user. I have been able to make a 3D surface plot using 'persp'. The surface is made by a grid of lines emanating perpendicularly from each of the x and y axes at regular intervals. I can get rid of that grid by setting 'border=NA'. Can anyone suggest some ways to replace the grid with contour lines, to create a 3-dimensional contour map? Thanks

Hi All, I'm using 'tck' option to *reduce* the length of tick marks but it is not working, can anyone please tell me where I'm going wrong... require(graphics) require(grDevices) x <- seq(-10, 10, length= 30) y <- x f <- function(x,y) { r <- sqrt(x^2+y^2); 10 * sin(r)/r } z <- outer(x, y, f) z[is.na(z)] <- 1 # 'bg' works but 'tck' is not showing

Dear all, I would like to put two perspective plots into one plot. The help page for ?persp shows how one can add points and lines but not another perspective plot. data(volcano) z <- 2 * volcano # Exaggerate the relief x <- 10 * (1:nrow(z)) # 10 meter spacing (S to N) y <- 10 * (1:ncol(z)) # 10 meter spacing (E to W) ## Don't draw the grid lines : border = NA persp(x,

Hello folks, i'm a R beginner and i want to put in a same plot both contour lines and persp plot. For example, fn<-function(x,y){sin(x)+2*y} #this looks like a corrugated tin roof x<-seq(from=1,to=100,by=1) #generates a list of x values to sample y<-seq(from=1,to=100,by=1) #generates a list of y values to sample z<-outer(x,y,FUN=fn) #applies the funct. across the combos of x

Dear Devels, I'm not sure if I'm going the right way by mailing to this list. I have two suggestions regarding the persp function from the graphics package: 1. In the documentation of the persp function, it is said, that the "col" argument will ignore transparent colors. However, this is not true, if you use strings as color code - at least it works with HEX-code (eg. :

Dear all; Does anyone know how to add grid lines to a persp plot? I've tried using lines(trans3d..) but the lines of course are superimposed into the actual 3d surface and what I need is something like the plot shown in the following link: http://thermal.gg.utah.edu/tutorials/matlab/matlab_tutorial.html I'll appreciate any ideas Thanks PM

Well, I started playing with fractals in R, and wrote a function to generate de Mandelbrot set, which might be of interest to some people ########################################################################### # Mandelbrot set ########################################################################### mandelbrot <- function(x = c(-3.0, 1.0), # x coordinates y =

Dear list members, I want to draw surfaces of Dirichlet distributions with different parameter settings. My code is the following: #<begin code> a1 <- a2 <- a3 <- 2 #a2 <- .5 #a3 <- .5 x1 <- x2 <- seq(0.01, .99, by=.01) f <- function(x1, x2){ term1 <- gamma(a1+a2+a3)/(gamma(a1)*gamma(a2)*gamma(a3)) term2 <- x1^(a1-1)*x2^(a2-1)*(1-x1-x2)^(a3-1)

Considering the R function/plot shown below, I wonder whether it is possible to do the following changes: (1) Change the color of each point to be picked from list of colors according to its z-value. (The range should be from blue (z=0) to red (z=1).) The grid should then be omitted. [I have seen "terrain.colors", but don't know how to use it for this purpose.] (2) Add two lines to

Hi guys, After all the discussion yesterday about persp and color, I decided to have a more closer look at demo(persp), and decided to write a function to generate 'topo-like' colours to plot perspectives (Thanks a lot to Uwe Ligges for his enlightning comments regarding the code in the demo). Here it goes, I believe that this function will be pretty useful to a lot of people:

Hi, I'm trying to use mars function in R to interpolate nonlinear multivariate functions. However, it seems that mars gives me a fit which uses only very few basis function and it underfits very badly. For example, I have tried the following code to test mars: require("mda") f <- function(x,y) { x^2-y^2 }; #f <- function(x,y) { x+2*y }; # Grid x <-

Hi folks, I'm experiencing some hassle to add some text on a persp plot... Here is the code to generate the persp plot: x <- seq(-1.95, 1.95, length = 30) y <- seq(-1.95, 1.95, length = 30) z <- outer(x, y, function(a,b) a*b^2) persp(x, y, z,xlim=c(-2,2),ylim=c(-2,2),zlim=c(-8,8), phi=30, theta=-30,nticks=8,ticktype="detailed") I've tried the text() function :

Hi, I'm unable to find a way to smooth data for a persp() graph. Example, suppose that I have data x,y,z like this: x <- 1:10 y <- 1:10 k <- 20 z <- outer(x, y, "*") + matrix( k*runif(100, -1, 1), 10, 10) persp(x, y, z, theta = 35, phi = 25) The graph is not very nice. Is there a way to smooth the z data so that at the end the graph would look more like something

Hi all, I need to plot two intersecting planes in a graph. Using persp() function (see code below) I am not able to add a second plane with good results. I tried with par(new=TRUE) but the output was very poor. I have found some previous discussions about this point but (apparently) without a clear solution. Any advice would be appreciated. Simone Tenan ################## library(popbio)

Dear WizaRds, I would like to create a 3d-plot with persp(). I sampled 17 points with xyz-coordinates (real-life example!), representing the peaks of the whole plane with "zero coordinates" x=3,y=3,z=3. My intention is to show which entries are above or below the "zero" level with persp() on a nicely created grid. I also tried scatterplot3d(), but, alas, I am unable

I have a set of data in the form x1, y1, z1 x1, y2, z2 ... x1, yN, zN x2, y1, z(N+1) x2, y2, z(N+2) ... x2, yN, z(2N) ...and so on... xM, yN, val(M*N) I have been trying to figure out how to get R to use this data in a persp plot. So far the only thing that I can figure out to do is to break the data file into three different files. The first file contains the x-coordinate data: x1 x2