mal.eq {Biodem}

Calculates the asymptotic generation for the Malecot model

Description

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.

Usage

mal.eq(S, P, N)

Arguments

S

is the Sistematic pressure matrix.

P

is the colum-stochastic migration matrix.

N

is the vector of effective population size.

Details

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.

Value

Returns one numeric value.

Note

This function has been coerced to use "only" six significant digits. ...

Author(s)

Federico C. F. Calboli [e-mail]

References

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

See Also

mal.phi for the function using the output of 'mal.eq'

Examples

# the data is originally from a paper by Swedlund et al. 1984.
data(S); data(P); data(N)
mal.eq(S,P,N)