library(foreign)
c = read.dta("C:/Users/BINH THANG/Dropbox/R - Learning/Tu/2018/NCKH3001T.dta")
c1=subset(c, select=c(tramcam, a2, a6, a11, a14,a16, a20, a25, a291, a292, a293, a33,a38, a41, a44, a45, ketquahoctap, bienco,hdtl, ruou, ngoi,chucap, hailongkqhoc, ppgiangday, apluchoc, dtdd, tghoatdong, klcthoc, hocnhom))
c1$a2=factor(c1$a2)
c1$a6=factor(c1$a6)
c1$a11=factor(c1$a11)
c1$a14=factor(c1$a14)
c1$a16=factor(c1$a16)
c1$a20=factor(c1$a20)
c1$a25=factor(c1$a25)
c1$a291=factor(c1$a291)
c1$a292=factor(c1$a292)
c1$a293=factor(c1$a293)
c1$a33=factor(c1$a33)
c1$a38=factor(c1$a38)
c1$a41=factor(c1$a41)
c1$a44=factor(c1$a44)
c1$ketquahoctap=factor(c1$ketquahoctap)
c1$bienco=factor(c1$bienco)
c1$hdtl=factor(c1$hdtl)
c1$ruou=factor(c1$ruou)
c1$ngoi=factor(c1$ngoi)
c1$hailongkqhoc=factor(c1$hailongkqhoc)
c1$ppgiangday=factor(c1$ppgiangday)
c1$apluchoc=factor(c1$apluchoc)
c1$dtdd=factor(c1$dtdd)
c1$tghoatdong=factor(c1$tghoatdong)
c1$klcthoc=factor(c1$klcthoc)
c1$hocnhom=factor(c1$hocnhom)
library(BMA)## Loading required package: survival## Loading required package: leaps## Loading required package: robustbase##
## Attaching package: 'robustbase'## The following object is masked from 'package:survival':
##
## heart## Loading required package: inline## Loading required package: rrcov## Scalable Robust Estimators with High Breakdown Point (version 1.4-3)yvar = c1[,1]
xvars = c1[,-1]
bma.search <- bic.glm(xvars, yvar, strict=F,
OR=20, glm.family="binomial")## Warning in bic.glm.data.frame(xvars, yvar, strict = F, OR = 20, glm.family
## = "binomial"): There were 76 records deleted due to NA'ssummary(bma.search)##
## Call:
## bic.glm.data.frame(x = xvars, y = yvar, glm.family = "binomial", strict = F, OR = 20)
##
##
## 13 models were selected
## Best 5 models (cumulative posterior probability = 0.8477 ):
##
## p!=0 EV SD model 1
## Intercept 100 -3.340579 0.85819 -3.0401
## a2.x 32.4
## .2 0.149486 0.24349 .
## a6.x 2.9
## .2 0.014445 0.09125 .
## .3 -0.019930 0.18466 .
## a11.x 0.0
## .2 0.000000 0.00000 .
## a14.x 0.0
## .2 0.000000 0.00000 .
## .3 0.000000 0.00000 .
## .4 0.000000 0.00000 .
## a16.x 0.0
## .2 0.000000 0.00000 .
## .3 0.000000 0.00000 .
## a20.x 100.0
## .2 2.011159 0.62484 2.0078
## a25.x 100.0
## .2 -0.564462 0.25342 -0.5740
## .3 0.620312 0.27725 0.6345
## .4 0.924669 0.26989 0.9200
## a291.x 0.0
## .thinh thoang 0.000000 0.00000 .
## .khong bao gio 0.000000 0.00000 .
## a292.x 0.0
## .thinh thoang 0.000000 0.00000 .
## .khong bao gio 0.000000 0.00000 .
## a293.x 0.0
## .thinh thoang 0.000000 0.00000 .
## .khong bao gio 0.000000 0.00000 .
## a33.x 0.0
## .2 0.000000 0.00000 .
## a38.x 100.0
## .2 1.332249 0.26278 1.3306
## a41.x 0.0
## .2 0.000000 0.00000 .
## a44.x 1.4
## .2 -0.007748 0.09978 .
## a45.x 3.7 0.038750 0.26352 .
## ketquahoctap.x 0.0
## .kha 0.000000 0.00000 .
## .trung binh 0.000000 0.00000 .
## bienco.x 65.3
## .khong -0.359843 0.30791 -0.5318
## hdtl.x 14.8
## .khong dat 0.055095 0.15298 .
## ruou.x 0.0
## .2 0.000000 0.00000 .
## .3 0.000000 0.00000 .
## .4 0.000000 0.00000 .
## ngoi.x 1.3
## .2 0.001992 0.02871 .
## chucap.x 0.0 0.000000 0.00000 .
## hailongkqhoc.x 0.0
## .binh thuong 0.000000 0.00000 .
## .khong hai long 0.000000 0.00000 .
## ppgiangday.x 0.0
## .binh thuong 0.000000 0.00000 .
## .khong hung thu 0.000000 0.00000 .
## apluchoc.x 0.0
## .binh thuong 0.000000 0.00000 .
## .nhe 0.000000 0.00000 .
## dtdd.x 0.0
## .2 0.000000 0.00000 .
## .3 0.000000 0.00000 .
## tghoatdong.x 0.0
## .thinh thoang 0.000000 0.00000 .
## .khong bao gio 0.000000 0.00000 .
## klcthoc.x 0.0
## .binh thuong 0.000000 0.00000 .
## .nhe 0.000000 0.00000 .
## hocnhom.x 0.0
## .thinh thoang 0.000000 0.00000 .
## .khong bao gio 0.000000 0.00000 .
##
## nVar 4
## BIC -3818.1450
## post prob 0.268
## model 2 model 3 model 4
## Intercept -3.3045 -3.3528 -3.6138
## a2.x
## .2 0.4791 . 0.4356
## a6.x
## .2 . . .
## .3 . . .
## a11.x
## .2 . . .
## a14.x
## .2 . . .
## .3 . . .
## .4 . . .
## a16.x
## .2 . . .
## .3 . . .
## a20.x
## .2 2.0631 1.9542 2.0000
## a25.x
## .2 -0.6280 -0.4923 -0.5370
## .3 0.5606 0.6774 0.6137
## .4 0.9290 0.9446 0.9566
## a291.x
## .thinh thoang . . .
## .khong bao gio . . .
## a292.x
## .thinh thoang . . .
## .khong bao gio . . .
## a293.x
## .thinh thoang . . .
## .khong bao gio . . .
## a33.x
## .2 . . .
## a38.x
## .2 1.3485 1.3208 1.3366
## a41.x
## .2 . . .
## a44.x
## .2 . . .
## a45.x . . .
## ketquahoctap.x
## .kha . . .
## .trung binh . . .
## bienco.x
## .khong -0.5709 . .
## hdtl.x
## .khong dat . . .
## ruou.x
## .2 . . .
## .3 . . .
## .4 . . .
## ngoi.x
## .2 . . .
## chucap.x . . .
## hailongkqhoc.x
## .binh thuong . . .
## .khong hai long . . .
## ppgiangday.x
## .binh thuong . . .
## .khong hung thu . . .
## apluchoc.x
## .binh thuong . . .
## .nhe . . .
## dtdd.x
## .2 . . .
## .3 . . .
## tghoatdong.x
## .thinh thoang . . .
## .khong bao gio . . .
## klcthoc.x
## .binh thuong . . .
## .nhe . . .
## hocnhom.x
## .thinh thoang . . .
## .khong bao gio . . .
##
## nVar 5 3 4
## BIC -3817.6634 -3817.4379 -3815.9984
## post prob 0.211 0.188 0.092
## model 5
## Intercept -3.2136
## a2.x
## .2 .
## a6.x
## .2 .
## .3 .
## a11.x
## .2 .
## a14.x
## .2 .
## .3 .
## .4 .
## a16.x
## .2 .
## .3 .
## a20.x
## .2 2.0279
## a25.x
## .2 -0.5914
## .3 0.5869
## .4 0.8566
## a291.x
## .thinh thoang .
## .khong bao gio .
## a292.x
## .thinh thoang .
## .khong bao gio .
## a293.x
## .thinh thoang .
## .khong bao gio .
## a33.x
## .2 .
## a38.x
## .2 1.3073
## a41.x
## .2 .
## a44.x
## .2 .
## a45.x .
## ketquahoctap.x
## .kha .
## .trung binh .
## bienco.x
## .khong -0.5742
## hdtl.x
## .khong dat 0.4040
## ruou.x
## .2 .
## .3 .
## .4 .
## ngoi.x
## .2 .
## chucap.x .
## hailongkqhoc.x
## .binh thuong .
## .khong hai long .
## ppgiangday.x
## .binh thuong .
## .khong hung thu .
## apluchoc.x
## .binh thuong .
## .nhe .
## dtdd.x
## .2 .
## .3 .
## tghoatdong.x
## .thinh thoang .
## .khong bao gio .
## klcthoc.x
## .binh thuong .
## .nhe .
## hocnhom.x
## .thinh thoang .
## .khong bao gio .
##
## nVar 5
## BIC -3815.9407
## post prob 0.089