similar to: Persp plotting of kernel density estimate.

Dear All, I am trying to visualize the surface of a bivariate kernel density estimate. I have a vector of bivariate observations(x,y), and a function which computes the kernel density estimate z corresponding to each observation. I cannot generate the (x,y) data in the ascending order needed by persp(x,y,z). I was wondering whether there is an R version of the S function interp. Would anybody

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

Dear R help, I am fitting a gam model with package GAM ( I use lo() to fit 2 non-parametric additive components). When I plot the GAM object, along with the graphics window displaying plots of each components, I read the following warnings at the prompt: 3: the condition has length 1 and only the first element will be used in: if (scale2 < scale) rep(mean(ylim), 2) + ((ylim - mean(ylim))

Hi, I am trying to plot densities given on a two dimensional grid. My data is in the an external file, and is arranged in three columns: x, y, density how may i get a plot of this? i would like to get (1) a three dimensional plot and (2) a color coded two dimensional plot. I have tried using image(x, y, density) but i am asked to put the data in ascending order. i am not sure how i may

Hello, I have a question about persp function: I made my classical matrix with x, y and z variables and I dont know why I obtain a 3D image with overestimate heights. Here is my script: microtopo = read.csv2("topo2.csv", dec=".") attach(microtopo) require(akima) # for persp or zz <-interp (x, y, z) persp(zz) See also my attach file Thanks, Best regards JS

Hi, This is probably documented, but I cannot find the right words or expression for a search. My attempts failed. I have a data frame of 3 vectors (x, y and z) and would like to transform this so that I could use persp. Presently I have y-level copies of each x level, and a z value for each x-y pair. I need 2 columns giving the possible levels of x and y, and then a transformation of

Hi, I am doing my way on R, with much experimenting. So, I am trying to plot a 3d graphic ussing "persp" and entering a data set (attached) of UTM coordinates as x,y and a pH values as z. However when I try an error message comes out telling that increasing x and y values are expected. Jus ordering the first vector does not help, and, of course, order the first two independently

Hello, I have a problem in drawing a 3d graphic. I simplified the program to show you the problem. I have a text file with hundreds of entries in 3 columns. I scaned these columns with matrix(scan). Then, I transformed some datas. Now I have the datas in h,x,y. I created a new matrix to use persp . But I got an error. The reason is probably the same length of x,y,z. But what can I do the draw

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

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

Dear All, When bootstrapping a statistics based on more than one vector, from a data.frame or a matrix object, it looks like I am not able to pass the data to R. What am I doing wrong? I use the library "bootstrap". Here is an example with a data.frame called "data" "boot2_bootstrap(data, theta, nboot) I get the following error message: Error in inherits(x,

Dear All, I'm trying to maximize a likelihood with respect one parameter using "optimize" on simulated data (without error component). I've iterated the maximization procedure 1000 times and I should always obtain the same estimate of the parameter (equal to the simulated one) but, instead, i obtain different results (the likelihood function shouldn't be flat). Does

Am trying to do some study on bivariate binomial distribution. Anyone knows if there is package in R that I can use to calculate the density function of bivariate binomial distribution and to generate random samples of it. Thanks, -- View this message in context: http://n4.nabble.com/How-to-calculate-density-function-of-Bivariate-binomial-distribution-tp992002p992002.html Sent from the R help

Hi, I would like to make a contour plot with the data below, x-axis= Day (increasing order, day number is out of 365 days of a year), y-axis= Depth (varies), and z=pH. I tried to use filled.contour function, but not sure the proper way to ordinate my z results into a matrix with different x and y length. After that, if I attempted to plot the data, the x and y values are not in ascending

What package/code could I use to create a 3-d surface graph of the predicted values over two of the explanatory variables' coefficients? Jim James Bang Department of Economics University of Illinois Well I AM missing the back of my head.you COULD cut me a little slack! -Homer Simpson

Dear R users: I have a dataset with two variables (>20000 observations, two samples from same subject) and I used "kernSur" from library(Genkern) to get a estimated bivariate density and corresponding plots as follows: new.data.normal<-data.normal[!is.na(data.normal[,2]),] x<-new.data.normal[,2] y<-new.data.normal[,3] op <- KernSur(x,y, xgridsize=50, ygridsize=50,

Dear R users, I have created a Loess surface in R, in which x is relative longitude by miles, y is relative latitude by miles, and z is population density at the neighborhood level. The purpose is to identify some population centers in the region. I'm wondering if there is a way to determine the coordinates (x,y) of each center, so I can know exactly where they are. Let me use the

Dear R users, I would like to plot the bivariate density surface with its marginal densities on the sides of the 3D box, just like in the picture I attach. I tried to found information about how to do it but did not find anything. Does anyone know how to do it? Thanks in advance, Eduardo.

dear R experts: I am trying to plot an empirical likelihood function in 3d. The values are not over a regular grid---I just searched the likelihood function to find the optimal value, and then computed a few values around it. (each point in the likelihood function takes a very long time to compute.) the likelihood values now sit in a csv file that has three columns, "mu",

Hello list! I have a problem with the function 'persp'. I have a data set with the variables X -> The X coordenate of the insect Y -> The Y coordenate of the insect deep-> how deep the insect buried into the soil. I would like to have a plot with a surface simulating the way the insect buried. Is it possible? If I type persp(x,y,deep) R gives me this error message Error in