箱线图以及小提琴图

生信喵

一、箱线图

1. ToothGrowth
   
2. boxplot(len ~ dose, data = ToothGrowth)
   
3. boxplot(len ~ supp, data = ToothGrowth)
   
4. boxplot(len ~ dose:supp, data = ToothGrowth)
   
5. # boxplot(len ~ dose+supp, data = ToothGrowth)
   
6. boxplot(len ~ supp+dose, data = ToothGrowth)
   
7. library(RColorBrewer)
   
8. brewer.pal(n = 3,name = 'Set2')
   
9. boxplot(len ~ supp+dose, data = ToothGrowth,
   
10.         col = c("#66C2A5","#FC8D62"))
    
11. boxplot(len ~ supp+dose, data = ToothGrowth,
    
12.         col = c("#66C2A5","#FC8D62"),notch=T)#画出凹口，可比较上下数据多少
    
13. boxplot(len ~ supp+dose, data = ToothGrowth,
    
14.         col = c("#66C2A5","#FC8D62"),horizontal = T)
    
15. t.test(ToothGrowth$len[ToothGrowth$supp == 'OJ'& ToothGrowth$dose == 0.5],
    
16.        ToothGrowth$len[ToothGrowth$supp == 'VC'& ToothGrowth$dose == 0.5])
    

复制代码

      
       带凹口的箱线图

二、小提琴图

1. library(vioplot)
   
2. vioplot(len ~ supp+dose, data = ToothGrowth,
   
3.         col = c("#66C2A5","#FC8D62"))
   
4. vioplot(len ~ supp+dose, data = ToothGrowth,
   
5.         col = c("#66C2A5","#FC8D62"),
   
6.         main = 'OJ vs VC',
   
7.         xlab = 'Supp and Dose',
   
8.         ylab = 'Tooth Length')
   
9. example(vioplot)
   
10. par(mfrow=c(1,1))#恢复案例后的绘图区域布局

复制代码