Utilizando a 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)
.