Latest
Open
RNA-Seq
ChIP-Seq
SNP
Assembly
Tutorials
Tools
Jobs
Forum
Planet
All »
Home
Log In
Question: QQ plot - how to add gris shadow showing confidence interval
0
gravatar for aleix.arnau1990
5.7 years ago by
aleix.arnau1990 • 10
United Kingdom
aleix.arnau1990 • 10 wrote:

 

Hi! I am trying to plot a QQplot with results from a gwas runned with plink. I would like to know how can I add the tipical grey shadow showing confidence interval (95%CI). I can not find it. Does someone know?

My plot looks like that:

assoc<-read.table("all_SNPs_assoc.txt",header=TRUE,dec=".")
png("QQplot.png")
a<-assoc$SNP
n<-length(a)
the<-(-1*log(pchisq(rchisq(n,df=1),df=1,lower.tail=F),base=10))
# pchisq gives de distribution function
# rchisq generates random deviates
act<-(-1*log(assoc$P[a],base=10))
qqplot(the,act,main="QQplot assoc",xlab="Expected -log p",ylab="Observed -log p")
qqline(probs=0.05)
abline(0,1)
dev.off()

Thanks!

 
qqplot R • 4.1k views
ADD COMMENT • link •
modified 5.7 years ago by PoGibas ♦ 4.8k • written 5.7 years ago by aleix.arnau1990 • 10
0
gravatar for Ahill
5.7 years ago by
Ahill • 1.9k
United States
Ahill • 1.9k wrote:

See the R snpStats package.  In the manual, page 29, there is a function qq.chisq() which will provide a QQ plot including a grey confidence band.  The confidence band is added using the polygon() function.

 
ADD COMMENT • link written 5.7 years ago by Ahill • 1.9k

 

Thanks Ahill! But now, with my code like this:

assoc<-read.table("all_SNPs_assoc.txt",header=TRUE,dec=".")
png("QQplot.png")
P<-assoc_allSNPs_mperm$P
qq.chisq(-2 * log(P), df = 2, conc=c(0.05, 0.95))
abline(0,1)
dev.off()

It provides a nice plot, but with x and y axis in chi-squared values and I need it in expected -log pvalues. In my data I just have p-values that are the values that I give to variable P. How can I have the x and y axis with the original p-values and not converted to chi-squared values? Or there is some kind of qq.pval() function? thanks.

 

ADD REPLY • link written 5.7 years ago by aleix.arnau1990 • 10

Don't worry! I found GWASTools package from bioconductor which solved my problem. Thanks for all.

ADD REPLY • link written 5.7 years ago by aleix.arnau1990 • 10
Please log in to add an answer.

Similar posts • Search »

Content
Help
Access
Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1656 users visited in the last hour
_