pm = singletons(c("V1", "V2"), sex = 1:2)
missing = paste0("M", 1:4)
am = list(nuclearPed(children = missing[1:3]),
nuclearPed(children = missing[4], sex = 2))
dvi = dviData(pm, am, missing)
generatePairings(dvi)
Run the code above in your browser using DataLab