[R] NLS: starting values

william Milton Karvajal Herradora wmkarherr at hotmail.com
Fri Dec 12 21:29:15 CET 2003


   Salud a los de R-proyect.

   Os envio, quizas un aporte, para el problema de valores iniciales  en
   regresion no lineal.

   En el Modelo Xenot  Y = m + exp(b0+b1*X)
   Quasi linealizando el modelo no lineal se transforma la variable
   dependiente acotando el parametro  m,  de donde
   (*)  log(Y-m) = b0+b1*X.    Aqui Y-m>0 solo si m<Y.
   Luego un valor inicial para m es m0 =min(Y)+eps, estimamos eps como el
   error standard: asi eps=sd(Y)/sqrt(n),  donde n es el tamaño muestral,
   sd(Y) la desviacion standar de Y.
   Luego de estimado m0, ajustamos el modelo lineal (*) para obtener
   aproximaciones iniciales a b0 y b1.

   Para el Modelo Logistico Y = M/(1+exp(b0+b1*X))
   despejamos  (*) log(M/Y-1) = b0+b1*X.           Aqui  M/Y-1>0 solo si
   M>Y.
   Asi un valor inicial a M es M0= max(Y) +eps, ...

   Adjunto documento con comandos de R para ejemplos de dietas, no
   comprimidos.
   Hasta otra...                 W. Milton K. H.
     _________________________________________________________________

   MSN Amor [1]Busca tu ½ naranja

References

   1. http://g.msn.com/8HMAES/2740??PS=


More information about the R-help mailing list