Statistics 多元回归中置信区间的Octave/Matlab实现

我需要在Octave/Matlab中实现多元回归系数的置信区间

任务以一种常见的方式定义：数据Y、设计矩阵X、系数β，因此Y=βX。β的代码简单如下：

beta = pinv(X)*Y

现在，作为一个愚蠢的物理学家，我对信心和预测区间有点迷茫。公式及其实现

注意：我知道有一个Matlab函数mvregresse，但我实际使用的倍频程中仍然缺少它

---

注2:这个问题是在交叉验证时提出的，并被标记为离题，因为它主要关注编程。

我想这就是您想要找到的：

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

其中bint是beta的置信区间

有关详情，请参阅