2010-09-14 76 views
0

假设我有一个指定的位置(经/纬度),高度列个人的R数据帧,性别:计算/绘制统计

x <- data.frame(
    lat=c(39.5,39.51,38,38.1,38.2), 
    long=c(86,86,87,87,87), 
    gender=c("M","F","F","M","F"), 
    height=c(72,60,61,70,80) 
) 

我想斌数据两个维度(例如,进入千米X千米正方形)和计算以下(然后在地图上显示):

  1. 什么个人的百分比在每个仓是女性
  2. 什么是男性的平均高度在每个仓

如果可能的话,我想用GGPLOT2。

+0

你应该看看SP和岩石圈包。也许有人向你展示如何处理这些问题,因为它们不是最容易使用的。 – 2010-09-14 14:49:41

回答

0

见切换功能的方法,从纬度/经度转换成箱(还看的帮助在削减其他选项也见部分)。这假定矩形网格足够好,如果区域足够大以便地球曲率影响东西,则需要查看空间包。

然后使用tapply或plyr包计算网格的每个单元内的摘要。

现在绘制不过你想要的结果。

+0

另请参阅reshape/plyr中的'round_any'函数。 – hadley 2010-09-15 03:09:14