Displaying 1 result from an estimated 1 matches for "cmerror".
Did you mean:
vmerror
2011 Apr 05
0
Changing parameter in local fdr R code
...lo0 & x < hi0]
x0 <- x[i0]
y0 <- l[i0]
if (nulltype == 3) {
X00 <- cbind((x0 - xmax)^2, pmax(x0 - xmax, 0)^2)
}
else {
X00 <- cbind(x0 - xmax, (x0 - xmax)^2)
}
lr <- lm(y0 ~ X00)
co <- lr$coef
if (nulltype == 3) {
cmerror = I(is.na(co[3]) | is.na(co[2]))
if (!cmerror)
cmerror = I(co[2] >= 0 | co[2] + co[3] >= 0)
}
else {
cmerror = is.na(co[3])
if (!cmerror)
cmerror = I(co[3] >= 0)
}
if (cmerror) {
if (nulltype == 3)
stop(&qu...