使用trees包绘制带有一些松鼠的树时出错

Posted 2023-03-31

【中文标题】使用trees包绘制带有一些松鼠的树时出错

【英文标题】：Error while plotting a tree with some squirrels using trees package

【发布时间】：2015-03-21 19:48:16

【问题描述】：

我正在使用@jbaums 找到的here 包trees 并解释了in this post。

我的数据如下：

这棵树是由 主干

    Trunk
[1] 13.60415

和分支机构

Tree
   TreeBranchLength TreeBranchID
1         10.004269            1
2          7.994269            2
3          9.028834           11
4         10.817401           12
5          8.551311          111
6         10.599798          112
7         11.073243          121
8         13.367392          122
9          9.625431         1111
10        10.793569         1112
11         9.896499        11121
12         8.687741        11122
13         7.791180         1211
14        12.506105         1212
15         6.768478         1221
16        10.441796         1222
17        10.751892         1121
18         9.458651         1122
19        10.768509        11221
20        10.150673        11222
21        12.377448       111211
22        12.235136       111212
23         9.074079        11211
24         9.996334        11212
25         9.807019       112221
26        10.895809       112222
27         6.741274      1122211
28        15.841272      1122212
29         5.753920     11222111
30         8.846389     11222112
31        11.925961       112111
32         9.780776       112112
33         8.207965        12221
34        10.079375        12222

50 只松鼠种群 -

Populations
   PopulationPositionOnBranch PopulationBranchID ID
1                  10.6321655             112111  1
2                   1.0644897                  1  2
3                   3.9315473                  1  3
4                   1.0310244                  0  4
5                   9.1768846                  0  5
6                  13.4267181                  0  6
7                   7.9461528                  0  7
8                   6.0533401                121  8
9                   2.1227425                121  9
10                  1.8256787                121 10
11                  4.7332588           11222112 11
12                  4.4837432           11222112 12
13                  4.6200834           11222112 13
14                  2.5622276               1221 14
15                  1.2446683               1221 15
16                  7.0674052                111 16
17                  1.3854674                111 17
18                  4.8735635                111 18
19                  9.5007998               1222 19
20                  6.6373468               1222 20
21                 12.6757728                122 21
22                  4.2685465                122 22
23                  3.9806540                  2 23
24                  3.1025403                  2 24
25                  3.9119065              11122 25
26                  1.5527653              11122 26
27                  1.6687957              11122 27
28                  8.0697456               1122 28
29                  6.7871391               1122 29
30                  9.8050713             111212 30
31                  8.5226920             111212 31
32                  3.6113379             111212 32
33                  7.3184965             111211 33
34                  8.6142984             111211 34
35                  1.3550870               1211 35
36                  8.3650639                 12 36
37                  4.6411446             112112 37
38                  3.2985541             112112 38
39                 12.2344148               1212 39
40                  9.0290776               1212 40
41                  1.3900249               1121 41
42                  0.9261425            1122212 42
43                 15.2522199            1122212 43
44                  4.0253771              12222 44
45                  8.7507678              11222 45
46                  4.6289841            1122211 46
47                  9.1799522                112 47
48                  5.1293838              12221 48
49                  1.1543080              12221 49
50                 10.1014837             112222 50

生成情节的代码

        g <- germinate(list(trunk.height=Trunk, 
                   branches=Tree$TreeBranchID,
                   lengths=Tree$TreeBranchLength), 
              left='1', right='2', angle=30))

    xy <- squirrels(g, Populations$PopulationBranchID, pos=Populations$PopulationPositionOnBranch, 
               left='1', right='2', pch=21, bg='white', cex=3, lwd=2)
    text(xy$x, xy$y, labels=seq_len(nrow(xy)), font=1)

，产生

正如您在下面的图上看到的那样，人口 43（蓝色箭头）超出了树的范围。图上的分支长度似乎与数据不对应。例如，人口 38 和 37 所在的分支（左侧绿色箭头）比人口 43 所在的分支（右侧绿色箭头）长，数据中并非如此。 我做错了什么？我是否正确理解如何使用trees？

【问题讨论】：

有史以来最好的标题 +1

大声笑。我不知道如何帮助你。但是因为您已经设法在涉及树上松鼠的编程问答中提出了一个完全合法的问题，先生您已经赢得了饼干。

:) 欣赏！希望这会让我的问题得到更多关注

感谢您的提醒 - 现在是 fixed。

【参考方案1】：

在研究 germinate 函数时，在我看来，您传递给它的 Tree 值需要在 TreeBranchId 字段中按升序排序。

您放置 43 的 BranchID: 1122212 不是实际的 1122212 分支。 由于您在树中输入值的顺序，该函数以某种方式弄乱了分支的位置。

我很好奇如果增加Branch ID:1122212的长度，会不会改变放置43的分支，你猜怎么着？它没有。实际显示长度增加的分支是您放置 37 和 38 的分支。

所以这个提示指出 germinate 函数有问题。在进一步调试时，我能够使用以下代码使其工作。

Tree<-read.csv("treeBranch.csv")
Tree<-Tree[order(Tree$TreeBranchID),] 
g <- germinate(list(trunk.height=15, 
                branches=Tree$TreeBranchID,
                lengths=Tree$TreeBranchLength), 
               left='1', right='2', angle=30)

xy <- squirrels(g, Populations$PopulationBranchID,pos=Populations$PopulationPositionOnBranch, 
            left='1', right='2', pch=21, bg='white', cex=3, lwd=2)
text(xy$x, xy$y, labels=seq_len(nrow(xy)), font=1)

【讨论】：

非常感谢！你让我开心:)

@Sulawesi：我很高兴它成功了。另外，谢谢你分享这个问题。我真的很喜欢这个树包。