Sorprendentemente, no pude encontrar una respuesta a la siguiente pregunta usando Google:
Tengo algunos datos biológicos de varias personas que muestran un comportamiento de crecimiento más o menos sigmoide en el tiempo. Por lo tanto, deseo modelarlo utilizando un crecimiento logístico estándar
P(t) = k*p0*exp(r*t) / (k+p0*(exp(r*t)-1))
siendo p0 el valor inicial en t = 0, k es el límite asintótico en t-> infinito yr es la velocidad de crecimiento. Hasta donde puedo ver, puedo modelar esto fácilmente usando nls (falta de comprensión de mi parte: ¿por qué no puedo modelar algo similar usando la regresión logit estándar al escalar el tiempo y los datos? EDITAR: Gracias Nick, aparentemente la gente lo hace, por ejemplo, por proporciones, pero rara vez http://www.stata-journal.com/article.html?article=st0147 . La siguiente pregunta sobre esta tangente sería si el modelo puede manejar valores atípicos> 1).
Ahora deseo permitir algunos efectos fijos (principalmente categóricos) y algunos aleatorios (una identificación individual y posiblemente también una identificación de estudio) sobre los tres parámetros k, p0 y r. ¿Es nlme la mejor manera de hacer esto? El modelo SSlogis parece sensato para lo que estoy tratando de hacer, ¿es correcto? ¿Es cualquiera de los siguientes un modelo sensato para empezar? Parece que no puedo obtener los valores iniciales correctamente y actualizar () solo parece funcionar para efectos aleatorios, no fijos, ¿alguna sugerencia?
nlme(y ~ k*p0*exp(r*t) / (k+p0*(exp(r*t)-1)), ## not working at all (bad numerical properties?)
data = data,
fixed = k + p0 + r ~ var1 + var2,
random = k + p0 + r ~ 1|UID,
start = c(p0=1, k=100, r=1))
nlme(y ~ SSlogis(t, Asym, xmid, scal), ## not working, as start= is inappropriate
data = data,
fixed = Asym + xmid + scal ~ var1 + var2, ## works fine with ~ 1
random = Asym + xmid + scal ~ 1|UID,
start = getInitial(y ~ SSlogis(Dauer, Asym, xmid, scal), data = data))
Como soy nuevo en los modelos mixtos no lineales en particular y en los modelos no lineales en general, agradecería algunas recomendaciones de lectura o enlaces a tutoriales / preguntas frecuentes con preguntas para novatos.