02.07.2013 Views

DMFitwebedition - Blog de ESPOL

DMFitwebedition - Blog de ESPOL

DMFitwebedition - Blog de ESPOL

SHOW MORE
SHOW LESS

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.

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!