(完整)R语言实战前三章图形初阶运行代码

1、 #管理 r 工作空间的函数getwd() #显示当前的工作目录setwd(mydirectory) #修改当前的的工作目录为 mydirectoryls() #列出当前工作空间中的对象rm(objectlist) #移除（删除）一个或多个对象help(options) #显示或设置当前选项history(#)#显示最近使用过的#个命令（默认值为 25）savehistory(mylife) #保存命令历史到文件 mylife 中（默认值为.rhistory）loadhistory(mylife) #载入一个命令历史文件save.image(mylife) #保存工作空间到文件 mylife

2、中（默认值为.rdata）save(objectlist,file=mylife) #保存指定对象到一个文件中load(mylife) #读取一个工作空间到当前会话中（默认值为.rdata）q() #退出 r，将会询问你是否保存工作空间#保存图形输出的函数bmp(filename.bmp)jpeg(filename.jpeg)pdf(filename.pdf)png(filename.png)postscript(filename.ps)svg(filename.svg)win.metafile(filename.wmf)libpaths() #显示库所在的位置library() #显示库中有

3、哪些包#数据的输入#1、使用键盘输入mydata-data.frame(age=numeric(0),gender=character(0),weight=numeric(0)mydata-edit(mydata)#等价于fix(mydata)#直接在程序中嵌入数据集mydatatxt-age gender weight25 m 16630 f 11518 f 120 mydata-read.table(header=t,text=mydatatxt)mydata#2、从带分隔符的文本文件导入数据#mydataframe-read.table(file,options)#函数 read.tab

4、le()的选项#header#sep表示是文件是否在第一行包含了变量名的逻辑型变量分开数据的分隔符。默认是 sep=，表示一个或多个空格、制表符、换行或回车。#sep=,，读取逗号分隔行内数据的文件。sep=t，读取制表符分隔行内数据的文件#s#s#用于指定一个或多个行标记符header=false 时，可以用 s 去指定一个包含变量名的字符向量。如果 header=false 以及 s 选项被省略，变量名为 v1,v2用于表示缺失值的字符向量。#na.strings#na.strings=c(-9,?)把-9 和？值在读取数据时

5、转换成 nacolclasses=c(numeric,“numeric”,character,null,numeric)用于对有特殊字符的字符串划定界限的字符串#colclasses#quote#skip#stringsasfactors 逻辑变量，标记处字符向量是否需要转化成因子。默认值时 true，除非它被 colclasses 所覆盖。#处理大型文本文件时，设置成 stringsasfactors=false 可以提升处理速度。#3、导入 excel 数据library(xlsx)workbook-c:/myworkbook.xlsxmydataframe-read.xlsx(workb

6、ook,1)#4、导入 xml 数据library(xml)#5、从网页抓取数据#readlines() 将网页上的文字下载到一个 r 的字符向量中grep()gsup()# rcurl# xml#6、导入 spss 数据#foreign 包中的函数 read.spss() install.packages(hmisc)library(hmisc)mydataframe-spss.get(mydata.sav,use.value.labels=true#让函数将带有值标签的变量导入为 r 中水平对应相同的因子)#7、导入 sas 数据#foreign 包中的 read.ssd()#hmisc

7、包中的 sas.get()安装了 sas 情况下的好的选择#sas7bdat 包中的 read.sas7bdat()#8、导入 stata 数据library(foreign)mydataframe-read.dta(mydata.dta)#9、导入 netcdf 数据#ncdf 包#ncdf4 包#10、导入 hdf5 数据#rhdf5 包#11、访问数据库管理系统install.packages(rodbc)library(rodbc)#rodbc 中的函数#odbcconnect(dsn,uid=,pwd=) #建立一个到 odbc 数据库的连接#sqlfetch(channel,sql

8、table)#sqlquery(channel,quuery)#读取 odbc 数据库中的某个表到数据框中#向 odbc 数据库提交一个查询并返回结果#sqlsave(channel,mydf,tablename=sqltable,append=false) # 将数据框写入或更新到odbc 数据库的某个表中#sqldrop(channel,sqtable)#删除 odbc 数据库中的某个表#close(channel)#关闭连接#处理数据对象的实用函数#length()#dim()#str() #显示某个对象的结构#class() #显示某个对象的类型 #mode() #显示对象的模式#na

