事后测试 你可以在内部受试者因子(这里是时间)的各个水平之间进行多次配对的t检验。p值使用Bonferroni多重检验校正方法进行调整。 # pairwise comparisons pwc <- selfesteem %>% pairwise_t_test( score ~ time, paired = TRUE, p.adjust.method = "bonferroni" ) pwc All the pairwise differences are stat...
1. t检验(pairwise.t.test):当数据满足正态性和方差齐性假设时,可以使用t检验来进行两两比较。该函数会对每对组进行t检验,计算出每对之间的差异显著性水平和置信区间。 ```R pairwise.t.test(data$group, data$value, p.adjust.method = "bonferroni") ``` 2. 方差分析(ANOVA):如果数据不满足t检验的...
函数pairwise.t.test()能够计算所有的两组比较。他能够针对多重检验做调整: > pairwise.t.test(folate,ventilation,p.adj="bonferroni") Pairwise comparisons using t tests with pooled SD data: folate and ventilation N2O+O2,24h N2O+O2,op N2O+O2,op 0.042 - O2,24h 0.464 1.000 P value adjustment ...
macs2 callpeak --nomodel --shift 0 --extsize 150 -g hs --keep-dup all -f BAM -t 2-read-align/HKE293-D210N-V5ChIP-Rep*.flt.bam -c 2-read-align/HKE293-D210N-Input-Rep*.flt.bam --outdir 5-peak-calling/narrow -n HKE293-D210 2> log/HKE293-D210.macs2.narrow.log & macs2...
pairwise.t.test(my_data$weight, my_data$group, p.adjust.method = "BH", pool.sd = FALSE) 2.3 检查正态性假设 2.3.1 残差的正态图。 在下面的图中,相对于正态分布的分位数绘制了残差的QQ图,并且还绘制了45度参考线。 残差的正态概率图用于检查残差呈正态分布的假设。如果它满足正态性,应该要大...
首先，Bonferroni法通过对检验水平进行调整来应对多组比较，尤其在比较次数小于10时效果较好。在R中，推荐使用pairwise.t.test函数进行多重比较，它能返回调整后的P值，避免混淆。以airquality数据为例，对比不同月份的航班数，结果显示5月与其他月份的差异显著。Dunnett检验则适用于实验设计中的特定比较，如...
在R中,使用pairwise.t.test()函数可以帮助我们完成多重T检验,该函数的原型如下: pairwise.t.test(x, g, p.adjust.method = p.adjust.methods, pool.sd = !paired, paired = FALSE, alternative = c("two.sided", "less", "greater"), ...) ...
pairwise.t.test(x,g,p.adjust.method="holm")——多重t检验,p.adjust.method是P值的调整方法,其方法由p.adjust()给出,默认值按Holm方法(”holm“)调整,若为”none“,表示P值不做任何调整。双因素交互作用时g=A:B shapiro.test(x)——数据的正态W检验 ...
bartlett.test(formala, data, subset,na.action…) formula是形如lhs一rhs的方差分析公式;data指明数据集:subset是可选项,可以用来指定观测值的一个子集用于分析:na.action表示遇到缺失值时应当采取的行为。 续上例: 代码语言:javascript 复制 >x=c(x1,x2,x3) ...
oneway.test(y~x,data=data) #p值为3.155e-06,说明在0.05水平下极显著, 说明四组有显著差异。#多重比较——进行两两t检验并不使用合并的标准差估计, 使用Holm方法进行p值调整以控制总错误率: with(data,pairwise.t.test(y,x,pool.sd=FALSE,p.adjust.method="holm")) ...