Functions for some unit-level small area estimators and their