OPTpsider

OPTpsider computes derivative of psi function (second derivative of rho function) for optimal weight function

expand all in page

Syntax

Description

example

psider =OPTpsider(u, c) Plot of derivative of psi function.

Examples

expand all

  • Plot of derivative of psi function.

    • x=-6:0.01:6;
psiOPTder=OPTpsider(x,1.2);
plot(x,psiOPTder)
xlabel('x','Interpreter','Latex')
ylabel('$\psi''(x)$','Interpreter','Latex')
title('Optimal')

    Input Arguments

    expand all

    u — scaled residuals or Mahalanobis distances. Vector.

    vector containing residuals or Mahalanobis distances for the n units of the sample

    Data Types: single| double

    c — tuning parameters. Scalar.

    Scalar greater than 0 which controls the robustness/efficiency of the estimator (beta in regression or mu in the location case ...)

    Data Types: single| double

    Output Arguments

    expand all

    psider —Values of the derivative of the optimal psi function associated to the residuals or Mahalanobis distances for the n units of the sample. vector

    More About

    expand all

    Additional Details

    Function OPTpsider transforms vector x as follows \[ \psi'(x) = \begin{cases} \frac{2.7692}{c^2} \qquad |x| \leq \frac{2}{3} c \\ -\frac{5.3834}{c^2} +\frac{3*43.0672 x^2}{c^4} -\frac{5*69.9840 x^4}{c^6} +\frac{7*32.3 x^6}{c^8} \qquad \frac{2}{3} c < |x| \leq c \\ 0 & \; \vert x \vert > c. \end{cases} \]

    Remark: Optimal psi-function is almost linear around u = 0 in accordance with Winsor's principle that all distributions are normal in the middle.

    This means that \psi (u)/u is approximately constant over the linear region of \psi, so the points in that region tend to get equal weight.

    References

    Maronna, R.A., Martin D. and Yohai V.J. (2006), "Robust Statistics, Theory and Methods", Wiley, New York.

    See Also

    | |

    This page has been automatically generated by our routine publishFS