很多朋友需要将多个时序图画在一起,方便比较。
[attachment=208]
代码如下:
<br />
<br />
library(IDPmisc)<br />
library(SwissAir) # data for the example<br />
<br />
## prepare the data<br />
<br />
st <- 6.5*30*48<br />
x.at <- seq(st,nrow(AirQual),48)<br />
<br />
## draw graph<br />
<br />
longtsPlot(<br />
y1 = AirQual[,c("ad.O3","ad.NOx")],<br />
y2 = AirQual$ad.T,<br />
names1 = c("O3","NOx"), names2 = "Temp",<br />
startP = st, upf = 3*48, fpp = 3,<br />
x.at = x.at, x.ann = substr(AirQual$start,1,6)[x.at],<br />
x.tick = seq(st,nrow(AirQual),12),<br />
y1.at = c(0,100), y1.tick = seq(0,150,50),<br />
y2.at = c(10,30), y2.tick = seq(10,30,10),<br />
y1lab="[ppb]", y2lab="[C]",<br />
y1lim = c(0,140), y2lim = c(10,30), xpd = TRUE,<br />
col2 = "red", type = "l",<br />
mar=c(1.2,3.2,1.2,3.2),main="2004-07-13 till 2004-07-22"<br />
)<br />
来自:
http://addictedtor.free.fr/graphiques/RGraphGallery.php?graph=160