It seems the code I've sent had typos... Here's a corrected version:
################################# x <- sapply(1:10, function(i)rnorm(1000)) f <- sapply(1:10, function(i)density(x[,i], from=-5,to=5)$y) grid <- density(x[,1], from=-5,to=5)$x win.graph() persp(grid, 1:10, f,theta=-50, phi=30, d=2) win.graph() opar <- par(mfrow=c(5,2), mar=c(2,2,1,1)) sapply(1:10, function(i)plot(grid, f[,i], ann=FALSE, type="l")) par(opar) ################################# Sorry for the mistake. Best regards, Eduardo On Mon, Feb 21, 2011 at 11:18 PM, Eduardo de Oliveira Horta <eduardo.oliveiraho...@gmail.com> wrote: > Hello, > > I'm willing to plot a sequence of densities on a 3d graph, something like > > ----------------------------------------------------------------- > x <- sapply(1:10, function(i)rnorm(1000)) > > f <- sapply(1:10, function(i)density(x[,i], from=-5,to=5)$y) > grid <- density(x[,1], from=-5,to=5)$x > > win.graph() > persp(grid1, 1:10, f,theta=-50, phi=30, d=2) > ----------------------------------------------------------------- > > That is, I simply want to stack the curves > > ----------------------------------------------------------------- > win.graph() > opar <- par(mfrow=c(5,2), mar=c(2,2,1,1)) > sapply(1:10, function(i)plot(grid, f[,i], ann=FALSE, type="l")) > par(opar) > ----------------------------------------------------------------- > > into a functional time series plot (unlike the example above, I'm > working with pdf's which are not independent). > > The result I get using persp() is 'almost' what I want, but I would > like to avoid the wireframes connecting the densities along 1:10, > because there's no continuity in that direction. In other words, I > wanted that the only lines appearing in the plot to be those > corresponding to the pdf's f[,1],...f[,10]. > > Thanks once again, and best regards, > > Eduardo > >> sessionInfo() > R version 2.11.1 (2010-05-31) > i386-pc-mingw32 > > locale: > [1] LC_COLLATE=Portuguese_Brazil.1252 LC_CTYPE=Portuguese_Brazil.1252 > [3] LC_MONETARY=Portuguese_Brazil.1252 LC_NUMERIC=C > [5] LC_TIME=Portuguese_Brazil.1252 > > attached base packages: > [1] stats graphics grDevices utils datasets methods base > > other attached packages: > [1] rgl_0.92.798 Revobase_4.2.0 RevoScaleR_1.1-1 lattice_0.19-13 > > loaded via a namespace (and not attached): > [1] grid_2.11.1 pkgXMLBuilder_1.0 revoIpe_1.0 tools_2.11.1 > [5] XML_3.1-0 > > ______________________________________________ > 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.
