Under gametophytic self-incompatibility (GSI), plants are heterozygous at the self-incompatibility locus (S-locus) and can only be fertilized by pollen with a different allele at that locus. The last century has seen a heated debate about the correct way of modeling the allele diversity in a GSI population that was never formally resolved. Starting from an individual-based model, we derive the deterministic dynamics as proposed by Fisher (1958), and compute the stationary S-allele frequency distribution. We find that the stationary distribution proposed by Wright (1964) is close to our theoretical prediction, in line with earlier numerical confirmation. Additionally, we approximate the invasion probability of a new S-allele, which scales inversely with the number of resident S-alleles. Lastly, we use the stationary allele frequency distribution to estimate the population size of a plant population from an empirically obtained allele frequency spectrum, which complements the existing estimator of the number of S-alleles. Our expression of the stationary distribution resolves the long-standing debate about the correct approximation of the number of S-alleles and paves the way to new statistical developments for the estimation of the plant population size based on S-allele frequencies.

Codes and datasets for generating the figures in the main text and the Supplementary Information of the manuscript entitled "Revisiting the number of self-incompatibility alleles: from old models to new results" published in the Journal of Evolutionary Biology.

Details are provided in the readme file.

Generally, the code is structured per figure in the manuscript, i.e., all code to reproduce Figure X in the manuscript is contained in the folder labeled "FigureX".

However, because Dryad requested me to do so, data and processing files need to be downloaded separately. Then the data needs to be copied into the respective folders of the software to run the codes. I am sorry for this complication, I do not like it either.