> print(round(cbind(betahat,s_betahat,t_betahat,p_betahat),4))
             s_betahat               
[1,] 28.3126   20.0806  1.4099 0.1732
[2,]  1.6965    0.3355  5.0572 0.0001
[3,]  3.9395    0.1351 29.1530 0.0000
> 
>
> 
> print(round(c(R_1_0,R_2_01),3))
[1]  2053.462 95478.329
> print(round(c(R_2_0,R_1_02),3))
[1] 94658.600  2873.191