Incorrect twins
Opened this issue · 0 comments
magnusdv commented
In the example below individuals 3 and 4 are MZ twins, but the plot incorrectly shows 3 and 5 as MZ twins.
x = kinship2::pedigree(
id = 1:7,
dadid = c(0, 0, 1, 1, 1, 3, 4),
momid = c(0, 0, 2, 2, 2, 5, 5),
sex = c(1, 2, 1, 1, 2, 1, 1),
relation = data.frame(id1 = 3, id2 = 4, code = 1)) # 3 and 4 MZ twins
plot(x)
#> Warning in autohint(ped): Unexpected result in autohint, please contact
#> developer
Created on 2023-02-15 with reprex v2.0.2