怎样生成这样的时间序列图？

1. rm(list=ls())
   
2. par(mfrow=c(3,2),mar=c(0,3.5,0,2), oma = c(3, 0, 2, 0))
   
3. data=read.csv("x1.csv",header=T)
   
4. attach(data)
   
5. plot(year,y1,xaxt="n", xlab = "", type = "l")
   
6. plot(year,y2,xaxt="n", xlab = "", type = "l")
   
7. plot(year,y3,xaxt="n", xlab = "", type = "l")
   
8. plot(year,y4,xaxt="n", xlab = "", type = "l")
   
9. plot(year,y5, xlab = "", type = "l")
   
10. plot(year,y7, xlab = "", type = "l")

复制代码

呵呵，不知是否满意

另外请问下版主，
1.怎么控制y轴标题到y轴的距离呢？用哪个参数？
2.如何控制x轴坐标名称显示？比如还想显示2001、2003等年份的坐标

第二个是不是使用mtext解决呢？

kk22boy 发表于 2012-3-14 17:05
呵呵，不知是否满意

太好了，但能否让y1等离纵轴近一些？或者让他们大一些？
因为看不清楚是哪个模型的了

多谢各位朋友！

个个是精英，受益匪浅

#读入数据
ds <- read.csv("c:/x1.csv", header = T)

#将宽数据变成长数据
my.melt <- melt(ds,id="year",measure=c("y1","y2","y3","y4","y5","y7"))

#作图
qplot(year,value,data=my.melt,geom="line")+facet_grid(variable~.,scales="free_y")

呵呵，还是ggplot2最漂亮

修改ylab的大小
data=read.csv('C:\\Users\\acer\\Desktop\\x1.csv',head=T)
par(mfrow=c(3,2),mar=c(0,5,0,2), oma = c(3, 0, 3, 0),cex.axis=1,cex.lab=2)
attach(data)
plot(y1,xaxt="n", xlab = "", type = "l")
plot(y2,xaxt="n", xlab = "", type = "l")
plot(y3,xaxt="n", xlab = "", type = "l")
plot(y4,xaxt="n", xlab = "", type = "l")
plot(y5, axes=F,type='l',xlab='year')
axis(1, 1:12, year-2000)
axis(2)
box()
plot(y7, axes = F, type = "l",xlab='year')
axis(1, 1:12, year-2000)
axis(2)
box()
text(-1,-1,'year')
mtext("Big Title", outer = TRUE, cex = 1,xlab='year')

吉林小王子 发表于 2012-3-14 21:44
修改ylab的大小
data=read.csv('C:\\Users\\acer\\Desktop\\x1.csv',head=T)
par(mfrow=c(3,2),mar=c(0 ...

怎样才能让y轴标题离y轴近一点呢
另外，小王子你发现没有x轴刻度10没有哦我运行了你的代码也没出现10，这是啥子问题？
呵呵

kk22boy 发表于 2012-3-14 22:08
怎样才能让y轴标题离y轴近一点呢
另外，小王子你发现没有x轴刻度10没有哦我运行了你的代码也没出现10，这 ...

你看上面那个图，稍微拉大一点就可以了，因为数据太多，图小无法显示完整，至于如何让ylab离的近一点，我还没有搞清楚，再等等[loveliness]