9、mes()#cbind() #按列合并对象rbind()#head() #列出对象的开始部分#ls() #显示当前的对象列表tail() #列出对象的最后部分#rm(object,object,.) #删除一个或多个对象乎所有对象rm(list = ls() #删除当前工作环境中的几#newobject-edit(object) #编辑对象并另存为 newobjectfix() #直接编辑对象#图形参数#pch=0 #绘图符号 025#cex=1#lty=1#绘图符号大小。1.5 放大为默认值的 1.5 倍#线条类型 16#lwd=1#线条宽度，默认为 1，=2 生成的线条宽度为默认宽度的 2

10、 倍#用于指定颜色的参数#col=c(red,blue) #lines pie 可以接受含有颜色值的向量并自动循环#col.axis #坐标轴刻度文字的颜色#col.lab坐标轴标签颜色#col.main#col.sub#fg #图形前景色#bg #图形的背景色#颜色向量函数：rainbow(),heat.colors(),terrain.colors(),topo.colors(),cm.colors().#rcolorbrewer 包library(rcolorbrewer)n-7mycolors-brewer.pal(n,set1)mycolorsbarplot(rep(1,n),col

11、=mycolors)#若要得到所有可选调色板的列表，输入：，或者输入：display.brewer.all()#多阶灰度色使用 gray()n-10mycolors-rainbow(n)mycolors pie(rep(1,n),labels=mycolors,col=mycolors)mygrays-gray(0:10/10)#生成 10 阶灰度色mygrayspie(rep(1,n),labels=mygrays,col=mygrays)#文本属性#指定文本大小的参数#cex#cex.axis 坐标轴刻度文字的缩放倍数cex.lab cex.main cex.

12、sub#指定字体族、字号、字样的参数#font 字体样式。1=常规，2=粗体，3=斜体，4=粗斜体，5=符号字体#font。axis font.lab font.main font.sub#ps 字体磅值。文本的最终大小为 ps*cex#family 绘制文本时使用的字体族。标准的取值为 serif（衬线），sans（无衬线），mono（等宽）#举例说明par(font.lab=3,cex.lab=1.5,font.main=4,cex.main=2)windowsfonts(a=windowsfont(arial black)b=windowsfont(bookkman old style)

13、c=windowsfont(comic sans ms)#图形尺寸和边界大小的参数#pin 以英寸表示的图形尺寸（宽和高）#mai#mar以数值向量标识的边界大小，顺序为“下，左，上，右”，单位为英寸以数值向量标识的边界大小，顺序为“下，左，上，右”，单位为英分par(pin=c(4,3),mai=c(1,.5,1,.2)#生成一幅 4 英寸宽，3 英寸高，上下边界为 1 英寸，左边界为 0.5，右边界为 0.2 英寸的图形dose-c(20,30,40,45,60)druga-c(16,20,27,40,60)drugb-c(15,18,25,31,40)#par()设定图形参数对所有图都有

14、效opar-par(no.readonly=true)par(pin=c(2,3)par(lwd=2,cex=1.5) par(cex.axis=.75,font.axis=3) #坐标轴刻度文本设置为斜体，缩小为默认大小的75%plot(dose,druga,type=b,pch=19,lty=2,col=red) #红色实心圆圈和虚线plot(dose,drugb,type=b,pch=23,lty=6,col=blue,bg=green) # 绿色填充的绿色菱形蓝色边框和蓝色虚线par(opar)#添加文本、自定义坐标轴和图例plot(dose,druga,type=b,pch=2,lt

15、y=2,col=red,lwd=2,main=clinicle trials for druga,sub=this is hypothetical data,xlab=dosage,ylab=drug response,xlim=c(0,60),ylim=c(0,70)#标题title(main=main title,sub=sub title,xlab = x-axis label,ylab = y-axis label)title(main=my title,col.main=red,sub=my subtitle,col.sub=blue,xlab=my x label,ylab=my

16、y label,col.lab=green,cex.lab=0.75)#坐标轴x-c(1:10)y-xz-10/xopar-par(no.readonly = true)#增加边界大小par(mar=c(5,4,4,8)+0.1)#绘制 x 对 y 的图形plot(x,y,type = b,pch=21,col=red,yaxt=n,lty=3,ann=false)#添加 x 对 1/x 的直线lines(x,z,type = b,pch=22,col=blue,lty=2)#绘制自己的坐标轴axis(2, #1=下，2=左，3=上，4=右at=x,labels=x, #字符型向量，表示刻度线

