Rで家系図

MikuHatsune2011-04-23

こちらから。
id
momid
dadid
sex(1は男、2は女、3は性別不明)
を指定すると勝手にがんばって描いてくれる。
オプションとして
affected(1は罹患なし、2は罹患あり、3は不明)
status(0は生存、1は死亡)
relarions(3行の行列指定。1、2列が兄妹、3列目が関係。1は一卵性双生児、2は二卵性双生児、3は不明、4は配偶者)

library(kinship)
test <- pedigree(
  id        = c(1,2,3,4,5,6,7,8,9,10,11,12,13,14),
  momid     = c(0,0,0,0,2,2,4,4,6, 2, 0, 0,12,13),
  dadid     = c(0,0,0,0,1,1,3,3,3, 7, 0, 0,11,10),
  sex       = c(0,1,0,1,0,1,0,0,3, 0, 0, 1, 1, 4),
  affected  = c(1,2,1,2,1,2,1,2,1, 2, 1, 2, 1, 2),
  status    = c(1,1,1,1,0,0,0,0,0, 0, 0, 0, 0, 1),
  relations = matrix(c(c(5,6,2),c(7,8,1)),nc=3,byrow=TRUE)
)
plot(test)