This paper suggests a methodology to estimate morbidity parameters in a three-state model, which includes the heterogeneity factor modeled by a gamma-distributed random variable. The cases of complete and incomplete information are considered and various risk estimation methods are discussed. We give estimates of morbidity radiation risks in four classes of diseases based on data from the Russian National Radiation Epidemiological Registry. The methodology can be used during the development of radiation protection systems.