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

[Javier Rubén Marcuzzi]

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]]