Question

我需要在R中绘制这个结构:

(1000) Diseases of the genitourinary system 
    (1580) Nephritis, nephrotic syndrome, and nephrosis 
        (580) Acute glomerulonephritis
            (580.9) Glomerulonephritis, acute, unspec.
        (581) Nephrotic syndrome
            (581.9) Nephrotic syndrome, unspec.
        (582) Chronic glomerulonephritis
            (582.9) Glomerulonephritis, chronic, unspec.
        (583) Nephritis and nephropathy, not specified as acute or chronic
        (584) Acute renal failure
            (584.5) Renal failure, acute w/ tubular necrosis

as a nice jpg/pdf/(or other) which will have nodes, connections and labels using R of the structure above. I looked at libraries which require GrafViz installed and had no luck so native solution (e.g., using ggplot2) would be best. I was not able to put together the code using igraph. I am also new to graphics in R with no textbook based foundation. Any hints or advice would be appreciated.

上述结构仅举一例。其它结构可能拥有50+的概念, 可以作为非常大的 PDF/ poser 绘制和打印。此图将是静态的( 没有使用鼠标与其互动 ) 。

Answer 1

这是对答案的一种努力。它接近于你需要的东西,最后的结果是这样的:

"https://i.sstatic.net/EBEEyq.png" alt="名词树"/>

我在igraph 中做了这个,代码使用您描述的数据类型的模拟。

library("igraph")

vertex.df <- read.table(text = "id    code  name
0   1000     Diseases of the genitourinary system  
1   1580     Nephritis, nephrotic syndrome, and nephrosis  
2   580  Acute glomerulonephritis 
3   580.9    Glomerulonephritis, acute, unspec. 
4   581  Nephrotic syndrome 
5   581.9    Nephrotic syndrome, unspec. 
6   582  Chronic glomerulonephritis 
7   582.9    Glomerulonephritis, chronic, unspec. 
8   583  Nephritis and nephropathy, not specified as acute or chronic 
9   584  Acute renal failure 
10  584.5    Renal failure, acute w/ tubular necrosis ",
                        header = TRUE,
                        stringsAsFactor = FALSE)

vertex.df$code <- as.character( vertex.df$code )

edge.df <- read.table(text = "from    to
0    1
1   2
1   4
1   6
1   8
1   9
2   3
4   5
6   7
9   10",
                      header = TRUE)

edges <- matrix(c(edge.df$from, edge.df$to), nc=2)

g <- graph.empty()
g <- add.vertices(g, nrow(vertex.df),
                  id=vertex.df$id, 
                  code=vertex.df$code, 
                  name=vertex.df$name)
g <- add.edges(g, t(edges))


plot(g, 
     layout = layout.kamada.kawai,
     vertex.label = V(g)$code,
     vertex.size = 35,
     vertex.color = "white",
     vertex.label.family = "sans")

我使用 ICD 代码作为顶点标签。这是因为在此刻度绘制时, 疾病名称的长文本看起来不整齐。

绘制时, 您可以将 < code> vertex. label 参数更改为 < code> V(g)$name , 如果您想要疾病名称, 而不是 ICD 代码。我怀疑, 如果您打印到一个大的 pdf 并删除顶点外观, 您可能会得到一棵好看的树。请查看 < code>? igraph. plotting , 了解您可以更改的参数的细节。

我希望这能让你的实验更进一步

友情链接