Как интегрировать, используя столбец в качестве верхней границы

Question

Я пытаюсь интегрировать функцию (y = a (1-d ^ -b) в R, где верхняя граница берется из столбца (Lay1), а затем прикрепить результаты интегрирования к фрейму данных. IЯ пытался сделать это с помощью функции mapply, но не могу заставить ее работать на 1 столбец, и при этом я не могу прикрепить ее к фрейму данных.

Код:

    SiteA<-read.csv("TransSiteA.csv",header=TRUE)
    output <- data.frame(SiteA)

    FunSiteA <-function(d){a*(1-(d^-b))}

    a <- -1.1441E-006   
    b <- 1.4346 

    KLay1 <- mapply(integrate, lower=SiteA$Lay1, upper=0.05, MoreArgs=list(f=FunSiteA))
    KLay1

    cbind(output, KLay1)

Пример данных:

ID  Date                      Lay1      Lay2    Lay3    Lay4
231.00  2011-04-10 9:30     0.020688    0.05    0.15    0.25
232.00  2011-04-10 9:45     0.019692    0.05    0.15    0.25
233.00  2011-04-10 10:00    0.019692    0.05    0.15    0.25
234.00  2011-04-10 10:15    0.021684    0.05    0.15    0.25
235.00  2011-04-10 10:30    0.019692    0.05    0.15    0.25
236.00  2011-04-10 10:45    0.019692    0.05    0.15    0.25
237.00  2011-04-10 11:00    0.019692    0.05    0.15    0.25
238.00  2011-04-10 11:15    0.019692    0.05    0.15    0.25
239.00  2011-04-10 11:30    0.021684    0.05    0.15    0.25
240.00  2011-04-10 11:45    0.02268     0.05    0.15    0.25
241.00  2011-04-10 12:00    0.023676    0.05    0.15    0.25
242.00  2011-04-10 12:15    0.016704    0.05    0.15    0.25
243.00  2011-04-10 12:30    0.0177      0.05    0.15    0.25
244.00  2011-04-10 12:45    0.018696    0.05    0.15    0.25

Answer 1

Как это?Ваше mapply заявление сработало потрясающе для меня, изменило его немного.SiteA теперь называется data кстати ..

mapply(integrate, lower= data$Lay1, upper=0.05, MoreArgs=list(f= FunSiteA), SIMPLIFY = FALSE) %>%
  lapply(`[`, "value") %>%
  unlist() %>%
  cbind(data, "integrated" = .)

               ID  Date     Lay1 Lay2 Lay3 Lay4   integrated
231.00 2011-04-10  9:30 0.020688 0.05 0.15 0.25 4.490543e-06
232.00 2011-04-10  9:45 0.019692 0.05 0.15 0.25 4.797246e-06
233.00 2011-04-10 10:00 0.019692 0.05 0.15 0.25 4.797246e-06
234.00 2011-04-10 10:15 0.021684 0.05 0.15 0.25 4.204397e-06
235.00 2011-04-10 10:30 0.019692 0.05 0.15 0.25 4.797246e-06
236.00 2011-04-10 10:45 0.019692 0.05 0.15 0.25 4.797246e-06
237.00 2011-04-10 11:00 0.019692 0.05 0.15 0.25 4.797246e-06
238.00 2011-04-10 11:15 0.019692 0.05 0.15 0.25 4.797246e-06
239.00 2011-04-10 11:30 0.021684 0.05 0.15 0.25 4.204397e-06
240.00 2011-04-10 11:45 0.022680 0.05 0.15 0.25 3.936582e-06
241.00 2011-04-10 12:00 0.023676 0.05 0.15 0.25 3.685198e-06
242.00 2011-04-10 12:15 0.016704 0.05 0.15 0.25 5.869608e-06
243.00 2011-04-10 12:30 0.017700 0.05 0.15 0.25 5.483335e-06
244.00 2011-04-10 12:45 0.018696 0.05 0.15 0.25 5.127134e-06