A maneira mais simples é usar uma regressão não-linear por mínimos quadrados (função nls
). Matematicamente, este é um método iterativo, que procura valores para os parâmetros para reduzir os resíduos abaixo de um valor limite.
Um exemplo de uso é o seguinte:
set.seed(0)
x <- rnorm(100)
y <- jitter(1234/(x + 1), 0.1)
dat <- data.frame(x = x, y = y)
nl <- nls(y ~ a/(x + b), data = dat, start = c(a = 1000, b = 2))
Os resultados podem ser observados com a função summary
:
> summary(nl)
# Formula: y ~ a/(x + b)
#
# Parameters:
# Estimate Std. Error t value Pr(>|t|)
# a 1132.5404 548.9356 2.063 0.0417 *
# b 1.9850 0.2013 9.863 2.4e-16 ***
# ---
# Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
#
# Residual standard error: 4639 on 98 degrees of freedom
#
# Number of iterations to convergence: 13
# Achieved convergence tolerance: 9.861e-06
Nem sempre a estimativa fica tão "boa", mas isso depende dos dados e da função, não do R.