Boa tarde colaboradores. Estou com dificuldade de estimar um parâmetro via programação em R para completar minha rotina de processamentos. A equação que estou utilizando é a descrita a seguir, onde tenho todos os valores, menos o valor da variável k
razao2 = (razao1)^((idade2/idade)^k)
onde, tenho os vetores
razao2<-c(1.0331279, 1.0416188, 1.0167619, 0.0000000, 0.9638364, 0.9668588)
razao1<-c(1.0151796, 1.0331279, 1.0416188, 1.0167619, 0.9822768, 0.9638364)
idade<-c(2.580833, 3.594167, 4.608333, 2.580833, 3.594167, 4.608333)
idade2<-c(3.594167, 4.608333, 5.625000, 3.594167, 4.608333, 5.625000)
razao1[1]
é respectivo a razao2[1]
, respectivo a idade[1]
e a idade2[1]
portanto para cada posição i
dos vetores , terei um valor calculado de k
que promove a igualdade
Tentei fazer via seguinte caminho, porém não obtive sucesso, estou com dificuldades de conseguir criar uma lógica q isole k para ser estimado, ou que de algum outro meio promova as estimativas de um valor de k para casa combinações de valores de uma dada posição i
dos vetores
kcalculado <- function(razao1, idade2, idade, k){
razao2 <- razao1^((idade2/idade)^(k))
razao2
}
optimise(kcalculado, razao2)
desde já agradeço a ajuda