r - 如何在循环中为 aov() 添加 Bonferroni 调整（在 R 中）

Question

Anger我有以下循环，可让我进行重复测量方差分析，以分析音乐曲目中 1 个平均评分（此处命名）与 8 个其他平均评分（列表中的 8 种不同情绪）之间是否存在显着差异。

由于它是重复测量方差分析，我知道必须进行p-value更正，因此我想进行Bonferroni调整。但是，我不确定如何为这个特定的代码编写它以及它应该放在循环中的什么位置。

代码：

trk <- fltr[fltr$Track == '001',]
lst <- c('Joy','Calm','Fear','Longing','Love','Power','Sad','Surprise')

for(k in 1:8){
    trgt <- trk[trk$Emotion == 'Anger',]  
    other <- trk[trk$Emotion == lst[k],]
    trgt <- rbind(trgt,other)  
    cat('Anger vs.', lst[k],'\n')  
    print(summary(aov(Rating ~ Emotion + Error(participant/Emotion),trgt))) 
    cat('\n')
}

任何帮助，将不胜感激。

谢谢你。