Dose painting by numbers (DPBN) refers to a voxel-wise prescription of radiation dose modelled from functional image characteristics, in contrast to dose painting by contours which requires delineations to define the target for dose escalation. The direct relation between functional imaging characteristics and DPBN implies that random variations in images may propagate into the dose distribution. The stability of MR-only prostate cancer treatment planning based on DPBN with respect to these variations is as yet unknown. We conducted a test-retest study to investigate the stability of DPBN for prostate cancer in a semi-automated MR-only treatment planning workflow. Twelve patients received a multiparametric MRI on two separate days prior to prostatectomy. The tumor probability (TP) within the prostate was derived from image features with a logistic regression model. Dose mapping functions were applied to acquire a DPBN prescription map that served to generate an intensity modulated radiation therapy (IMRT) treatment plan. Dose calculations were done on a pseudo-CT derived from the MRI. The TP and DPBN map and the IMRT dose distribution were compared between both MRI sessions, using the intraclass correlation coefficient (ICC) to quantify repeatability of the planning pipeline. The quality of each treatment plan was measured with a quality factor (QF). Median ICC values for the TP and DPBN map and the IMRT dose distribution were 0.82, 0.82 and 0.88, respectively, for linear dose mapping and 0.82, 0.84 and 0.94 for square root dose mapping. A median QF of 3.4% was found among all treatment plans. We demonstrated the stability of DPBN radiotherapy treatment planning in prostate cancer, with excellent overall repeatability and acceptable treatment plan quality. Using validated tumor probability modelling and simple dose mapping techniques it was shown that despite day-to-day variations in imaging data still consistent treatment plans were obtained.