Блог

У меня есть данные с двумя переменными. Каждая переменная имеет 10 наблюдений. Я хотел бы найти 10 интегралов, используя каждое наблюдение в переменной one в качестве нижней границы. Я попробовал код следующим образом,

    data = data.table(x=x, y=y)
     for(i in 1:n){
     F <- ecdf(data$x)
     f <- function(h){((1-F(h))/(1-F(data[i,1])))^((1-0.1*(theta-lambda))/(0.1*lambda))}
    m_i <-data$x[i]
   data$k <- integrate(f, lower=m_i, upper=1)$value
    }
 

Но я получил «Ошибку в статистике :: integrate (…): значение не конечной функции», любая помощь будет оценена. Спасибо за ваше время и внимание.