mal.eq {Biodem} | R Documentation |
Mal.eq calculates the Malecot model iteratively, stopping when one more cycle adds 0 to every value of the matrix obtained by the model. Once equilibrium is reached, Mal.eq returns the number of cycles ("generations") needed to reach it.
mal.eq(S, P, N)
S |
is the Sistematic pressure matrix. |
P |
is the colum-stochastic migration matrix. |
N |
is the vector of effective population size. |
The use of mal.eq is necessary before the calculation of the Malecot model proper because the value returned by Mal.eq is one of the arguments of the Malecot model function Phi.
Returns one numeric value.
This function has been coerced to use "only" six significant digits. ...
Federico C. F. Calboli f.calboli@gmail.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.
Swedlund, A. C., L. B. Jorde and J. H. Mielke. 1984. Population structure in the Connecticut valley. I. Marital migration. American Journal of Physical Anthropology 65: 61-70
mal.phi
for the function using the output of 'mal.eq'
# the data is originally from a paper by Swedlund et al. 1984. data(S); data(P); data(N) mal.eq(S,P,N)