Использование agregate в R с пользовательской функцией, которая зависит от двух столбцов

Question

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

year Concentration Precipitaion
2000 19.01 23.7251908396947
2000 3.38 25.8842975206612
2000 7.78 12.9504950495049
2000 50.33 22.6153846153846
2000 0.04 19.063829787234
2000 1.4 67.8202247191011
2000 5.11 88.4383561643836
2000 4.53 42.7234042553192
2001 9.57 33
2001 5.25 21.3023255813953
2001 2.33 28.9491525423729
2001 9.29 42.7428571428571
2001 4.01 16.8813559322034
2001 0.39 125.093525179856
2001 1.14 50
2001 6.1 51.0909090909091
2001 1.19 25.0833333333333
2001 4.09 35.921568627451
2001 1.89 127.396226415094
2002 1.28 100.266666666667
2002 5.96 29.5922330097087
2002 2.36 49.0526315789474
2002 5.47 121.756097560976
2002 13.03 53.6978417266187
2002 6.57 23.7575757575758
2002 5.11 74.4375
2002 0.65 29.3592233009709
2002 0.39 180.512195121951
2002 3.35 20.5423728813559
2002 12.92 53.5789473684211
2002 10.01 24.5274725274725
2002 4.66 39.6363636363636
2002 2.25 13.6901408450704
2002 1.31 96.24
2002 1.13 13.1428571428571
2002 5.45 19.8347107438017
2002 6.4 57.375
2002 1.06 186
2002 3.09 59.2142857142857

И у меня есть своя собственная функция, которая зависит от двух столбцов - Концентрация и Осадки

user_function <- function(Concentration, Precipitaion){
  return(mean(rnorm(10000, 
                    mean = mean(Concentration),
                    sd = sd(Precipitaion))))}

Я пытался использовать эту функцию

aggregate(data_sulfat, by=list(data_sulfat$year), 
          FUN = user_function(data_sulfat$Concentration,
                              data_sulfat$Precipitaion))

И у меня ошибка.Скажите пожалуйста, как правильно я использую агрегат с функцией, которая зависит от двух или более столбцов

Answer 1

У вас есть несколько вариантов

d <- structure(list(
    year = c(2000L, 2000L, 2000L, 2000L, 2000L, 2000L, 2000L, 
        2000L, 2001L, 2001L, 2001L, 2001L, 2001L, 2001L, 2001L, 2001L, 2001L, 
        2001L, 2001L, 2002L, 2002L, 2002L, 2002L, 2002L, 2002L, 2002L, 2002L, 
        2002L, 2002L, 2002L, 2002L, 2002L, 2002L, 2002L, 2002L, 2002L, 2002L, 
        2002L, 2002L), 
    Concentration = c(19.01, 3.38, 7.78, 50.33, 0.04, 1.4, 5.11, 4.53, 9.57, 
        5.25, 2.33, 9.29, 4.01, 0.39, 1.14, 6.1, 1.19, 4.09, 1.89, 1.28, 
        5.96, 2.36, 5.47, 13.03, 6.57, 5.11, 0.65, 0.39, 3.35, 12.92, 
        10.01, 4.66, 2.25, 1.31, 1.13, 5.45, 6.4, 1.06, 3.09), 
    Precipitaion = c(23.7251908396947, 25.8842975206612, 12.9504950495049, 
        22.6153846153846, 19.063829787234, 67.8202247191011, 88.4383561643836, 
        42.7234042553192, 33, 21.3023255813953, 28.9491525423729, 
        42.7428571428571, 16.8813559322034, 125.093525179856, 50, 
        51.0909090909091, 25.0833333333333, 35.921568627451, 127.396226415094, 
        100.266666666667, 29.5922330097087, 49.0526315789474, 121.756097560976, 
        53.6978417266187, 23.7575757575758, 74.4375, 29.3592233009709, 
        180.512195121951, 20.5423728813559, 53.5789473684211, 24.5274725274725, 
        39.6363636363636, 13.6901408450704, 96.24, 13.1428571428571, 
        19.8347107438017, 57.375, 186, 59.2142857142857)), 
    class = "data.frame", row.names = c(NA, -39L))

. Простое решение - использовать split и lapply

user_function <- function(d) {
    mean(rnorm(10000, mean = mean(d$Concentration), sd = sd(d$Precipitaion)))
}

lapply(split(d, d$year), user_function)