Как рассчитать доверительный интервал нормы расхода, рассчитанной с использованием функции nls в R?

Question

Я вычислил норму запасов для TV GRP, используя функцию nls R. Я использовал приведенный ниже код для ее вычисления:

# Adstock function
adstock <- function(inp_x, rate=0){
  return(as.numeric(stats::filter(x= inp_x, filter=rate, method="recursive")))
}


rate_TV = nls((Total_Traffic)~b0+b1*adstock(TV_GRP, rate), data = loans, 
               start=c(b0=0, b1=1,rate=0.1), trace = TRUE, nls.control(maxiter = 100))
summary(rate_TV)

Здесь переменная y - Total_Traffic, а переменная x - TV_GRP.Теперь я хочу рассчитать доверительный интервал ставки.Есть ли способ, которым я могу рассчитать это.Пожалуйста, помогите.

Answer 1

было бы здорово иметь воспроизводимый пример, но функция confint должна сработать, как указано в документе nls.Ниже с lm

x <- c(1:100)
y <- rnorm(100)
mod <- lm(y~x)
confint(mod)

                   2.5 %      97.5 %
(Intercept) -0.759751731 -0.03777529
x           -0.002136328  0.01027562

Надеюсь, это поможет.