-------------
      name:  <unnamed>
       log:  C:\Users\lenovo\Desktop\jay thesis proposal\Quality of Family Planning care as enabling factor of modern contraceptive use in Indonesia.smcl
  log type:  smcl
 opened on:   9 May 2021, 00:39:44

. do "C:\Users\lenovo\Desktop\jay thesis proposal\Quality of Family Planning Care as enabling factor of Modern Contraceptive Use in Indonesia.do"

. use IDR2_HH_NearestFPDistance_16Apr2020_Dataset, clear

.
. *OUTCOME VARIABLE
.
. gen currently_use_modern_method = current_user
(942 missing values generated)

. label var currently_use_modern_method "Current use of modern contraceptive"

. recode currently_use_modern_method (-99=.)
(currently_use_modern_method: 44 changes made)

. recode currently_use_modern_method (1 = 0) if current_methodnum==13 | current_methodnum==14 | current_methodnum==30| current_methodnum==31| current_methodnum==39
(currently_use_modern_method: 178 changes made)

. recode currently_use_modern_method (1 = .) if current_methodnum==-99
(currently_use_modern_method: 16 changes made)

. label define modern_method_label 0 "Non-user or Traditional method user" 1 "Modern contraceptive user"

. label values currently_use_modern_method modern_method_label

.
. gen modern = current_methodnum
(6,696 missing values generated)

. recode modern (1=1) (2=2) (3=3) (4=4) (5 6 = 5) (7 8 = 6) (9 10 = 7) (else=.)
(modern: 1534 changes made)

. label define modern_label 1 "1. female_sterilization" 2 "2. male_sterilization" ///
> 3 "3. implants" 4 "4. IUD" 5 "5. injectables" 6 "6. pills" 7 "7. condoms"

. label values modern modern_label

.
. *PUT QUALITY OF FP CARE VARIABLES INTO MAIN DATASET
.
. drop if currently_use_modern_method==.
(1,002 observations deleted)

. merge m:1 EA_ID facility_ID using SDP_qoc_data

    Result                           # of obs.
    -----------------------------------------
    not matched                         1,143
        from master                       598  (_merge==1)
        from using                        545  (_merge==2)

    matched                             9,612  (_merge==3)
    -----------------------------------------

.
. *PUT COMMUNITY-LEVEL VARIABLES INTO MAIN DATASET
.
. drop if currently_use_modern_method==.
(545 observations deleted)

. drop _merge

. merge m:1 EA_ID using community_level_data

    Result                           # of obs.
    -----------------------------------------
    not matched                             0
    matched                            10,210  (_merge==3)
    -----------------------------------------

.
. *RECODING OTHER COVARIATES/CONFOUNDERS
.
. gen age_group = FQ_age

. label var age_group "Age Group"

. recode age_group (15/19 = 0) (20/24 = 1)(25/29 = 2)(30/34 = 3) (35/39 = 4) (40/44 = 5) (45/49 = 6)
(age_group: 10210 changes made)

. label define age_label 0 "15-19" 1 "20-24" 2 "25-29" 3 "30-34" 4 "35-39" 5 "40-44" 6 "45-49"

. label values age_group age_label

.
. gen marital = FQmarital_status

. label var marital "Marital Status"

. recode marital (5 = 0) (1/2 = 1)(3/4 = 2)
(marital: 2937 changes made)

. label define marital_label 0 "0. Never Married" 1 "1. Married/Living with partner" 2 "2. Divorced/Separated/Widow"

. label values marital marital_label

.
. gen education = school

. label var education "Education"

. recode education (0/1 = 0)(2/3 = 1) (4/5 = 2)
(education: 10035 changes made)

. label define education_label 0 "0. None/Primary School" 1 "1. High School" 2 "2. Academy/University"

. label values education education_label

.
. gen parity = children_born

. label var parity "Parity"

. recode parity (-99 = .) (0 = 0) (1 = 1)(2 = 2) (3 = 3) (4/11 = 4)
(parity: 520 changes made)

. label define parity_label 0 "No children" 1 "1 child" 2 "2 children" 3 "3 children" 4 "more than 3 children"

. label values parity parity_label

.
. gen womans_FPdecision_autonomy1 = fp_final_decision
(5,437 missing values generated)

