Wykreślanie dwóch obiektów xts

Question

Używam xtsExtra do wykreślania dwóch obiektów XTS.

Rozważmy następujące wywołanie plot.xts:

plot.xts(merge(a,b),screens=c(1,2)) 

który służy do działki XT obiektów a i b w dwóch odrębnych paneli.

Jak kontrolować odstępy osi y? W szczególności mam problem, w którym etykiety osi y są zbyt blisko lub nawet zachodzą na siebie.

Idealnie, chciałbym określić minimalne wypełnienie, które ma być utrzymywane pomiędzy dwiema etykietami osi Y. Każda pomoc jest doceniana!

EDIT: Powtarzalny przykład:

#install if needed 
#install.packages("xtsExtra", repos="http://R-Forge.R-project.org") 
library(xtsExtra) 

ab=structure(c(-1, 0.579760106421202, -0.693649703427259, 0.0960078627769613, 
0.829770469089809, -0.804276208608663, 0.72574639798749, 0.977165659135716, 
-0.880178529686181, -0.662078620277974, -1, 2.35268982675599, 
-0.673979231663719, 0.0673890875594205, 1.46584597734824, 0.38403707067242, 
-1.53638088345349, 0.868743976582955, -1.8394614923913, 0.246736581314485 
), .Dim = c(10L, 2L), .Dimnames = list(NULL, c("a", "b")), index = structure(c(1354683600, 
1354770000, 1354856400, 1354942800, 1355029200, 1355115600, 1355202000, 
1355288400, 1355374800, 1355461200), tzone = "", tclass = "Date"), class = c("xts", 
"zoo"), .indexCLASS = "Date", .indexTZ = "", tclass = "Date", tzone = "") 

plot.xts(ab,screens=c(1,2)) 

która produkuje:

Answer 1

Niestety to trwało tak długo. Próbowałem ustalić, dlaczego mój wykres zaczyna się w dniu 04 grudnia 2012 r. I kończy się 13 grudnia 2012 r., Kiedy zaczyna się od 05 grudnia 2012 r., A kończy się 14 grudnia 2012 r. Czy możesz sprawdzić, czy opublikowany powyżej ab jest taki sam ab użyty do wykreślenia wykresu?

Ponadto użyłem biblioteki xts zamiast xtsExtra. Czy istnieje powód do korzystania z xtsExtra?

Oto kod:

library(xts) 

ab=structure(c(-1, 0.579760106421202, -0.693649703427259, 0.0960078627769613, 
      0.829770469089809, -0.804276208608663, 0.72574639798749, 0.977165659135716, 
      -0.880178529686181, -0.662078620277974, -1, 2.35268982675599, 
      -0.673979231663719, 0.0673890875594205, 1.46584597734824, 0.38403707067242, 
      -1.53638088345349, 0.868743976582955, -1.8394614923913, 0.246736581314485), .Dim = c(10L, 2L), .Dimnames = list(NULL, c("a", "b")), index = structure(c(1354683600, 
      1354770000, 1354856400, 1354942800, 1355029200, 1355115600, 1355202000, 
      1355288400, 1355374800, 1355461200), tzone = "", tclass = "Date"), class = c("xts", 
      "zoo"), .indexCLASS = "Date", .indexTZ = "", tclass = "Date", tzone = "") 

#Set up the plot area so that multiple graphs can be crammed together 
#In the "par()" statement below, the "mar=c(0.3, 0, 0, 0)" part is used to change 
#the spacing between the graphs. "mar=c(0, 0, 0, 0)" is zero spacing. 
par(pty="m", plt=c(0.1, 0.9, 0.1, 0.9), omd=c(0.1, 0.9, 0.2, 0.9), mar=c(0.3, 0, 0, 0)) 

#Set the area up for 2 plots 
par(mfrow = c(2, 1)) 

#Build the x values so that plot() can be used, allowing more control over the format 
xval <- index(ab) 

#Plot the top graph with nothing in it ========================= 
plot(x=xval, y=ab$a, type="n", xaxt="n", yaxt="n", main="", xlab="", ylab="") 
mtext(text="ab", side=3, font=2, line=0.5, cex=1.5) 

#Store the x-axis data of the top plot so it can be used on the other graphs 
pardat <- par() 

#Layout the x axis tick marks 
xaxisdat <- index(ab) 

#If you want the default plot tick mark locations, un-comment the following calculation 
#xaxisdat <- seq(pardat$xaxp[1], pardat$xaxp[2], (pardat$xaxp[2]-pardat$xaxp[1])/pardat$xaxp[3]) 

#Get the y-axis data and add the lines and label 
yaxisdat <- seq(pardat$yaxp[1], pardat$yaxp[2], (pardat$yaxp[2]-pardat$yaxp[1])/pardat$yaxp[3]) 
axis(side=2, at=yaxisdat, las=2, padj=0.5, cex.axis=0.8, hadj=0.5, tcl=-0.3) 
abline(v=xaxisdat, col="lightgray") 
abline(h=yaxisdat, col="lightgray") 
mtext(text="ab$a", side=2, line=2.3) 
lines(x=xval, y=ab$a, col="red") 
box() #Draw an outline to make sure that any overlapping abline(v)'s or abline(h)'s are covered 

#Plot the 2nd graph with nothing in it ================================ 
plot(x=xval, y=ab$b, type="n", xaxt="n", yaxt="n", main="", xlab="", ylab="") 

#Get the y-axis data and add the lines and label 
pardat <- par() 
yaxisdat <- seq(pardat$yaxp[1], pardat$yaxp[2], (pardat$yaxp[2]-pardat$yaxp[1])/pardat$yaxp[3]) 
axis(side=2, at=yaxisdat, las=2, padj=0.5, cex.axis=0.8, hadj=0.5, tcl=-0.3) 
abline(v=xaxisdat, col="lightgray") 
abline(h=yaxisdat, col="lightgray") 
mtext(text="ab$b", side=2, line=2.3) 
lines(x=xval, y=ab$b, col="blue") 
box() #Draw an outline to make sure that any overlapping abline(v)'s or abline(h)'s are covered 

#Plot the X axis ================================================= 
axis(side=1, label=format(as.Date(xaxisdat), "%b %d\n%Y\n") , at=xaxisdat, padj=0.4, cex.axis=0.8, hadj=0.5, tcl=-0.3) 
mtext(text="Date", side=1, line=2.5) 

Answer 2

Gram z niektórych parametrów

plot.xts(ab, bty = "n", las = 1,  cex.axis = 0.5) 

Answer 3

dałem się wykreślenie z XTS, ze względu na słabą formatowania w oś y. Dokładnie ten sam problem, o którym mowa w tym pytaniu.

autoplot.zoo było miłą alternatywą.

library(ggplot2) 
autoplot.zoo(ab) + theme(panel.background = element_blank(), 
         panel.grid.major = element_blank(), 
         panel.grid.minor = element_blank()) + 
         theme_bw()