Description Usage Arguments Examples
Function which computes the weights for the iteratively reweighted LS estimate of the HBR fit.
1 |
x |
|
y |
|
percent |
|
ehat0 |
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 | ##---- Should be DIRECTLY executable !! ----
##-- ==> Define data, use random,
##-- or do help(data=index) for the standard data sets.
## The function is currently defined as
function (x, y, percent = 0.95, ehat0 = ltsreg(x, y)$residuals)
{
x = as.matrix(x)
n = dim(x)[1]
p = dim(x)[2]
robdis2 <- robdist.hbrfit(x)
cut = qchisq(percent, p)
sigma = mad(ehat0)
m = psi(cut/robdis2)
a = ehat0/(sigma * m)
c = (median(a) + 3 * mad(a))^2
h = sqrt(c)/a
hbrrewted = psi(abs(h))
mat <- cbind(hbrrewted, ehat0, robdis2)
return(mat)
}
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.