Plotting positions close to the exact unbiased solution: application to the Pozzuoli's bradeysism earthquake data