Retinal microsurgery is a high-precision surgery performed on a delicate tissue requiring the skill of highly trained surgeons. Given the restricted range of instrument motion in the confined intraocular space, snake-like robots may prove to be a promising technology to provide surgeons with greater flexibility, dexterity, and positioning accuracy during retinal procedures such as retinal vein cannulation and epiretinal membrane peeling. Kinematics modeling of these robots is an essential step toward accurate position control. Unlike conventional manipulators, modeling these robots does not follow a straightforward method due to their complex mechanical structure and actuation mechanisms. The hysteresis problem can especially impact the positioning accuracy significantly in wire-driven snake-like robots. In this paper, we propose a data-driven kinematics model using a probabilistic Gaussian mixture model (GMM) and Gaussian mixture regression (GMR) approach with a hysteresis compensation algorithm. Experimental results on the two-degree-of-freedom (DOF) integrated robotic intraocular snake (I2RIS) show that the proposed model with the hysteresis compensation can predict the snake tip bending angle for pitch and yaw with 0.45° and 0.39° root mean square error (RMSE), respectively. This results in overall 60% and 70% improvements of accuracy for yaw and pitch over the same model without the hysteresis compensation.