rel.phi {Biodem}

Calculates an 'a priori' kinship matrix from isonymy data

Description

"rel.phi" calculates an 'a priori' kinship matrix starting from isonymy data.

Usage

rel.phi(x,R, method="A")

Arguments

x

is a square Unbiased Random Isonymy matrix, possibly obtained using the "uri" function on the raw surname data

R

is an unbiased estimate of Regional Random Isonymy, calculated by the function "rri"

method

a character string specifying the method to be used in the calculation of the coefficients. The available options are "A" and "B". Both the methods give similar results. The "A" method is given as the default option

Details

The function implements Relethford's method to calculate kinship coefficients starting from surname data.

Value

Returns a square symmetric 'a priori' kinship matrix.

Note

The term 'a priori kinship' refers to kinship relative to a founding population

Author(s)

Federico C. F. Calboli and Alessio Boattini [e-mail]

References

Relethford, J. H. 1988. Estimation of kinship and genetic distance from surnames. Human Biology, 60(3): 475-492.

See Also

uri to calculate Unbiased Random Isonymy starting from tables of surname frequencies, rri to calculate an an unbiased estimate of Regional Random Isonymy, rel.cond to calculate a conditional kinship matrix from isonymy data

Examples

# starting from a raw marriage records dataset:
data(valley)
tot <- sur.freq(valley,valley$PAR,valley$SURM,valley$SURF)
tot # a frequency table calculated above all the surnames
iso.matrix <- uri(tot)
iso.matrix # an unbiased random isonymy matrix
reg <- rri(tot)
reg # a coefficient of unbiased Regional Random Isonymy
kin <- rel.phi(iso.matrix,reg)
kin # an 'a priori' kinship matrix

#starting from a generic surname frequency table
data(surnames)
surnames # a made-up dataset
iso.matrix <- uri(surnames)
iso.matrix # an unbiased random isonymy matrix
reg <- rri(surnames)
reg # a coefficient of unbiased Regional Random Isonymy
kin <- rel.phi(iso.matrix,reg)
kin # an 'a priori' kinship matrix