## ----preamble,warning=F-------------------------------------------------- library(lfa) library(gcatest) dim(sim_geno) length(sim_trait) ## ----lfa----------------------------------------------------------------- LF = lfa(sim_geno, 3) dim(LF) ## ----gcat---------------------------------------------------------------- gcat_p = gcat(sim_geno, LF, sim_trait) ## ----gcat2--------------------------------------------------------------- gcat_p[1:5] ## ----fig1, fig.height=3-------------------------------------------------- library(ggplot2) dat = data.frame(gcat_p[6:10000]) colnames(dat) = "gcat_p" ggplot(dat, aes(gcat_p, ..density..)) + geom_histogram(binwidth=1/20) + theme_bw()