dat <- read.csv("x1.csv", header = T, sep = ",")
attach(dat)
x <- cbind(c(1,2,3), c(4,5,6))
layout(x, c(1,1), c(4,3,4))

par(mar=c(0,5,4,2), oma = c(0, 0, 3, 0))
plot(year, y1, xaxt = "n", xlab = "", type = "l")
par(mar=c(0,5,0,2))
plot(year, y2, xaxt = "n", xlab = "", type = "l")
par(mar=c(5,5,0,2))
plot(year, y3, xlab = "Time", type = "l")
par(mar=c(0,5,4,2))
plot(year, y4, xaxt = "n", xlab = "", type = "l")
par(mar=c(0,5,0,2))
plot(year, y5, xaxt = "n", xlab = "", type = "l")
par(mar=c(5,5,0,2))
plot(year, y7, xlab = "Time", type = "l")
mtext("Big Title", outer = TRUE, cex = 1)

直接运行的错误是：
> plot(year, y1, xaxt = "n", xlab = "", type = "l")
错误于as.POSIXlt.numeric(x) : 'origin' must be supplied
> par(mar=c(0,5,0,2))
> plot(year, y2, xaxt = "n", xlab = "", type = "l")
错误于as.POSIXlt.numeric(x) : 'origin' must be supplied
> par(mar=c(5,5,0,2))
> plot(year, y3, xlab = "Time", type = "l")
错误于as.POSIXlt.numeric(x) : 'origin' must be supplied
> par(mar=c(0,5,4,2))
> plot(year, y4, xaxt = "n", xlab = "", type = "l")
错误于as.POSIXlt.numeric(x) : 'origin' must be supplied
> par(mar=c(0,5,0,2))
> plot(year, y5, xaxt = "n", xlab = "", type = "l")
错误于as.POSIXlt.numeric(x) : 'origin' must be supplied
> par(mar=c(5,5,0,2))
> plot(year, y7, xlab = "Time", type = "l")
错误于as.POSIXlt.numeric(x) : 'origin' must be supplied
> mtext("Big Title", outer = TRUE, cex = 1)
错误于mtext("Big Title", outer = TRUE, cex = 1) : 还没有调用plot.new

修改程序：用dat$year代替了year,
新的提示是：
> par(mar=c(0,5,4,2), oma = c(0, 0, 3, 0))
> plot(dat$year, y1, xaxt = "n", xlab = "", type = "l")
错误于xy.coords(x, y, xlabel, ylabel, log) :
  'x' and 'y' lengths differ
> par(mar=c(0,5,0,2))
> plot(dat$year, y2, xaxt = "n", xlab = "", type = "l")
错误于xy.coords(x, y, xlabel, ylabel, log) :
  'x' and 'y' lengths differ
> par(mar=c(5,5,0,2))
> plot(dat$year, y3, xlab = "Time", type = "l")
错误于xy.coords(x, y, xlabel, ylabel, log) :
  'x' and 'y' lengths differ
> par(mar=c(0,5,4,2))
> plot(dat$year, y4, xaxt = "n", xlab = "", type = "l")
错误于xy.coords(x, y, xlabel, ylabel, log) :
  'x' and 'y' lengths differ
> par(mar=c(0,5,0,2))
> plot(dat$year, y5, xaxt = "n", xlab = "", type = "l")
错误于xy.coords(x, y, xlabel, ylabel, log) :
  'x' and 'y' lengths differ
> par(mar=c(5,5,0,2))
> plot(dat$year, y7, xlab = "Time", type = "l")
错误于xy.coords(x, y, xlabel, ylabel, log) :
  'x' and 'y' lengths differ
> mtext("Big Title", outer = TRUE, cex = 1)
错误于mtext("Big Title", outer = TRUE, cex = 1) : 还没有调用plot.new
>

麻烦版主了！

呵呵，不知是否满意

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

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

太好了，但能否让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')

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

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

data=read.csv('C:\\Users\\acer\\Desktop\\x1.csv',head=T)
par(mfrow=c(3,2),mar=c(0,8,0,2), oma = c(3, 0, 3, 0),cex.axis=1,cex.lab=2,mex=0.5)
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)

调节mar和mex即可，但是具体调多少，还是你自己尝试

调节mar和mex即可，但是具体调多少，还是你自己尝试
data=read.csv('C:\\Users\\acer\\Desktop\\x1.csv',head=T)
par(mfrow=c(3,2),mar=c(0,8,0,2), oma = c(3, 0, 3, 0),cex.axis=1,cex.lab=2,mex=0.5)
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)

> my.melt <- melt(ds,id="year",measure=c("y1","y2","y3","y4","y5","y7"))
错误: 没有"melt"这个函数

这是为什么？

melt是ggplot2里面的函数。

reshape包里面也有melt函数
a <- array(1:24, c(2,3,4))
melt(a)

学习了

我又重新安装的ggplot2,仍然和上次一样，提示找不到my.melt。
我怀疑是出在了参数：id="year"上。因为这是一个时间序列，应该用ts()生成，但我却没有用，不会用，只是用excel硬生生地输入了year。
版主，是否是这个问题呢？

可以啊，不知道你的是怎么回事
library(ggplot2)
data=read.csv('C:\\Users\\acer\\Desktop\\x1.csv',head=T)
my.melt <- melt(data,id="year")

找到原因了：我用最新的R.2.14.2做，才出了问题。
刚才用2.14做，就没有问题了。
不过更不解了，最新的R2.14.2有毛病？？
多谢版主和各位朋友热情相助！

请问在R中如何改变线条的颜色。谢谢！