重回帰分析

重回帰分析

> x
[1] 1 2 3 4 5 6
> q1<-cbind(rep(1,6),x)
> q1
       x
[1,] 1 1
[2,] 1 2
[3,] 1 3
[4,] 1 4
[5,] 1 5
[6,] 1 6
> solve(t(q1)%*%q1)
                       x
   0.8666667 -0.20000000
x -0.2000000  0.05714286

>|sh|
[,1]
[1,] 0.0616738
>
> q2<-cbind(x-mean(x))
> 0.258^2*solve(t(q2)%*%q2)
[,1]
[1,] 0.003803657
|