Deriving p values for gene correlations by permutation analysis · Biostar

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

Deriving p values for gene correlations by permutation analysis

Hi

I am calculating gene correlations for all patients gene expression data against clinical parameters. Then I want to derive p values by permutation. Is it OK to randomly permute the sample columns and then re calculate the correlation for every gene N times. Then ask how many times are correlations observed equal to or greater than a genes correlation by chance?

E.g. If my correlation is 0.3

Then the p value would be, n obs > 0.3/ total obs

So I have one distribution per gene?

Thanks,

Robert

rna-seq next-gen

• 3,150 views • link

updated 4.1 years ago by Ram

• written 10.5 years ago by chris86

It sounds about right to me. Have a look also at the boot package in R (and this FAQ http://www.ats.ucla.edu/stat/r/faq/boot.htm). Once you have the p-values, I guess you should correct them for the multiple tests you performed (i.e. one test per gene).

• 0 views • link

updated 4.1 years ago by Ram

• written 10.5 years ago by dariober

Yeah I hadn't done multiple testing correction, that is a good point

• 0 views • link

written 10.5 years ago by chris86

0 answers

No answers yet.

Log in to answer this question.

More posts like this

Permutation test for convergent gene expression changes
written 2 months ago by Iulia •

Hello, I have 4 independent species pairs (i.e., 4 x Species A vs. Species B contrasts). For each species, I have 4 biological replicates with …
Multiple correlations RNAseq data - correction for multiple testing?
written 4.5 years ago by psm

Hello all, I have performed bulk RNAsequencing of tissue (e.g. kidney) from patients with a disease (e.g. diabetes), and I would like to identify genes …
Correlation analyses on TCGA data
written 5.2 years ago by mario.keller.1988 •

Good morning, I also made the same post in the Bioconductor Forum. From our experiments in cell lines we identified genes and exons differentially regulated …
Statistical test for QQ-plot comparison of GWAS summary statistics
written 5.3 years ago by f-rasmussen •

I haver performed a GWAS and I want to know if there are significantly more significant p-values than I would expect by random chance. Normally …
Tutorial: Finding the significance of the overlap between 2 or more gene sets using simulation in R.
written 5.9 years ago by rpolicastro

TLDR: Example R function to calculate significance of overlap of 2 or more gene sets. `genes_all` is a vector that contains all genes, and `gene_sets` …
Correlating Correlations: Conceptual statistics questions behind correlating multi-dimensional data
written 6.6 years ago by sako242 •

Hi all, I have an RNAseq and a proteomics dataset for patients along some disease spectrum vs control samples (n = 50+). My lab has …
Statistical test for finding driver mutations.
written 9.4 years ago by Gene_MMP8

One approach to identify driver mutations that drive cancer progression is to look for recurrent mutations. Now BMR (background mutation rate) gives the frequency of …
Calculating Expression Conservation Score Between Two Gene Co-Expression Networks
written 12.9 years ago by Bioinfo_2006

<p>Suppose these two matrices </p> <pre><code> set.seed(1) x <- matrix(rnorm(1000), nrow=100, ncol=10) y <- matrix(rnorm(1000), nrow=100, ncol=10) </code></pre> <p>represent the genes (as rows) and genotypes …
Bootstrap On Correlation Coefficients In R
written 12.9 years ago by suganyasivagurunathan •

<p>Bootstrapping on a correlation is useful because we know that the distribution of correlations is not normal since it's bounded between -1 and 1. In …
Permutation Method In Cancer Data Set (Same As Mutsic)
written 13.0 years ago by imagineyd

<p>I have one question about permutation method that used in TCGA paper (<strong>Somatic mutations affect key pathways in lung adenocarcinoma</strong>, 2008, Nature) and MutSic method.</p> …