. recode womans_FPdecision_autonomy1 (-99 = -99) (2 3 5 6 = 0) (1 4 = 1)
(womans_FPdecision_autonomy1: 1999 changes made)

. gen womans_FPdecision_autonomy2 = why_not_decision
(4,955 missing values generated)

. recode womans_FPdecision_autonomy2 (-99 = -99) (2 3 96 = 0) (1 = 1)
(womans_FPdecision_autonomy2: 2070 changes made)

. recode womans_FPdecision_autonomy1 (mis = 0)
(womans_FPdecision_autonomy1: 5437 changes made)

. recode womans_FPdecision_autonomy2 (mis = 0)
(womans_FPdecision_autonomy2: 4955 changes made)

. gen womans_FPdecision_autonomy = womans_FPdecision_autonomy1 + womans_FPdecision_autonomy2

. label var womans_FPdecision_autonomy "Woman's Autonomy in Family Planning Decision"

. recode womans_FPdecision_autonomy (0 = 0) (1 2 = 1) (else = .)
(womans_FPdecision_autonomy: 303 changes made)

. label define womans_FPdecision_autonomy_label 0 "0. Low" 1 "1. High"

. label values womans_FPdecision_autonomy womans_FPdecision_autonomy_label

.
. gen region = province

. label var region "Region"

. recode region (11 12 13 14 21 15 16 17 18 19 = 0) (31 32 33 34 35 36 51 = 1) ///
> (52 53 = 2) (61 62 63 64 = 3) (71 72 75 73 74 76 = 4) (81 82 91 94 = 5)
(region: 10210 changes made)

. label define region_label 0 "0. Sumatera" 1 "1. Java and Bali" ///
> 2 "2. Nusa Tenggara" 3 "3. Kalimantan" 4 "4. Sulawesi" 5 "5. Maluku and Papua"

. label values region region_label

.
. gen sdp_type = facility_type

. label var sdp_type "SDP type"

. recode sdp_type (1 2 = 0) (3 4 5 = 1) (6 7 = 2) (8 10 = 3) (11 96 = 4)
(sdp_type: 10210 changes made)

. label define sdp_type_label 0 "0. Hospital" 1 "1. Health Center/clinic" 2 "2. Village health post/Delivery post" 3 "3. Private doctor/midwife" 4 "4. Others"

. label values sdp_type sdp_type_label

.
.
. *DESCRIPTIVE ANALYSIS
.
. *OUTCOME VARIABLE
. tab currently_use_modern_method

Current use of modern contraceptive |      Freq.     Percent        Cum.
------------------------------------+-----------------------------------
Non-user or Traditional method user |      5,888       57.67       57.67
          Modern contraceptive user |      4,322       42.33      100.00
------------------------------------+-----------------------------------
                              Total |     10,210      100.00

. tab modern

                 modern |      Freq.     Percent        Cum.
------------------------+-----------------------------------
1. female_sterilization |        266        6.15        6.15
  2. male_sterilization |          5        0.12        6.27
            3. implants |        403        9.32       15.59
                 4. IUD |        336        7.77       23.37
         5. injectables |      2,305       53.33       76.70
               6. pills |        907       20.99       97.69
             7. condoms |        100        2.31      100.00
------------------------+-----------------------------------
                  Total |      4,322      100.00

.
. *INDIVIDUAL LEVEL VARIABLES
. tab age_group

  Age Group |      Freq.     Percent        Cum.
------------+-----------------------------------
      15-19 |      1,477       14.47       14.47
      20-24 |      1,303       12.76       27.23
      25-29 |      1,292       12.65       39.88
      30-34 |      1,532       15.00       54.89
      35-39 |      1,679       16.44       71.33
      40-44 |      1,586       15.53       86.87
      45-49 |      1,341       13.13      100.00
------------+-----------------------------------
      Total |     10,210      100.00

. tab education

             Education |      Freq.     Percent        Cum.
-----------------------+-----------------------------------
0. None/Primary School |      2,868       28.09       28.09
        1. High School |      5,852       57.32       85.41
 2. Academy/University |      1,490       14.59      100.00
-----------------------+-----------------------------------
                 Total |     10,210      100.00

. tab parity

              Parity |      Freq.     Percent        Cum.
---------------------+-----------------------------------
         No children |      2,896       28.39       28.39
             1 child |      1,823       17.87       46.27
          2 children |      2,728       26.75       73.02
          3 children |      1,581       15.50       88.52
