R 中的几个成对 t 检验(使用 Bonferroni 校正)

Several pairwise t-tests in R (with Bonferroni correction)

提问人:oden_g 提问时间:11/14/2023 最后编辑:JayPeerachaioden_g 更新时间:11/22/2023 访问量:23

问:

我正在进行一项实验,参与者必须在不同的维度上对几种不同的气味进行评分。下面是我的数据示例。换句话说,参与者在四个维度上对气味 A 进行评分,然后在相同的维度上评估气味 B,然后是气味 C。


A_pleasant <-c(4,5,4,6,7,5,5,4,4,1,5,6)
A_complex <-c(1,1,2,3,2,2,1,2,1,4,2,3)
A_familiar <-c(7,7,7,6,7,5,7,5,7,5,5,2)
A_relaxing <-c(4,5,4,6,3,4,2,4,5,2,5,5)
B_pleasant <-c(4,5,7,6,7,5,5,7,7,2,5,6)
B_complex <-c(3,3,3,2,3,5,1,3,2,3,3,3)
B_familiar <-c(6,6,6,4,3,4,3,4,7,6,5,6)
B_relaxing <-c(4,4,4,3,3,3,4,3,5,2,5,2)
C_pleasant <-c(1,1,1,2,2,1,2,1,4,1,2,1)
C_complex <-c(2,2,2,2,2,3,1,1,1,2,1,3)
C_familiar <-c(1,1,2,1,4,3,4,2,1,2,1,2)
C_relaxing <-c(2,2,2,2,2,1,1,1,1,1,2,2)

df <-data.frame(A_pleasant,A_complex,A_familiar,A_relaxing,B_pleasant,B_complex,B_familiar,B_relaxing,C_pleasant,C_complex,C_familiar,C_relaxing)
df

我有兴趣进行受试者内 t 检验,对多重比较进行 Bonferroni 校正,看看这些气味之间的差异在哪里。例如,我想知道哪些气味在愉悦感方面有显着差异。

如何在 R 上运行它?我的数据设置是否正确,还是需要以某种方式对变量进行分组?

谢谢!

我在 R 上使用过 t.test,但我不确定如何使用多个样本运行此测试,也不确定如何校正多个比较。

R 统计量 t 检验 Bonferroni

评论


答: 暂无答案