Dolor de cabeza de predicción de estabilidadHola,
Tengo un problema de análisis sobre el que me gustaría recibir algún consejo. A continuación se muestra una tabla con algunos datos de estabilidad reales, de los cuales intento extraer una predicción de estabilidad razonable.
Normalmente uso el generador de perfiles de degradación para esto, y hasta ahora lo he hecho sin ningún problema en muchos otros proyectos, pero este conjunto de datos me está causando dolor de cabeza; Creo que es porque tuvimos un T0 anormalmente bajo que sesgó la tendencia haciendo las predicciones... Optimistas.
Tabla vinculada: con los límites de aceptación asociados que tenemos.
Notas:
1-Este es un ejemplo del mundo real, por lo que es un sistema completamente comprendido.
2-Sabemos, por tanto, que es imposible que este material adquiera actividad. es decir, la tendencia alcista se debe al ruido del ensayo.
3-También sabemos por datos y experimentos anteriores que este material se degrada de forma aproximadamente lineal, y que las caducidades (es decir, cuando los niveles de confianza cruzan el límite inferior) en la región de 730 a 1095 días no son inusuales.
Mi pregunta entonces es cómo analizar esto. Me sentiría cómodo usando algo como un intervalo de confianza más amplio y expandiéndolo como predicción, pero no sabría cómo hacerlo.
El próximo punto de recopilación de datos es dentro de 3 meses, por lo que, si bien estoy seguro de que la tendencia "real" surgirá con el tiempo y el ruido inicial del ensayo se promediará, necesito hacer algún tipo de predicción realista ahora. ¡Cualquier sugerencia/orientación (a pesar de regañar al científico que hizo el T0 
) sería genial!
jon
Análisis y modelado de datos básicosModelado predictivo y aprendizaje automáticoIngeniería de Calidad y ProcesosRe: dolor de cabeza de predicción de estabilidadHola @Jon_Armer : Algunas observaciones/pensamientos/preguntas... sin ningún orden en particular.
1. ¿Puede la evaporación influir aquí? Veo que las unidades de respuesta son mg/mL. ¿Podría la evaporación hacer que parezca que la tendencia está aumentando?
2. T5 y la segunda repetición de T32 también son bajas...
3. ¿Por qué "necesitas hacer algún tipo de predicción realista ahora"? ¿Está (1) intentando obtener una lectura temprana de este lote en relación con las especificaciones, o (2) está intentando obtener una predicción precisa? Si (1), es posible que pueda decir simplemente que, según estos datos , no hay evidencia de que la especificación esté en peligro. Si (2)... bueno, según sus suposiciones (de sus "Notas"), ese es un animal diferente.
4. ¿Podría ser que haya algo diferente en este lote que haga que se comporte de manera diferente a los "datos anteriores"?
Re: dolor de cabeza de predicción de estabilidadHola @MRB3855
Buena pregunta:
1-no, están sellados para evitar esto mismo
2-De acuerdo. Con los sistemas biológicos no es raro ver cierto grado de variación en los estudios longitudinales, especialmente teniendo en cuenta los ensayos, etc., pero estoy de acuerdo en que la variación es amplia.
3-Predicción de la vida útil de este lote. Necesitamos hacer una evaluación al respecto y poner una 'línea en la arena' para que podamos completar nuestro cronograma anual.
4-Sí, potencialmente, es un sistema complejo, por lo que existe un margen de variación que sería difícil de capturar/identificar.
Re: dolor de cabeza de predicción de estabilidadHola @Jon_Armer : Mmmm. Dice que su próximo punto de atracción es a los 3 meses. ¿Cuál es su vida útil objetivo (¿registrada?)?
Re: dolor de cabeza de predicción de estabilidadAh, ligera corrección (disculpas si engañé): el próximo momento es después de 3 meses adicionales ; que en realidad será el punto de tiempo de 6 meses.
El lote anterior tenía una estabilidad prevista de 18 meses a partir de datos de 8 meses (confirmados a los 12 meses).
Re: dolor de cabeza de predicción de estabilidadA menudo pienso que el término "simple" se utiliza en estos contextos como una forma de no asustar a los tímidos estudiantes de doctorado para que no se dediquen a la estadística....
OK entonces
1.Correcto. Este es un enfoque bastante típico utilizado en la industria; los primeros puntos de tiempo están muy juntos y luego se van dispersando cada vez más a medida que avanzas.
Ahora bien, a menudo se realizaría un estudio de alcance para identificar cualquier punto de inflexión en las tendencias y luego los puntos de tiempo del estudio se diseñarían de manera que aumentaran la densidad de datos alrededor de ese punto. A veces, sin embargo, especialmente al comienzo de un nuevo proyecto o ciclo de desarrollo, no tienes eso, por lo que estás un poco ciego (y tiendes a recurrir a enfoques "probados y probados").
Desafortunadamente, estamos "atascados" con la propagación aquí, al menos para este lote.
2. Los lotes anteriores tuvieron una caída generalmente lineal, con una tendencia que no se hizo evidente hasta más tarde.
3.Sí, esa es una opción, aunque no resuelve mi problema inmediato (de tener que asignar una fecha de vencimiento preliminar para este lote)
4-Creo que esta es la forma en que me inclino. Ese tipo de caída sería consistente con lotes anteriores (generalmente es bastante lineal), solo lleva tiempo.
5, sí de nuevo, aunque solo tengo otro lote en este 'conjunto'; la tendencia fue similar pero 'ligeramente a la baja', lo que me permitió estimar un cruce de nivel inferior.
6-Estamos de acuerdo mucho aquí, lo cual siempre es bueno (¡en Internet precisamente!). Lo que creo que estamos viendo aquí es una variación del ensayo superpuesta al material (relativamente) estable. Los dioses de la variación se alinearon con el lote anterior para presentarme una tendencia ligeramente negativa que permitió una predicción de estabilidad relativamente sencilla (con las obvias advertencias habituales). Esta vez, de alguna manera he enojado a los dioses de la variabilidad y no se ofrece tal tendencia, lo que me provoca dolor de cabeza.
El ruido en estos datos también es más amplio que en el lote anterior.
Entonces, aquí es donde está mi cabeza actualmente:
Tengo datos adicionales limitados a los que recurrir aquí: el lote anterior (y el único otro) dio una tendencia ligeramente negativa, lo que me permite asignar una estabilidad inicial de 18 meses en este punto (t3m). Obviamente somos conscientes de los peligros de una extrapolación excesiva en este caso, pero a veces deben intervenir consideraciones prácticas.
Esta extensión en este lote no me permite hacer eso. Por lo tanto, la única opción que actualmente veo que sería defendible sería volver a la estabilidad de los lotes anteriores sobre la base de que este nuevo lote, hasta el momento, no me ha ofrecido ningún indicio de que sea peor . Como tal, se podría argumentar que el uso de la estabilidad de los lotes anteriores es el peor de los casos y, por lo tanto, es seguro/defendible.
Luego superviso y actualizo según corresponda a medida que llegan nuevos datos. A menos que haya alguna forma inteligente de utilizar la variabilidad agregada en estos datos para ampliar los intervalos de confianza (y así "forzar" la predicción a cruzar el nivel inferior)
¿Pensamientos?
Re: dolor de cabeza de predicción de estabilidadHola @Jon_Armer : Creo que tu pensamiento tiene sentido aquí; Según estos datos, no hay evidencia de que la especificación esté en peligro o de que sea peor que el lote anterior... y continúe monitoreando.
Y, si es cierto que los lotes anteriores tuvieron una caída generalmente lineal, con una tendencia que no se hizo evidente hasta más tarde... entonces me inclinaría, si es posible, a tirar más de aquí a el siguiente punto de extracción programado para que se pueda estimar el "nudo" (el momento en que la tendencia comienza a bajar). Luego podría modelar toda la "curva" introduciendo otro término en su modelo:
max(0, timepoint - t)
Donde t es el momento en el que comienza la tendencia a la baja. Entonces, tendrías Y = b0 + b1*Punto de tiempo + b2*max(0, Punto de tiempo - t). La pendiente inicial es b1, la segunda pendiente es b1+b2. Puede usar la plataforma no lineal para estimar t (y b0, b1 y b2), o puede observar t y usar Ajustar modelo para estimar b0, b1 y b2. Estrictamente hablando, la plataforma no lineal es lo correcto... aunque, a simple vista, puede ser adecuada para su propósito.
Re: dolor de cabeza de predicción de estabilidadGenial, ¡gracias por la discusión/consejo!
