If you want to use abline on a multivariate plot it must be issued within a panel function like this:
# same set up as in prior email p <- function(x, y, ...) { points(x, y); lines(x, y); abline(v = time(s)[s > 0], col = "green") } plot(cbind(x, s), panel = p) There are further examples of panel functions in the examples section of ?plot.zoo in case what you want is a variation of the above. On Tue, Mar 3, 2009 at 9:04 AM, Sergey Goriatchev <serg...@gmail.com> wrote: > Gabor, yes, I want to color portions of EACH plot of the MULTIPLE plot > done with plot.zoo() > I tried to do: > plot(multivariate zoo object) > abline(v=...) > > but that does not work. > > I will check your suggestions of the examples. > Thank you for your help, as always! > > Best, > Sergey > > On Tue, Mar 3, 2009 at 14:44, Gabor Grothendieck > <ggrothendi...@gmail.com> wrote: >> Do you mean you want to shade a portion of the plot? >> >> There are two examples of that in the examples section of ?plot.zoo >> and a further example using xyplot.zoo in the examples section of >> ?xyplot.zoo >> >> On Tue, Mar 3, 2009 at 8:37 AM, Sergey Goriatchev <serg...@gmail.com> wrote: >>> Hi, Gabor >>> >>> No, what I am trying to do is similar to: >>> abline(v=time(spread)[spread[,"Indicator"]==(-1)], col="yellow"), >>> >>> where spread is the multivariate zoo object (say, 5 timeseries). >>> >>> That is, I want to color parts of the plots where indicator==(-1), but >>> do the coloring >>> without using layout() and then repeating plot() and abline() for each >>> of the 5 timeseries. >>> >>> Best, >>> Sergey >>> >>> >>> On Tue, Mar 3, 2009 at 14:22, Gabor Grothendieck >>> <ggrothendi...@gmail.com> wrote: >>>> Are you trying to color the points themselves? This plots the >>>> first two series in frame 1 (they are the same but one is plotted >>>> as points and the other as a line) and the third series is shown >>>> in frame 2 and for the series of points it colors them green or red. >>>> The lines are all colored black: >>>> >>>> library(zoo) >>>> set.seed(1) >>>> x <- zoo(rnorm(10)) >>>> s <- sign(x) >>>> >>>> plot(cbind(x, x, s), screen = c(1, 1, 2), type = c("p", "l", "l"), >>>> col = list(ifelse(s > 0, "green", "red"), 1, 1), pch = 20) >>>> >>>> >>>> >>>> On Tue, Mar 3, 2009 at 7:48 AM, Sergey Goriatchev <serg...@gmail.com> >>>> wrote: >>>>> Hello, everyone >>>>> >>>>> I have a zoo object containing several time series of daily frequency. >>>>> One of these timeseries is an indicator function with value (-1) at >>>>> certain times, and (+1) at the other. >>>>> I do a plot of several of the timeseries in one go (a multiple plot). >>>>> I wonder if I can automatically in EACH plot color the area where >>>>> indicator variable is (-1)? >>>>> >>>>> Of course, I could do it simply with layout() and then for each >>>>> timeseries do plot and a color overlay, but I wonder if with >>>>> plot.zoo() simething similar is possible to do automatically. >>>>> >>>>> Thanks in advance for help! >>>>> >>>>> Best, >>>>> Sergey >>>>> >>>>> ______________________________________________ >>>>> 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. >>>>> >>>> >>> >>> >>> >>> -- >>> I'm not young enough to know everything. /Oscar Wilde >>> Experience is one thing you can't get for nothing. /Oscar Wilde >>> When you are finished changing, you're finished. /Benjamin Franklin >>> Tell me and I forget, teach me and I remember, involve me and I learn. >>> /Benjamin Franklin >>> Luck is where preparation meets opportunity. /George Patten >>> >> > > > > -- > I'm not young enough to know everything. /Oscar Wilde > Experience is one thing you can't get for nothing. /Oscar Wilde > When you are finished changing, you're finished. /Benjamin Franklin > Tell me and I forget, teach me and I remember, involve me and I learn. > /Benjamin Franklin > Luck is where preparation meets opportunity. /George Patten > ______________________________________________ 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.