Question: QQ plot - how to add gris shadow showing confidence interval
0
gravatar for aleix.arnau1990
5.7 years ago by
United Kingdom
aleix.arnau199010 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 COMMENTlink modified 5.7 years ago by PoGibas4.8k • written 5.7 years ago by aleix.arnau199010
0
gravatar for Ahill
5.7 years ago by
Ahill1.9k
United States
Ahill1.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 COMMENTlink written 5.7 years ago by Ahill1.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 REPLYlink written 5.7 years ago by aleix.arnau199010

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

ADD REPLYlink written 5.7 years ago by aleix.arnau199010
Please log in to add an answer.

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
_