Hyperglycemia is a common metabolic problem in premature, low-birth-weight infants. Blood glucose homeostasis in this group is often disturbed by immaturity of endogenous regulatory systems and the stress of their condition in intensive care. A dynamic model capturing the fundamental dynamics of the glucose regulatory system provides a measure of insulin sensitivity (S(I)). Forecasting the most probable future S(I) can significantly enhance real-time glucose control by providing a clinically validated/proven level of confidence on the outcome of an intervention, and thus, increased safety against hypoglycemia. A 2-D kernel model of S(I) is fitted to 3567 h of identified, time-varying S(I) from retrospective clinical data of 25 neonatal patients with birth gestational age 23 to 28.9 weeks. Conditional probability estimates are used to determine S(I) probability intervals. A lag-2 stochastic model and adjustments of the variance estimator are used to explore the bias-variance tradeoff in the hour-to-hour variation of S(I). The model captured 62.6% and 93.4% of in-sample S(I) predictions within the (25th-75th) and (5th-95th) probability forecast intervals. This overconservative result is also present on the cross-validation cohorts and in the lag-2 model. Adjustments to the variance estimator found a reduction to 10%-50% of the original value provided optimal coverage with 54.7% and 90.9% in the (25th-75th) and (5th-95th) intervals. A stochastic model of S(I) provided conservative forecasts, which can add a layer of safety to real-time control. Adjusting the variance estimator provides a more accurate, cohort-specific stochastic model of S(I) dynamics in the neonate.