[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