I would like to perform post-hoc analysis on adjusted variances rather than on adjusted means (lsmeans or emmeans). Is it possible to calculate adjusted variances? and to compare them?
I have two different linear models:
model1 = lm((Phenotype)~ Batch + Tray + Genotype, data=data)
model2 = glm((Phenotype)~ Batch + Tray + Genotype, data=data, family = binomial())
I expect to obtain adjusted variances from linear models.
Thanks for your help,
LM