[R-es] Resolución de ecuación

jmcontreras en ugr.es jmcontreras en ugr.es
Vie Oct 23 12:33:16 CEST 2015


Hola a todos

Tengo una duda existencial, como siempre, jejeje

Quiero resolver una ecuación (F=0) que depende de un parámetro z, pero 
este viene en función de una suma que depende de un vector de tiempos ti 
(1:25) y unos valores de un vector x

Lo he intentado resolver mediante una función para luego utilizar 
solve, pero me he atrancado.

La función es esta:

ti<-1:25
x<-sample(1:10, size=25, replace=T)
beta<-0.1

F<-function(z) {
suma<-0
for(i in 2:n){
sum<-sum+exp(1)^(-beta*ti[i])/(z-beta)^2*(log(x[i])-exp(1)^(-beta)*log(x[i-1]))*(exp(1)^(ti[i]*(beta-z))*(1-ti[i]*(beta-z))-exp(1)^(-ti[i-1]*(beta-z))*(1-ti[i-1]*(beta-z)))
}
}

Quiero calcular el valor de z para que la suma valga cero. ¿He metido 
mucho la pata?

Gracias de antemano

Saludos



Más información sobre la lista de distribución R-help-es