R中的直方图跟踪库存水平答案

【问题标题】：Histogram tracking inventory levels in RR中的直方图跟踪库存水平
【发布时间】：2016-12-13 15:54:34
【问题描述】：

我正在寻找一种方法来可视化全天的库存。数据集如下所示，最后两列的摘要如下：

                 Time  Price Inventory Duration
1 9/1/2016 9:25:06 AM 13.960    318        0
2 9/1/2016 9:36:42 AM 13.980    106      696
3 9/1/2016 9:40:52 AM 13.990   -599      250
4 9/1/2016 9:52:54 AM 14.015     68      722
5 9/1/2016 9:52:54 AM 14.015    321        0
6 9/1/2016 9:54:17 AM 14.010     74       83

库存

 Min.  1st Qu.   Median     Mean  3rd Qu.     Max. 
-1120.00   -98.75     9.00     0.00   100.00  1988.00

持续时间

Min.   1st Qu.    Median      Mean   3rd Qu.      Max. 
 0.00     40.25    205.50   2100.00    529.00 272700.00

我想通过显示在不同库存水平上花费的时间来可视化数据。你会推荐什么作为这个功能？到目前为止，我只发现了基于频率的直方图，而不是时间。我的预期结果类似于：

https://postimg.org/image/z074waij1/

提前致谢

【问题讨论】：

你可以试试这个：barplot(prop.table(table(df$Inventroy,df$Duration)),beside=T)

标签： r histogram

【解决方案1】：

我为我的需要编写了以下函数。希望对你有帮助

inv.barplot.IDs <- function(inv.list, IDs = 1:1620)
  {
    # Subset according to the IDs
    myinvs <- as.data.frame(matrix(nrow = 0, ncol = 14))
    names(myinvs) <- inv.names
    Volume <- Duration <- vector("numeric")

    for (i in IDs)
    {
      #myinvs <- rbind(myinvs, inv.list[[i]])
      Volume <- c(Volume, as.numeric(inv.list[[i]]$Volume))
      Duration <- c(Duration, as.numeric(inv.list[[i]]$Duration))
    }

    # Design a sequence of skatules
    minimum <- min(Volume)
    maximum <- max(Volume)
    width <- (maximum + abs(minimum)) / 18
    width <- round(width, -1)
    seq.pos <- seq(width, maximum + width, by = width)
    seq.neg <- - seq(0, abs(minimum) + width, by = width)
    seq <- c(rev(seq.neg), seq.pos)

    # Categorize the dataframe (new column)
    Skatule <- numeric(length = length(Volume))
    for (i in 1:length(Volume))
    {
      Skatule[i] <- seq[head(which(seq > Volume[i]), 1) - 1]
    }

    barplot.data <- tapply(Duration, Skatule, sum)

    # Save the barplot
    #jpeg(filename = file.barplot, width = 480 * (16/9))
    inv.barplot <- barplot(barplot.data, border = NA, ylim = c(0, max(barplot.data)), main = "Total time spent on various inventory levels", xlab = "Inventory", ylab = "Log of Hours")
    #print(inv.barplot)
    #dev.off()
  }

【讨论】：