similar to: contour & persp

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

Hi, I draw a surface corresponding to bivariate density of independent variables (rho=0) using persp(). Then I add a contour line (i.e., circle in my case) at a particular density. Below is a minimal example of what I have so far. # Bivariate density dnorm2d <- function(x, y, rho = 0) { xoy = (x^2 - 2 * rho * x * y + y^2)/(2 * (1 - rho^2)) density = exp(-xoy)/(2 * pi * sqrt(1 - rho^2))

Dear Uwe Ligges , I still can't finish it. *> aa* #my data x1 x2 y 5 0.05 6 4.4180 1 0.50 3 2.6979 4 0.50 9 2.9000 7 0.95 6 2.6230 8 0.95 6 2.9078 9 0.95 6 2.6727 3 1.40 3 2.4203 2 1.40 9 2.5329 6 1.85 6 2.4867 *> attach(aa)* *> persp(x1,x2,y* error in persp.default(x1, x2, y) : increasing 'x' and 'y' values expected

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

I've been thinking of adding the possibility of including "z" among the axes to be logged in image, contour, and persp. In the first two, it would only affect where the breaks were set if they are calculated automatically; it would have a bigger effect in persp. For example, image(x, y, z, log="z") would set 12 colours evenly spaced on a log scale of the z values.

R-sters: I'm interested in keeping data plotted in persp to preserve the wireframe look, I'd just like to change one of the lines drawn (in either the x or y direction) into a different color so that it stands out. Or is there some way to add a line (say, via lines(), or abline()) to a persp() plot at the designated x or y that would follow the z surface contour? I could add a line

Greetings, I have what I hope is a simple question. I would like to change my contour interval on the vis.gam( plot.type="contour") in the mgcv package. Is this a situation where I need to modify the function or is there a default value I can change? Thanks

The following plot is a first approximation to what I need: *********************************** mu1 <- 0 mu2 <- 5 s <- 1 x <- seq(-2.5, 7.5, length = 41) y <- seq(-2.5, 2.5, length = 41) f <- function(x,y){ term1 <- 1/(2*pi*sqrt(s*s)) term2 <- -1/2 term3 <- (x - mu1)^2/s term4 <- (y - mu1)^2/s term5 <- (x - mu2)^2/s term1*(.5 * exp(term2*(term3 + term4)) + .5 *

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 R-sters, I'm trying to add several lines to a response surface that I've plotted using persp(). I've tried lines() using the "trans3d" function but I've been unsuccessful in getting it to work (R v2.0.1). Essentially, I'm trying to highlight one or more of the surface wireframe lines in a bolder (or different) color. Any tips from those of you who have some

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)

someone sent in a question earlier about doing something in 3D so i took a stab at it purely for educational purposes ( i'm not even sure that I understood the question actually ). Unfortunately, persp gives me an error that I don't understand because it says "object y not found". I'm sending y in as a parameter to persp similar to what ?persp shows in one of oits examples

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

The following matrix of normed likelihoods should give a smooth surface but instead gives a series of spikes in both persp and contour (the dim labels are the axes values). I know an algorithm cannot be infallible but it would be nice to have some parameters to control the smoothing. (It takes close to an hour to produce this matrix on a Pentium II 300mh. It was after that that it crashed with the

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 R-help! I am trying to draw series of lines in 3d, and I am not sure how to get the depth set up properly - lines that are in front of other lines appear to be behind and vice versa. I've been doing this by first generating an empty persp plot, then adding lines via the data into trans3d followed by lines (this is probably a bad approach, but it was the first one I could get working).

Dear R help contributors, I'd like to plot ground temperature with time on X-axis and depth on Y-axis on this datasets ( http://r.789695.n4.nabble.com/file/n3301033/NEdaily.csv NEdaily.csv ), and to do so I use the following commands: library(RSEIS) xNE <- seq(1, as.numeric(as.Date(max(NEdaily[[1]])) - as.Date(min(NEdaily[[1]]))), 1) yNE <- rev(c(-0.3, -0.5, -0.7, -0.9, -1.1,

Dear all, As a newbie in R I would like to do the following (simple?) thing: to plot a filled.contour plot over a map showing country boundaries (e.g. for Europe) What i do is: map('worldHires',xlim=c(-10,40),ylim=c(35,70),boundary = TRUE,border=0.1) map.axes() filled.contour(mslp, zlim=c(1000,1020),color.palette = colorRampPalette(c("blue", "white",

Hello, I want to make a contour plot. I have a matrix (nrow=4,ncol=58) . and two lists, one with 4 values (the depths) and one with 58 values (time in days)( Depth (should be later on the y-axis),and Time(should be later on the x-axis)) In the matrix there is the temperature for a specified depth and time. there are missing values (NAs). I also put ranges for x-,y-,zlim. I

> From: Martin Maechler <maechler at stat.math.ethz.ch> > To: royle at bearmtn-e0.cgd.ucar.edu > CC: lgygax at access.unizh.ch, r-help at stat.math.ethz.ch > Subject: contour labelling [was "Re: [R] image ()"] > Mime-Version: 1.0 (generated by tm-edit 7.106) > > >>>>> "Jeffrey" == Jeffrey A Royle <royle at bearmtn-e0.cgd.ucar.edu>