Logo Questions Linux Laravel Mysql Ubuntu Git Menu
 

How to plot a set of densities in 3D using R?

Tags:

plot

r

3d

I need to plot, in 3D, a set of densities associated to a time series. More precisely, I would like to be able in R to build an image close to this example alt text

This image is taken from [1]. The transparency plays an important role as let us see the trajectory of the "measures" in the x-y plane.

Any help will be greatly appreciated.

[1]: Juban and Kariniotakis, "Uncertainty Estimation of Wind Power Forecasts", presentation at EWEC 2008 - 01 April - Brussels, Belgium. (I can't post the link, google will help interested readers).

like image 212
José Luis Avatar asked Jul 21 '10 08:07

José Luis


2 Answers

In 1996 I wrote a paper (published in JCGS) with a figure very similar to that but without the transparency. See http://robjhyndman.com/papers/estimating-and-visualizing-conditional-densities/ for the details. The plotting function is implemented in the R package hdrcde available on CRAN. The package contains some examples in the help files. You should be able to adapt my code to add the transparency.

like image 142
Rob Hyndman Avatar answered Sep 21 '22 23:09

Rob Hyndman


This is how far I got thanks to Rob's hint. I used persp() to create an empty plot and added polygons and lines to it:

alt text

However, it is not as pretty as the original one... :(

like image 23
José Luis Avatar answered Sep 18 '22 23:09

José Luis