如何计算在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

希望对您有所帮助。