The purpose of this article was to extend the model used to predict selection response with selfed progeny from 2 alleles per locus to a model which is general for number and frequency of alleles at loci. To accomplish this, 4 areas had to be dealt with: 1) simplification of the derivation and calculation of the condensed coefficients of identity; 2) presentation of the genetic variances expressed among and within selfed progenies as linear function of 5 population parameters; 3) presentation of selection response equations for selfed progenies as functions of these 5 population parameters; and 4) to identify a set of progeny to evaluate, such that one might be able to estimate these 5 population parameters.The five population parameters used in predicting gains were the additive genetic variance, the dominance variance, the covariance of additive and homozygous dominance deviations, the variance of the homozygous dominance deviations and a squared inbreeding depression term.