FIRST STEPS: Load package and load the data | STEPS 2: Fitting HB Model | STEP 3 Extract mean estimation | Subarea Estimation | Area Estimatio | Coefficient Estimation | Random effect variance estimation | STEP 3 : Extract CV and MSE Subarea | STEP 4 Extract CV and MSE of Area | STEP 5 : You can also compare the CV between subarea direct estimator and HB Twofold estimator