R语言--图形初阶
第一个例子
dev.off()
pdf('my_firstgraph')
attach(mtcars)
plot(wt,mpg) #plot后面第一个参数是横坐标
abline(lm(mpg~wt))
title('Regression of MPG on Weight')
detach(mtacrs)
image-20210417081643587
第二个例子
rm(list=ls())
dose <- c(20, 30, 40, 45, 60)
drugA <- c(16, 20, 27, 40, 60)
drugB <- c(15, 18, 25, 31, 40)
par(mfrow=c(1,2))
plot(dose,drugA)
title("relations between dose and drug")
plot(dose,drugA,type='b')
lines(dose,drugB,type='b')#再添加一条线
title("relations between dose and drug")
图形参数
method
#方法一
rm(list=ls())
dose <- c(20, 30, 40, 45, 60)
drugA <- c(16, 20, 27, 40, 60)
drugB <- c(15, 18, 25, 31, 40)
opar <- par(no.readonly = T)
par(lty=2,pch=17) #lty改变线类型 pch改变符号类型
plot(dose,drugA,type='b')
par(opar)
#方法二
plot(dose,drugA,type='b',lyt=3,pch=20)
符号参数
指定符号和线条类型的参数 P48
| 参数 | 描述 |
|---|---|
| pch | 改变符号类型 |
| cex | 指定符号的大小 |
| lty | 指定线条的类型 1~6 |
| lwd | 指定线条的宽度 #默认为1 |
par(mfrow=c(1,3))
plot(dose,drugA,type='b',lyt=3,pch=20,cex=1.5,lwd=4)
plot(dose,drugA,type='b',lyt=3,pch=20,cex=0.5,lwd=2)
plot(dose,drugA,type='b',lyt=3,pch=20,cex=0.85,lwd=2.5)
image-20210417095525432
颜色
col 默认的绘图颜色（就是图形颜色，比如上面的符号与线条）
col.axis 坐标轴刻度文字的颜色（一般就是那些数字）
col.lab 坐标轴标题的颜色
col.main 标题颜色
col.sub 副标题颜色
fg 图形前景色（轴和框的颜色）
bg 图形背景色
#可以直接用基本颜色英文表示，也支持十六进制，RGB等等
col="red"
colors() #可返回所有可用颜色的名称，下面的网站可以查看具体所对应的颜色
# http://www.docin.com/p-1491680339.html
rainbow(3) #可以返回指定数目的“彩虹型”颜色
#[1] "#FF0000FF" "#00FF00FF" "#0000FFFF"
#R包RColorBrewer可以创建出吸引人的颜色配对
library(RColorBrewer)
mycolors <- brewer.pal(5, "Set1") #从Set1中抽取5中颜色
barplot(rep(1,10), col=mycolors) #循环重复
image-20210417100941239
文本参数
文本大小参数
cex.axis 为数值，表倍数
cex.lab
cex.main
cex.sub
#### 字体参数
font.axis 1=常规，2=粗体，3=斜体，4=粗斜体，5=符号字体
font.lab
font.main
font.sub
尺寸
pin 以英寸表示的图形尺寸，宽和高
mai 表图的四周的宽大小，顺序为下，左，上，右；单位为英寸
mar 同上，单位为英分。1英分=1/12英寸=0.21厘米
终合实例
opar <- par(no.readonly = T)
par(pin=c(2,3))
par(lwd=2,cex=1.5)
par(cex.axis=0.75,font.axis=2)
par(mfrow=c(1,2))
plot(dose,drugA,type='b',pch=0,lty=2,col='red',bg='blue',fg='green')
plot(dose,drugA,type='b',pch=23,lty=6,col='red',bg='blue',fg='yellow')
par(opar)
image-20210417103517579
添加文本、自定义坐标轴和图例
par(pin=c(3,3))
plot(dose,drugA,type='b',
col='red',lwd=2,lty=2,pch=2,
main='Clinical Trials for Drug A',col.main='red',
sub = 'This is hypothetical data',col.sub='green',
xlab='Dosage',ylab='Drug Reponse',col.lab='blue',
xlim = c(0,60),ylim=c(0,70))
image-20210417105258546
#添加标题
#title()函数
title(main="主标题",sub="副标题",xlab="x轴标题",ylab="y轴标题")
#同时可以参考上述参数进行细节修改
#添加次要刻度线
library(Hmisc)
minor.tick(nx=2,ny=3,tick.ratio=0.5)
#意为在x轴的的每两条主刻度线间添加一条次要刻度线（分成2段）
#y添加两条，分成3段。大小为主刻度的0.5倍。
添加参考线
abline(h=c(1,5,7)) #添加水平参考线
abline(v=c(2,5),lty=2,col="blue") #添加蓝色的垂直虚线
##添加图例
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,drugB,type='b',
pch=17,lty=2,col='blue')
abline(h=c(30),lwd=1.5,lty=2,col='gray')
legend("topleft",inset=.02,
title = "Drug Type",
c("A","B"),cex=0.7,bty="n", lty = c(2,6),pch = c(19,23),col=c("red","blue"),text.width=0.5
) #去掉边框bty="n" #调整图例大小cex=0.7
par(opar)
image
image-20210417114339723
文本标注
通过函数text（）和mtext（）将文本添加到图形上
text（）向绘图区域内部添加文本
mtext（）向图形的四个边界之一添加文本
函数text（）和mtext（）的选项
| 选项 | 描述 |
|---|---|
| location | 文本的位置参数 |
| pos | 文本相对于未知参数的方位 |
| side | 指定用来放置文本的边 1=下，2=左，3=上，4=右 |
attach(mtcars)
colnames(mtcars)
plot(wt,mpg,
main='Mileage vs. Car Weight',
xlab = 'Weight',ylab = 'Mileage',
pch=16,col='blue')
text(wt,mpg,
row.names(mtcars),
cex = 0.6,pos=4,col ='red')
detach(mtcars)
image-20210417203702726
多幅图显示在一起/拼图
#mfrow=c() 函数，默认 mfrow=c(1,1)
opar <- par(no.readonly=TRUE)
par(mfrow=c(2,3)) #按行填充，行数为2，列数为3的图形矩阵
par(mfcol=c(2,3)) #按列填充
#layout() 函数，可以设置“品”字形的图形矩阵
layout(matrix(c(1,1,2,3),2,2,byrow=TRUE),widths=c(3,1),heights=c(1,2))
par(opar)
#第一幅图被置于第一行，另两幅图被放在第二行。
#第一行图形高度是第二行的1/2，右下图宽度是左下宽度的1/3
图形布局的精细调整
opar <- par(no.readonly=TRUE)
par(fig=c(0,0.8,0,0.8))
plot(mtcars$wt,mtcars$mpg,col='blue',
xlab="Miles Per Gallon",
ylab="Car Weight",col.lab='red')
par(fig=c(0,0.8,0.55,1),new=TRUE)
# new=TRUE 很关键的参数，表示在原有图上加图
hist(mtcars$wt,main = "",xlab = "",col='green',col.lab='yellow',labels = T,ylim=c(0, 10))
par(fig=c(0.55,1,0,0.8),new=TRUE)
boxplot(mtcars$mpg,axes=F,col = 'orange')
mtext("Enhanced Scatterplot",side=3,outer = TRUE,line=-3,col='red')
par(opar)
dev.off()
image-20210417205748430
网友评论