projects/01-gait1/R/plotqqman.R

library(qqman)

load("resultsF11.Rdata")


dataP <- A$snpf[,c("SNP", "pSNP", "pos", "chr")]
names(dataP)[2] <- "P"
names(dataP)[3] <- "BP"
dataP$CHR <- as.numeric(dataP$chr)



pos <- which(is.na(dataP$CHR))

dataPP<- dataP[-pos,]

png("plotF11.png")
manhattan(dataPP)
dev.off()

png("qqF11.png")
qq(dataPP$P)
dev.off()
ugcd/solarius documentation built on May 3, 2019, 2:22 p.m.