Given probability weighted moments \(\beta_0, \beta_1, \beta_2, \beta_3\),
the first four sample L-moments are:
\(l_1 = \beta_0\)
\(l_2 = 2\beta_1 - \beta_0\)
\(l_3 = 6\beta_2 - 6\beta_1 + \beta_0\)
\(l_4 = 20\beta_3 - 30\beta_2 + 12\beta_1 - \beta_0\)
Then, the sample L-skewness is \(t_3 = l_3 / l_2\) and the sample L-kurtosis
is \(t_4 = l_4 / l_2\).