My question is the same blood than this post C: DESeq2: can I correct for relatedness when using data from multiplex families?
In your solution, the model takes care only if the samples belong to the same family. I would like to know how we can add the level of relatedness such as MZ twin (100%) and DZ twin (50%). Currently, I have a variable zygosity with 2 levels (MZ and DZ).
Do you thus suggest to use the following model?
~ Covariate + Technical.Confounder + familyID + zygosity + condition
and how we do the same thing if we are also also other members of family.