将现有直方图导入R

Question

我已经使用ROOT分析了一个庞大的数据集，并在文本文件中以“（double）bin position”“（int）bin value”的形式保存了直方图。 E.g：

0.1 0
0.2 1
0.3 4
0.4 2
0.5 1
...

我已将txt文件导入R并且现在有两个向量。到目前为止，我无法在R中找到允许导入已存在的直方图的函数。这样的功能应该允许设置直方图的各个区间的区间位置和区间值（=频率）。如果你能帮助我，那就太好了。

干杯

Answer 1

你可以尝试barplot（可能是最简单的方法）。

或者查看hist中的代码：

# ....
r <- structure(list(breaks = breaks, counts = counts, density = dens, 
    mids = mids, xname = xname, equidist = equidist), class = "histogram")
if (plot) {
    plot(r, freq = freq1, col = col, border = border, angle = angle, 
        density = density, main = main, xlim = xlim, ylim = ylim, 
        xlab = xlab, ylab = ylab, axes = axes, labels = labels, 
        ...)
# ....

所以组建一个类似上面r的结构，为它分配类＆＃34;直方图＆＃34;并且只需使用plot。或者，如果您不想为其分配课程，请使用graphics:::plot.histogram。

Answer 2

这应该有效：

x <- c(0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7)
y <- c(0, 1, 4, 2, 1, 0, 3)

library(MASS)
truehist(rep(x,y), prob = F, h = 0.1)