用python计算residuals

经常涉及到这些小操作，每次都没整理，下次用的时候又得百度，麻烦，干脆整理下，下次使用直接翻笔记了！

采用numpy库，如何计算residuals(x,y)？记得这里的x,y都是numpy.ndarray，假设x.shape=(3L,3L), y.shape=(3L,)

    >>> xarray([[1, 2, 3],[4, 5, 6],[7, 8, 9]])>>> yarray([2, 5, 8])

那怎么计算x与y之间的residuals呢？用到一些np中的简单操作，比如函数原型：numpy.vstack(tup)，等价于：np.concatenate(tup, axis=0） if tup contains arrays thatare at least 2-dimensional.

A=np.vstack([x,np.ones(x.shape[-1])]).T
>>> A
array([[1., 4., 7., 1.],[2., 5., 8., 1.],[3., 6., 9., 1.]])

直观的从例子可以看出A就是x另外加上全1的行，两个连接起来，然后转置。
4. 接下来使用lstsq(ndarray)来计算AB=y的最小二乘解:

B=np.linalg.lstsq(A,y)[0]
>>> B
array([ 2.57894737,  1.        , -0.57894737, -0.52631579])

不好理解了？再深入点，要拟合y=mx+c这条线了，所以要求解下AB=y的最小二乘解，实在不懂，参考numpy帮助，看看那个解释，偶就不想再琢磨了，哈哈，懒一点哇。找到m,c的值：

m=B[:-1]
>>> m 注意这是除常数的其它系数了
array([ 2.57894737,  1.        , -0.57894737])
c=B[-1]
>>> c  注意这是常数项系数
-0.5263157894736847

那显然可以知道，拟合出的线是啥（为了简单，我们就当成直线哇，和帮助吻合，看起来更容易懂）？pre=np.sum(m*x.T,axis=1)+c
原始的线是y, So:

pre=np.sum(m*x.T,axis=1)+c
>>> pre  注意这里m.shape=(3L,),x.shape=(3L,3L),这俩乘积是用m的每个元素分别乘以x的每一列，形成新的(3L,3L)矩阵，要注意理解背后的原理，别被忽悠了！
array([2., 5., 8.])

那最后两者之间的residuals就可以直接如下计算：
res=y-pre

用python计算residuals相关推荐

用python计算邮费考虑是否加急,用python计算residuals
经常涉及到这些小操作,每次都没整理,下次用的时候又得百度,麻烦,干脆整理下,下次使用直接翻笔记了! 采用numpy库,如何计算residuals(x,y)?记得这里的x,y都是numpy.ndarra ...
使用OpenCV和Python计算图像的“彩色度”
使用OpenCV和Python计算图像"彩色度" 1. 效果图 2. 炫彩度量方法是什么? 3. 源代码参考你是否尝试过计算每个图像的炫彩值,并根据炫彩值对自己的图像数据集进行 ...
Python计算训练数据集（测试集）中某个分类变量阴性(阳性)标签样本的不同水平（level）或者分类值的统计个数以及比例
Python计算训练数据集(测试集)中某个分类变量阴性(阳性)标签样本的不同水平(level)或者分类值的统计个数以及比例目录
Python计算两个numpy数组的交集（Intersection）实战：两个输入数组的交集并排序、获取交集元素及其索引、如果输入数组不是一维的，它们将被展平(flatten)，然后计算交集
Python计算两个numpy数组的交集(Intersection)实战:两个输入数组的交集并排序.获取交集元素及其索引.如果输入数组不是一维的,它们将被展平(flatten),然后计算交集目录
Python使用datetime中的timedelta模块实现时间增减：python计算100天后是哪年那月那日？
Python使用datetime中的timedelta模块实现时间增减:python计算100天后是哪年那月那日? 目录
python计算特征的统计值并文本输出
python计算特征的统计值并文本输出 # 输出统计分位数 df.describe(percentiles=[0.05, 0.25, 0.5, 0.75, 0.95]) featname:A => ...
python计算特征与目标的相关性并可视化
python计算特征与目标的相关性并可视化 pandas计算相关性 # 相关性 tips.corr() Out[2]: total_bill tip size total_bill 1.000000 ...
python计算时间差
python计算时间差 # 使用datetime 和 timedelta import datetime old_time = datetime.datetime.now() print(old_ti ...
python计算两个字典的相同点(从key、value、item角度)
python计算两个字典的相同点(从key.value.item角度) 字典是由一个一个的 key-value 构成的,key 是找到数据的关键,Python 对字典的操作都是通过 key 来完成的. ...

用python计算residuals

用python计算residuals相关推荐

最新文章

热门文章