| # QQ plot using GenABEL estlambda function | |
| png(paste(target, "_QQplot.png", sep = ""), width = 500, height = 500) | |
| lambda <- estlambda(GWASout$t.value**2, plot = T, method = "median") | |
| dev.off() |
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment