Displaying 1 result from an estimated 1 matches for "robdist".
Did you mean:
probdist
2011 Dec 06
2
Why can't I figure this out? :S
...= quantile[2], nsamp = dist_nsamp)
rd <- sqrt(mahalanobis(X, cv$center, cv$cov))
rd_nok <- rd > critval[2]
## ROBUST results
nok <- rr_nok & rd_nok
NameC_rob <- lm(formula, data[!nok,])
rval <- list(ols = NameC, lts = NameC_lts, robust = NameC_rob,
cov.rob = cv, robresid = rr, robdist = rd,
high_residuals = rr[rr_nok], high_leverage = rd[rd_nok],
bad_leverage = nok, psamp = psamp, method = method,
nsamp = list(lts = nsamp, dist = dist_nsamp),
quantile = list(lts = quantile[1], dist = quantile[2]))
return(rval)
}
Error: object of type 'closure' is not subsettable
Any id...