多元迴歸中置信區間的Octave/Matlab實現

Question

我需要在Octave/Matlab中實現多元迴歸係數的置信區間。

該任務是以一種常見的方式定義的：數據Y，設計矩陣X，係數β，使得Y =βX。對於β的代碼是後來乾脆：

beta = pinv(X)*Y 

現在，作爲一個愚蠢的物理學家，我有點信心和預測區間丟失。公式以及它們的實現。

注意：我知道有一個Matlab函數mvregress，但它仍然從我正在使用的Octave中丟失。

注2：在CrossValidated上詢問此問題並將其標記爲關閉主題，因爲它側重於編程。

Answer 1

我想這是你想要查找：

[b, bint, r, rint, stats] = regress (y, X, [alpha]). 

其中賓特是測試版的置信區間。

詳情請參考https://octave.sourceforge.io/statistics/function/regress.html。