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()
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.