Displaying 1 result from an estimated 1 matches for "xxlo".
Did you mean:
xxl
2011 Apr 05
0
Changing parameter in local fdr R code
...p0theo * f00)/f, 1)
f0p <- p0 * f0
if (nulltype == 0)
f0p <- p0theo * f00
F0l <- cumsum(f0p)
F0r <- cumsum(rev(f0p))
Fdrl <- F0l/Fl
Fdrr <- rev(F0r/Fr)
Int <- (1 - fdr) * f * (fdr < 0.9)
if (sum(x <= xmax & fdr == 1) > 0)
xxlo <- min(x[x <= xmax & fdr == 1])
else xxlo = xmax
if (sum(x >= xmax & fdr == 1) > 0)
xxhi <- max(x[x >= xmax & fdr == 1])
else xxhi = xmax
if (sum(x >= xxlo & x <= xxhi) > 0)
fdr[x >= xxlo & x <= xxhi] <- 1
i...