Displaying 1 result from an estimated 1 matches for "p0theo".
2011 Apr 05
0
Changing parameter in local fdr R code
...}
if (nulltype < 2) {
delhat = xmax = xmaxx = mlests[1]
sighat = mlests[2]
p0 = mlests[3]
f0 = dnorm(x, delhat, sighat)
f0 = (sum(f) * f0)/sum(f0)
}
fdr = pmin((p0 * f0)/f, 1)
f00 <- exp(-x^2/2)
f00 <- (f00 * sum(f))/sum(f00)
p0theo <- sum(f[i0])/sum(f00[i0])
fp0["thest", 3] = p0theo
fdr0 <- pmin((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 <...