V Matlab nápovědě k funkci „polyfit“ se píše:
„If the data Y are random, an estimate of the covariance matrix of P (Rinv*R
inv‘)*normr^2/df, where Rinv is the inverse of R.“
Přepsání výrazu má vypadat takto:
ste = sqrt(diag(inv(s.R)*inv(s.R‘)).*s.normr.^2./s.df)