Optimize lambda0 and pi0, using Poisson version of analytical approach in Efron's "Microarrays, Empirical Bayes, and the Two-Groups Model"
In the paper, he derives formulas for arbitrary support. We use support [0, 1] for estimating the null. We initialize the search with lambda=0.05 and pi=0.95.
1 | AnalyticalOptim(r0, r1, N)
|
r0 |
The number of 0's in the data. |
r1 |
The number of 1's in the data. |
N |
The total number of points in the data. |
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.