smoothArea function

Small area estimation via basic area level model