Calculates a kinship matrix using the Malecot Migration Model
Calculates a kinship matrix using the Malecot Migration Model, in the form described by L. B. Jorde 1982.
the sistematic pressure matrix, where the diagonal elements are 1-sk, with sk the sistematic pressure for the k-th population, and the non diagonal elements are 0
the column stochastic migration matrix, possibly obtained using col.sto on the "raw" migration matrix
the vector of effective populations, where each element is the population size for all the n populations divided by 3
the number of iterations needed to reach the equilibrium, calculated by the function Mal.eq
The Malecot model is simply an iterative markow-chain-like process that gives rise to an asymptotic growth curve, so that an equilibrium is reached after a number of iterations.
Returns a square and symmetrical matrix.
Federico C. F. Calboli email@example.com
Imaizumi, Y., N. E. Morton and D. E. Harris. 1970. Isolation by distance in artificial populations. Genetics 66: 569-582.
Jorde, L. B. 1982. The genetic structure of the Utah mormons: migration analysis. Human Biology 54(3): 583-597.
mal.eq for the function generating the number of cycles needed to reach the asymptotic value
1 2 3 4 5
Want to suggest features or report bugs for rdrr.io? Use the GitHub issue tracker.