I performed an experiment where I raised different families coming from two different source populations, where each family was split up into a different treatments. After the experiment I measured several traits on each individual. Now I would like to have an overall and hence multivariate statistic that tests for an effect of either treatment or source as well as their interaction but correcting for family effects.
So, basically I would like to perform a MANOVA with a random effect, which would translate to something like this:
However, the standard R MANOVA function does not support random effects and I was unable to find another function that would do the trick. So I would like to ask if anyone has a suggestion that may help.