Closed yoonsucho closed 3 years ago
library(simulateGP) library(tidyverse)
map <- tibble(snp=1:100000, af=runif(100000, 0.01, 0.99)) ss <- generate_gwas_params(map, h2=0.5, S=1)
d1 <- generate_gwas_ss(ss, nid=100000) d2 <- generate_gwas_ss(ss, nid=10000)
library(simulateGP) library(tidyverse)
map <- tibble(snp=1:100000, af=runif(100000, 0.01, 0.99)) ss <- generate_gwas_params(map, h2=0.5, S=1)
d1 <- generate_gwas_ss(ss, nid=100000) d2 <- generate_gwas_ss(ss, nid=10000)