Entering edit mode
6.8 years ago
bha
▴
80
Same question is asked at stackoverflow, but other members asking to put it at this forum. sorry for repetition.
I want to generate some phenotype data, example trait could be human height. The expression i want to use should be as:
P corresponds to phenotype
P ~ N (∑XiBi, δp)
i couldn't directly pasted above equation here,
(∑) varies from i=1 to L, is number of loci contribute to trait, could be 40
(Bi) is beta effect size
heritability should be 40%
(δp) standard deviation for phenotype
Do apology if i make this complicated. Can i generate this phenotype in R? if so, then what's procedure should be?
What is your intended use for this? - why do you want to do it in R? There is minimal context provided in yuor question.
Apologies for late reply. Actually i wanted to simulate the phenotype, but now i found some programmes which can do that, for instances, gcta, dissect etc. Many thanks for your message1