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

[Carlos J. Gil Bellosta]

Hola, ¿qué tal?

Tu función F está mal definida. Supongo que en el bucle quieres escribir
suma <- suma + ...

Tampoco tengo claro qué devuelve tu función. Si quieres que devuelva suma,
hazlo explícitamente.

Un saludo,

Carlos J. Gil Bellosta
http://www.datanalytics.com



El 23 de octubre de 2015, 12:32, José Miguel Contreras García <
jmcontre...@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@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [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 <
jmcontre...@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@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@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es



[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

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

[Carlos Ortega]

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 <
jmcontre...@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@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@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

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

[Olivier Nuñez]

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"  
Para: "r-help-es"  
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@r-project.org 
https://stat.ethz.ch/mailman/listinfo/r-help-es 


[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es