R绘图系统2重点.docx

2.1par(mfrow=c(2, 2), cex=0.6, mar=c(4, 4, 1, 1))y - rnorm(20)plot(y, type=p)plot(y, type=l)plot(y, type=b)plot(y, type=h)2.2par(mfrow=c(2, 2), cex=0.6, mar=c(4, 4, 4, 2), mex=0.8)plot(lm.SR - lm(sr ~ pop15 + pop75 + dpi + ddpi, data = LifeCycleSavings),id.n=1, cex.caption=0.8, which=1:4,panel=function(...) { panel.smooth(..., col.smooth=gray) })2.3library(cluster)subset - sample(1:150, 20)cS - as.character(Sp - iris$Species[subset])cS[Sp == setosa] - ScS[Sp == versicolor] - VcS[Sp == virginica] - gai - agnes(iris[subset, 1:4])par(mfrow=c(2, 1), cex=0.5, pty=s)plot(ai, which=1, col=c(gray90, gray), labels = cS)plot(ai, which=2, labels = cS)2.4plottitle - function(plotfun, funarg, outer=FALSE, cex=.7, line=1) {ncp - nchar(plotfun)nca - nchar(funarg)mtext(paste(plotfun, (, paste(rep( , nca), collapse=), ), sep=),family=mono, cex=cex, line=line, font=2, outer=outer)mtext(paste(paste(rep( , ncp + 1), collapse=),funarg, , sep=),family=mono, col=gray60, cex=cex, line=line, font=2, outer=outer)}plot2title - function(plotfun, funarg, extrafn, extraarg, outer=FALSE, cex=.7, line=.5) {ncp - nchar(plotfun)nca - nchar(funarg)ncep - nchar(extrafn)ncea - nchar(extraarg)mtext(paste(plotfun, (, paste(rep( , nca), collapse=), )\n, extrafn, (,paste(rep( , ncea), collapse=), ), sep=),family=mono, cex=cex, line=line, font=2, outer=outer)mtext(paste(paste(rep( , ncp + 1), collapse=),funarg, \n, paste(rep( , ncep + 1), collapse=),extraarg, , sep=),family=mono, col=gray60, cex=cex, line=line, font=2, outer=outer)}dohplot - function(plotfn, ..., funarg, extrafn=NULL, extraarg=NULL, namefudge=FALSE, main=, xlab=, ylab=, axes=FALSE, box=TRUE) {if (is.null(xlab) || is.null(ylab)) {do.call(plotfn,list(..., main=)) } else if (is.null(axes)) {do.call(plotfn,list(..., main=, xlab=, ylab=)) } else {do.call(plotfn,list(..., main=, axes=FALSE, xlab=, ylab=)) }if (is.null(extrafn)) {plottitle(plotfn, funarg)