Simulate a phenotype

NQTN = 10
h2 = 0.8
alpha = 0.6
distribution = "pseudo_geom"
set.seed(1337) # Fix randomization
sim_pheno <- G2P(genotypes, h2, alpha, NQTN, distribution) # Simulate a phenotype

Investigate properties of the simulated phenotype

## [1] "GWAS correlation = 1"