R语言笔记PART1

R语言笔记画图函数列表数列函数seq()seq()是一个输入数列的函数：seq(from,to,step)e.g.seq(1,10,2)输出：1,3,5,7,9格式函数par()e.g.opar-par(no.readonly=TRUE)---复制当前画图所有参数par(pin=c(3,7),mai=c(1,2,3,4))--改变图片参数：尺寸（宽=3高=4）及位置（下=1左=2上=3右=4，单位英寸（英分为mar））par(lty=3,pch=17)--改变图片参数：线型=3点型=17par(font.lab=2)—改变坐标轴标签的字体plot(x,y)–用更改后的参数画图par(opar)–恢复原始参数（即第一句复制的参数）画图函数plot()参数：type:图形类型（b,s,l,m等）col:主图颜色（线型、点型）lty:线型pch:点型lwd:线型宽度cex:缩放倍数（线型、点型）main:主标题sub:副标题xlab(ylab):x(y)轴标签xlim(ylim):x(y)轴最大刻度范围cex.:某项指标的缩放倍数cex.axis:坐标轴宽度缩放倍数cex.xlab:x轴标签缩放倍数cex.main:主标题缩放倍数cex.sub:副标题缩放倍数Font.:某项指标的字体font.axis:坐标轴刻度的字体font.xlab:x轴标签字体font.main:主标题字体font.sub:副标题字体e.gplot(x=c(1:10),y=2*x+3,type=”b”,col=”red”,lty=2,pch=2,lwd=2,main=”Testmaintitle”,sub=”testsubtitle”,xlab=”level”,ylab=”score”,xlim=c(1:12),ylim=c(1:30),col.main=”pink”,cex.main=2,cex.sub=1.5)线型图形lines()lines()语句可以为一副现有的图形添加新的图形元素。e.g.x-c(1:10)y-xz-10/xplot(x,y,type=”b”,pch=22,lty=2,col=”red”,ann=FALSE)lines(x,z,type=”b”,pch=21,col=”blue”,lty=2)设置坐标轴格式axis()e.g(接上)axis(2,at=x,labels=x,col.axis=”red”,las=2”)axis(4,at=z,labels=round(z,digits=2),col.axis=”blue”,lax=2,cex.axis=0.7,tck=-.01)添加绘图区文本函数text()/mtext()text()可向绘图区域内部添加文本。text(location,”texttoplace”,pos=n,…)mtext用于在图形的边界添加文本。mtext(“testtoplace”,side=n,line=n,…)e.g.(接上)mtext(“y=1/x”,side=4,line=3,cex.lab=1,las=2,col=”blue”)e.g.1attach(mtcars)plot(wt,mpg,main=Mlieagevs.CarWeight,xlab=Weight,ylab=Mileage,pch=8,col=blue,cex=0.5)text(wt,mpg,row.names(mtcars),cex=0.3,pos=4,col=red)detach(mtcars)e.g.2opar-par(no.readonly=TRUE)par(cex=1.5)plot(1:7,1:7,type=n)text(4,4,family=mono,Exampleofmono-spacedtext)text(5,5,family=serif,Exampleofmono-spacedtext)text(5,5,family=serif,Exampleofseriftext)text(5,7,family=serif,Exampleofseriftext)text(3,3,Exampleofdefaulttext)par(opar)次要刻度线（Hmisc包）加载包：Hmisc函数：minor.tick()Install.packages(Hmisc”)library(Hmisc)minor.tick(nx=n,ny=n,tick.ratio=n)添加次要刻度线。其中nx、ny笔试次要刻度线现对于主刻度线的大小比例。当前的主刻度线长度可以使用par（”tck”）获取。e.g.minor.tick(nx=2,ny=3,tick.ratio=0.5)表示在x轴的每两条主刻度线之间添加一条次要刻度线，并在y轴的每两条主刻度线之间添加2条次要刻度线。次要刻度线的长度将是主刻度线的一半。参考线abline()用来为图形添加参考线。abline(h=yvalues,v=xvalues)e.gabline(h=c(1,5,7))在y为1,5,7位置添加了水平实线（默认线型为实线lty=1）。abline(v=seq(1,10,2),lty=2,col=”blue”)则在x为1,3,5,7,9的位置添加了垂直的蓝色虚线。图例1-1opar-par(no.readonly=TRUE)plot(dose,drugA,lty=1,pch=15,type=b,ylim=c(0,60),xlim=c(20,60),main=DrugAvs.DrugB,xlab=DrugDosage,ylab=DrugResponse,col=red)lines(dose,drugB,lty=2,pch=17,type=b,col=blue)abline(h=c(30),lwd=1.5,lty=2,col=grey)install.packagse(“Hmisc”)library(Hmisc)minor.tick(nx=3,ny=3,tick.ratio=0.5)legend(topleft,inset=0.05,title=DrugType,c(A,B),lty=c(1,2),pch=c(15,17),col=c(red,blue))图形组合图形组合参数mfrow()mfrow(nrows,ncols)用来创建按行填充的、任意行列数的图形矩阵高级绘图函数hist()得到柱状图。hist()函数包含了一个默认标题（见图例2-2），可以用main=””来禁用它，或者用ann=FALSE来禁用所有的标题和标签。图形组合函数layout()layout(mat),其中mat为一个矩阵，指定了所要组合的多个图形的所在位置。参见图例2-3。mat元素的数量决定了一个outputdevice被等分成几份，为了方便我把一份叫做一个格子。这样mat内的每个元素根据他们的行列序号对应一个格子。而元素本身的值代表它属于第几个figure。举例来看。layout(matrix(c(1,2,3,0,2,3,0,0,3),nr=3))matrix有9个元素，具有这样的形式：[,1][,2][,3][1,]100[2,]220[3,]333把这个矩阵传入layout函数，我们就能得到这样的outputdevice如此，figure1占据了左上角的一个格子，第二行的前两个格子属于figure2，figure3占满最下一行的三个格子。为了醒目，figure1，2，3分别标记了黄绿红颜色。在输出figure时，会按照先后顺序，将figure绘制在与其顺序相同的区域内。在我的这个例子内，就是按照黄色区域，绿色区域，红色区域的顺序。当然你可以通过更改matrix，使得各个figure按照你需要出现在不同区域，不一定按照从上到下或从左到右的传统顺序。范围参数fig()fig=c(x1,x2,y1,y2)。其中x1,x2表示图形横向占据范围fromx1tox2,;y1,y2表示图形纵向占据范围fromy1toy2。图例2-1attach(mtcars)opar-par(no.readonly=TRUE)par(mfrow=c(2,2))---建立一个数组元素为图形的2*2矩阵plot(wt,mpg,main=Scatterploofwtvs.mpg)plot(wt,disp,main=Scatterploofwtvs.disp)hist(wt,main=Histogramofwt)boxplot(wt,main=Boxplotogwt)par(opar)detach(mtcars)图例2-2opar-par(no.readonly=TRUE)hist(wt)hist(mpg)hist(disp)par(opar)detach(mtcars)图例2-3attach(mtcars)layout(matrix(c(1,1,2,3),2,2,byrow=TRUE),widths=c(3,1),heights=c(1,2))--matrix(c(1,1,2,3),2,2,byrow=TRUE)表示2*2矩阵，矩阵元素为（1,1,2,3）—widths=c(3,1)第一二行的宽度分别为(3,1):表示宽度比例为3:1—hights=c(1,2)第一二行的高度分别为(3,1):表示一个比例为1:2hist(mpg)hist(wt)hist(disp)detach(mtcars)图例2-4opar-par(no.readonly=TRUE)par(fig=c(0,0.8,0,0.8))plot(mtcars$wt,mtcars$mpg,xlab=MilesperGallon,ylab=CarWeight)par(fig=c(0,0.8,0.55,1),new=TRUE)boxplot(mtcar$wt,horizontal=TRUE,axes=FALSE)错误于boxplot(mtcar$wt,horizontal=TRUE,axes=FALSE):找不到对象'mtcar'boxplot(mtcars$wt,horizontal=TRUE,axes=FALSE)par(fig=c(0.65,1,0,0.8),new=TRUE)boxplot(mtcars$mpg,axes=FALSE)mtext(EnhancedScatterlpot,side=3,outer=True,line=-3)错误于mtext(EnhancedScatterlpot,side=3,outer=True,line=-3):找不到对象'True'mtext(EnhancedScatterlpot,side=3,outer=TRUE,line=-3)par(opar)基本数据管理函数within()wthin()与with作用类似，不同的是within允许修改数据框。举例参见下面的“变量重编码”变量重编码更新现有值的过程。格式为：变量名[conditinon]-变量新值仅在condition为TRUE时执行赋值。e.g.创建数据集leadershipmanager-c(1,2,3,4,5)date-c(10/24/08,10/28/08,10/1/08,10/12/08,5/1/09)country-c(US,US,UK,UK,UK)gender-c(M,F,F,M,F)age-c(32,45,25,39,99)q1-c(5,3,3,3,2)q2-c(4,5,5,3,2)q3-c(5,2,5,4,1)q4-c(5,5,25,NA,2)q5-c(5,5,2,NA,1)leadership-data.frame(manager,date,country,gender,age,q1,q2,q3,q4,stringAsFactors=FALSE)leadershipmanagerda