R语言箱式图的制作代码

Python016

R语言箱式图的制作代码,第1张

#使用并列箱线图进行跨组比较，如图6-12#公式y~A*B则将为类别型变量A和B所有水平的两两组合生成数值型变量y的箱线图， 例子的mpg~cyl表示为四缸、六缸、八缸发动机对每加仑汽油行驶的英里数的影响。boxplot(mpg~cyl,data=mtcars,main="Car Mileage Data",xlab="Number of Cylinders",ylab="Miles Per Gallon")

#添加notch=TRUE，可以得到含凹槽的箱线图， 若两个箱的凹槽互不重叠， 则表明它们的中位数有显著差异； 而

varwidth=TRUE则使箱线图的宽度与它们各自的样本大小成正比。>boxplot(mpg~cyl,data=mtcars,notch=TRUE,varwidth=TRUE,col="red",main="Car Mileage Data",xlab="Number of Cylinders",ylab="Miles Per Gallon")

输出结果如下

最终出图如下

这里自动做统计检验的函数是 stat_compare_means()

读入数据

作图

这个函数来自于ggpubr这个包，只需要指定根据那一列来分组就可以了

默认的是Wilcoxon Rank Sum and Signed Rank Tests，如果要用t检验指定method参数

如果想把P值改成星号，直接加label=“p.signif”参数

这里如果不显著会在图上显示ns，如果不想要ns，可以加 hide.ns = TRUE 参数

星号的位置可以手动指定，用 label.y = c(26,31) 参数

使用到的是 ggsignif 这个包

小明的数据分析笔记本