similar to: size of subplots with par() / layout()

Hello list. I very often need 3d scatterplots, and use scatterplot3D quite a lot. I am trying to modify persp() to plot scatterplots, and make use of the theta and phi arguments that persp() offers. I am having some difficulty passing the correct arguments to persp(). Here is my function so far. Much of it is copied from the persp() manpage. points3d <- function(x,y,z,

Hello list. I very often need 3d scatterplots, and use scatterplot3D quite a lot. I am trying to modify persp() to plot scatterplots, and make use of the theta and phi arguments that persp() offers. I am having some difficulty passing the correct arguments to persp(). Here is my function so far. Much of it is copied from the persp() manpage. points3d <- function(x,y,z,

Hi the manpage for persp() has a wonderful section where a the trans3d () function is used with points() and lines() to add red dots and a green sinusoid to the Mexican hat surface. Does anyone have a way to tell what distance a point is from the eye/ camera? Take the following line: lines (trans3d(x, y=10, z= 6 + sin(x), pm = res), col = 3) Is there a function like trans3d() that returns

Hello I am trying to draw a basic black and white map of two European countries. After searching some key words in google and reading many pages I arrived to the conclusion that persp() could be used to draw that map. I have prepared three small example files, which are supposed to be the files required for running that function. xvector is a vector with the longitudes yvector is a vector with

Hi, I'm new to R and am learning the ropes of r-help and programming. I'm attempting to plot a 2-D mesh in 3-D using the persp function. I've positioned the graph using phi and theta, but I'm having trouble overlaying color onto the actual surface according to the surface's "z" values. I've tried the "col" feature but it doesn't put a flowing

Jacobi's theta functions crop up here and there in various branches of mathematics, such as unimodular and elliptic functions. They have Taylor expansions, but the powers increase as the square of n, as in 1 + z + z^4 + z^9 + z^16 + z^25 + . . . so they converge very quickly, provided |z|<1 The following toy function shows how I'm implementing these objects. I just add terms

Hello One of my packages, untb_1.3-2, passes R CMD check under MacOSX (and apparently the systems used in the package check summary page on CRAN) but fails with the following message on R-2.6.0.tgz compiled last night on my (home) linux box. I hasten to add that I have never seen this error before on home-compiled pre-releases of R-2.6.0. Can anyone help me understand what is going on?

Hi I want to write a little function that takes a matrix X of size m-by-n, and a list L of length "m", whose elements are matrices all of which have the same number of columns but possibly a different number of rows. I then want to get a sort of dumbed-down kronecker product in which X[i,j] is replaced by X[i,j]*L[[j]] where L[[j]] is the j-th of the "m" matrices. For

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, I'm new to R (and math ;) Would somebody please be so kind as to direct me in plotting a 3D sphere? I tried something in the lines of: #### y <- x <- seq(-pi, pi, length=pi*10) f <- function(x,y) { z <- sqrt(pi - x^2 - y^2) #z[is.na(z)] <- 0 z } z <- outer(x, y, f) persp(x, y, z, theta = 120, phi = 30) #### I've also tried: .... make.surface.grid(...) ..

Dear R-devel mailing list, I am wondering whether the "theta" and "phi" parameters of the "persp" function in the graphics commands are named in error (the names seem to have been swapped). Also, in the documentation to "persp", reference is made to "colatitude" when the effect in the image is actually that of "latitude". Thank you,

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

Hi I'm having difficulty with a multiple plot. What I want is 12 plots, all of which are the same size and shape, differing only in colour. Each one is a square, so there is an asp=1 in the plot call. I'm working in an Sweave environment so I am free to choose the height and width of the plot. I want the columns to be labelled at the head (here t=1,2,3) and the rows to be labelled at

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 :

Hello List does anyone have an R function for the Lambert W function? I need complex arguments. [the Lamert W function W(z) satisfies W(z)*exp(W(z)) = z but I could'nt even figure out how to use uniroot() for complex z] -- Robin Hankin Uncertainty Analyst Southampton Oceanography Centre SO14 3ZH tel +44(0)23-8059-7743 initialDOTsurname at soc.soton.ac.uk (edit in obvious way; spam

Hi all, I'm making an animation of a rotating persp plot of some data. It's easy, really: I just plot the persp in a loop over the theta parameter, saving each plot, and then chain them together with ImageMagick (as suggested on Paul Johnson's R-tips page). But... when you run the animation, because each graph is individually re-scaled, the image is not the same size in each

R'istas: I am trying to plot a marginal distribution on the plane behind a persp() plot. My existing code is: library(MASS) X <- mvrnorm(1000,mu=c(0,0),Sigma=matrix(c(1,0,0,1),2)) X.kde <- kde2d(X[,1],X[,2],n=25) # X.kde is list: $x 1*n, $y 1*n, $z n*n persp(X.kde,phi=30,theta=60,xlab="x_b",ylab="x_a",zlab="f") ->res Any suggestions are very

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,

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,

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