Simulation of the evolution of inbreeding load under partial selfing

Starting with allele frequencies corresponding to the mutation-selection balance, this program simulates the evolution of independent loci under mutation and selection in populations with N monoecious individuals. These populations are run for t generations (generaciones) assuming a given selfing rate. The inbreeding load is computed as B=Σ2dpq (the addition is over loci; see article), and its value at the Mutation-Selection-Drift-Inbreeding balance must be obtained by averaging B over a number of generations well after it reached a plateau. Values were averaged over a number of simulation replicates (repeticiones).