powered by
Compute concave function values
compute_g(z, cfun, s, delta=0.0001)
Concave function values
vector nonnegative values from dfun, e.g., u^2/2
dfun
integer from 1-8, concave function as in irglm_fit
irglm_fit
a numeric value, see details in irglmreg_fit
irglmreg_fit
a positive small value, see details in irglmreg_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_g(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