qq-plot not start from the 0 · Biostar

This is a test version of Biostars. For the public version, visit https://www.biostars.org.

qq-plot not start from the 0

Hi everyone, I did the GWAS analysis and got the result from TASSEL with GLM. Then I tried to make a qq-plot based on the p-value from GWAS output result. However, the qq-plot always started from y-axis 1 but not from 0 no matter how I modified the calculation. Following is the linkage for one of my qq-plot, please check it.

https://www.dropbox.com/s/7b1eo4flefb3q10/QQ_plot_Leaf_GWAS%20_%20miR396e-5p%20_%20length_specific_boxcox_transform_5genetic_2Latent.jpg?dl=0

Please, any suggestion?

Best, Shu-Yun Chen

qq-plot

• 2,618 views • link

updated 8.5 years ago by theobroma22

• written 8.5 years ago by ubs717 •

1 answer

It makes sense to me - guessing from your plot it shows that your pvalues range from 0.15 to 1E-7, right? If so, -log of 0.15 = 0.82 and -log of 1E-7 = 7. If you thought it should start from 0 means you have pvalues equal 1, cause -log of 1 = 0. Do you have pvalues equal to 1??

• 0 views • link

written 8.5 years ago by theobroma22

Log in to answer this question.

More posts like this

Which GWAS tool to use that can work on Czech cases and British controls, on a binary trait?
written 1.2 years ago by anita.szabo08 •

Hi, I am running a GWAS on 482 Czech cases and 1928 British controls for a binary trait (1:4 case - control ratio) (Attached is …
Tutorial: Single-cell RNA-seq: Annotation: Celltype auto annotation with SCSA
written 2.6 years ago by Julia Ma •

Part-2 here: https://www.biostars.org/p/9584987/ Part-3 here: https://www.biostars.org/p/9584989/ Full article lifted from: https://omicverse.readthedocs.io/en/latest/Tutorials-single/t_cellanno/ --- Single-cell transcriptomics allows the analysis of thousands of cells in a single experiment …
Tutorial: Unsupervised clustering on gene expression data
written 4.8 years ago by Hamid Ghaedi

Clustering is a data mining method to identify unknown possible groups of items solely based on intrinsic features and no external variables. Basically, clustering includes …
How to covert Plot_ly codes to ggplot2 codes?
written 5.3 years ago by anikcropscience •

Hello, I am using the following codes for producing an overlapping QQ-plot from two GWAS result files. # first model a<-read.table("GWAS_result_1.txt" , header=T) o = …
PCA with PLINK and SmartPCA: same imput file but different results in my plot. Why?
written 5.5 years ago by francesca94.minnai •

Hi everyone, and thank you in advance for any kind of help! I'm trying to perform a **Principal Component Analysis** both with **PLINK** and **smartPCA**. …
Manhattan plot with Bayes factor
written 7.1 years ago by biosol

Hi everyone, I'm trying to do a Manhattan plot with my GWAS SNPs for which I've calculated Bayes factor. I've seen that people usually transform …
Differences between RTCGA and TCGAbiolinks data
written 8.2 years ago by atakanekiz

Hello Biostars community, I'm trying to analyze miRNASeq data of TCGA melanoma (SKCM) samples. This is a prototype analysis in which I picked `hsa-mir-155` and …
Correcting for population structure in GWAS?
written 8.4 years ago by maya123z •

I'm running a GWAS for the first time and modeling my analysis based on [this paper][1], since the nature of our data is similar. The …
Clustering Million Points : Transcript Start And End Sites
written 14.3 years ago by Abhi

<p>Hey Guys</p> <p>I have about million points per chromosome that are basically putative start and end regions of transcripts. </p> <p>sample data here : <a …
Testing For Significance Of Observed Synteny
written 14.6 years ago by Eric Normandeau

<p>Hi,</p> <p>We have built a strongly supported map from RAD markers in a fish species and wanted to assess whether our map showed strong synteny …