[版]
我终于用第三种方式获得了一些东西。我用过:
ddply()( plyr 包)获取每组点的质心(组由因子变量确定),
< / LI>table()计算另一个因子变量的每个级别的频率,并从中创建单独的数值变量(用于馅饼的切片),并最终将数据集简化为唯一的行每个小组,
add.pie()( mapplots 包)在 for 循环中显示所有这些馅饼。结果接近于phyllogeography研究的结果,除了数据完全不同:
Here's the result(这里有假数据)。无法在此处直接发布任何图像,我需要在StackOverflow上获得更多声望。
我认为可以通过改进数据排列方式来完成某些事情。它需要进一步的工作,但我现在遇到的问题已经解决了。
[原帖]
我希望在地图上绘制全欧洲的GPS数据。这些点是分批收集的,即。,当记录GPS点时,很有可能在附近保存了5或6个其他GPS点(比方说&lt; 2 km)。当在如此大的区域的地图上绘制时,这些批处理点彼此太靠近而无法区分,并且我无论如何都不需要高精度的地图,但我希望看到所有的点。我试过抖动点,但是,即使这个特定地图的精度不是强制性的,我也不习惯抖动GPS坐标。
显示所有点的另一种方法是例如向日葵表示,或者饼图在这里更合适,因为我可以在每个饼上显示不同颜色的切片,一种颜色与一个因子变量水平相关联对于每个数据点。但是然后将我需要的附近点分组:
舍入坐标,以便关闭点实际上具有 完全相同的坐标(但找到最佳点来围绕数据 会很棘手),
在点周围设置缓冲区,以便在有多个缓冲区时 叠加,然后相关的点全部绘制在 平均坐标。也许是距离上的聚类方法 从坐标计算的矩阵也可以起作用。
计算每组的质心(我可以使用现有的 因子变量用于按区域分组点)在那里绘制饼图,但随后是 饼图切片仍然需要从个人收集数据 分。
第二种和第三种解决方案似乎对我来说最强大(但我不确定,请告诉我),但我完全不知道如何做到这一点。有任何想法吗?对于解决方案3,我计算了每组点的质心(组=区域,因子变量),以便数据框中的每一行还包含其对应组的质心坐标。然后我可以在这些坐标处绘制饼图(每组一个饼图),但如果我希望饼图片代表从每个组中的各个行中获取的值,我不知道如何继续。
我已经使用 maps 包的worldHires的更新版本创建了地图。我的数据点目前在地图上单独显示。我对饼图没有经验,但我想我可以在这里找到一些东西,或者在第一个关于分组的问题得到解决时再问一些问题。
非常感谢您的帮助和/或建议!