Non-orthogonal multiple access (NOMA) has been envisioned as a promising multiple access technique for 5G and beyond wireless networks due to its significant enhancement of spectral efficiency. In this paper, we investigate a robust energy efficiency design for multiuser multiple-input singleoutput (MISO) NOMA systems where imperfect channel state information is available at the base station. A clustering algorithm is applied to group the users into different clusters, and then NOMA technique is employed to share the available resources fairly among the users in each cluster. To remove the interference between clusters, two different types of zero-forcing (ZF) designs, namely, hybrid-ZF and full-ZF are employed at the BS. The full-ZF scheme completely removes the interference leakage at the cost of more number of antennas and the hybrid-ZF scheme partially mitigates the interference leakage. To solve the problem, the Dinkelbach's algorithm is employed to convert the non-linear fractional programming problem into a simple subtractive form. Finally, simulation results reveal that hybrid-ZF outperforms the full-ZF scheme with a few clusters, while full-ZF shows a better performance with the higher number of clusters. The numerical results confirm that our proposed robust scheme outperforms the non-robust scheme in terms of the rate-satisfaction ratio at each user.