Calculates a kinship matrix using the Malecot Migration Model, in the form described by L. B. Jorde 1982.

1 |

`S` |
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 |

`P` |
the column stochastic migration matrix, possibly obtained using col.sto on the "raw" migration matrix |

`N` |
the vector of effective populations, where each element is the population size for all the n populations divided by 3 |

`n` |
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 federico.calboli@helsinki.fi

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 |

Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.

All documentation is copyright its authors; we didn't write any of that.