原始数据框如下所示:
status[1:15,]
time dt_5
1 2014-07-31 17:00:01.0 0
2 2014-07-31 17:00:01.2 0
3 2014-07-31 17:00:01.4 1
4 2014-07-31 17:00:01.5 1
5 2014-07-31 17:00:02.2 1
6 2014-07-31 17:00:02.5 1
7 2014-07-31 17:00:03.0 1
8 2014-07-31 17:00:04.0 1
9 2014-07-31 17:00:04.0 1
10 2014-07-31 17:00:05.0 1
11 2014-07-31 17:00:09.0 1
12 2014-07-31 17:00:10.0 1
13 2014-07-31 17:00:10.2 1
14 2014-07-31 17:00:11.2 1
15 2014-07-31 17:00:11.5 1
为了将它分成3分钟间隔组,我使用了George Dontas的解决方案using plyr to group并计算了值的百分比" 1"在每个小组中。
percentage.occupied <- function(x) (NROW(subset(x,dt_5==0)))/(NROW(x))
splitbytime <- dlply(selectstatus327, .(cut(time,"3 mins")),percentage.occupied)
然后我得到了结果:
> splitbytime
$`2014-07-31 17:00:00`
[1] 0.2857143
$`2014-07-31 17:03:00`
[1] 0.2835821
$`2014-07-31 17:06:00`
[1] 0.6044304
$`2014-07-31 17:09:00`
[1] 0.5735294
$`2014-07-31 17:12:00`
[1] 0.8811475
$`2014-07-31 17:15:00`
[1] 0.599359
$`2014-07-31 17:18:00`
[1] 0.8021201
$`2014-07-31 17:21:00`
[1] 0.7956522
$`2014-07-31 17:24:00`
[1] 0.7022059
$`2014-07-31 17:27:00`
[1] 0.6568266
$`2014-07-31 17:30:00`
[1] 0.8390805
$`2014-07-31 17:33:00`
[1] 0.3838863
$`2014-07-31 17:36:00`
[1] 0.7686567
$`2014-07-31 17:39:00`
[1] 0.6580311
$`2014-07-31 17:42:00`
[1] 0.6923077
$`2014-07-31 17:45:00`
[1] 0.5813008
$`2014-07-31 17:48:00`
[1] 0.9066148
$`2014-07-31 17:51:00`
[1] 0.59375
$`2014-07-31 17:54:00`
[1] 0.9487179
$`2014-07-31 17:57:00`
[1] 0.7325581
我的目标是绘制价值百分比&#34; 1&#34;按时间(每3分钟),但我不知道如何使用奇怪的数据框进行绘图(x_axis是时间,y_axis是百分比)。我是一个新手(= ^ ^ =)
答案 0 :(得分:1)
您现在有一个列表但不是数据帧。 (如果你想要一个数据帧,那么使用ddply
应该已成功。)你可以恢复&#34;次&#34;至少作为names(splitbytime)
的角色向量。原始日期时间值可以重构为:
as.POSIXct( names(splitbytime) )
可以使用以下内容将值提取为数字向量:
unlist( splitbytime )
所以试试:
plot( x= as.POSIXct( names(splitbytime) ),
y= unlist( splitbytime ) )