In this article, a weighted empirical likelihood technique for constructing the empirical likelihood confidence regions is applied to study the heteroscedastic varying coefficient partially non-linear models with missing response data. We first provide an estimator of the error variance based on the Nadaraya-Watson kernel estimation method. Then a weighted empirical log-likelihood ratio of the unknown parameter is constructed based on the inverse probability weighted technique. The maximum empirical likelihood (MEL) estimator of the unknown parameter is obtained. Further, a weighted empirical log-likelihood ratio of the varying coefficient function is introduced based on the MEL estimator and the inverse probability weighted method. The limiting distributions of the resulting statistics for both the unknown parameter and varying coefficient function are shown to have the standard chi-squared distribution. A simulation study and a real data set example are undertaken to investigate the finite sample performance of the proposed methods.