如果两个样本具有方差齐性,那么做独立样本t检验时,直接套用t检验的公式,计算t值,,查表的自由度为n1+n2-2,然后用函数pt( t value, n1+n2-2)给出p值,小于0.05即为显著。
如果方差不齐,需要计算校正后的自由度,
最终出图如下
这里自动做统计检验的函数是 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 这个包
小明的数据分析笔记本
形式如下:
输出文件如下