Beschreibung

Lizenz

Ich, der Urheberrechtsinhaber dieses Werkes, veröffentliche es als gemeinfrei. Dies gilt weltweit. In manchen Staaten könnte dies rechtlich nicht möglich sein. Sofern dies der Fall ist: Ich gewähre jedem das bedingungslose Recht, dieses Werk für jedweden Zweck zu nutzen, es sei denn, Bedingungen sind gesetzlich erforderlich.

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)")

 

Dieses Chart wurde mit R erstellt.