重回帰分析
> 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
|