На этот вопрос уже есть ответ:

У меня есть фрейм данных:

station     date        classification
 1    June - 01/16          A
 2    June - 03/16          B
 1    June - 01/16          A
 7    June - 01/16          C
 1    June - 03/16          A
 2    June - 03/16          B
 2    June - 03/16          B

Я хочу получить общее количество вхождений A, B и C, агрегированных по станции № и дате:

Например, на станции 1 01 июня есть 2 As, а на станции 2 3 июня - 3 Bs.

Я пытался

aggregate(x = list(data_frame$classification), by = list(station=data_frame$station, Date=data_frame$date), function(x) length(unique(x))