Question

我有像这样名为PACKS的数据：

error    CENTRE_B
    -13  1104
    -13  1303
    -13  1303
    2    1204
    2    1403
    2    1403
    2    1403
    2    1502
    3    1503

我的目标是比较分布。我想为CENTRE_B的每个值绘制误差的直方图。问题是直方图必须是相同的比例。

我试过了：

new.par <- par(mfrow=c(5, 5))
histograms = aggregate(error ~ CENTRE_B, PACKS, hist)

绘制直方图。但是，我不知道如何在聚合中为hist提供额外的参数（breaks = c（-80,80））。

另一个问题是直方图似乎是尺寸为2x45的data.table，因此它不包含直方图。所以我不知道如何自动更改参数。

提前感谢任何建议。

Answer 1

忘掉aggregate。只需使用split并执行简单的for循环。然后，您可以将所有参数输入hist。

new=split(PACKS$error,PACKS$CENTRE_B)
for(i in 1:length(new))
hist(new[[i]],main="Some fancy title", sub="don't forget you can use paste to change the title between loops")

聚合并绘制具有相同比例和轴的直方图

1 个答案: