Question: RegioneR library in R
gravatar for panacotaforcota
7 months ago by
panacotaforcota0 wrote:


I try to use regioneR library in R for permutation test. I get result, where z-score is equal to 79 or -90 and p-value is equal to 0.002 for both z-scores. Number of permutations is 500. I don't understand how those z-score and p-value are calculated.


permutation test R • 218 views
ADD COMMENTlink modified 7 months ago by ATpoint41k • written 7 months ago by panacotaforcota0

Did you check their publication?

regioneR: an R/Bioconductor package for the association analysis of genomic regions based on permutation tests

ADD REPLYlink written 7 months ago by venu6.7k

Yes, and there are no explanation why for different z-score there are the same p-value.

ADD REPLYlink modified 7 months ago • written 7 months ago by panacotaforcota0

They are not the same z-score.

ADD REPLYlink written 7 months ago by Kevin Blighe67k
gravatar for ATpoint
7 months ago by
ATpoint41k wrote:

These are bootstraping p-values. The second formula on this page here gives a reference where this comes from.

From the source code of the regioneR package (permTest function) they implement it as:

pval <- (sum(orig.ev >= rand.ev, na.rm = TRUE) + 1) / (num.valid.values + 1)

So the p-value, here testing that the overlap (or whatever you test) is greater for the observed events compared to the randomly-permutated events, is the number of times that the observed event is larger or equal divided by the number of permutations.

The Z-score is the deviation of the observed events by the mean of the permutated events (corrected for standard deviation) so basically a significant event should have a large Z-score since the observed event should be quite different from the mean of the random events.

zscore <- round((orig.ev - mean(rand.ev, na.rm = TRUE))/stats::sd(rand.ev,  na.rm = TRUE), 4)
ADD COMMENTlink written 7 months ago by ATpoint41k
Please log in to add an answer.


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