more than 3 children |      1,171       11.48      100.00
---------------------+-----------------------------------
               Total |     10,199      100.00

. tab marital

                Marital Status |      Freq.     Percent        Cum.
-------------------------------+-----------------------------------
              0. Never Married |      2,433       23.83       23.83
1. Married/Living with partner |      7,338       71.87       95.70
   2. Divorced/Separated/Widow |        439        4.30      100.00
-------------------------------+-----------------------------------
                         Total |     10,210      100.00

. tab womans_FPdecision_autonomy

    Woman's |
Autonomy in |
     Family |
   Planning |
   Decision |      Freq.     Percent        Cum.
------------+-----------------------------------
     0. Low |      4,061       40.99       40.99
    1. High |      5,846       59.01      100.00
------------+-----------------------------------
      Total |      9,907      100.00

.
. *HOUSEHOLD LEVEL VARIABLES
. tab wealth

Household's |
     Wealth |      Freq.     Percent        Cum.
------------+-----------------------------------
 0. Poorest |      1,675       16.41       16.41
    1. Poor |      1,791       17.54       33.95
  2. Middle |      1,873       18.34       52.29
    3. Rich |      2,591       25.38       77.67
 4. Richest |      2,280       22.33      100.00
------------+-----------------------------------
      Total |     10,210      100.00

. tab insurance

                       Health Insurance |      Freq.     Percent        Cum.
----------------------------------------+-----------------------------------
                 0. No health insurance |      4,004       39.22       39.22
1. Free National/District Health Insura |      2,165       21.20       60.42
  2. Non-Free National Health Insurance |      1,975       19.34       79.76
              3. Other health insurance |      2,066       20.24      100.00
----------------------------------------+-----------------------------------
                                  Total |     10,210      100.00

.
. *COMMUNITY LEVEL VARIABLES
. tab residence

  Residence |      Freq.     Percent        Cum.
------------+-----------------------------------
   0. Rural |      4,572       44.78       44.78
   1. Urban |      5,638       55.22      100.00
------------+-----------------------------------
      Total |     10,210      100.00

. tab region

             Region |      Freq.     Percent        Cum.
--------------------+-----------------------------------
        0. Sumatera |      1,870       18.32       18.32
   1. Java and Bali |      4,974       48.72       67.03
   2. Nusa Tenggara |        386        3.78       70.81
      3. Kalimantan |        489        4.79       75.60
        4. Sulawesi |      2,271       22.24       97.85
5. Maluku and Papua |        220        2.15      100.00
--------------------+-----------------------------------
              Total |     10,210      100.00

. sum mean_ideal_number_children

    Variable |        Obs        Mean    Std. Dev.       Min        Max
-------------+---------------------------------------------------------
mean_ideal~n |     10,133    2.662062    .4799439   1.772727   5.333333

. sum percent_dissatisfied_with_FP

    Variable |        Obs        Mean    Std. Dev.       Min        Max
-------------+---------------------------------------------------------
percent_di~P |     10,210    15.04939    12.52461          0     111.11

. sum percent_visited_by_CHW

    Variable |        Obs        Mean    Std. Dev.       Min        Max
-------------+---------------------------------------------------------
percent_vi~W |     10,210    9.488956    12.59335          0      96.43

.
. *SDP LEVEL VARIABLES
. tab sdp_type

                            SDP type |      Freq.     Percent        Cum.
-------------------------------------+-----------------------------------
                         0. Hospital |        249        2.44        2.44
             1. Health Center/clinic |      5,924       58.02       60.46
2. Village health post/Delivery post |      1,923       18.83       79.29
           3. Private doctor/midwife |      1,539       15.07       94.37
                           4. Others |        575        5.63      100.00
-------------------------------------+-----------------------------------
                               Total |     10,210      100.00

. tab overall_quality

    Overall |
 Quality of |
     Family |
   Planning |
       Care |      Freq.     Percent        Cum.
------------+-----------------------------------
     0. Low |      2,733       28.43       28.43
  1. Medium |      5,636       58.64       87.07
    2. High |      1,243       12.93      100.00
------------+-----------------------------------
      Total |      9,612      100.00

