Utilizando a fórmula final da [resposta do Marcelo Shiniti Uchimura](https://pt.stackoverflow.com/questions/363208/estimar-vari%c3%a1vel-de-dif%c3%adcil-isolamento/363214#363214), 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)`.