mayoverse/kinship2

Incorrect twins

Opened this issue · 0 comments

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