# NOT RUN {
data(valley)
valley #a subset of a real marriage data base
mar <- sur.freq(valley,valley$PAR,valley$SURM,valley$SURF,freq.table="marriages")
mar # frequency tables of the observed pairs of surnames in each population
iso <- mar.iso(mar)
iso # a data frame containing Pt and Pr values for each (sub)population
# }
Run the code above in your browser using DataLab