h******1 发帖数: 94 | 1 If I run
> qobj = qvalue(p, lambda=seq(0, 0.90, 0.30), pi0.method="bootstrap")
[1] "ERROR: The estimated pi0 <= 0. Check that you have valid p-values or
use another lambda method."
> qobj = qvalue(p, lambda=seq(0, 0.90, 0.05), pi0.method="bootstrap")
[1] "ERROR: The estimated pi0 <= 0. Check that you have valid p-values or
use another lambda method."
> qobj = qvalue(p, lambda=seq(0, 0.90, 0.05), pi0.method="smoother")
[1] "ERROR: The estimated pi0 <= 0. Check that you have valid p-values or
use |
g********r 发帖数: 8017 | 2 这些p-value都非常小啊.qvalue是假定你有一个mixture.
一部分pvalue应当是uniform(0,1)
【在 h******1 的大作中提到】![](/moin_static193/solenoid/img/up.png) : If I run : > qobj = qvalue(p, lambda=seq(0, 0.90, 0.30), pi0.method="bootstrap") : [1] "ERROR: The estimated pi0 <= 0. Check that you have valid p-values or : use another lambda method." : > qobj = qvalue(p, lambda=seq(0, 0.90, 0.05), pi0.method="bootstrap") : [1] "ERROR: The estimated pi0 <= 0. Check that you have valid p-values or : use another lambda method." : > qobj = qvalue(p, lambda=seq(0, 0.90, 0.05), pi0.method="smoother") : [1] "ERROR: The estimated pi0 <= 0. Check that you have valid p-values or : use
|
h******1 发帖数: 94 | 3 You are right. If I add some big p-values to the vector, it runs correctly.
How to perform FDR correction to get q-values for small p-values?
【在 g********r 的大作中提到】![](/moin_static193/solenoid/img/up.png) : 这些p-value都非常小啊.qvalue是假定你有一个mixture. : 一部分pvalue应当是uniform(0,1)
|
g********r 发帖数: 8017 | 4 If your problem doesn't generate a mixture, then q-value doesn't help you.
The whole business of FDR build on the assumption that there are many null
in
the large number of tests you conduct.
If it does, just feed it all p-values without pre-selection.
correctly.
【在 h******1 的大作中提到】![](/moin_static193/solenoid/img/up.png) : You are right. If I add some big p-values to the vector, it runs correctly. : How to perform FDR correction to get q-values for small p-values?
|
g********r 发帖数: 8017 | 5 If the pvalues you listed are all you got, then maybe try p.adjust(, method=
"BH")
correctly.
【在 h******1 的大作中提到】![](/moin_static193/solenoid/img/up.png) : You are right. If I add some big p-values to the vector, it runs correctly. : How to perform FDR correction to get q-values for small p-values?
|
h******1 发帖数: 94 | 6 Yes, p.adjust is what I used before. But my supervisor asked me to use
pvalue.
I did preselection by fold change and this is why the p-values are so
small.
So I have to either qvalue() before filtering by fold change or p.adjust()
after filtering.
Thank you so much.
method=
【在 g********r 的大作中提到】![](/moin_static193/solenoid/img/up.png) : If the pvalues you listed are all you got, then maybe try p.adjust(, method= : "BH") : : correctly.
|
s*r 发帖数: 2757 | 7 you just have a problem with q0 estimation.
just manually insert a range of q0 |
h******1 发帖数: 94 | 8 How?
【在 s*r 的大作中提到】![](/moin_static193/solenoid/img/up.png) : you just have a problem with q0 estimation. : just manually insert a range of q0
|
s*r 发帖数: 2757 | 9 read the manual?
i do not understand your question and i forgot the detail of the package.
【在 h******1 的大作中提到】![](/moin_static193/solenoid/img/up.png) : How?
|
h******1 发帖数: 94 | 10 manual does not say any thing about presetting p0.
It only says
D. Why am I getting an error that says my pi_0 estimate is less than zero?
11
Usually this indicates that the p-values were incorrectly calculated.
Sometimes, π0
is very small, in which case this error can be avoided by using the
bootstrap
option rather than the smoother option.
【在 s*r 的大作中提到】![](/moin_static193/solenoid/img/up.png) : read the manual? : i do not understand your question and i forgot the detail of the package.
|
|
|
g********r 发帖数: 8017 | 11 Forget pre-selection and use qvalue on all the p-values.
If you did pre-selection, then all FDR generated from p.adjust are wrong.
()
【在 h******1 的大作中提到】![](/moin_static193/solenoid/img/up.png) : manual does not say any thing about presetting p0. : It only says : D. Why am I getting an error that says my pi_0 estimate is less than zero? : 11 : Usually this indicates that the p-values were incorrectly calculated. : Sometimes, π0 : is very small, in which case this error can be avoided by using the : bootstrap : option rather than the smoother option.
|
h******1 发帖数: 94 | 12 How come?
Then how do you know which tests are included in the multitest? How do you
define a group to perform a multi-test?
Could you provide any links or papers for your theory?
【在 g********r 的大作中提到】![](/moin_static193/solenoid/img/up.png) : Forget pre-selection and use qvalue on all the p-values. : If you did pre-selection, then all FDR generated from p.adjust are wrong. : : ()
|
g********r 发帖数: 8017 | 13 I don't have a theory. Just intuition. B-H procedure requires low dependency
.
I think by imposing fold change, you are selecting a subset with high
dependency.
you
【在 h******1 的大作中提到】![](/moin_static193/solenoid/img/up.png) : How come? : Then how do you know which tests are included in the multitest? How do you : define a group to perform a multi-test? : Could you provide any links or papers for your theory?
|
h******1 发帖数: 94 | 14 I cannot find the low dependency assumption for B-H method on internet.
I need to show it to my boss.
dependency
【在 g********r 的大作中提到】![](/moin_static193/solenoid/img/up.png) : I don't have a theory. Just intuition. B-H procedure requires low dependency : . : I think by imposing fold change, you are selecting a subset with high : dependency. : : you
|
s*r 发帖数: 2757 | 15 he said it is intuition
【在 h******1 的大作中提到】![](/moin_static193/solenoid/img/up.png) : I cannot find the low dependency assumption for B-H method on internet. : I need to show it to my boss. : : dependency
|
h******1 发帖数: 94 | 16 If I include all genes in the study, all of the q-values turn out to be
very large.
dependency
【在 g********r 的大作中提到】![](/moin_static193/solenoid/img/up.png) : I don't have a theory. Just intuition. B-H procedure requires low dependency : . : I think by imposing fold change, you are selecting a subset with high : dependency. : : you
|
s*r 发帖数: 2757 | 17 haha
【在 h******1 的大作中提到】![](/moin_static193/solenoid/img/up.png) : If I include all genes in the study, all of the q-values turn out to be : very large. : : dependency
|
g********r 发帖数: 8017 | 18 If you do a histogram of all p-values, does it look like uniform ? If so,
the
q-values should all be large. Claiming any significant gene with a twisted
FDR is 忽悠.
The BH method assumes independent tests. It is in their paper.
【在 h******1 的大作中提到】![](/moin_static193/solenoid/img/up.png) : If I include all genes in the study, all of the q-values turn out to be : very large. : : dependency
|
h******1 发帖数: 94 | 19 The histogram looks like the attached file.
But the genes are not independent after filtering with fold change.
What can I do?
【在 g********r 的大作中提到】![](/moin_static193/solenoid/img/up.png) : If you do a histogram of all p-values, does it look like uniform ? If so, : the : q-values should all be large. Claiming any significant gene with a twisted : FDR is 忽悠. : The BH method assumes independent tests. It is in their paper.
|
h******1 发帖数: 94 | 20 It cannot be attached.
It is quite uniform with a small peak on the front middle.
【在 h******1 的大作中提到】![](/moin_static193/solenoid/img/up.png) : The histogram looks like the attached file. : But the genes are not independent after filtering with fold change. : What can I do?
|
|
|
g********r 发帖数: 8017 | 21 If you look at the qvalue paper by Storey, it is clear in what case
you can get small qvalue - when the histogram of p-value has
a big peak near zero.
【在 h******1 的大作中提到】![](/moin_static193/solenoid/img/up.png) : It cannot be attached. : It is quite uniform with a small peak on the front middle.
|
h******1 发帖数: 94 | 22 Yes, I saw it.
But I don't have the ideal histogram for my data, so I probably cannot use
Storey's method.
In order to use B-H method, I have to filter the genes by fold change, in
order to get small q-values, but the preselection could give only dependent
genes in the study.
Do I have any other options?
【在 g********r 的大作中提到】![](/moin_static193/solenoid/img/up.png) : If you look at the qvalue paper by Storey, it is clear in what case : you can get small qvalue - when the histogram of p-value has : a big peak near zero.
|
g********r 发帖数: 8017 | 23 not every experiment generates expected results you know. What if there's no
real biological difference?
Using a stringent p-value cutoff + a fold change cutoff might work, because
p-value or q-value don't cover absolute amount of change.....and in real
biology each condition may not be a single population blah blah blah ....
use
in
dependent
【在 h******1 的大作中提到】![](/moin_static193/solenoid/img/up.png) : Yes, I saw it. : But I don't have the ideal histogram for my data, so I probably cannot use : Storey's method. : In order to use B-H method, I have to filter the genes by fold change, in : order to get small q-values, but the preselection could give only dependent : genes in the study. : Do I have any other options?
|
g********r 发帖数: 8017 | 24 Could consider trying more liberal methods, such as local FDR. |
h******1 发帖数: 94 | 25 Thank you so much goldmember.
I think you are very knowledgeable.
【在 g********r 的大作中提到】![](/moin_static193/solenoid/img/up.png) : Could consider trying more liberal methods, such as local FDR.
|