[R-es] resolución de ecuación
Olivier Nuñez
onunez en unex.es
Vie Oct 23 13:16:26 CEST 2015
José,
varios comentarios:
1) exp(1)^a =exp(a)
2) solve es para ecuaciones lineales. La tuya es nolineal
3) ¿cuales el rango de valores posibles de z?
4) Prueba la función uniroot.all del paquete "rootSolve".
Un saludo. Olivier
----- Mensaje original -----
De: "José Miguel Contreras García" <jmcontreras en ugr.es>
Para: "r-help-es" <r-help-es en r-project.org>
Enviados: Viernes, 23 de Octubre 2015 10:58:06
Asunto: [R-es] resolución de ecuación
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
_______________________________________________
R-help-es mailing list
R-help-es en r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es
[[alternative HTML version deleted]]
Más información sobre la lista de distribución R-help-es