Nothing
plot.WGassociation <- function(x, ...){
if (!inherits(x, "WGassociation"))
stop("x must be an object of class 'WGassociation'")
xx <- data.frame(SNP=rownames(x), data.frame(x)[,2:6])
names(xx)[6] <- "additive"
dat <- tidyr::gather(xx, key="model", value="p.value", -"SNP")
dat$model <- factor(dat$model,
levels=c("codominant", "dominant",
"recessive", "overdominant",
"additive"))
plt <- ggplot(dat, aes(x=dat$SNP, y=-log10(dat$p.value))) +
geom_point() +
xlab("SNPs") + ylab(expression(-log[10]("p-value"))) +
theme(axis.text.x = element_text(angle = 90, hjust = 1)) +
geom_hline(aes(yintercept = -log10(0.05),
linetype = "Nominal"),
colour="blue") +
geom_hline(aes(yintercept = -log10(0.05/nrow(xx)),
linetype = "Bonferroni"),
colour="red") +
scale_linetype_manual(name = "Significance", values = c(2, 2),
guide = guide_legend(override.aes = list(color = c("red", "blue"))))
plt + facet_wrap( ~ model, ncol=1) +
theme(strip.text.x = element_text(face="bold")) +
theme(legend.position="top", legend.direction="horizontal")
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.