where L is the likelihood and edf is the equivalent degrees of
freedom (i.e., the number of free parameters for usual parametric
models) of fit. The backward stepwise algorithm starts with the full
model, then sequentially drops the variables that least hurt the fit,
as measured by AIC. Forward and backward stepwise procedures
usually result in similar models [31].