Se ha borrado un adjunto en formato HTML... URL: <https://stat.ethz.ch/pipermail/r-help-es/attachments/20151023/300f7387/attachment.html>
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]]
Hola, Otra alternativia que tienes, por la complejidad de la ecuación que la tienes expresada como un sumatorio (sin desarrollar) es utilizar "R" unido al paquete de matemática simbólica "Yacas", mediante "RYacas": https://cran.r-project.org/web/packages/Ryacas/index.html Saludos, Carlos Ortega www.qualityexcellence.es El 23 de octubre de 2015, 10:58, José Miguel Contreras García < jmcontreras en ugr.es> escribió:> 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 >-- Saludos, Carlos Ortega www.qualityexcellence.es [[alternative HTML version deleted]]
Estimados Posiblemente algo de lo expuesto en https://cran.r-project.org/web/views/Optimization.html pueda ser útil, aunque esa ecuación escapa a mi experiencia. Javier Rubén Marcuzzi Técnico en Industrias Lácteas Veterinario De: Carlos Ortega Enviado: viernes, 23 de octubre de 2015 8:52 Para: José Miguel Contreras García CC: r-help-es Asunto: Re: [R-es] resolución de ecuación Hola, Otra alternativia que tienes, por la complejidad de la ecuación que la tienes expresada como un sumatorio (sin desarrollar) es utilizar "R" unido al paquete de matemática simbólica "Yacas", mediante "RYacas": https://cran.r-project.org/web/packages/Ryacas/index.html Saludos, Carlos Ortega www.qualityexcellence.es El 23 de octubre de 2015, 10:58, José Miguel Contreras García < jmcontreras en ugr.es> escribió:> 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 >-- Saludos, Carlos Ortega www.qualityexcellence.es [[alternative HTML version deleted]] _______________________________________________ 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]]