title: plotmeans legends formatting
body: <p>I am plotting means of grouped data and I'm having trouble getting the legends to be right. The text is so large that one can only see the names of two groups, not all four. I have spent a long time trying to use <code>cex-like</code> commands to change the size, but it doesn't work. I have tried rotating them with <code>las=3</code>, but it doesn't work.</p>
  
  <p><img src="https://i.stack.imgur.com/kEQP7.png" alt="enter image description here"></p>
  
  <p>I cannot share the data, but the code is here:</p>
  
  <pre><code>plot.question = function(number){
    #which question to plot? get ID
    question = names(sorted.by.n)[number]
    #the formula
    form = paste0("DF.scored.g.scale ~ ",question)
    #fit it to data
    fit = lm(form, DF.merged.g)
    #get ANOVA results
    fit.anova = anova(fit)
    #get ANOVA p value
    p.value = round(fit.anova[[5]][2],4) #p value
    #plot it
    plotmeans(as.formula(form), DF.merged.g,
              ylab = "4 g-items sumscore",
              xlab = "Answer",
              main = paste0(questions.unique[question,"text"],"\nANOVA p=",p.value),
              cex.main = .8,
              cex.axis = .8,
              cex.lab = .8,
              cex.sub = .8,
              las=3,) #size of main title
  }
  </code></pre>
  
  <p>Preferably, I'd like to simply make the text smaller, so it can fit. Alternatively, I'd like to rotate it so it can fit (perhaps along with a margin change). If not what else?</p>
  
  <p>One can suppress the legends with <code>xaxt="n"</code>, but then one has to add them some other way. Can it really not be done within the <code>plotmeans()</code> function?</p>