R help es - Oct 2015 - 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
------------ próxima parte ------------
A non-text attachment was scrubbed...
Name: smime.p7s
Type: application/pkcs7-signature
Size: 1089 bytes
Desc: no disponible
URL:
<https://stat.ethz.ch/pipermail/r-help-es/attachments/20151023/c68955ed/attachment.bin>

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 <
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
>
	[[alternative HTML version deleted]]