On Tue, 27 May 2008, Stefano Balietti wrote: Too complicated, IMO. A case for do-it-yourself. After mle: matrix stats = $coeff ./ $stderr stats = stats ~ (2 * pvalue(z, abs(stats))) colnames(stats, "t-ratio p-value") printf "%12.5f\n", stats Glad you asked, because it made me realize that the use of pvalue() with a matrix argument was broken by a recent "fix" -- but this now works properly in CVS. Allin.