我正在尝试像书中那样使线性回归与无信息的先验相适应，我的数字不匹配，我也不明白为什么。
我正在尝试使用qr分解或仅使用linalg.solve来找到线性回归后NIG的均值，因此我得到的数字与本书示例中的数字相同。
链接到数据：http://www.ceremade.dauphine.fr/~xian/BCS/caterpillar
我的代码是：
#mass - data and y - target(last columns of data)
w_=scipy.linalg.inv(mass.T@mass)@mass.T@y
and I get w :
0 8.243695
1 -0.002875
2 -0.035048
3 0.027882
4 -0.466332
5 0.103799
6 0.075530
7 -0.260833
8 0.023532
9 -0.886833
10 -0.332745
与书中相似但不相同