用R语言绘制显示taxa特征值的散点图

    在系统发育分析时，经常会对树上taxa的性状进行描述性统计。对于离散型的变量，R中各种包（APE，phytools，etc）自带的函数可以较好解决问题，但对于连续型变量（如体长、体重等），则没有对应的函数进行显示。

   其实，可以通过Hmisc包中的dotchart3函数绘制“克利夫兰点图”，再将其拼接到系统发育树上解决此问题。效果如下图：

       图中，系统发育树右侧的dot plots是用dotchart3()绘制的taxa某个连续变量性状的值。其中，每个点的颜色、图标也可根据需要进行设定，用于同步显示与连续变量相关的离散变量的值。在此图中，点图表示的是各皿蛛的足长比，点颜色表示气管类型，点图标代表体色类型，点的位置和taxa在树上的位置对应。这样一图多用，是一种系统发育树上综合显示taxa各性状的好方法。

下面简单说一下实现方法：

1. R语言中安装Hmisc包

2.使用dotchart3()函数绘图。

   dotchart3 (x=[x],label=[tiplabel],
             groups = [len_x]:1,gdata =[x],gpch=[gpch],gcolor = [gcolor],
             color ="white",lcolor = "white")

其中࿱