Как добавить корректировку Бонферрони для aov () в цикле (в R) - PullRequest
Новая
       104

Как добавить корректировку Бонферрони для aov () в цикле (в R)

0 голосов
/ 19 октября 2018

У меня есть следующий цикл, который позволяет мне выполнять повторные измерения ANOVA, чтобы проанализировать, есть ли значительная разница между 1 средним рейтингом (здесь он назван Anger) против 8 других средних рейтингов (8 различных эмоций в списке) вмузыкальный трек.

Поскольку это ANOVA с повторными измерениями, я знаю, что должна быть сделана коррекция 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')
}

Любая помощь будет принята.

Спасибо.

...