Совокупность (подсчет) вхождений значений за произвольный период времени

Мы можем прочитать его, используя read.csv, преобразовать первый столбец в дату и время, разделенные на 6-минутные интервалы, и добавить фиктивный столбец из единиц. Затем перечитайте его, используя read.zoo разбиение по типу и агрегирование по фиктивному столбцу:

# test data

Lines <- 'date,type
"Sep 22, 2011 12:54:53.081240000","2"
"Sep 22, 2011 12:54:53.083493000","2"
"Sep 22, 2011 12:54:53.084025000","2"
"Sep 22, 2011 12:54:53.086493000","2"
"Sep 22, 2011 12:54:53.081240000","3"
"Sep 22, 2011 12:54:53.083493000","3"
"Sep 22, 2011 12:54:53.084025000","3"
"Sep 22, 2011 12:54:53.086493000","4"'

library(zoo)
library(chron)

# convert to chron and bin into 6 minute bins using trunc
# Also add a dummy column of 1's 
# and remove any leading space (removing space not needed if there is none)

DF <- read.csv(textConnection(Lines), as.is = TRUE)
fmt <- '%b %d, %Y %H:%M:%S'
DF <- transform(DF, dummy = 1,
         date = trunc(as.chron(sub("^ *", "", date), format = fmt), "00:06:00"))

# split and aggregate

z <- read.zoo(DF, split = 2, aggregate = length)

С приведенными выше тестовыми данными решение выглядит так:

> z
                    2 3 4
(09/22/11 12:54:00) 4 3 1

Обратите внимание, что вышеприведенное было сделано в широкой форме, поскольку эта форма представляет собой временной ряд, а длинная форма — нет. Для каждого типа есть одна колонка. В наших тестовых данных у нас были типы 2, 3 и 4, поэтому есть три столбца.

(Здесь мы использовали chron, так как его метод trunc хорошо сочетается с объединением в 6-минутные группы. chron не поддерживает часовые пояса, что может быть преимуществом, поскольку вы не можете сделать одну из многих возможных ошибок часового пояса, но если вы все равно хотите POSIXct преобразовать его в конце, например, time(z) <- as.POSIXct(paste(as.Date.dates(time(z)), times(time(z)) %% 1)) . Это выражение показано в таблице в одной из статей R News 4/1, за исключением того, что мы использовали as.Date.dates вместо as.Date, чтобы обойти ошибку, которая, кажется, появилась с тех пор. Мы могли бы также используйте time(z) <- as.POSIXct(time(z)), но это приведет к другому часовому поясу.)

РЕДАКТИРОВАТЬ:

Исходное решение было разбито на даты, но впоследствии я заметил, что вы хотите разбить на 6-минутные периоды, поэтому решение было пересмотрено.

РЕДАКТИРОВАТЬ:

Отредактировано на основе комментария.

Вы почти все там. Все, что вам нужно сделать сейчас, это создать зоо-версию этих данных и сопоставить ее с кодомaggregate.zoo. Поскольку вы хотите классифицировать как по времени, так и по типу, ваш второй аргумент дляaggregate.zoo должен быть немного сложнее, и вам нужны подсчеты, а не средства, поэтому вы должны использовать length(). Я не думаю, что count является базовой функцией R или зоопарком, и единственная функция count, которую я вижу в своем рабочем пространстве, исходит от pkg:plyr, поэтому я не знаю, насколько хорошо она будет работать с агрегатом. length работает так, как большинство людей ожидает от векторов, но часто удивляет людей при работе с data.frames. Если вы не получаете то, что хотите, с length, тогда вы должны посмотреть, работает ли вместо этого NROW (и с вашим макетом данных они оба преуспевают): с новым объектом данных необходимо сначала поставить аргумент типа. И получается, что агрегат/зоопарк обрабатывает только классификаторы одной категории, поэтому вам нужно добавить as.vector, чтобы удалить его зоопарковство:

with(cured_data, 
     aggregate(as.vector(x), list(type = type, 
                                   interval=as.factor(time(x) - as.numeric(time(x)) %% 360)),
                             FUN=NROW) 
 )

#  interval            x 
#1 2011-09-22 09:24:00 12
#2 2011-09-22 09:24:00 11

Это пример, измененный из того места, где вы получили код (пример SO WizaRd Dirk): Совокупность (количество) вхождений значений за произвольный период времени

tt <- seq(10, 2000, 10)
x <- zoo(tt, structure(tt, class = c("POSIXt", "POSIXct")))
aggregate(as.vector(x), by=list(cat=as.factor(x), 
     tms = as.factor(index(x) - as.numeric(index(x)) %% 600)), length)

   cat                 tms  x
1    1 1969-12-31 19:00:00 26
2    2 1969-12-31 19:00:00 22
3    3 1969-12-31 19:00:00 11
4    1 1969-12-31 19:10:00 17
5    2 1969-12-31 19:10:00 28
6    3 1969-12-31 19:10:00 15
7    1 1969-12-31 19:20:00 17
8    2 1969-12-31 19:20:00 16
9    3 1969-12-31 19:20:00 27
10   1 1969-12-31 19:30:00  8
11   2 1969-12-31 19:30:00  4
12   3 1969-12-31 19:30:00  9