R中的函数“效果”有什么作用?


17

我不理解effects()R的帮助文件中的解释:

对于用lm或拟合的线性模型aov,其影响是通过在拟合过程中将数据投影到QR分解生成的连续正交子空间上而获得的不相关的单自由度值。

谁能解释这是什么意思?

一维子空间所涉及的正交子空间是否被QR分解的Q部分的列所覆盖(因此彼此正交)?还是应该与其他东西正交?

Answers:


3

给定响应向量,解释变量矩阵X及其QR分解X = Q R,返回的效果是向量Q T yÿXX=[RRŤÿ

这是数值示例,证实了上述内容:

> set.seed(1001)
> x<-rnorm(100)
> y<-1+2*x+rnorm(100)
> mod<-lm(y~x)
> xqr<-qr(cbind(1,x))
> sum(abs(qr.qty(xqr,y)-effects(mod)))
[1] 0
By using our site, you acknowledge that you have read and understand our Cookie Policy and Privacy Policy.
Licensed under cc by-sa 3.0 with attribution required.