```r
ms <-ggs(output) 
mt<-filter(ms,grepl("Difbeta",Parameter))
ggs_caterpillar(mt) +geom_vline(xintercept = 0,col="red")
```

<img src="004a_anova_files/figure-html/unnamed-chunk-35-1.png" width="672" />