17、旁边的文字标签（如为null，则将直接使用 at中的值）col.axis=red,#数值型向量，表示需要绘制刻度线的位置las=2)#标签是否平行于（=0）或垂直于（=2）坐标轴axis(4,at=z,labels = round(z,digits = 2), col.axis=blue,las=2,cex.axis=0.7,tck=-.01 #刻度线的长度，以相对于绘图区域大小的分数表示（负值在图形外侧，正值在图形内侧#0 表示禁用刻度，1 表示绘制网格线）默认值为-0.01)#添加标题和文本mtext(y=1/x,side=4,line=3,cex.lab=1,las=2,col=blue

18、)title(an example of creative axes,xlab=x values,ylab=y=x)par(opar)#创建次要刻度线#需要使用 hmisc 包中的 minor.tick()函数library(hmisc)#在 x 轴的每两条主刻度线之间添加 1 条次要刻度线，#并在 y 轴的每两条刻度线之间添加 2 条次要刻度线，#次要刻度线时主要刻度线的一半minor.tick(nx=2,ny=3,tick.ratio=0.5)#在 y 轴的 1,5,7 位置添加水平实线abline(h=c(1,5,7)#在 x 为 1，3，5，7，9 的位置添加垂直的蓝色虚线abline

19、(v=seq(1,10,2),lty=2,col=blue)#图例dose-c(20,30,40,45,60)druga-c(16,20,27,40,60)drugb-c(15,18,25,31,40)#par()设定图形参数对所有图都有效opar-par(no.readonly=true)par(lwd=2,cex=1.5,font.lab=2)plot(dose,druga,type=b,pch=15,lty=1,col=red,ylim=c(0,60),main=drug a vs.drug b,xlab=drug dosage,ylab=drug response)lines(dose

20、,drugb,type=b,pch=17,lty=2,col=blue) abline(h=c(30),lwd=1.5,lty=2,col=gray)library(hmisc)minor.tick(nx=3,ny=3,tick.ratio=0.5)legend(topleft, #bottom,bottomleft,left,topleft,top,topright,right,center,inset=.05, #指定图例向图形内侧移动的大小（以绘图区域大小的分数表示）title=drug type,c(a,b),lty=c(1,2),pch=c(15,17),col=c(red,blue

21、)par(opar)#文本标注text(location,text to place,pos,.)mtext(text to place,side,line=n,.)#location 文本的位置参数，可为一对 x,y 坐标，#pos 文本相对于位置参数的方位：1234，下左上右。指定了 pos 就可以同时指定参数 offset=作为位置偏移量，以相对于单个字符宽度的比例表示#side 指定用来放置文本的边。可以指定参数 line=来内移或外移文本，随着值的增加，文本将外移。也可以使用 adj=0 将文本向左下对齐，=1 向右上对齐attach(mtcars)plot(wt,mpg,main=

22、mileage vs. car weight,xlab=weight,ylab=mileage,pch=18,col=blue)text(wt,mpg,s(mtcars),cex=0.6,pos=4, col=red)detach(mtcars)#展示不同字体族opar-par(no.readonly=true)par(cex=1.5)plot(1:7,1:7,type=n)text(3,3,example of default text)text(4,4,family=mono,example of mono-spaced text)text(5,5,family=serif

23、,example of serif text)par(opar) #图形的组合#1 2 行 2 列按行填充的图形排列attach(mtcars)opar-par(no.readonly=true)par(mfrow=c(2,2)plot(wt,mpg,main=scatterplot of et vs. mpg)plot(wt,disp,main=scatterplot of wt vs. disp)hist(wt,main=histogram of wt)boxplot(wt,main=boxplot of wt)par(opar)detach(mtcars)#2 三行一列图形排列attach(mtcars)opar-par(no.readonly=true)par(mfrow=c(3,1)hist(wt)#hist()为高级作图函数，包含一个默认的标题，使用 main=