STEP 1 Load package and load the data | STEP 2 Data Exploration | STEPS 3 Fitting Twofold HB Beta Model | STEP 4 Check Convergence via Plot MCMC | STEP 5 Extract Result for Area | Area Mean Estimation | Area Random Effect | Calculate Area Relative Standard Error (RSE) or CV and Mean Squared Error (MSE) | STEP 6 Extract Result for Subarea | Subarea Mean Estimation | Subarea Random Effect | Calculate Subarea Relative Standard Error (RSE) or CV and Mean Squared Error (MSE) | Extract Coefficient Estimation $\beta$ | Extract Area Random Effect Variance $\sigma_u^2$ and Subarea Random Effect Variance $\sigma_v^2$ | STEP 7 Visualize The Result