> print(cbind(X,Y,Yhat,e))
      [,1] [,2] [,3] [,4]     [,5]        [,6]
 [1,]    1 43.7 19.0  177 177.3019  -0.3018776
 [2,]    1 43.7 43.0  279 271.8505   7.1495370
 [3,]    1 43.7 56.0  346 323.0643  22.9357198
 [4,]    1 54.6 13.0  160 172.1569 -12.1568874
 [5,]    1 54.6 19.0  193 195.7940  -2.7940338
 [6,]    1 54.6 43.0  280 290.3426 -10.3426192
 [7,]    1 54.6 56.0  335 341.5564  -6.5564364
 [8,]    1 55.7 13.0  169 174.0231  -5.0230683
 [9,]    1 55.7 26.0  212 225.2369 -13.2368854
[10,]    1 55.7 34.5  244 258.7228 -14.7228428
[11,]    1 55.7 43.0  285 292.2088  -7.2088001
[12,]    1 58.8 13.0  181 179.2823   1.7176946
[13,]    1 58.8 43.0  298 297.4680   0.5319628
[14,]    1 60.5 19.0  212 205.8035   6.1964505
[15,]    1 60.5 43.0  317 300.3521  16.6478650
[16,]    1 60.5 56.0  347 351.5660  -4.5659521
[17,]    1 61.9 13.0  186 184.5415   1.4584576
[18,]    1 61.9 19.0  216 208.1787   7.8213112
[19,]    1 61.9 34.5  265 269.2413  -4.2413169
[20,]    1 61.9 43.0  306 302.7273   3.2727257
[21,]    1 61.9 56.0  348 353.9411  -5.9410914
[22,]    1 66.7 13.0  209 192.6849  16.3151227
[23,]    1 66.7 43.0  324 310.8706  13.1293909
[24,]    1 66.7 56.0  352 362.0844 -10.0844263
> XXinv
             [,1]          [,2]          [,3]
[1,]  3.589339357 -5.801127e-02 -6.157377e-03
[2,] -0.058011273  1.001746e-03  1.084358e-05
[3,] -0.006157377  1.084358e-05  1.625484e-04
> betahat
          [,1]
[1,] 28.312636
[2,]  1.696528
[3,]  3.939524
> s2
         [,1]
[1,] 112.3413
> V_betahat
            [,1]         [,2]         [,3]
[1,] 403.2310528 -6.517061844 -0.691727735
[2,]  -6.5170618  0.112537458  0.001218182
[3,]  -0.6917277  0.001218182  0.018260902