R中的构建图（Plot）：通过时间间隔的中位数价格

Question

我有数据框架的价格和一些拍卖的结束日期。我想检查何时出现，例如，具有最小和最大价格（也是中位数）的销售，取决于一天的时间。

更确切地说，我有数据框mtest：

> str(mtest)
'data.frame':   9144 obs. of  2 variables:
$ Price      : num  178 188 228 305 202 ...
$ EndDateTime: POSIXct, format: "2015-05-25 05:00:59" "2015-05-23 00:06:01"  ...

我想建立图表（图表），在X轴上有30分钟时间内部（00：00-00：30,00：31-01：00等），并且中位数（最大，最小价格）在Y轴。
另一个想法是为每个时间间隔绘制一个简单的直方图，如hist(mtest$Price, breaks=10, col="red")

我怎样才能以最好的方式做到这一点？

Answer 1

试试这个：

cutt=seq(from=min(mtest$EndDateTime),to=max(mtest$EndDateTime), by=30*60)
if (max(mtest$EndDateTime)>max(cutt))cutt[length(cutt)+1]=max(cutt)+30*60
mtest$tint=cut(mtest$EndDateTime,cutt)
stats=do.call(rbind,tapply(mtest[,"Price"],mtest[,"tint"],
                     function(p)c(min=min(p),median=median(p),max=max(p))))
bp=boxplot(mtest[,"Price"]~mtest[,"tint"],xaxt="n",
           col=1:length(levels(mtest$tint)))
axis(1,at=1:length(levels(mtest$tint)),labels=format.Date(levels(mtest$tint),"%Y-%m-%d %H:%M"),
     las=2,cex.axis=.5)
stats

or wilth plot

plot(NA,ylim=range(stats),xlim=c(1,lint),type="n",xaxt="n",xlab="",ylab="")
sapply(1:3,function(z)points(stats[,z]~c(1:lint),col=z))
axis(1,at=1:lint,labels=format.Date(levels(mtest$tint),"%Y-%m-%d %H:%M"),
     las=2,cex.axis=.5)

你会有这样的事情：