[pdl-stats] A PDL::Stats change I'd love to see

[<dw...@em...>]

Thank you Maggie.



What I'd like to be able to do, with things like the logistic function (any curve fitter
in PDL::Stats, actually) is something like this:


my %m = $y->logistic( COV => 1 );
print $m{cov}->corr_table();


I don't think logistic should return the covariance matrix by default, but I would like
to be able to get to it.


A single line patch in GLM.pm, after %ret is declared, around line 2000 or so..


$ret{cov} = $cov if ( $opt{COV} );


would do the trick for logistic, I believe. I don't think $cov is out of scope by that point.


Sincerely,


David Myers.