==Family-based Association Exercise==
cd exercises/cordell/FASTLMM_GCTA
ls -l
head(res1)
source("qqmanHJCupdated.R")
jpeg("mh1.jpeg", height=500, width=800)
manhattan(res1, pch=20, suggestiveline=F, genomewideline=F, ymin=2, cex.x.axis=0.65, colors=c("black","dodgerblue"), cex=0.5)
dev.off()
jpeg("qqplot1.jpeg", height=500, width=500)
qq(res1$P)
dev.off()
chi<-(qchisq(1-res1$P,1))
lambda=median(chi)/0.456
names(new)<-c("SNP", "CHR", "BP", "P")
head(new)
jpeg("qqplot2.jpeg", height=500, width=500)
qq(new$P)
dev.off()
jpeg("mh2.jpeg", height=500, width=800)
manhattan(new, pch=20, suggestiveline=F, genomewideline=F, ymin=2, cex.x.axis=0.65, colors=c("black","dodgerblue"), cex=0.5)
dev.off()