如何使用R中的plot()函数使基本图更美观

时间:2012-08-01 23:48:04

标签: r plot

我正在创建一个R脚本来生成由多行组成的图表。我已经制作了一个基本的图表,并设法平滑每条绘制的线条。但是,我想通过以下方式美化图表:

  1. 目前,每条线都被绘制为平滑线。我想使用其他一些线型绘制平滑线,例如点划线,破折号点等。

  2. 我想在图表中添加一个图例,并在图例中显示一条线(与主图表中代表的线颜色相同)。

  3. 指定主图表标题的字体大小,x和y

  4. 这是我到目前为止所做的:

    dat  <- read.csv(filename, sep=',')
    
    xvals <- dat$xvals
    y1  <- dat$y1
    y2  <- dat$y2
    y3  <- dat$y3
    y4  <- dat$y4
    y5  <- dat$y5
    y6  <- dat$y6
    
    lo1 <- loess(y1~xvals)
    lo2 <- loess(y2~xvals)
    lo3 <- loess(y3~xvals)
    lo4 <- loess(y4~xvals)
    lo5 <- loess(y5~xvals)
    lo6 <- loess(y6~xvals)
    
    plot(xvals,y1, xlab='X label', ylab='Y label', type='n')
    xl <- seq(min(xvals),max(xvals), (max(xvals) - min(xvals))/1000)
    lines(xl, predict(lo1,xl), col='gray', lwd=1)
    lines(xl, predict(lo2,xl), col='pink', lwd=1)
    lines(xl, predict(lo3,xl), col='red', lwd=1)
    lines(xl, predict(lo4,xl), col='cyan', lwd=1)
    lines(xl, predict(lo5,xl), col='black', lwd=1)
    lines(xl, predict(lo6,xl), col='green', lwd=1)
    
    legend("topright", "(x,y)", pch=1, lty=c(1,1), # gives the legend appropriate symbols (lines)
           , lwd=c(1,1),col=c("blue","red"), inset = .02)
    

    如何修改上述代码以实现上述要求1?

    顺便说一句,我知道ggplot。我想我会坚持使用情节() - 我觉得ggplot太神秘,令人费解和令人愤怒。

1 个答案:

答案 0 :(得分:5)

因为我已经开始研究答案了:

制作一些数据:

dat <- data.frame(xvals=1:100,
                  y=matrix(rnorm(600),ncol=6))

创建公式,应用loesspredict,将结果合并到矩阵中:

formulae <- lapply(names(dat)[-1],
                   function(yvar) as.formula(paste(yvar,"xvals",sep="~")))
loessfits <- lapply(formulae,loess,data=dat)

xl <- with(dat,seq(min(xvals),max(xvals), length.out=1000))
pvals <- do.call(cbind,lapply(loessfits,predict,newdata=data.frame(xvals=xl)))

设置参数(总是值得回去,花一些时间阅读?par):

## set box type (open), horizontal tick labels, axis label and title size
par(bty="l",las=1,cex.main=2,cex.lab=2)

设置colo(u)rs和线型:

cvec <- c("gray","pink","red","cyan","black","green")
lvec <- 1:6

创建情节:

matplot(xl,pvals,type="l",col=cvec,lty=lvec,
        xlab='X label', ylab='Y label')    
legend("topright", paste0("L",1:6), pch=1, lty=lvec,col=cvec, inset = .02)

enter image description here