如何计算在R中使用nls函数计算的广告费率的置信区间？

Question

我已经使用R的nls函数计算了电视GRP的广告库存率。我使用以下代码进行了计算：

# 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

希望有帮助。