根据遗传数据使用饼图创建地图

Question

我正在进行群体遗传学研究并处理橡树物种的叶绿体数据。我进行了 DAPC，发现有 5 个叶绿体基因型簇。我想制作一张地图，根据我通过饼图或类似的东西采样的 26 个人口来显示每个人的比率。 我用以下列创建了一个 sf 对象：单个样本（Row.names 和 TreeNum）、它们所属的群体（代码）、它们的坐标（纬度、经度和几何=个体之间的平均值，以便给我一个单一的坐标来表示每个人口）和表示它们属于 5 个聚类中哪一个的因子值 (grp.grp)。

  <I<chr>>  <chr>   <chr>     <dbl>    <dbl> <fct>           <POINT [°]>
1 A01       A01     AZH       -123.     38.0 2       (-122.6256 37.9655)
2 A02       A02     AZH       -123.     38.0 4       (-122.6256 37.9655)
3 A04       A04     AZH       -123.     38.0 4       (-122.6256 37.9655)
4 A07       A07     AZH       -123.     38.0 2       (-122.6256 37.9655)
5 A08       A08     AZH       -123.     38.0 4       (-122.6256 37.9655)
6 A09       A09     AZH       -123.     38.0 2       (-122.6256 37.9655)

我已经能够使用 R 中的 tmap 包生成一个地图，该地图按颜色显示每个单独样本的因子值，但这些开始相互堆叠并且难以阅读。 我想生成一个带有饼图的地图，该饼图显示每个总体中因子值的比率，并记录该总体中的样本数量。有些群体有 10 个样本，其他群体只有 4 个样本，所以最好展示一下。 我不会将其用于出版物，只是为了探索数据。

我一直在考虑解决这个问题的方法是改变当前的 sf 对象并创建新的列，这些新列为我提供了我可以使用的每个因子值的频率百分比。如果这样做了，我将需要生成一个新的 sf 对象，该对象被修改为只为每个总体生成一行，然后生成描述我想要使用的值的列。还不清楚我如何将描述每个总体因子值百分比的每一列放入单个列中。 我觉得必须有一种更简单的方法来做到这一点。有人可以给一些好的建议吗？如果这是解决这个问题的最直接的方法，您将如何编码、格式化和创建这个 sf 对象？ 另外，从那里，您建议如何构建这些映射的饼图？我在网上查看并发现了一些有用的帖子，但如何将 tmaps（我当前使用的包）与 ggplot & 结合起来？ scatterpie 似乎不像其他示例那样直接处理我的数据。

任何帮助将不胜感激。谢谢！