with(mdata, pairwise.t.test(score, mood, p.adjust.method="none", paired=T))
lin <- -3*y$control-y$semantic+y$lexical+3*y$phobic
phobdiff <- y$semantic+y$lexical-2*y$phobic
smdiff <- y$semantic-y$lexical

with(mdata, pairwise.t.test(score, mood, p.adjust.method="none", paired=T))
with(mdata, pairwise.t.test(score, mood, p.adjust.method="bonferroni", paired=T))

library(MHTdiscrete)
Sidak.p.adjust(c(0.132,0.004,0.106,0.0000092,0.0044,0.1212))
