平时用lm()用predict()习惯了，没想到里面这么多复杂的东西要考虑，算法本身是倒腾点矩阵乘法配置个迭代就完了，前前后后的涉及到好多材料啊。前段时间看@谢益辉 SAMA在微博上说“得设计阵者，得半边（传统）统计天下”，诚不我欺。
