- You don't need a panel function since there is only one.
- There should not be more than 2 components on lty and lwd
since tab has two columns and that is the cause of the
warning message you are seeing.
- with the recent version of zoo plot.type uses screen so
you can write the shorter screen = 1 instead of the lengthier
plot.type = "single"
Redone with an abline at the end to give the grid lines gives:
plot(tab, screen = 1, lty = c(1,2), lwd = c(8,2),
col = c("darkorange", "red2"),
xlab = " ", ylab = "euro/MWh", xaxt = "n")
fmt <- "%b-%Y" # format for axis labels
axis(1, at = time(tab), lab = FALSE)
ix <- seq(1, length(tab), 3)
labs <- format(time(tab), fmt)
axis(1, at = time(tab)[ix], lab = labs[ix], tcl = -0.7, cex.axis = 0.9)
abline(v = time(tab)[ix], lty = 2, col = "grey")
On Dec 9, 2007 12:33 PM, vittorio <vdemart1 at tin.it>
wrote:> I have the following code
>
> ####################################################################
>
> library(zoo)
>
> miedate <- yearmon((2006)+seq(0,23)/12)
> tab <- zoo(cbind(A = c(79.47, 89.13, 84.86, 75.68, 72.82, 78.87, 93.46,
> 78.18, 82.46, 77.25, 80.95, 84.39, 81.7, 74.76, 65.29, 60.3,
> 66.59, 73.19, 92.39, 65.76, 77.45, 74.22, 101.36, 100.01), B = c(77.95,
> 76.73, 51.2, 51.86, 51.29, 49.45, 53.88, 47.96, 55.07, 45.34,
> 37.07, 37.53, 47.79, 37.5, 30.35, 37.78, 34.13, 39.14, 39.89,
> 35.46, 36.54, 38.39, 47.33, 45.34)),miedate)
>
> par(mai=c(0.75, 0.75, 0.1, 0.1))
> my.panel <- function(...) {
> fmt <- "%b-%Y" # format for axis labels
> lines(...)
> panel.number <- parent.frame()$panel.number
> # if bottom panel
> if (!length(panel.number) || panel.number == NCOL(tab)) {
> # next line only if non-labelled ticks wanted for each point
> axis(1, at = time(tab), lab = FALSE)
> ix <- seq(1, length(tab), 3)
> labs <- format(time(tab), fmt)
> axis(1, at = time(tab)[ix], lab = labs[ix], tcl = -0.7, cex.axis =
0.9)
> }
> }
>
> plot.zoo(tab, plot.type="single", panel=my.panel,
lty=c(1,2,3,4,5,1),
> lwd=c(8,2,2,2,2,8),
> col=c("darkorange",
"red2","darkblue","darkblue","green","midnightblue"),
> xlab=" ",ylab =
"euro/MWh",xaxt="n")
>
> ####################################################################
>
> This code produces a graph with a non-decimal x-axis but by far more
(human)
> readable because ticks corresponds to months exactly (and not to a weird
> decimal ticking!).
> I would like to add a vertical grid with the same criterium.
>
> How can I obtain this?
>
> Ciao
> Vittorio
>
> ______________________________________________
> R-help at 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.
>