Displaying 1 result from an estimated 1 matches for "namec_lts".
2011 Dec 06
2
Why can't I figure this out? :S
...default: choose psamp = quantile
n <- length(residuals(NameC))
k <- length(coef(NameC))
if(is.null(quantile)) quantile <- c(floor((n + k + 1)/2),
floor((n + k)/2))
quantile <- 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...