similar to: To add cut-off points in surface response with lattice

Hello list! I am trying to project the fitted surface to a 3d plot of the data, similar to figures 13.7 or 6.5 in Deepayan Sarkar's "Lattice, Multivariate Data Visualization with R", but replace the contour/map lines with "levelplot". Problem is I can't get the color regions to line up after the coordinate transformation. Is there a simple solution my geometry

Hello List! I asked this before (with no solution), but maybe this time... I'm trying to project a surface to the XY under a 3d cloud using lattice. I can project contour lines following the code for fig 13.7 in Deepayan Sarkar's "Lattice, Multivariate Data Visualization with R", but it fails when I try to "color them in" using panel.levelplot. ?utilities.3d says there

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'd like to compare SST data for year 2000 with SST for year 2001. I managed to get filled contour plots showing monthly SST for both years (12 plots for each year, 24 plots in total). In order to compare year 2000 and year 2001, however, I'd need to have the same legend scale (same color bar) in both years plots. Thus, my question is: is there a way to set the legend scale in

Hi all, To my knowledge, the current version of contour.default() does not handle the 'asp' parameter. This can be embarassing when displaying eg geographical maps, etc... Submitted to the opinion of more experienced R programmers, contour.defaut() function should be changed according to the followings: line 7: add = FALSE,asp=NA,...) line 33: plot.window(xlim, ylim,

Hi, I'm trying to follow the suggestions given by Deepayan Sarkar in this message: http://tolstoy.newcastle.edu.au/R/help/05/11/16135.html to plot 3-D points on a wireframe plot. The problem is that I keep getting a partly formed plot- with the colored lattice visible but no axis labels or additional points- with the error message "error using packet 1, 'x' and 'units'

Hello all, I’m trying to add some data points to a wireframe. X an Y axis are independent variables, Z axis is predicted probability. I’d like to add the original data points on which the predicted probabilities are based to the wireframe. I’ve followed some of the previous post on this but get stuck : ## first part yields the wireframe## setwd("C:/Temp") rnp <-

Hi All, Does anyone know of a function to plot a regression surface for two predictors? RSiteSearch()s and findFn()s have not turned up what I was looking for. I was thinking something along the lines of: http://mallit.fr.umn.edu/fr5218/reg_refresh/images/fig9.gif I like the rgl package because showing it from different angles is nice for demonstrations. I started to write my own, but it has

Rers: I can't seem to locate the command to suppress what appears to be a faint internal grid when running the following command to make a filled contour plot of some data I have (x,y,z being the inputs):

Dear expeRts, I would like to add two points to a wireframe plot. The points have (x,y,z) coordinates where z is determined to be on the wireframe [same z-value]. Now something strange happens. One point is perfectly plotted, the other isn't shown at all. It only appears if I move it upwards in z-direction by adding a positive number. So somehow it disappears in the wireframe-surface

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,

I'm trying to make a set of contour plots of bivariate kernel density estimates, showing three such plots overlaid, similar to this plot http://euclid.psych.yorku.ca/SCS/Private/Test/ridge-boot2.pdf except that I would like to have the contours *filled* (using transparent colors). To make this reproducible, I've saved the results of KernSmooth::bkde2D() in the following file:

Hello, Is it possible to draw level lines on the surface of a bivariate function? In the following example, to draw surface and levels lines for a multivariate normal law, I use persp, trans3d, contourLines and lines, but if the lines are correctly drawn, some parts of them are, of course, visible even if they are drawn on a non visible "face". Any suggestion to avoid this problem

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 :

I have looked through the documentation and have not been able to find a way of using an xlim, ylim, or zlim type option on rgl.surface. I know that persp3d has the option, but seems to only be able to expand the axes not reduce them. Is there anyone who has an idea of how to do this? Thank you for your time. [[alternative HTML version deleted]]

Hello list, I just started R today and tried something quite simple. I wanted to create a colored plot and eventually after hours of fiddling around got it working. However, my solution seems very suboptimal and I'd really appreciate your hints on how to improve. I believe that R already offers many functions I coded (e.g. distance between two vectors, vector length, vector normalization and

Greetings all! I'm wanting to plot contours of a function, and I don't want to evaluate it at a dense grid of points (i.e. I don't want a huge array of values). Say I have a vector of x-values such as x <- 0.1*c(1:10), and the same for y <- 0.1*(0:10). I then evaluate a 10*10 matrix z of values of z = f(x,y). so I can then do CL <- contourLines(x,y,z) and get a list of

I am working on plotting a response surface using wireframe(). The default style/orientation is z | | y | \ | \ | \ | \ | \ | \ | \ | \|________________x 0 Now what I want the orientation of axes is: z | | | |

Hi, I want to draw a figure similar to http://dsarkar.fhcrc.org/lattice/book/images/Figure_13_07_stdBW.png from http://lmdvr.r-forge.r-project.org/figures/figures.html (figure 13.7) . However instead of using a contour plot as a panel for the wireframe I want to use a levelplot. Can somebody help me with this? Thanks in advance for any help, Kind regards, Anton Bossenbroek

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 *