Yes, thanks Giovanni. Numpy worked.
On Jul 21, 1:51 am, "giovanni.marche...@ds.unifi.it"
<giovanni.m.marche...@gmail.com> wrote:
> I read in the listing of the function
>
> warnings.warn("""\
> scipy.stats.corrcoef is deprecated; please update your code to use
> numpy.corrcoef.
> Please note that:
> - numpy.corrcoef rowvar argument defaults to true, not false
> - numpy.corrcoef bias argument defaults to false, not true
> """, DeprecationWarning)
> if y is not None:
> x = np.transpose([x,y])
> y = None
> c = cov(x, y, rowvar=rowvar, bias=bias)
> d = np.diag(c)
> return c/np.sqrt(np.multiply.outer(d,d))
>
> Anyway
> bias is used only for the covariance computation.
> In this case if bias = True the denominator is set to
> n-1.
>
> As the correlation coefficient is cov(x,y)/ sqrt( cov(x,x)* cov(y,y))
> the denominator choice is irrelevant as it cancels out.
>
> -- Giovanni Marchetti
--~--~---------~--~----~------------~-------~--~----~
To post to this group, send email to sage-support@googlegroups.com
To unsubscribe from this group, send email to
sage-support-unsubscr...@googlegroups.com
For more options, visit this group at
http://groups.google.com/group/sage-support
URLs: http://www.sagemath.org
-~----------~----~----~----~------~----~------~--~---
