statistics - 多元回归中置信区间的 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]).

其中 bint 是 beta 的置信区间。

详情请引用 https://octave.sourceforge.io/statistics/function/regress.html .