発生している問題
pythonで重回帰分析のプログラムを書いているのですが、なかなか思った結果になりません。結果的に0よりも小さな値が出ることを望んでいるのですが、非常に大きい値が出てしまいます。計算に使っている配列の要素はすべて0以下なので、このように値が大きくなるとは思えません。
該当のソースコード
python
1sigma = np.cov(sample[:,2:4], rowvar=0, bias=0) 2print("sigma",sigma) 3print("X2",X2[:,1:3]) 4h = (X2[:,1:3] @ np.linalg.inv(sigma) @ X2[:,1:3].T) 5print(h)
出力
sigma [[5.31666667e-01 4.38888889e-04]
[4.38888889e-04 9.88888889e-07]]
X2 [[ 9.00e-04 1.15e+00]
[ 1.90e-03 4.50e-01]
[ 9.00e-04 -5.00e-02]
[-1.00e-04 1.15e+00]
[-1.00e-04 2.50e-01]
[-1.00e-04 -6.50e-01]
[-1.10e-03 -5.50e-01]
[-1.00e-04 -3.50e-01]
[-1.10e-03 -8.50e-01]
[-1.10e-03 -5.50e-01]]
h [[ 2110635.52147809 825898.51117836 -91768.18464405 2110637.0365399
458834.25694829 -1192968.52264332 -1009433.36540466 -642367.59611278
-1560034.2919352 -1009433.36540466]
[ 825898.51117836 323176.76066518 -35909.18768022 825899.10402404
179543.32988011 -466812.44426383 -394994.54317994 -251360.51954919
-610446.46789458 -394994.54317994]
[ -91768.18464405 -35909.18768022 3989.98293674 -91768.25051909
-19949.6248335 51869.00085209 43889.08767865 27929.45895689
67828.62957385 43889.08767865]
[ 2110637.0365399 825899.10402404 -91768.25051909 2110638.55160468
458834.58631044 -1192969.37898381 -1009434.08999745 -642368.05721906
-1560035.4117622 -1009434.08999745]
[ 458834.25694829 179543.32988011 -19949.6248335 458834.58631044
99746.67497409 -259341.23636226 -219442.25018497 -139645.26591681
-339138.22063042 -219442.25018497]
[-1192968.52264332 -466812.44426383 51869.00085209 -1192969.37898381
-259341.23636226 674286.9062593 570549.58962751 363077.52538544
881758.97050137 570549.58962751]
[-1009433.36540466 -394994.54317994 43889.08767865 -1009434.08999745
-219442.25018497 570549.58962751 482771.99394445 307218.97635669
746102.60721527 482771.99394445]
[ -642367.59611278 -251360.51954919 27929.45895689 -642368.05721906
-139645.26591681 363077.52538544 307218.97635669 195503.26161803
474793.2401241 307218.97635669]
[-1560034.2919352 -610446.46789458 67828.62957385 -1560035.4117622
-339138.22063042 881758.97050137 746102.60721527 474793.2401241
1153068.33759254 746102.60721527]
[-1009433.36540466 -394994.54317994 43889.08767865 -1009434.08999745
-219442.25018497 570549.58962751 482771.99394445 307218.97635669
746102.60721527 482771.99394445]]
補足情報(FW/ツールのバージョンなど)
ここにより詳細な情報を記載してください。
回答2件
あなたの回答
tips
プレビュー
バッドをするには、ログインかつ
こちらの条件を満たす必要があります。