File:Mest weightfunc.jpg

Summary[edit]

Licensing[edit]

I, the copyright holder of this work, release this work into the public domain. This applies worldwide. In some countries this may not be legally possible; if so: I grant anyone the right to use this work for any purpose, without any conditions, unless such conditions are required by law.

z <- (-1000:1000)/100
# Hubert k
k  <- 1.339
wh <- 1*(abs(z)<=k)+
      k*sign(z)/z*(abs(z)>k)

# Hampel
a   <- 1.7
b   <- 3.4
c   <- 8.5
wha <- 1*(abs(z)<=a)+
       a*sign(z)/z*((abs(z)>a) & (abs(z)<=b))+

       a*(c-abs(z))/(c-b)*sign(z)/z*((abs(z)>b) & (abs(z)<=c))
# Andrews wave
a   <- 1.34*pi
wav <- a/(pi*z)*sin(pi*z/a)*(abs(z)<=a)

# Tukey's biweight
a   <- 4.685
wtb <- (1-z^2/a^2)*(1-z^2/a^2)*(abs(z)<=a)

plot(z, wh, type="l", main="Hubert (k=1.339)")
plot(z, wha, type="l", main="Hampel (a=1.7, b=3.4, c=8.5)")
plot(z, wav, type="l", main="Andrews wave (a=1.34*pi)")
plot(z, wtb, type="l", main="Tukey's biweight (a=4.685)")

 

This chart was created with R.