DMFitwebedition - Blog de ESPOL
DMFitwebedition - Blog de ESPOL
DMFitwebedition - Blog de ESPOL
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
DMFit web edition<br />
Cabeza-Herrera, E.A.<br />
Departamento <strong>de</strong> Microbiología, Facultad <strong>de</strong> Ciencias Básicas, Universidad <strong>de</strong><br />
Pamplona. Campus Universitario, Km. 1, vía a Bucaramanga. Pamplona, Colombia.<br />
1. Introducción<br />
enalcahe@unipamplona.edu.co<br />
La edición web <strong>de</strong>l DMFit es una aplicación on-line aplicada para ajustar curvas bacterianas<br />
don<strong>de</strong> la fase lineal (fase exponencial) está precedida <strong>de</strong> una fase <strong>de</strong> adaptación y una posterior<br />
estacionaria. La edición web <strong>de</strong>l DMFit ha sido elaborada con financiamiento <strong>de</strong>l UK Food<br />
Standards Agency (Agencia <strong>de</strong> Normas Alimentarias <strong>de</strong>l Reino Unido).<br />
La versión <strong>de</strong> escritorio <strong>de</strong>l DMFit (<strong>de</strong>scargable <strong>de</strong> http://www.ifr.ac.uk/safety/DMFit) es parte <strong>de</strong>l<br />
sistema utilizado en el Institute of Food Research (Instituto <strong>de</strong> Investigación Alimentaria) para<br />
mo<strong>de</strong>lar la variación <strong>de</strong>l logaritmo <strong>de</strong>l número <strong>de</strong> células bacterianas <strong>de</strong> un cultivo en función <strong>de</strong>l<br />
tiempo. (DM: mo<strong>de</strong>lización dinámica). Los mo<strong>de</strong>los <strong>de</strong>l ComBase predictor, al cual pue<strong>de</strong><br />
acce<strong>de</strong>rse a través <strong>de</strong>l sitio web ComBase (http://www.combase.cc) fueron <strong>de</strong>sarrollados<br />
utilizando DMFit.<br />
Esta edición <strong>de</strong>l DMFit permite al usuario:<br />
• Ver una representación gráfica <strong>de</strong> los datos <strong>de</strong> crecimiento/supervivencia microbiana.<br />
• Ajustar los datos un mo<strong>de</strong>lo <strong>de</strong> crecimiento/supervivencia para obtener estimaciones <strong>de</strong><br />
parámetros cinéticos como:<br />
o Velocidad Máxima <strong>de</strong> crecimiento/muerte<br />
o Fase <strong>de</strong> latencia (o interfase)<br />
o Número inicial <strong>de</strong> células<br />
o Número final <strong>de</strong> células<br />
o Estimación <strong>de</strong> errores típicos o estándar <strong>de</strong> estos parámetros.<br />
2. Descripción <strong>de</strong> la interfaz <strong>de</strong>l usuario<br />
La pantalla se compone <strong>de</strong> cuatro áreas fundamentales (ver figura 1):<br />
1. Un espacio <strong>de</strong>dicado al ingreso <strong>de</strong>l conjunto <strong>de</strong> datos, el cual se ubica en la parte<br />
superior izquierda <strong>de</strong> la pantalla.<br />
2. El área <strong>de</strong> selección por parte <strong>de</strong>l usuario <strong>de</strong>l tipo <strong>de</strong> mo<strong>de</strong>lo que se utilizará, esta área<br />
se ubica en la parte superior <strong>de</strong>recha <strong>de</strong> la pantalla.<br />
3. A continuación en la tercera área ubicada en el recuadro inferior <strong>de</strong>recho se muestra los<br />
'Resultados' con la estimación <strong>de</strong> los parámetros <strong>de</strong>l mo<strong>de</strong>lo y la bondad <strong>de</strong>l ajuste.<br />
4. Finalmente, en la parte inferior izquierda <strong>de</strong> la pantalla pue<strong>de</strong> visualizarse <strong>de</strong> forma<br />
gráfica los resultados.
2.1. Cargando datos<br />
El tiempo vs conteo logarítmico –log10– (o conteo LN) <strong>de</strong>be ser cargado en la "entrada <strong>de</strong> datos"<br />
<strong>de</strong> texto, en el lado izquierdo <strong>de</strong> la pantalla. El conjunto <strong>de</strong> datos <strong>de</strong>be ser escrito directamente<br />
en el cuadro <strong>de</strong> texto como sigue: Valor <strong>de</strong> tiempo, espacio, valor log <strong>de</strong>l conteo bacteriano,<br />
enter, (repetir para cada línea). Alternativamente, los datos pue<strong>de</strong>n ser copiados <strong>de</strong>s<strong>de</strong> otra<br />
aplicación (por ejemplo, hoja <strong>de</strong> cálculo <strong>de</strong> Excel, fila <strong>de</strong> texto) y pegar directamente en la<br />
'entrada <strong>de</strong> datos "<strong>de</strong> texto.<br />
Figura 1. Interfaz principal (pantalla) <strong>de</strong>l DMFit, con la ubicación <strong>de</strong> las áreas principales.<br />
1 2<br />
3 4<br />
Cuando las concentraciones <strong>de</strong> bacterias son más bajas que el umbral <strong>de</strong> <strong>de</strong>tección, pue<strong>de</strong> ser<br />
usado "-0,01" como código para "No Detectado". Para cargar este punto <strong>de</strong> datos, simplemente<br />
escriba: el tiempo <strong>de</strong> medición, espacio, -0,01, enter. Aunque se muestra en el gráfico, este<br />
punto no se utiliza en el ajuste.
2.2. ¿Cómo introducir los datos el conteo vs tiempo?<br />
Los datos <strong>de</strong>ben ser insertados en el espacio “entrada <strong>de</strong> datos” (data input – área 1) en la forma<br />
como se muestra a continuación:<br />
0 -0.01<br />
1 3.5<br />
2 4.1<br />
4 5.95<br />
5 6.35<br />
8 8.21<br />
10 9.1<br />
11 8.6<br />
• Solo se aceptan caracteres numéricos (0-9), tabulación, espacio, enter, “-“ están<br />
permitidos en el cuadro <strong>de</strong> texto.<br />
• Para cada grupo <strong>de</strong> datos, el tiempo <strong>de</strong>be ser ingresado antes que el conteo log<br />
bacteriano.<br />
• Debe usarse el símbolo “punto” como separador <strong>de</strong> <strong>de</strong>cimales.<br />
• Tenga en cuenta que “-0.01” sólo <strong>de</strong>be utilizarse como un código para N/D (no<br />
<strong>de</strong>tectado).<br />
2.3. Mediciones por <strong>de</strong>bajo <strong>de</strong>l umbral <strong>de</strong> <strong>de</strong>tección<br />
Cuando la concentración <strong>de</strong> microorganismos es menor al nivel <strong>de</strong> <strong>de</strong>tección, “-0.01” pue<strong>de</strong> ser<br />
usado como un código para N/D (No <strong>de</strong>tectado).<br />
Tenga en cuenta que “-0.01” se utilizará tan sólo como un código, y aunque estos datos se<br />
muestren en el gráfico, no se utilizan durante el ajuste.<br />
2.4. Otros requisitos<br />
• Los recuentos bacterianos <strong>de</strong>ben ser ingresados como logaritmo, ya sea log10 ufc/ml o<br />
ln ufc/ml (o log10 ufc/g o ln ufc/g).<br />
• Debe haber un mínimo <strong>de</strong> 2 datos y máximo 100 datos (log recuento vs tiempo) para<br />
cada grupo <strong>de</strong> datos.<br />
3. Mo<strong>de</strong>los<br />
3.1. Mo<strong>de</strong>los disponibles y la importancia <strong>de</strong> los parámetros<br />
Una <strong>de</strong>scripción <strong>de</strong> los parámetros que <strong>de</strong>finen los mo<strong>de</strong>los disponibles en DMFit se<br />
presenta a continuación:<br />
• Valor inicial: Describe el logaritmo inicial <strong>de</strong> la <strong>de</strong>nsidad <strong>de</strong> células bacterianas (en log10<br />
ufc/g ó Ln ufc/g).<br />
• Fase Lag/interfase: El tiempo <strong>de</strong> latencia (o interfase en el caso <strong>de</strong> una curva <strong>de</strong><br />
supervivencia). El tiempo <strong>de</strong> latencia se <strong>de</strong>fine generalmente como la intersección entre
la tangente <strong>de</strong> la fase <strong>de</strong> crecimiento exponencial y el valor inicial. Las unida<strong>de</strong>s son las<br />
mismas que las unida<strong>de</strong>s utilizadas por los datos <strong>de</strong> tiempo.<br />
• Velocidad o Tasa máxima: Hace relación a la máxima velocidad o tasa <strong>de</strong> crecimiento<br />
(o la máxima tasa <strong>de</strong> mortalidad en caso <strong>de</strong> una curva <strong>de</strong> supervivencia). Las unida<strong>de</strong>s<br />
<strong>de</strong> las estimaciones están relacionada con las unida<strong>de</strong>s <strong>de</strong> los datos. Por ejemplo, si el<br />
conteo bacteriano está expresado en log10 ufc/g y el tiempo en horas, la máxima<br />
velocidad <strong>de</strong> crecimiento/muerte es log10 ufc/g/h.<br />
• Valor final: Se refiere al logaritmo <strong>de</strong> la <strong>de</strong>nsidad celular final, en las mismas unida<strong>de</strong>s<br />
que los datos.<br />
Las curvas bacterianas pue<strong>de</strong>n ser ajustadas a dos tipos diferentes <strong>de</strong> mo<strong>de</strong>los (o <strong>de</strong> forma<br />
parcial a estos mo<strong>de</strong>los):<br />
• El mo<strong>de</strong>lo <strong>de</strong> Baranyi y Roberts<br />
• El mo<strong>de</strong>lo trilineal, mo<strong>de</strong>los bifásicos y mo<strong>de</strong>los lineales<br />
3.2. Mo<strong>de</strong>lo <strong>de</strong> Baranyi y Roberts<br />
El mo<strong>de</strong>lo <strong>de</strong> Baranyi y Roberts (1994) <strong>de</strong>scribe una curva bacteriana sigmoi<strong>de</strong>a. La diferencia<br />
principal entre este mo<strong>de</strong>lo y otras curvas sigmoidales como la función modificada <strong>de</strong> Gompertz,<br />
el mo<strong>de</strong>lo logístico, etc., es que las interfases son mas lineales que los <strong>de</strong> aquellas curvas<br />
sigmoi<strong>de</strong>as clásicas, las cuales tienen una mayor pronunciación en la curvatura <strong>de</strong> estas<br />
interfases. El mo<strong>de</strong>lo <strong>de</strong> Baranyi y Roberts tiene 4 parámetros principales (valor inicial, fase <strong>de</strong><br />
latencia o retraso/interfase, tasa máxima, valor final) y 2 parámetros <strong>de</strong> curvatura: mCurv y<br />
nCurv, los cuales <strong>de</strong>scriben la curvatura <strong>de</strong> la curva sigmoi<strong>de</strong>a respectivamente al principio y al<br />
final <strong>de</strong> la fase <strong>de</strong> crecimiento. En la versión <strong>de</strong> escritorio <strong>de</strong> DMFit, estos parámetros, aunque no<br />
estén optimizados, se pue<strong>de</strong>n omitir uno <strong>de</strong> ellos o ambos, por el programa si así se requiere. En<br />
esta versión, los valores <strong>de</strong> mCurv y nCurv <strong>de</strong>pen<strong>de</strong>n <strong>de</strong>l mo<strong>de</strong>lo seleccionado por el usuario:<br />
Cuando se selecciona “Mo<strong>de</strong>lo <strong>de</strong> Baranyi y Roberts – no lag”, el parámetro mCurv se<br />
establece como cero, es <strong>de</strong>cir, el mo<strong>de</strong>lo <strong>de</strong>scribe solo el crecimiento o la muerte y la fase<br />
estacionaria.<br />
Figura 2. Curva <strong>de</strong> crecimiento ajustada al Mo<strong>de</strong>lo <strong>de</strong> Baranyi y Roberts sin fase <strong>de</strong> latencia.
Al seleccionar “Mo<strong>de</strong>lo <strong>de</strong> Baranyi y Roberts – no asymtot ", el parámetro nCurv se establece<br />
en cero. En consecuencia, el mo<strong>de</strong>lo <strong>de</strong>scribe sólo la fase <strong>de</strong> latencia o retraso/interfase y la<br />
fase <strong>de</strong> crecimiento o la fase <strong>de</strong> muerte.<br />
Figura 3. Curva <strong>de</strong> crecimiento ajustada al Mo<strong>de</strong>lo <strong>de</strong> Baranyi y Roberts sin interfase <strong>de</strong><br />
exponencial a estacionaria y sin fase estacionaria.<br />
Al seleccionar el “Mo<strong>de</strong>lo <strong>de</strong> Baranyi y Roberts – mo<strong>de</strong>lo completo", se utilizan para mCurv y<br />
nCurv los siguientes valores por <strong>de</strong>fecto:<br />
• mCurv=10<br />
• nCurv=1<br />
Figura 4. Curva <strong>de</strong> crecimiento ajustada al Mo<strong>de</strong>lo <strong>de</strong> Baranyi y Roberts completo: fase <strong>de</strong><br />
latencia/interfase, fase exponencial/interfase, fase estacionaria.
3.3. Mo<strong>de</strong>lo trilineal, mo<strong>de</strong>los bifásicos o mo<strong>de</strong>los lineales<br />
Mo<strong>de</strong>lo trilineal: Como su nombre lo indica, el mo<strong>de</strong>lo trilineal <strong>de</strong>scribe la curva <strong>de</strong> crecimiento<br />
bacteriana con tres líneas: la fase <strong>de</strong> latencia y la fase estacionaria se <strong>de</strong>scriben por dos rectas<br />
horizontales. La pendiente <strong>de</strong> la tercera línea recta <strong>de</strong>scribe la fase <strong>de</strong> crecimiento / muerte y es<br />
llamada “velocidad o tasa máxima”.<br />
Figura 5. Curva <strong>de</strong> crecimiento ajustada al Mo<strong>de</strong>lo <strong>de</strong> Baranyi y Roberts trilineal.<br />
Mo<strong>de</strong>los bifásicos: Cuando la curva bacteriana no presenta fase <strong>de</strong> latencia / interfase o<br />
ninguna fase estacionaria, los mo<strong>de</strong>los bifásicos sin fase <strong>de</strong> latencia (figura 6) o estacionaria<br />
(figura 7) pue<strong>de</strong>n preferirse al mo<strong>de</strong>lo trilineal.<br />
Figura 6. Curva <strong>de</strong> crecimiento ajustada a un Mo<strong>de</strong>lo bifásico sin fase <strong>de</strong> latencia/interfase<br />
(mo<strong>de</strong>lo no-lag)
Figura 7. Curva <strong>de</strong> crecimiento ajustada a un Mo<strong>de</strong>lo bifásico sin interfase/fase estacionaria<br />
(mo<strong>de</strong>lo no-asymtot)<br />
Mo<strong>de</strong>lo lineal: Cuando el conteo bacteriano <strong>de</strong>scribe sólo la fase <strong>de</strong> crecimiento / fase <strong>de</strong><br />
muerte, los datos pue<strong>de</strong> ser ajustados a un mo<strong>de</strong>lo lineal.<br />
Figura 8. Curva <strong>de</strong> crecimiento ajustada a un Mo<strong>de</strong>lo lineal (sin fases <strong>de</strong> latencia y estacionaria<br />
ni interfases))<br />
3.4. Ajuste <strong>de</strong>l Mo<strong>de</strong>lo<br />
Después <strong>de</strong> haber cargado los datos en la "Entrada <strong>de</strong> datos – Data input" <strong>de</strong> texto y<br />
seleccionado el mo<strong>de</strong>lo más a<strong>de</strong>cuado, los datos se pue<strong>de</strong>n ajustar pulsando el botón "Ajustar -<br />
Fit".<br />
Cuando se usa el mo<strong>de</strong>lo lineal, los coeficientes <strong>de</strong> regresión son calculados con base a una<br />
regresión lineal normal.
Para cualquier otro mo<strong>de</strong>lo, haga clic en el botón 'Fit' el cual hace una primera estimación <strong>de</strong> los<br />
parámetros <strong>de</strong> crecimiento y, a continuación, el programa utiliza este cálculo como punto <strong>de</strong><br />
partida para la estimación <strong>de</strong> los parámetros no lineales por un método <strong>de</strong> mínimos cuadrados.<br />
Si el algoritmo <strong>de</strong> ajuste no converge a una solución, aparece el siguiente mensaje <strong>de</strong> error<br />
(mostrado en la figura 9): " Datos no compatibles con el mo<strong>de</strong>lo seleccionado – Data not<br />
consistent with mo<strong>de</strong>l selected". Usted <strong>de</strong>be entonces tratar <strong>de</strong> ajustar los datos con otro<br />
mo<strong>de</strong>lo.<br />
Figura 9. Mensaje <strong>de</strong> error <strong>de</strong>bido a la incompatibilidad <strong>de</strong>l mo<strong>de</strong>lo seleccionado (mo<strong>de</strong>lo<br />
completo) con respecto a los datos.<br />
Cuando el programa encuentra las estimaciones para los parámetros, la estimación <strong>de</strong> los<br />
parámetros y sus errores estándar se muestran en el panel <strong>de</strong> “Resultados", junto con la grafica<br />
<strong>de</strong> ajuste <strong>de</strong> datos (figura 10):<br />
"R-square (R-cuadrado): R-cuadrado <strong>de</strong>l ajuste<br />
'SE of Fit' (Desviación estándar): Error estándar <strong>de</strong>l ajuste
Figura 10. Resultado <strong>de</strong>l ajuste <strong>de</strong> los datos con un mo<strong>de</strong>lo compatible, en esta caso ajustado al<br />
mo<strong>de</strong>lo <strong>de</strong> Baranyi y Roberts sin fase <strong>de</strong> latencia.<br />
3.5. Unida<strong>de</strong>s<br />
Las unida<strong>de</strong>s <strong>de</strong> las estimaciones se refieren a las unida<strong>de</strong>s <strong>de</strong> los datos en la "Entrada <strong>de</strong><br />
datos" <strong>de</strong> texto. Por ejemplo, si el conteo bacteriano está expresado en log10 ufc/g y el tiempo en<br />
horas, la máxima velocidad <strong>de</strong> crecimiento/velocidad <strong>de</strong> muerte se expresará en log10 ufc/g/h. Si<br />
el conteo se encuentra expresado en Ln ufc/g y el tiempo en horas, la máxima velocidad se<br />
expresará en Ln ufc/g/h (la velocidad expresada en estas unida<strong>de</strong>s también se <strong>de</strong>nomina<br />
“Máxima velocidad específica”).<br />
De manera similar, si el tiempo está en horas, la fase <strong>de</strong> latencia/interfase será expresado en<br />
horas. Las estimaciones para el valor inicial y el valor final también se encuentran en las mismas<br />
unida<strong>de</strong>s que los datos (normalmente log10 ufc/g ó Ln ufc/g).
4. Referencias Bibliográficas<br />
Baranyi, J., Roberts, T.A. (1994). A dynamic approach to predicting bacterial growth in food. Int.<br />
J. Food Microbiol. 23, 277 – 294.<br />
DMFit. (2009). Software <strong>de</strong> mo<strong>de</strong>lamiento dinámico edición on-line. Disponible en:<br />
http://ifrsvwww<strong>de</strong>v.ifrn.bbsrc.ac.uk/CombasePMP/GP/DMFit.aspx.