Re: [R] Fwd: Re: How to create following chart for visualizing multivariate time series

From: jim holtman <jholtman_at_gmail.com>
Date: Sat, 01 Mar 2008 19:21:43 -0500

If you want color, then a slight addition to Henrique's solution will do it:

x <- y <- seq(-4*pi, 4*pi, len=27)
r <- sqrt(outer(x^2, y^2, "+"))
library(lattice)
colkey <- colorRampPalette(c('red','yellow','green'))(32) levelplot(r, colorkey=list(col=colkey),

col.regions=(col=colkey))

On Sat, Mar 1, 2008 at 6:08 PM, Henrique Dallazuanna <wwwhsd_at_gmail.com> wrote:
> This works for me:
>
> x <- y <- seq(-4*pi, 4*pi, len=27)
> r <- sqrt(outer(x^2, y^2, "+"))
> library(lattice)
> levelplot(r, colorkey=list(col=gray((0:32)/32)),
> col.regions=(col=gray((0:32)/32)))
>
> 'r' is a matrix for you?
>
> On 01/03/2008, David Winsemius <dwinsemius_at_comcast.net> wrote:
> > "Henrique Dallazuanna" <wwwhsd_at_gmail.com> wrote in
> >
> > library(lattice)
> > levelplot(r, colorkey=list(col=gray((0:32)/32)),
> > col.regions=(col=gray((0:32)/32)))
> >
> > When I try that example, I get an error, even after updating lattice.
> >
> > > levelplot(r, colorkey=list(col=gray((0:32)/32)),
> > + col.regions=(col=gray((0:32)/32)))
> > Error in UseMethod("levelplot") : no applicable method for "levelplot"
> >
> > If I simply change colorkey=FALSE to colorkey=TRUE in the first levelplot
> > help page example, I have what looks to me as "success".
> >
> > levelplot(z~x*y, grid, cuts = 50, scales=list(log="e"), xlab="",
> > ylab="", main="Weird Function", sub="with log scales",
> > colorkey = TRUE,
> > region = TRUE)
> >
> > --
> > David Winsemius
> >
> > >
> > > On 29/02/2008, Megh Dal <megh700004_at_yahoo.com> wrote:
> > >> Hi Jim, i think you could not get my point. I did not want to put
> > >> red-blue color there. I want to put a pallet which will describe
> > >> the values of r. please have a look on following :
> > > http://bp0.blogger.com/_k3l6qPzizGs/RvDVglPknRI/AAAAAAAAAKo/itlWOvuuO
> > > tI/s1600-h/pairwise_kl_window60.png. Please see how a color pallate
> > > is added on the right side of this plot describing the value of red
> > > color, value of blue color etc.
> > >>
> > >> Is there any solution?
> > >>
> > >> Regards,
> > >>
> > >>
> > >> jim holtman <jholtman_at_gmail.com> wrote:
> > >> Try something like this:
> > >>
> > >> require(grDevices) # for colours
> > >> x <- y <- seq(-4*pi, 4*pi, len=27)
> > >> r <- sqrt(outer(x^2, y^2, "+"))
> > >> image(x, y, r, col=gray((0:32)/32))
> > >> colors <- colorRampPalette(c('red', 'yellow', 'blue')) # create
> > >> you color spectrum
> > >> image(x,y,r, col=colors(100))
> > >>
> > >>
> > >>
> > >> On Thu, Feb 28, 2008 at 9:28 PM, Megh Dal wrote:
> > >> > I used ?image function to do that, like below :
> > >> >
> > >> > require(grDevices) # for colours
> > >> > x <- y <- seq(-4*pi, 4*pi, len=27)
> > >> > r <- sqrt(outer(x^2, y^2, "+"))
> > >> > image(x, y, r, col=gray((0:32)/32))
> > >> >
> > >> > However my next problem to add a color pallet for color
> > >> > description [as shown in following link]. If anyone here tell
> > >>>> me how to do that, it will be good for me.
> > >> >
> > >>
> > >> > Megh Dal wrote: Hi all,
> > >>
> > >> >
> > >> > Can anyone here please tell me whether is it possible to produce
> > >> > a chart displayed in http://www.datawolf.blogspot.com/ in R for
> > > visualizing multivariate time series? If possible how?
> > >> >
> > >> >
> >
> > ______________________________________________
> > R-help_at_r-project.org mailing list
> > https://stat.ethz.ch/mailman/listinfo/r-help
> > and provide commented, minimal, self-contained, reproducible code.
> >
>
>
> --
> Henrique Dallazuanna
> Curitiba-Paraná-Brasil
> 25° 25' 40" S 49° 16' 22" O
>
> ______________________________________________
> R-help_at_r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> and provide commented, minimal, self-contained, reproducible code.
>

```--
Jim Holtman
Cincinnati, OH
+1 513 646 9390

What is the problem you are trying to solve?  Tell me what you want to
do, not how you want to do it.

______________________________________________
R-help_at_r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help