2012-07-01 6 views
12

私はファミリーツリーをプロットする方法について調べてきましたが、再現できるものは見つけられませんでした。私はハドレーのggplotに関する本を見てきましたが、同じことをしています。ファミリーツリーをプロットする方法R

私はソースこれに似たデータフレームとして持つ家系図をプロットしたい:

dput(head(familyTree)) 
structure(
    list(
    id = 1:6, 
    cnp = c("11", NA, "22", NA, NA, "33"), 
    last_name = c("B", "B", "B", NA, NA, "M"), 
    last_name_alyas = c(NA, NA, NA, NA, NA, "M"), 
    middle_name = c("C", NA, NA, NA, NA, NA), 
    first_name = c("Me", "P", "A", NA, NA, "S"), 
    first_name_alyas = c(NA, NA, NA, NA, NA, "F"), 
    maiden_name = c(NA, NA, "M", NA, NA, NA), 
    id_father = c(2L, 4L, 6L, NA, NA, 8L), 
    id_mother = c(3L, 5L, 7L, NA, NA, 9L), 
    birth_date = c("1986-01-01", "1963-01-01", "1964-01-01", NA, NA, "1936-01-01"), 
    birth_place = c("City", "Village", "Village", NA, NA, "Village"), 
    death_date = c("0000-00-00", NA, NA, NA, NA, "2007-12-23"), 
    death_reason = c(NA, NA, NA, NA, NA, "stroke"), 
    nr_brothers = c(NA, 1L, NA, NA, NA, NA), 
brothers_names = c(NA, "M", NA, NA, NA, NA), 
    nr_sisters = c(1L, NA, 1L, NA, NA, 2L), 
    sisters_names = c("A", NA, "E", NA, NA, NA), 
    school = c(NA, "", "", NA, NA, ""), 
    occupation = c(NA, "", "", NA, NA, ""), 
    diseases = c(NA_character_, NA_character_, NA_character_, NA_character_, NA_character_, NA_character_), 
    comments = c(NA_character_, NA_character_, NA_character_, NA_character_, NA_character_, NA_character_) 
), 
    .Names = c("id", "cnp", "last_name", "last_name_alyas", "middle_name", "first_name", "first_name_alyas", "maiden_name", "id_father", "id_mother", "birth_date", "birth_place", "death_date", "death_reason", "nr_brothers", "brothers_names", "nr_sisters", "sisters_names", "school", "occupation", "diseases", "comments"), 
    row.names = c(NA, 6L), 
    class = "data.frame" 
) 

は私がggplotで家系図をプロットすることができます任意の方法はありますか?もしそうでなければ、どうすれば別のパッケージを使ってプロットすることができますか?

主キーは 'id'で、 "id_father"と "id_mother"を使用して他のメンバーに接続します。

+4

あなたは何を試してみましたか? Rとそのグラフィックスの使い方を知っていますか?あなたのデータとそのグラフィックスメソッドを表現するためにigraphを使ってみましたか?ファミリツリーは一種のグラフなので、igraphと全体のグラフィカルモデルのタスクビューが適しています。あなたはまだそれをお読みですか? – Spacedman

+0

おそらくggplot2の拡張である[ggphylo](https://github.com/gjuggler/ggphylo)のプロットを改変して、系統樹のプロットを作ることができます。 – sckott

答えて

8

コメントに記されているように、igraphを試してください。ここではクイックスタートは、次のとおりです。

require(igraph) 
mothers=familyTree[,c('id','id_mother','first_name', 'last_name')] 
fathers=familyTree[,c('id','id_father','first_name', 'last_name')] 
mothers$name=paste(mothers$first_name,mothers$last_name) 
fathers$name=paste(fathers$first_name,fathers$last_name) 
names(mothers)=c('parent','id','first_name','last_name','name') 
names(fathers)=c('parent','id','first_name','last_name','name') 
links=rbind(mothers,fathers) 
links=links[!is.na(links$id),] 
g=graph.data.frame(links) 
co=layout.reingold.tilford(g, flip.y=F) 
plot(g,layout=co) 

enter image description here

は、任意の名前はありません、そして矢印は間違った方向に行っているが、あなたはそこから行くことができるはずです。

+1

ありがとうございます。私はigraphを見渡しましたが、どこから始めるべきか分かりませんでした。あなたの例は良いスタートです。 –

+0

私は別の質問があります。私はツリーを180度回転させ、頂点上の数字の代わりに名前を配置しました。問題は、ツリーが混雑しすぎるということです。 x軸上の頂点間にスペースを入れるためにレイアウトを指定する方法はありますか? –

+0

正直言って私は数日前に 'igraph'を学んだだけです。彼らは、彼らのウェブサイト上にいくつか素晴らしいドキュメント(http://igraph.sourceforge.net/)を持っている、私はそれが可能であると確信しています。 – nograpes

0

kinship2パッケージを試しましたか?

library(kinship2) 
df <- data.frame(id = c(1,2,3,4,5,6), sex = c(1,2,1,2,2,2), dadid = c(0,0,0,0,1,3), momid = c(0,0,0,0,2,4), famid = 1) 
relation1 <- matrix(c(2,3,4,1), nrow = 1) 
foo <- pedigree(id = df$id, dadid = df$dadid, momid = df$momid, sex = df$sex, relation = relation1, famid = df$famid) 
ped <- foo['1'] 
plot(ped) 

you can see the resulting plot here

関連する問題