.
.
. *MULTILEVEL NULL MODEL
.
. qui logit currently_use_modern_method i.age_group i.education i.parity i.marital i.wealth i.residence i.region ///
> i.womans_FPdecision_autonomy i.insurance ///
> mean_ideal_number_children percent_dissatisfied_with_FP /// *community social norms
> percent_visited_by_CHW /// *community demand generation
> i.sdp_type i.overall_quality, nolog or

. est store fixedeffect

.
. melogit currently_use_modern_method || EA_ID: , nolog or

Mixed-effects logistic regression               Number of obs     =     10,210
Group variable:           EA_ID                 Number of groups  =        372

                                                Obs per group:
                                                              min =          4
                                                              avg =       27.4
                                                              max =         50

Integration method:  mvaghermite                Integration pts.  =          7

                                                Wald chi2(0)      =          .
Log likelihood = -6703.5872                     Prob > chi2       =          .
----------------------------+----------------------------------------------------------------
 currently_use_modern_method | Odds Ratio    Std. Err.       z    P>|z|      [95% Con f. Interval]
----------------------------+----------------------------------------------------------------
                      _cons |    .7297399    .0308998     -7.44    0.000      .6716224     .7928865
----------------------------+----------------------------------------------------------------
 EA_ID                       |
                  var(_cons)|    .4895969     .053002                    .3959959     .6053224
----------------------------+----------------------------------------------------------------
LR test vs. logistic model: chibar2(01) = 505.75      Prob >= chibar2 =  0.0000

. estat ic

Akaike's information criterion and Bayesian information criterion

-------------+---------------------------------------------------------------
       Model |        Obs  ll(null)  ll(model)      df         AIC        BIC
-------------+---------------------------------------------------------------
.            |     10,210         .  -6703.587       2    13411.17   13425.64
-------------+---------------------------------------------------------------


Note: N=Obs used in calculating BIC; see [R] BIC note.


. est store nullmodel

. lrtest nullmodel fixedeffect, force

Likelihood-ratio test                                 LR chi2(36) =   5719.56
(Assumption: nullmodel nested in fixedeffect)         Prob > chi2 =    0.0000

. est tab fixedeffect nullmodel

-------------+--------------------------
Variable   | fixedeff~t     nullmodel
-------------+--------------------------
currently_~d |
   age_group |
          1  | -.08721707
          2  | -.78981231
          3  | -1.1789268
          4  | -1.3124408
          5  | -1.6383167
          6  | -2.427054
             |
   education |
          1  | .08644427
          2  | -.07884459
             |
      parity |
          1  | 3.3635968
          2  | 4.3129757
          3  | 4.6493781
          4  | 4.4963579
             |
     marital |
          1  | 3.8536024
          2  | .10320168
             |
      wealth |
          1  | .05548454
          2  | -.04606226
          3  | -.01106332
          4  | -.01457891
             |
   residence |
          1  | -.31030223
             |
      region |
          1  | .21751668
          2  | -.4555481
          3  | .33300759
          4  | -.35951313
          5  | .04777623
             |
womans_FPd~y |
          1  | .94920099
             |
   insurance |
          1  | .11890512
          2  | .09400938
          3  | .09818314
             |
mean_ideal~n | -.56219593
percent_di~P | -.03466091
percent_vi~W | .00611647
             |
    sdp_type |
          1  | .04513865
          2  | .22450067
          3  | .15012871
          4  | .21790404
             |
overall_qu~y |
          1  | .00234549
          2  | .31118687
             |
       _cons | -4.8098862     -.31506707
