Skip to main content
2 de 2
Texto corrigido, tornado mais exato.
Rui Barradas
  • 19mil
  • 2
  • 14
  • 28

Utilizando a fórmula anterior à fórmula final da resposta do Marcelo Shiniti Uchimura, podemos ajustar um modelo linear.

#ln(ln(razao2) = ln(ln(razao1)) + ln(idade2/idade1)*k

log.log.razao1 <- log(log(razao1))
log.log.razao2 <- log(log(razao2))
log.id2.id <- log(idade2/idade)

modelo <- lm(log.log.razao2 ~ log.id2.id)

r1 <- exp(exp(coef(modelo)[1]))
k <- exp(coef(modelo)[2])

Neste resultado, r1 é o valor calculado de razao1, que não está muito afastado do valor médio dos dados.

mean(razao1)
#[1] 1.0088

r1
#(Intercept) 
#   1.009239

E o valor de k será o seguinte.

k
#log.id2.id 
#  73.33412 

Se estes nomes forem incomodativos, pode-se fazer unname(r1) e unname(k).

Rui Barradas
  • 19mil
  • 2
  • 14
  • 28