I am trying to carry out heritability analysis with gene expression data.
After going through several tutorials and studies, I understood that a saturated model must be fitted for the data as a first step. The data must be tested for equal mean and variances by twin order and also by zygosity in subsequent steps. And all the three models must be compared to check if any differences exist between these models. The examples which I came across did not have significant differences between these models. So, my question is: What would be the appropriate way to proceed if there exists a significant difference between these models ?
Any response would be highly appreciated.