试用Scientific模块极不爽

niexining

UID: 23543
帖子: 2
积分: 4
在线时间: 10 分钟

1^# niexining 发表于 2007-06-05 18:30

试用Scientific模块极不爽

代码：

[Copy to clipboard] [ - ]

CODE:

# 试用Scientific模块
# 功能: 求解相关系数

from Scientific.Statistics import correlation
import numarray
data1=[1, 2, 3, 4, 5, 6, 4, 3]
data1=numarray.array(data1)
data2=2*data1
c=correlation(data1,data2)
print "相关系数为", c
raw_input()

运行结果：

QUOTE:

>>>
相关系数为 0.948683298051

众所周知，线性相关的两向量之相关系数为1
而模块给出如此不爽的答案。

结论为：
python 在科学计算领域非常差劲，装个样子而已，很不实用；
要做数学工作还是返回matlab的好。

niexining

UID: 23543
帖子: 2
积分: 4
在线时间: 10 分钟

2^# niexining 发表于 2007-06-08 23:01

用了numpy.corrcoef(data1, data2), 算得跟matlab一样好。

是谁说numarray 是要取代 numpy的？
我以亲身经历证明，无论 numarray, scientific 都不如 numpy 好用。
不仅仅是这个例子，还有 array()函数的使用，还有一个determinant()函数，这是我所用过的。