On the non-Gaussian errors in high-z supernovae type Ia data