Standard operating procedures for ANOVA

I hear you, but I have to decide whether I want to incorporate any changes, as it implements user preferences as it stands now…

> options(digits = 3)
> emmeans(fit, pairwise ~ Trt)
NOTE: Results may be misleading due to involvement in interactions
$`emmeans`
 Trt emmean lower.HPD upper.HPD
 0     1.80      1.57      2.02
 1     1.45      1.23      1.68

Results are given on the log (not the response) scale. 
HPD interval probability: 0.95 

$contrasts
 contrast estimate lower.HPD upper.HPD
 0 - 1        0.35    0.0361     0.676

Results are given on the log (not the response) scale. 
HPD interval probability: 0.95 

On the other hand, I do limit precision for some frequentist summaries, e.g. t rations and P values; and I have been known to write letters to newspaper editors saying things like “86.234862 percent of people who report results to exaggerated precision are just trying to impress you.”

2 Likes