---

In this paper, we first introduce semi-parametric heteroscedastic hierarchical models. Then, we define a new version of the empirical likelihood function (Restricted Joint Empirical likelihood) and use it to obtain the shrinkage estimators of the models' parameters in these models. Under different assumptions, a simulation study investigates the better performance of the restricted joint empirical likelihood function in the analysis of semi-parametric heterogeneity hierarchical models. Furthermore, we analyze an actual data set using the RJEL method.