Is this what you are trying to do? x <- rnorm(25) y <- rnorm(25) xd <- density(x, n=25, from=-3.5, to=3.5) yd <- density(x, n=25, from=-3.5, to=3.5) z <- outer(xd$y, yd$y) persp(xd$x, yd$x, z)
The density function gives x values at which the density is estimated, so that gives you the axes for the plot, but the z value must be computed from the y values that density returns. This example multiplies the two y values to get z. ---------------------------------------------- David L Carlson Associate Professor of Anthropology Texas A&M University College Station, TX 77843-4352 > -----Original Message----- > From: r-help-boun...@r-project.org [mailto:r-help-bounces@r- > project.org] On Behalf Of Maximilian Lklweryc > Sent: Wednesday, November 28, 2012 10:07 AM > To: r-help@r-project.org > Subject: [R] Plot 3d density > > I want to create a 3d plot with densities. > > I use the function density to first create a 2d dimensional plot for > specific x values, the function then creates the density and puts them > into a y variable. Now I have a second set of x values and put it > again into the density function and I get a second set of y variables > and so on.... I want to put those sets into a 3d plot, I hope you know > what I mean. So I have a surface of densities.... > > E.g. I have: > > x1<-c(1:10) > x2<-c(2:11) > y1<-c(1,1,2,1,3,4,2,3,2,2) > y2<-c(1,2,3,1,3,6,2,8,2,2) > . > . > . > . > > Now I want to put on the x axis for the value 1 the first set that > means for x=1 on the y axis of the 3d plot the x values from the first > set and on the z axis the densities. So I have a "disk" for x=1, for > x=2 I have the second "disk" and so on, so I get a density "mountain". > > I hope I am understandable, if you have a better idea to realize it > then you are welcome! I want to do it with the persp function, would > be nice if you make an example with that function, > > The function kde2d does not help. > > Thanks a lot for your help. > > ______________________________________________ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting- > guide.html > and provide commented, minimal, self-contained, reproducible code. ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
