Hi zcurve developer I am curious about how to get the false discover

False discovery rate about zcurve HOT 1 OPEN

Yefeng0920 commented on June 26, 2024

False discovery rate

from zcurve.

Comments (1)

FBartos commented on June 26, 2024

Hi Yefeng,

You can obtain the false discovery risk (Soric's FDR) estimate by adding all = TRUE argument to the summary function.
(I recognize that this should be a bit better documented)

> library(zcurve)
> fit <- zcurve(OSC.z)
> summary(fit, all = TRUE)
Call:
zcurve(z = OSC.z)

model: EM via EM

              Estimate  l.CI   u.CI
ERR              0.615 0.458  0.746
EDR              0.383 0.072  0.691
Soric FDR        0.085 0.023  0.675
File Drawer R    1.608 0.446 12.826
Expected N         222   123   1175
Missing N          132    33   1085

Model converged in 53 + 62 iterations
Fitted using 73 z-values. 90 supplied, 85 significant (ODR = 0.94, 95% CI [0.87, 0.98]).
Q = -60.61, 95% CI[-71.04, -47.86]

Thanks for the suggestion for a comparison with the swfdr---it would make for a very nice vignette! I will keep it in mind when having some free time over the weekends :)

Cheers,
Frantisek

from zcurve.

Recommend Projects