Я хотел бы вычислить интеграл из двух столбцов в моей базе данных, когда интеграл основан на функции; 1 / x, мне нужно, чтобы я написал цикл. Можете ли вы помочь в написании цикла и функции интеграции?
Это примерный кадр данных;
upper_concentration<-c(1:200, 1)
lower_concentration<-upper_concentration*0.9
df = data.frame(upper_concentration,lower_concentration)for (i in 1:(length(df))){
integral <- function (x) {1/df$upper_concentration}
result <- integrate(integral, lower = df$upper_concentration, upper =upper_concentration*0.9)
}
Встроенная функция
integrate
выполняет интеграл для полного диапазона.Нет, если вы хотите, чтобы интеграл для каждой части кадра данных то необходим цикл.
result <- vector("list", nrow(df))