Hi all, I have thousands samples from different families. Now I want to find differential expressed gene using Limma. I have family correlation matrix from kinship, Is there any way I can include these correlation in lmFit? Thanks, Jack [[alternative HTML version deleted]]