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=