A maximum likelihood (ML) method similar to that proposed
by McLeod and Hipel (1978) is used here to estimate
the scaling coefficient H of the equivalent Normal variates.
First, the data is detrended by using Sen’s non-parametric
trend estimator (Sen, 1968) which is given by