R/mestse.R

mestse <-
function(x,bend=1.28,op=2){
#
#   Estimate the standard error of M-estimator using Huber's Psi
#   using estimate of influence function
#
n<-length(x)
mestse<-sqrt(sum((ifmest(x,bend,op=2)^2))/(n*(n-1)))
mestse
}
musto101/wilcox_R documentation built on May 23, 2019, 10:52 a.m.