General linear modeling with multiple responses (MANCOVA). An overall p-value for each model term is calculated by the 50-50 MANOVA method by Langsrud (2002) , which handles collinear responses. Rotation testing, described by Langsrud (2005) , is used to compute adjusted single response p-values according to familywise error rates and false discovery rates (FDR). The approach to FDR is described in the appendix of Moen et al. (2005) . Unbalanced designs are handled by Type II sums of squares as argued in Langsrud (2003) . Furthermore, the Type II philosophy is extended to continuous design variables as described in Langsrud et al. (2007) . This means that the method is invariant to scale changes and that common pitfalls are avoided.