Displaying 1 result from an estimated 1 matches for "rr_nok".
Did you mean:
rd_nok
2011 Dec 06
2
Why can't I figure this out? :S
...- rep(quantile, length.out = 2)
if(is.null(psamp)) psamp <- quantile[1]
## LTS results with robust residuals
NameC_lts <-
lqs(gpanew~female+female:lastinit+agenew+canadian+mom_ed+yearstudy,
quantile = quantile[1], psamp = psamp, nsamp = nsamp)
rr <- residuals(NameC_lts)/NameC_lts$scale[2]
rr_nok <- abs(rr) > critval[1]
## robust leverage via MCD (or MVE)
X <- model.matrix(NameC)[,-1]
cv <- cov.rob(X, method = method,
quantile = quantile[2], nsamp = dist_nsamp)
rd <- sqrt(mahalanobis(X, cv$center, cv$cov))
rd_nok <- rd > critval[2]
## ROBUST results
nok <- rr_nok &am...