信心regino in matlab polyval

时间:2012-11-13 13:42:54

标签: matlab linear-regression

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

  1. [E,R] = qr(V,0) different from [Q,R]=qr(V);
  2. 怎么样?
  3. 如何使用E(=V/R)来计算置信区间(下面的delta)?

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

    E = V / R;

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

    ...

    delta = normr / sqrt(df)* e;

  4. 非常感谢!

    大辅

1 个答案:

答案 0 :(得分:0)

我在optimization tips and tricks doc。

中详细讨论了这个问题

对于qr(V,0)与qr(V)的差异,这就是所谓的经济qr。阅读qr。

的帮助文档