R2D2 image reconstruction with model uncertainty quantification in radio astronomy