similar to: Dates in 'persp' plots

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,

Hi All. I have created two, 3-D, color graphs using persp(). Z values range from 1 to 100 on the first plot, and 0.5 to 50 on the second plot. I would like to keep the color range constant between the two graphs, rather than each graph showing the full range of colors. So the first graph should go from yellow to red and the second graph should go from yellow to orange. Below I present a

Dear R-users, I'd like to plot a three-dimensional surface and at the meantime I'm using an array. I would like to have the values of my first matrix in the heights of the plot and the colors of the single facet taking into account the second matrix. I hope that the next code will help all of you to understand better my issue, Thanks in advance, Giancarlo ############################ ##

Hi R-listers, I am having trouble with persp plotting and I hope some knowledgeable person can help me. I have searched the help files to no avail. I am sure there is a way of achieving this simple task. I thought in par maybe but that seems not to work for 3d plots. Maybe I have missed something. I read in ?persp that it expects increasing values for both my x and y axes, and a matrix

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 am trying to specify four ticks (at 0,1,2,3 for the y axis) in a persp plot but the defaults overrule my specification and I obtain seven of them. Is it possible to gain full control over them in such a plot? Here is my code: matlines=matrix(c(1:7,3:9,3:9,2:8),nrow=7,ncol=4)

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

Full_Name: Ramin Nakisa Version: 1.1.0 OS: i386 linux Submission from: (NULL) (192.127.94.7) Labels created with expression seem to disappear, e.g. the x axis label in x <- -6:16 contour(outer(x, x), method = "edge", vfont = c("sans serif", "plain"), xlab=expression(theta[0]), ylab=expression(theta[1])) and all three labels in x <- -6:16 persp(outer(x, x),

Dear R list, I want to plot the yield curve in a 3D graph: the x-axis refers to time; the y-axis refers to the maturities and the z-axis are the yields. The following code works fine: persp(y, x, as.matrix(ypper), xlab="Last 50 periods", ylab="Maturities (months)",zlab="Yields", zlim=c(3,6),theta=130, phi=15, col="Seagreen", box=T,

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

Thank you for your replies so far. Sorry for bothering you again, but I'm still not able to get what I need as I don't understand all parts of the replies (just using R for easy things....). Is there a code for plane3d() like some of you sent me for points3d()? I was not able to get that out of the scatterplot3d package... What I can do is to get the x,y and z-range for the xlim,ylim and

Hi, I have been using the persp() function to plot species response surfaces based on GAMs (using package mgcv). Everything works fine, except that I have a slight problem when it comes to producing a figure with multiple (4 - 6) persp() plots on it. The axis labels remain the same size, no matter how large or small the plot is, causing some overlap when the multiplot figure is produced. Is

I am having trouble getting my facets colored with "concentric" colors by z level. This code give a fair representation of my problems dealing a real dataset cmtx <- matrix(, nrow=10, ncol=10) cc <- (row(cmtx)-5)^2 + (col(cmtx)-5)^2 cc #---------- [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [1,] 32 25 20 17 16 17 20 25 32 41 [2,]

I am using the mgcv package to develop vis.gam plots and having trouble figuring out how to relabel the z-axis (image attached). It is currently labeled as "linear predictor," but I would like to change it to a different name. Currently I am using this code: vis.gam(model1,theta=320,ticktype="detailed",color="gray",nCol=12, zlab="BCS") However, when run

Dear R users, Could somebody tell me why the grey color ramp is repeated in this graph, eventhough the ramp values go from 0 to 1? I must be missing something obvious, but I can not see it! z <- c(0.064329041,0.117243316,0.161565116,0.19923015,0.231642175,0.259835539,0.284571226, 0.038507288,0.094184749,0.140959431,0.180803984,0.215159105,0.245096084,0.271412845,

Dear All, I am very new to R projects.May be i am wrong in some steps.I have given the code which i tried for drawing 3d surface using persp.I need to label the axes with scales z <- array(topnew2$V2, dim=c(600,2)) x <- 10 * (1:nrow(z)) y <- (1:ncol(z)) persp(x, y, z, theta = 30, phi = 30, expand = 0.5, col = "lightblue", xlab ="fluidlevel", ylab

Hello Masters, wishing you all a great 2011 I was also going to ask if anyone knows a quick and efficient way to plot a regression plane (z~x*y). I have tried the regr2.plot{HH} function but it is only an educational tool and has poor graphical properties. I also tried to run the following script on a fictitious longitudinal problem, with poor results set.seed(1234)

Hello Experts I have the following codes and data for 2 interpolation plots. http://www.nabble.com/file/7206/3d_plot_data.txt 3d_plot_data.txt data<-read.table("3d_plot_data.txt", header=T) attach(data) par(mfrow=c(1,2)) library(akima) interpolation<-interp(rr,veg_r,predict) persp(interpolation,theta = -45, phi = 30, ticktype = "detailed", nticks=4, cex=0.8,

Hi, I have thought a lot about my problem and also posted it on stackexchange, the post can be found here: http://stats.stackexchange.com/questions/44276/how-to-plot-3d-gbm but I have got no useful answer, that's why I trry the last possibility, which is to ask you professional guys. I want to recreate the following picture with my own code, the picture can be found in my other post:

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)