Итак, изначально у меня был следующий объект:

> head(gs)
  year disturbance lek_id  complex tot_male
1 2006           N     3T  Diamond        3
2 2007           N     3T  Diamond       17
3 1981           N   bare 3corners        4
4 1982           N   bare 3corners        7
5 1983           N   bare 3corners        2
6 1985           N   bare 3corners        5

После этого я вычислил общую статистику min, max, mean и sd of tot_male для года в комплексе. Я использовал функции разделения данных R и назначил логические имена столбцов там, где это казалось подходящим, и в конечном итоге сделал их различными объектами.

> tyc_min = aggregate(gs$tot_male, by=list(gs$year, gs$complex), FUN=min)
> names(tyc_min) = c("year", "complex", "tot_male_min")
> tyc_max = aggregate(gs$tot_male, by=list(gs$year, gs$complex), FUN=max)
> names(tyc_max) = c("year", "complex", "tot_male_max")
> tyc_mean = aggregate(gs$tot_male, by=list(gs$year, gs$complex), FUN=mean)
> names(tyc_mean) = c("year", "complex", "tot_male_mean")
> tyc_sd = aggregate(gs$tot_male, by=list(gs$year, gs$complex), FUN=sd)
> names(tyc_sd) = c("year", "complex", "tot_male_sd")

Пример вывода (2-й объект - Tyc_max):

year  complex tot_male_max
1 2003                     0
2 1970 3corners           26
3 1971 3corners           22
4 1972 3corners           26
5 1973 3corners           32
6 1974 3corners           18

Теперь мне нужно добавить количество образцов в год / сложную комбинацию. Затем мне нужно объединить их в один объект данных и экспортировать в виде файла .csv

Я знаю, что мне нужно использовать функцию merge () вместе с all.y, но я не знаю, как обработать эту ошибку:

Error in fix.by(by.x, x) : 
  'by' must specify one or more columns as numbers, names or logical

Или .. добавить количество образцов в год и сложный. Какие-либо предложения?