#To get information (variance-covariance) matrix of beta1 and beta2 where
#beta1 and 2 are regression coefficients from a multiple regression model.
dat=dat1
omat=cor(dat)[1:3,1:3]
#omat
#1.0000000 0.1958636 0.1970060
#0.1958636 1.0000000 0.9981003
#0.1970060 0.9981003 1.0000000
nv=length(dat$V1)
output=olkin_beta_inf(omat,nv)
output
#output$info (2x2 information (variance-covariance) matrix)
#0.2531406 -0.2526212
#-0.2526212 0.2530269
#output$var1 (variance of beta1)
#0.2531406
#output$var2 (variance of beta2)
#0.2530269
#output$var1_2 (variance of difference between beta1 and beta2)
#1.01141
Run the code above in your browser using DataLab