Displaying 1 result from an estimated 1 matches for "efdr".
Did you mean:
efd
2011 Apr 05
0
Changing parameter in local fdr R code
...f (nulltype == 1) {
fdr[x >= mlests[1] - mlests[2] & x <= mlests[1] + mlests[2]] = 1
fdr0[x >= mlests[1] - mlests[2] & x <= mlests[1] + mlests[2]] = 1
}
p1 <- sum((1 - fdr) * f)/N
p1theo <- sum((1 - fdr0) * f)/N
fall <- f + (yall - y)
Efdr <- sum((1 - fdr) * fdr * fall)/sum((1 - fdr) * fall)
Efdrtheo <- sum((1 - fdr0) * fdr0 * fall)/sum((1 - fdr0) *
fall)
iup <- (1:K)[x >= xmax]
ido <- (1:K)[x <= xmax]
Eleft <- sum((1 - fdr[ido]) * fdr[ido] * fall[ido])/sum((1 -
fdr[ido]) * fall[id...