© 2019 Author(s). In this work, analysis and simulation of all experimentally observed switching modes in hafnium oxide based resistive random access memories are carried out using a simplified electrical conduction model. To achieve switching mode variation, two metal-insulator-metal cells with identical stack combination, but varying oxygen stoichiometry of the hafnia layer, namely, stoichiometric vs highly deficient, are considered. To access the individual switching modes, the devices were subjected to a variety of cycling conditions comprising different voltage and current ranges. For modeling the device behavior, a single or two antiserially connected memdiodes (diode with memory) were utilized. In this way, successful compact simulation of unipolar, bipolar, threshold, and complementary resistive switching modes is accomplished confirming the coexistence of two switching mechanisms of opposite polarity as the basis for all observable switching phenomena in this material. We show that only calibration of the outer current-voltage loops with the memdiode model is necessary for predicting the device behavior in the defined region revealing additional information on the switching process. The correspondence of each memdiode device with the conduction characteristics of the individual top and bottom metal-oxide contacts allows one to assess the role played by each interface in the switching process separately. This identification paves the path for a future improvement of the device performance and functionality by means of appropriate interface engineering.