R语言绘图系统
基础绘图包 ⚠️
Lattice （语法复杂）
ggplot2家族 ⚠️
其他，比如：sjplot; plots; plotrix...

1. 基础绘图包的绘图函数

高级绘图函数

低级绘图函数见R语言可视化及作图2--低级绘图函数

2. par函数（进行绘图参数的设置）

bty=c('o','l','7','c',']','u') #将bty所有格式设置为 一个向量
par(mfrow=c(3,2)) #切分画布
for (i in 1:6){
  par(bty=bty[I])
  plot(1:5,main = paste('the bty is',bty[i],sep = ':'))
}

bty参数

plot(rep(1:5,times=5),rep(1:5,each=5),pch=1:25,cex=1.5,bty='l',xlim=c(1,5.4))
text(rep(1:5,time=5)+0.15,rep(1:5,each=5),labels=1:25,col='red')

pch参数。R语言中默认有25种点的形式

3. 散点图 plot函数

type=c('p','b','l','s','o','n')
par(mfrow=c(2,3))
for (i in 1:6){plot(1:10,type=type[i],main=paste('The plot type is: ',type[i]))} 

4. box plot盒形图

par(mfrow=c(1,3)) #mfrow设置横轴可以画三张图
set.seed(100)
data=rnorm(40)
boxplot(data,range=0.5,xlab='range is 0.5') #range determines how far the plot whiskers extend out from the box. 
boxplot(data,range=1,xlab='range is 1')
boxplot(data,range=1.5,xlab='range is 1.5')

boxplot(len~dose,data = ToothGrowth,#根据dose对len画图
        boxwex=0.25,at=1:3-0.2,
        #boxwex设置盒子宽度，at参数定义了盒形图中盒子横坐标的位置，传入一个向量，也就是1-0.2=0.8，2-0.2=1.8，3-0.2=2.8)
        subset = supp=='VC',col='orange',
        names=c('dose=0.5','dose=1.0','dose=2.0'),
        main="Guinea Pigs' Tooth Growth",
        xlab='Vitamin C dose mg',
        ylab='Tooth Length',
        xlim=c(0.5,3.5),ylim=c(0,35),yaxe='i') #xlim和ylim设置xy轴的取值范围，yaxe设置y轴样式

boxplot(len~dose,data = ToothGrowth, add = TRUE, #在原有图上增加新图
        boxwex=0.25,at=1:3+0.2,
        names=c('dose=0.5','dose=1.0','dose=2.0'),
        subset = supp=='OJ',col='yellow')

#准备数据 
names=c(rep('Maestro',20),rep('Presto',20),
        rep('Nerak',20),rep('Eskimo',20),rep('Nairobi',20),rep('Artiko',20))
value=c(sample(3:10,20,replace = T),sample(2:5,20,replace = T),sample(6:10,20,replace = T)
        ,sample(6:10,20,replace = T),sample(1:7,20,replace = T),sample(3:10,20,replace = T))
data=data.frame(names,value)

boxplot(data$value~data$names,
        col=ifelse(levels(data$names)=='Nairobi',rgb(0.1,0.1,0.7,0.5),
                   ifelse(levels(data$names)=='Eskimo',rgb(0.8,0.1,0.3,0.6),'grey90')),
        ylab='disease',xlab='-variety-') #没有颜色不知道为什么
legend('bottomleft',legend = c('Positive control','Negative control'),
       col=c(rgb(0.1,0.1,0.7,0.5),rgb(0.8,0.1,0.3,0.6)),bty='n',pch=20,pt.cex = 3,)

order_name <- with(data,reorder(names,value,median,na.rm=TRUE))
boxplot(data$value~order_name,col='darkorchid',ylab='disease',xlab='-variety-') #按均值排列

5. bar plot 条形图（离散型数值）

data <- sample(c(50:80),5)
barplot(data,col=heat.colors(5))

my_matrix <- matrix(data=sample(10:40,9),nrow = 3,dimnames = list(c('A','B','C'),paste('dose',1:3)))
my_matrix
barplot(t(my_matrix),beside = TRUE,col = rainbow(3))

par(mfrow=c(2,2))
average_gdp <- c(6500,8000,13000,9200)
country <- c('China','Korea','Japan','Singapore')
barplot(average_gdp,names.arg=country)
barplot(average_gdp,names.arg=country,horiz=TRUE)
barplot(average_gdp,names.arg=country,horiz=FALSE,width=0.5)
barplot(average_gdp,names.arg=country,horiz=FALSE,space=1.2)

#加error bar(基础绘图包中不能直接添加，需要手动添加，ggplot2可以直接添加)
data(iris)
data_mean <- apply(iris[,1:3],2,mean)#生成均值
data_sd <- apply(iris[,1:3],2,sd)#生成标准差
barcenters <- barplot(data_mean,name.arg=names(data_mean),ylim=c(0,8),main='ugly')
barcenters
#将barplot()存入一个对象中，这个对象保存的是条柱中心位置的横坐标。这一步的目的是找到三个bar的横坐标最中心，用于确定error bar的位置
segments(barcenters,data_mean-data_sd,barcenters,data_mean+data_sd,lty=1.2) #添加中间的线段
arrows(barcenters,data_mean-data_sd,barcenters,data_mean+data_sd,code = 3,angle = 90) #添加上下的横杠

6. hist 直方图（连续型数值）

op <- par(mfrow=c(2,3))
data <- rnorm(100,10,5)
hist(data,col = 'light green')
hist(data,col='sky blue',breaks = 10) #分成10组（break传入数值表示分成多少组）
hist(data,col = 'orange',breaks = seq(-5,25,1)) #自定义组距（break也可通过seq()设置组距）
hist(data,col = 'pink',breaks = seq(-5,25,1),density = TRUE) #加shading line
hist(data,col = 'pink',breaks = seq(-5,25,1),freq = FALSE) #概率密度图
lines(density(data),col='blue',lty=1,lwd=2)
#添加概率密度曲线。注意：添加概率密度曲线前，一定要把freq设置成FALSE。
par(op) #释放参数


#同向叠加直方图
titanic <- carData::TitanicSurvival
#数据读入
hist(titanic$age,main = 'Passenger Survival by Age',xlab = 'Age',col = 'sky blue',breaks = seq(0,80,2))
hist(titanic$age[which(titanic$survived=='no')],col = 'orange',add=TRUE,breaks = seq(0,80,2)) #在原图上添加一个死亡人数的直方图


#背靠背直方图
df=data.frame(x=rnorm(100),x2=rnorm(100,mean = 2)) #创建数据
h1=hist(df$x,plot=FALSE)
h2=hist(df$x2,plot = FALSE)
#绘制两个直方图，储存在h1和h2这两个对象中
h2$counts=-h2$counts
#将h2的值反转过来
hmax=max(h1$counts)
hmin=min(h2$counts)
#这两步是为了找到y轴的取值范围
X=c(h1$breaks,h2$breaks)
xmax=max(X)
xmin=min(X) #这三步是为了找到x轴的取值范围
plot(h1,ylim = c(hmin,hmax),col = 'green',xlim = c(xmin,xmax))
lines(h2,col='blue')#注意，用敌机会图函数lines将h2添加上去