【问题标题】:Change x-variable bin counts of a histogram更改直方图的 x 变量 bin 计数
【发布时间】:2017-06-09 03:58:18
【问题描述】:

我有 0.25 公顷的地块,我需要将我的数据显示为 1 公顷。我正在尝试制作下图,但将计数乘以 4(所以我有整公顷而不是四分之一)。但是,所有帖子似乎都在处理更改轴标题、值等,但我需要更改实际的直方图频率计数。

Histogram x-variable in size classes plotted by factor variable

ggplot(liveTrees, aes(diam1DBH)) +
geom_histogram(binwidth =10) +
facet_wrap(~site) +
ggtitle("Stems/0.25ha by Size Class") +
ylab("Stems/0.25ha") +
xlab("Diameter Class")

liveTrees = my data 
diam1DBH  = diameter (numeric, continuous) 
site = plot location (factor)

原代码: 我试过的:`

for (i in 1:length(unique(liveTrees$site))) { 
    test<-hist(liveTrees[liveTrees$site== unique(liveTrees$site)[i], "diam1DBH"],   plot = F) 
    b <- barchart(test$counts*4, width = 10, xlim=c(0,350), cex.axis = 0.85) 
    axis(side = 1, at = "b", cex.axis = 0.85)
}

但我不断得到

轴错误(side = 1, at = "b", cex.axis = 0.85):没有位置 有限另外:警告消息:在轴(边= 1,在=“b”, cex.axis = 0.85) : 强制引入的 NAs

所以,有了这个我可以得到计数,但数字不正确,而且它们不是有用的格式。

我的数据是一个data.frame,例如:data example

我需要的是每个直径等级的总和,每个 bin 频率量,乘以 4。我一直在尝试这样做,但无法让它工作,感谢任何帮助!

【问题讨论】:

    标签: r count histogram


    【解决方案1】:

    如果将频率乘以 4,值会发生变化,但图形看起来仍然相同,因此有两种选择,一种是简单地更改轴值标签,另一种更简单的方法是添加数据4次。例如:

    ggplot(rbind(data, data,data,data), aes(variable_X)) + geom_histogram(binwidth =10)
    

    这样数据会成倍增加,不会产生新的data.frame,这可能会混淆以后的分析。

    【讨论】:

      猜你喜欢
      • 2013-10-08
      • 2016-12-25
      • 1970-01-01
      • 1970-01-01
      • 1970-01-01
      • 1970-01-01
      • 1970-01-01
      • 1970-01-01
      • 2018-06-10
      相关资源
      最近更新 更多