The HPSEVERITY Procedure

Example 9.7 Estimating Parameters Using Cramér-von Mises Estimator

PROC HPSEVERITY enables you to estimate model parameters by minimizing your own objective function. This example illustrates how you can use PROC HPSEVERITY to implement the Cramér-von Mises estimator. Let $F(y_ i; \Theta )$ denote the estimate of CDF at $y_ i$ for a distribution with parameters $\Theta $, and let $F_ n(y_ i)$ denote the empirical estimate of CDF (EDF) at $y_ i$ that is computed from a sample ${y_ i}$, $1 \leq i \leq N$. Then, the Cramér-von Mises estimator of the parameters is defined as

\[  \hat{\Theta } = \arg \min _{\Theta } \sum _{i=1}^{N} (F(y_ i; \Theta ) - F_ n(y_ i))^2  \]

This estimator belongs to the class of minimum distance estimators. It attempts to estimate the parameters such that the squared distance between the CDF and EDF estimates is minimized.

The following PROC HPSEVERITY step uses the Cramér-von Mises estimator to fit four candidate distribution models, including the LOGNGPD mixed-tail distribution model that was defined in Defining a Model for Mixed-Tail Distributions. The input sample is the same as is used in that example.

/*--- Set the search path for functions defined with PROC FCMP ---*/
options cmplib=(work.sevexmpl);

/*-------- Fit LOGNGPD model with PROC HPSEVERITY by using -------
  -------- the Cramer-von Mises minimum distance estimator -------*/
proc hpseverity data=testmixdist obj=cvmobj print=all;
   loss y;
   dist logngpd burr logn gpd;

   * Cramer-von Mises estimator (minimizes the distance *
   * between parametric and nonparametric estimates)    *;
   cvmobj = _cdf_(y);
   cvmobj = (cvmobj -_edf_(y))**2;
run;

The OBJ= option in the PROC HPSEVERITY statement specifies that the objective function cvmobj should be minimized. The programming statements compute the contribution of each observation in the input data set to the objective function cvmobj. The use of keyword functions _CDF_ and _EDF_ makes the program applicable to all the distributions.

Some of the key results prepared by PROC HPSEVERITY are shown in Output 9.7.1. The "Model Selection" table indicates that all models converged. When you specify a custom objective function, the default selection criterion is the value of the custom objective function. The "All Fit Statistics" table indicates that LOGNGPD is the best distribution according to all the statistics of fit. Comparing the fit statistics of Output 9.7.1 with those of Output 9.3.1 indicates that the use of the Cramér-von Mises estimator has resulted in smaller values for all the EDF-based statistics of fit for all the models, which is expected from a minimum distance estimator.

Output 9.7.1: Summary of Cramér-von Mises Estimation

The HPSEVERITY Procedure

Input Data Set
Name WORK.TESTMIXDIST
Label Lognormal Body-GPD Tail Sample

Model Selection
Distribution Converged cvmobj Selected
logngpd Yes 0.02694 Yes
Burr Yes 0.03325 No
Logn Yes 0.03633 No
Gpd Yes 2.96090 No

All Fit Statistics
Distribution cvmobj -2 Log
Likelihood
AIC AICC BIC KS AD CvM
logngpd 0.02694 * 419.49635 * 429.49635 * 430.13464 * 442.52220 * 0.51332 * 0.21563 * 0.03030 *
Burr 0.03325   436.58823   442.58823   442.83823   450.40374   0.53084   0.82875   0.03807  
Logn 0.03633   491.88659   495.88659   496.01030   501.09693   0.52469   2.08312   0.04173  
Gpd 2.96090   560.35409   564.35409   564.47780   569.56443   2.99095   15.51378   2.97806  
Note: The asterisk (*) marks the best model according to each column's criterion.