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

Question

我有以下循环，使我可以进行重复测量方差分析，以分析1个平均评分（在此处称为Anger与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')
}

任何帮助将不胜感激。

谢谢。