昔、家系図を作成したのだが、某業界で有名な一家の家系図をRで作成したいと思った。
ここやここを参考にデータ入力し、作図した。
もっと色づけとかやりたかったけど保留。
こんな家系の家族歴問診は激ムズだろう。
library(kinship) school_days_cast <- seq(59) Mother <- c(6, 9, 11, 0, 0, 0, 0, 0, 5, 5, 0, 0, 0, 7, 8, 10, 10, 10, 11, 11, 12, 12, 13, 0, 30, 14, 14, 14, 14, 9, 0, 9, 17, 18, 25, 21, 21, 21, 24, 23, 23, 0, 30, 30, 30, 42, 31, 32, 45, 49, 49, 0, 52, 52, 50, 50, 55, 43, 30) Father <- c(4, 1, 1, 0, 0, 0, 0, 0, 4, 4, 0, 0, 0, 1, 1, 2, 2, 1, 1, 1, 1, 2, 1, 0, 1, 2, 2, 2, 1, 1, 0, 1, 2, 2, 20, 20, 20, 1, 3, 3, 3, 0, 1, 1, 1, 2, 2, 3, 1, 1, 1, 0, 1, 1, 1, 1, 20, 3, 2) Sex <- c(1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 1, 2, 2, 2, 2, 1, 2, 1, 2, 2, 2, 1, 1, 2, 1, 2, 2, 2, 1, 2, 3, 1, 1, 3, 3, 2, 2, 2, 2, 2, 2, 2, 1, 2, 2, 2, 3, 2, 2, 1, 2, 3, 3, 2, 1) Affected <- rep(1, length(school_days_cast)) Affected[1:3] <- 2 #沢越止、伊能歩、間瞬 Status <- rep(0, length(school_days_cast)) Status[c(48, 54, 58)] <- 1 #伊藤誠、西園寺世界、清澄刹那 cbind(school_days_cast, Mother, Father) Makoto <- pedigree( id = school_days_cast, momid = Mother, dadid = Father, sex = Sex, affected = Affected, status = Status ) plot(Makoto) plot(Makoto, col=2)