I want to generate a 1600 samples (800 cases and 800 controls) and only 2 snps data for interaction analysis.
The condition that I know is MAF=0.1 and genetic heritability h^2 = 0.03, prevalence p(D)=0.1. The epistasis model is a multiplicative model, with alpha=0.09989 and theta=3.4481.
penetrance for each genotype will be as follow:
P(D|AABB)=0.0909,
P(D|AABb)=0.0909,
P(D|AAbb)=0.0909,
P(D|AaBB)=0.0909,
P(D|AaBb)=0.3079,
P(D|Aabb)=0.6645,
P(D|aaBB)=0.0909,
P(D|aaBb)=0.6645,
P(D|aabb)=0.9751.
Does anyone have any ideas?
Maybe it is a little difficult understand...
谢谢~
