我需要了解如何计算线性回归的置信区间。我自己获得的价值与使用Matlab的价值不同。所以,我一直试图理解它是如何在Matlab中完成的。在polyval.m中,我不明白下面代码的一部分。似乎E等于[E,R] = qr(V,0)。但不确定它究竟是什么。 我的问题是:
[E,R] = qr(V,0) different from [Q,R]=qr(V);
如何使用E(=V/R)
来计算置信区间(下面的delta)?
%S是一个包含三个元素的结构:原始X的Vandermonde矩阵的三角因子,自由度和残差的范数。
E = V / R;
e = sqrt(1+sum(E.*E,2));
...
delta = normr / sqrt(df)* e;
非常感谢!
大辅
答案 0 :(得分:0)
我在optimization tips and tricks doc。
中详细讨论了这个问题对于qr(V,0)与qr(V)的差异,这就是所谓的经济qr。阅读qr。
的帮助文档