2012-11-13 45 views
0

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

  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)); 
    

    ...

    ,δ= normr/SQRT(DF)* E;

非常感谢!

大辅

回答

0

我在我的文档optimization tips and tricks一些长度讨论这个问题。

至于在QR(V,0)到QR(V)的区别,也就是所谓的经济QR。阅读qr的帮助文档。

+0

您的文档非常有帮助。谢谢! –