matlab - matlab polyval中的置信度

Question

我需要了解如何计算线性回归的置信区间。我自己得到的值与使用 Matlab 得到的值不同。所以，我一直试图了解它是如何在 Matlab 中完成的。在 polyval.m 中，我不明白下面的代码部分。似乎 E 等价于 [E,R] = qr(V,0)。但不确定它到底是什么。我的问题是：

1. 怎么[E,R] = qr(V,0) different from [Q,R]=qr(V);

2. 如何E(=V/R)用于计算置信区间（下面的增量）？

   % S 是一个包含三个元素的结构：原始 X 的 Vandermonde 矩阵的三角因子、自由度和残差范数。

   E = V/R；

   e = sqrt(1+sum(E.*E,2));
   

   ...

   delta = normr/sqrt(df)*e;

非常感谢！

大辅

Answer 1

我在我的优化提示和技巧文档中详细讨论了这一点。

至于qr(V,0)与qr(V)的差异，就是所谓的经济qr。阅读 qr 的帮助文档。