-------------+--------------------------
var(_cons[~) |
       _cons |      .48959694
-------------+--------------------------

.
. *MULTILEVEL MIXED EFFECT MODEL OF MODERN CONTRACEPTIVE USE
.
. collin age_group education parity marital wealth residence region ///
> womans_FPdecision_autonomy insurance ///
> mean_ideal_number_children percent_dissatisfied_with_FP /// *community social norms
> percent_visited_by_CHW /// *community demand generation
> sdp_type overall_quality
(obs=9,246)

  Collinearity Diagnostics

                        SQRT                   R-
  Variable      VIF     VIF    Tolerance    Squared
----------------------------------------------------
 age_group      2.13    1.46    0.4700      0.5300
 education      1.33    1.15    0.7520      0.2480
    parity      2.20    1.48    0.4547      0.5453
   marital      1.79    1.34    0.5593      0.4407
    wealth      1.39    1.18    0.7212      0.2788
 residence      1.31    1.15    0.7624      0.2376
    region      1.08    1.04    0.9238      0.0762
womans_FPdecision_autonomy      1.03    1.02    0.9672      0.0328
 insurance      1.03    1.01    0.9715      0.0285
mean_ideal_number_children      1.12    1.06    0.8934      0.1066
percent_dissatisfied_with_FP      1.11    1.06    0.8984      0.1016
percent_visited_by_CHW      1.05    1.02    0.9533      0.0467
  sdp_type      1.13    1.06    0.8865      0.1135
overall_quality      1.13    1.06    0.8878      0.1122
----------------------------------------------------
  Mean VIF      1.34

                           Cond
        Eigenval          Index
---------------------------------
    1    10.3245          1.0000
    2     0.8497          3.4858
    3     0.6624          3.9481
    4     0.5007          4.5411
    5     0.4776          4.6494
    6     0.4407          4.8402
    7     0.4041          5.0544
    8     0.3322          5.5751
    9     0.3098          5.7732
    10     0.2244          6.7826
    11     0.1659          7.8900
    12     0.1195          9.2940
    13     0.0979         10.2670
    14     0.0786         11.4583
    15     0.0120         29.2898
---------------------------------
 Condition Number        29.2898
 Eigenvalues & Cond Index computed from scaled raw sscp (w/ intercept)
 Det(correlation matrix)    0.1238

.
. melogit currently_use_modern_method i.age_group i.education i.parity i.marital i.wealth i.residence i.region ///
> i.womans_FPdecision_autonomy i.insurance ///
> mean_ideal_number_children percent_dissatisfied_with_FP /// *community social norms
> percent_visited_by_CHW /// *community demand generation
> i.sdp_type i.overall_quality || EA_ID: womans_FPdecision_autonomy , nolog or

Mixed-effects logistic regression               Number of obs     =      9,246
Group variable:           EA_ID                 Number of groups  =        350

                                                Obs per group:
                                                              min =          1
                                                              avg =       26.4
                                                              max =         50

Integration method:  mvaghermite                Integration pts.  =          7

                                                Wald chi2(37)     =    1277.24
Log likelihood = -3822.5568                     Prob > chi2       =     0.0000
---------------------------------------------+----------------------------------------------------------------
                  currently_use_modern_method | Odds Ratio    Std. Err.       z    P>|z|      [95% Con f. Interval]
---------------------------------------------+----------------------------------------------------------------
                                             |
                                   age_group |
                                      20-24  |    .9545924     .264652     -0.17    0.867      .5544087     1.643637
                                      25-29  |     .464037    .1268434     -2.81    0.005       .271568     .7929149
                                      30-34  |    .3150639    .0866476     -4.20    0.000      .1837827     .5401231
                                      35-39  |    .2727417    .0753895     -4.70    0.000      .1586601     .4688514
                                      40-44  |    .1949543    .0542241     -5.88    0.000      .1130268     .3362668
                                      45-49  |    .0860547    .0242501     -8.70    0.000      .0495344     .1495001
                                             |
                                   education |
                             1. High School  |    1.130217     .080603      1.72    0.086      .9827825      1.29977
                      2. Academy/University  |    .9519344    .1065822     -0.44    0.660      .7643693     1.185525
                                             |
                                      parity |
                                    1 child  |     30.0488    7.956834     12.85    0.000      17.88264     50.49203
                                 2 children  |    80.87231    21.92482     16.20    0.000      47.53741     137.5828
                                 3 children  |    115.6331    32.05492     17.14    0.000      67.16126     199.0883
                       more than 3 children  |    99.54237    28.18637     16.25    0.000      57.14552     173.3939
                                             |
                                     marital |
             1. Married/Living with partner  |    55.07519    26.99821      8.18    0.000      21.07131     143.9529
                2. Divorced/Separated/Widow  |    1.168606    .6388872      0.28    0.776      .4002291     3.412145
                                             |
                                      wealth |
                                    1. Poor  |     1.08527    .1121344      0.79    0.428      .8863151     1.328885
                                  2. Middle  |    .9938977    .1057032     -0.06    0.954      .8068904     1.224247
                                    3. Rich  |    1.034326    .1114763      0.31    0.754      .8373709     1.277607
                                 4. Richest  |    1.033583    .1217344      0.28    0.779      .8205243     1.301964
                                             |
                                   residence |
                                   1. Urban  |     .707886    .0598716     -4.08    0.000      .5997501     .8355189
                                             |
                                      region |
                           1. Java and Bali  |    1.300387    .1354195      2.52    0.012      1.060304     1.594832
                           2. Nusa Tenggara  |       .6738    .1363051     -1.95    0.051      .4532511     1.001667
                              3. Kalimantan  |    1.391368    .2613431      1.76    0.079      .9628515     2.010595
                                4. Sulawesi  |    .6782971    .0855121     -3.08    0.002      .5297976     .8684203
                        5. Maluku and Papua  |    1.118248     .317063      0.39    0.693      .6414935     1.949324
                                             |
                  womans_FPdecision_autonomy |
                                    1. High  |    2.862215     .202992     14.83    0.000      2.490772     3.289051
                                             |
                                   insurance |
 1. Free National/District Health Insurance  |     1.14685    .0931595      1.69    0.092      .9780538     1.344777
      2. Non-Free National Health Insurance  |    1.081947    .0944766      0.90    0.367      .9117554     1.283907
                  3. Other health insurance  |    1.111559    .0918218      1.28    0.200      .9454047     1.306914
                                             |
                  mean_ideal_number_children |    .5745655    .0484788     -6.57    0.000      .4869897     .6778902
                percent_dissatisfied_with_FP |     .964244    .0032778    -10.71    0.000       .957841     .9706899
                      percent_visited_by_CHW |    1.006922    .0030622      2.27    0.023      1.000938     1.012941
                                             |
                                    sdp_type |
                    1. Health Center/clinic  |      1.0522    .2697614      0.20    0.843      .6366026     1.739114
       2. Village health post/Delivery post  |    1.300415    .3454849      0.99    0.323      .7725755     2.188886
                  3. Private doctor/midwife  |    1.220911    .3265393      0.75    0.455       .722812     2.062255
                                  4. Others  |    1.400215    .4184036      1.13    0.260      .7795497     2.515043
                                             |
                             overall_quality |
                                  1. Medium  |    1.024491    .0883063      0.28    0.779       .865244     1.213047
                                    2. High  |    1.422621    .1802727      2.78    0.005      1.109752     1.823696
                                             |
                                       _cons |    .0055349    .0031592     -9.10    0.000      .0018083     .0169417
---------------------------------------------+----------------------------------------------------------------
 EA_ID                                        |
              var(womans_FPdecision_autonomy)|    .2651738    .0789757                    .1479178     .4753798
                                   var(_cons)|      .09696    .0339624                    .0488023     .1926394
---------------------------------------------+----------------------------------------------------------------
LR test vs. logistic model: chi2(2) = 42.50               Prob > chi2 = 0.0000


Note: LR test is conservative and provided only for reference.


. estat ic

Akaike's information criterion and Bayesian information criterion

-------------+---------------------------------------------------------------
       Model |        Obs  ll(null)  ll(model)      df         AIC        BIC
-------------+---------------------------------------------------------------
.            |      9,246         .  -3822.557      40    7725.114   8010.391
-------------+---------------------------------------------------------------


Note: N=Obs used in calculating BIC; see [R] BIC note.


.
. est store mixedeffect

. lrtest mixedeffect fixedeffect, force

Likelihood-ratio test                                 LR chi2(2)  =     42.50
(Assumption: fixedeffect nested in mixedeffect)       Prob > chi2 =    0.0000

. lrtest mixedeffect nullmodel, force

Likelihood-ratio test                                 LR chi2(38) =   5762.06
(Assumption: nullmodel nested in mixedeffect)         Prob > chi2 =    0.0000

. est tab fixedeffect nullmodel mixedeffect

-------------+---------------------------------------
Variable   | fixedeff~t     nullmodel     mixedeff~t
-------------+---------------------------------------
currently_~d |
   age_group |
      20-24  | -.08721707          -.04647088
      25-29  | -.78981231          -.767791
      30-34  | -1.1789268          -1.1549797
      35-39  | -1.3124408          -1.2992301
      40-44  | -1.6383167          -1.6349903
      45-49  | -2.427054          -2.4527726
             |
   education |
1. High S..  | .08644427          .12241007
2. Academ..  | -.07884459          -.04925918
             |
      parity |
    1 child  | 3.3635968          3.4028228
 2 children  | 4.3129757          4.3928715
 3 children  | 4.6493781          4.7504225
more than..  | 4.4963579          4.6005834
             |
     marital |
1. Marrie..  | 3.8536024          4.0086993
2. Divorc..  | .10320168          .15581146
             |
      wealth |
    1. Poor  | .05548454          .08182859
  2. Middle  | -.04606226          -.00612096
    3. Rich  | -.01106332          .03375033
 4. Richest  | -.01457891          .033031
             |
   residence |
   1. Urban  | -.31030223          -.34547225
             |
      region |
1. Java a..  | .21751668          .26266204
2. Nusa T..  | -.4555481          -.39482194
3. Kalima..  | .33300759          .33028738
4. Sulawesi  | -.35951313          -.38816983
5. Maluku..  | .04777623          .11176311
             |
womans_FPd~y |
    1. High  | .94920099          1.0515959
             |
   insurance |
1. Free N..  | .11890512          .1370187
2. Non-Fr..  | .09400938          .07876224
3. Other ..  | .09818314          .10576323
             |
mean_ideal~n | -.56219593          -.55414114
percent_di~P | -.03466091          -.03641086
percent_vi~W | .00611647          .00689771
             |
    sdp_type |
1. Health..  | .04513865          .05088302
2. Villag..  | .22450067          .2626835
3. Privat..  | .15012871          .19959704
  4. Others  | .21790404          .33662559
             |
overall_qu~y |
  1. Medium  | .00234549          .024196
    2. High  | .31118687          .35250098
             |
       _cons | -4.8098862     -.31506707     -5.1966846
-------------+---------------------------------------
var(_cons[~) |
       _cons |      .48959694     .09696005
-------------+---------------------------------------
var(         |
       _cons |           .2651738
-------------+---------------------------------------

.
. predict random_intercept random_slope, reffects
(calculating posterior means of random effects)
(using 7 quadrature points)

. predict yhat, fitted
(predictions based on fixed effects and posterior means of random effects)
(using 7 quadrature points)
(964 missing values generated)

. tabstat random_intercept random_slope yhat if EA_ID<=8060, by(EA_ID)

Summary statistics: mean
  by categories of: EA_ID (EA number (random))

EA_ID   |  random~t  random~e      yhat
---------+------------------------------
8004   |  .0988832 -.1001581 -1.153787
8005   |  -.221927   .193685  .1360555
8010   |  .1746993  .0563468 -1.174551
8013   |  .1111906 -.0356388  -2.36823
8020   |  .0556005 -.1672615 -1.394754
8022   |  .0447763  .0849933 -1.601934
8026   |  .2886707   .139508 -.4382866
8029   |  .0687178  .0704818 -2.735673
8030   | -.1345847 -.0135205 -1.213742
8032   |   .329122 -.2348595 -1.190232
8036   |  .8365927  .1238277 -.6093433
8038   | -.3235621 -.1444452 -3.780792
8043   | -.1422481  .2877121 -2.634701
8044   |  .4785458 -.1025905 -.4938952
8047   |  .0817869  .0299051 -1.002193
8048   |  8.91e-18  .0496482 -.0493169
8053   | -.0603664 -.0815166 -2.335384
8060   | -.1111453  .1692882  .4291258
---------+------------------------------
Total   |  .1066271  .0201799 -1.441358
---------+------------------------------

.
. graph hbar (mean) random_intercept if EA_ID<=8060, over(EA_ID) ytitle("Random effect for intercepts of modern contraceptive use model") name(graph1, replace)

. graph save graph1 Random_Intercepts.gph, replace
(file Random_Intercepts.gph saved)

.
. graph hbar (mean) random_slope if EA_ID<=8060, over(EA_ID) ytitle("Random effect for slopes of modern contraceptive use model") name(graph2, replace)

. graph save graph2 Random_Slopes.gph, replace
(file Random_Slopes.gph saved)

.
.
end of do-file

. log close
      name:  <unnamed>
       log:  C:\Users\lenovo\Desktop\jay thesis proposal\Quality of Family Planning care as enabling factor of modern contraceptive use in Indonesia.smcl
  log type:  smcl
 closed on:   9 May 2021, 00:40:53
-------------