Hello,
I think it might be a very silly question, therefore sorry about my ignorance.
Imagine a RNA-seq data set with n samples, and I want to determine the correlation between each other.
After count normalization, we can perform a correlation between each sample, using Pearson, Spearman, or other one. In a gene to gene correlation I do can understand how these tests works.
But if I have a huge data set (n samples) with m genes, how the correlation test works? Can anybody answer me please?
Regards,
1 answer
Correlation is a pairwise measure. You can calculate correlation between two samples (by considering each gene as an independent observation) but not between n samples. To get a sense of the overall concordance of your dataset, I would calculate all pairwise correlations, which would generate a symmetric nxn correlation matrix, and should identify any outlier samples. The corrplot function in R provides a nice heatmap-style visualization of this.
•
written
by
ethan.kaufman
Log in to answer this question.
-
PCA, Pearson correlation
written by Irving •
-
RNA seq normalization and gene gene correlation
written by Fish •
-
Can I calculate the correlation of hazard ratios of different features for comparative purposes?
written by joker33 •
-
Correlation in RNAseq
written by feng920308 •
-
How to correlate phenotype values with gene expression levels
written by afli
-
Comparing similarity of RNAseq datasets
written by m.laarman •
-
How to correlate chip seq peaks?
written by Zee_S
-
FastQTL-Finding a candidate QTL per phenotype
written by ATCG
-
How I can compute Pearson or Spearman correlation in R?
written by shahryary •
-
EDASeq - Problem "longer object length is not a multiple of shorter object length"
written by kaihami •
