[Matplotlib-users] matplotlib.mlab.corrcoef() not working?

[Sebastian S. <st...@tu...>]

How is corrcoef() supposed to work? When I do the following:

---------------------------------------
from numarray import *
from numarray import random_array
from matplotlib.mlab import corrcoef
randn = random_array.standard_normal

x = randn((10,4))
r = corrcoef(x)
---------------------------------------

it gives:

corrcoef_test.py:4: DeprecationWarning: Non-ASCII character '\xc3' in 
file /usr/lib/python2.3/site-packages/matplotlib/mlab.py on line 1143, but no 
encoding declared; see http://www.python.org/peps/pep-0263.html for details
  from matplotlib.mlab import corrcoef
Numeric import failed... trying numarray.
Traceback (most recent call last):
  File "corrcoef_test.py", line 10, in ?
    r = corrcoef(x)
  File "/usr/lib/python2.3/site-packages/matplotlib/mlab.py", line 303, in 
corrcoef
    r = divide(C,denom)
  File "/usr/lib/python2.3/site-packages/numarray/generic.py", line 585, in 
_dualbroadcast
    raise ValueError("Arrays have incompatible shapes");
ValueError: Arrays have incompatible shapes


Any ideas, help?


(Just to note: There's another bug in mlab.py, in line 294 I think it should 
read "elif len(args)==1:" and not "elif len(args==1):". I changed that, 
otherwise corrcoef() would have bailed out way earlier.)


Thanks,
Sebastian

-- 
Sebastian Stark -- http://www.kyb.tuebingen.mpg.de/~stark
Max Planck Institute for Biological Cybernetics
Spemannstr. 38, 72076 Tuebingen
Phone: +49 7071 601 555 -- Fax: +49 7071 601 552