powered by
Weight value from concave function
compute_wt(z, weights, cfun, s, delta=0.0001)
vector nonnegative values from dfun, e.g., u^2/2
dfun
optional numeric vector of weights.
integer from 1-8, concave function as in irglm_fit
irglm_fit
a numeric value, see details in irglm_fit
a positive small value, see details in irglm_fit
Zhu Wang <zwang145@uthsc.edu>
Zhu Wang (2024) Unified Robust Estimation, Australian & New Zealand Journal of Statistics. 66(1):77-102.
irglmreg
u <- seq(0, 4, by=0.01) z <- u^2/2 ### this is dfun res <- compute_wt(z, cfun=1, s=1) plot(z, res, ylab="Weight", type="l", lwd=2, main=expression(paste("hcave", " (", sigma, "=1)", )))
Run the code above in your browser using DataLab