MutSigCV output has p-value of 0
1
0
Entering edit mode
7.8 years ago
varsha619 ▴ 90

Hi, I am running MutSigCV 1.0 using the public server and noticed that the top gene in my list of significantly mutated genes has a p-value of 0. I was wondering if I can get the exact p-value with decimals instead of rounded off to 0. Please advise.

mutsigcv • 2.6k views
ADD COMMENT
3
Entering edit mode
7.8 years ago
poisonAlien ★ 3.2k

A p-value of 0 means the value is really low or close enough to zero. MutSig always produce 0 pvalue for overly significant genes such as TP53 in solid tumors. If you really want some value instead of zero, replace them with lowest machine epsilon values.

In r you could do this by

mutsig = read.delim('sig_genes.txt', sep = '\t', stringsAsFactors = FALSE, header = TRUE)
mach.epsi = .Machine$double.eps
mutsig$q = ifelse(test = mutsig$q == 0, yes = mach.epsi, no = mutsig$q)
ADD COMMENT
0
Entering edit mode

Got it, thank you so much!

ADD REPLY

Login before adding your answer.

Traffic: 787 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6