python类,使用numpy矩阵求逆实现线性模型的一般最小二乘拟合
linearfit的Python项目详细描述
linear fit
=
**python类使用numpy矩阵求逆实现线性模型的一般最小二乘拟合。**
所有输入都必须是numpy矩阵。
math基于press'
numeric receipes p661:第15.2节将数据拟合到直线
numeric receipes p671:第15.4节general linear least squares
code基于damien segransan的早期yorick实现
(日内瓦大学)
请参阅test_linearfit.py中的测试,以获取示例用法。
因变量
>y=NP.数组([5.9,5.4,5.4,4.4,4.4,4.6,3.5,3.5,3.7,2.7,2.8,2.8,2.4,4.4,4.4,4.4,4.4,4.br/>
>因变量权重
y=NP.数组([1,1.8,4,8,8,20,20,70,70,70,70,70,70,70,70>
>准备矩阵
>m=NP mat(y);
y>m=NP因变量的对角协方差矩阵自变量中,只有一个变量
c=np.mat(np.vstack(np.vstack(len(x))(np.one(len(x))(x))
;初始化对象
res=linearfit.linearfit(m,s,c);
;做拟合
>res fit
>res displa结果
`
文档
,我们所有的类和方法/函数都是
包括基本文件。
安装说明
----
此软件包是在python 2.7环境中开发的,但也使用python 3.5进行了成功的测试。
需要以下python软件包:
*numpy
``pip install linearfit``
或获取源文件,例如:
``git clone https://github.com/johannes sahlmann/linearfit````
install pygacs:
``cd linearfit````
``python setup.py install--user```
运行测试,做:
`` python test_linearfit.py```
许可证
----
linearfit是开源和自由软件:您可以根据
自由软件基金会发布的GNU Lesser General Public许可证的条款重新分发和/或修改
,许可证的第3版,或者(在您的选择)
任何以后的版本。
BR/>此程序被分发,希望它是有用的,但没有任何
保修;甚至没有隐含的保证适销性或适合于
特定用途。有关详细信息,请参阅GNU通用公共许可证。
您应该已经收到GNU通用公共许可证的副本以及此程序。如果没有,请参见http://www.gnu.org/licenses/
。:变更日志:
0.1(2017年5月)
++++++++++++++