Análisis Gráfico y Estadístico del Ajuste Polinomial

¿Qué polinomio debemos escoger para ajustar datos?

Muchos profesionales y estudiantes de ingeniería, cuando quieren representar matemáticamente sus datos experimentales, después de graficarlos en Excel, optan por agregar una línea de tendencia de ajuste de un polinomio de cuarto orden o superior, pensando que mientras más términos tenga la ecuación (mayor grado del polinomio), mayor será su precisión con la que represente a los datos. Este criterio puede llevarnos a cometer errores en procesos y operaciones en las que la precisión de los datos juega un papel fundamental, o también, hacernos perder el tiempo con ecuaciones larguísimas que no justifican su robustez.

Para tener la certeza de que la ecuación seleccionada es la correcta debemos hacer uso de herramientas gráficas y estadísticas. A continuación detallamos un ejemplo de correlación de datos para explicar sobre dichas herramientas y ver cómo se emplean para seleccionar el mejor modelo.

Ajuste Polinomial de los datos de la Presión de Vapor del Benceno

Tenemos los siguientes datos:

TABLA 1

Mediante el ajuste polinomial se obtendrá la ecuación que permita calcular la presión de vapor del benceno (P) en función de la temperatura (T). Esta ecuación será de la siguiente forma:

donde “n” es el grado del polinomio y las “a” son sus coeficientes.

Con el programa de regresión polinomial del software POLYMATH obtenemos los coeficientes de las ecuaciones de primer orden hasta el quinto, vistos en la siguiente tabla:

TABLA 2: Ecuaciones de ajuste obtenidas con POLYMATH

Ahora nos toca analizar y escoger la ecuación que mejor represente a los datos de la presión de vapor del benceno.

Análisis Estadístico

a) Coeficiente de correlación (R^2)
El Coeficiente de Correlación es usado para juzgar si el modelo representa correctamente a los datos. Si su valor se aproxima a 1 entonces el modelo es correcto. Observamos en la tabla 2 que casi todos los R^2 de los polinomios se aproximan dramáticamente a 1, a excepción del de primer orden (el cual da la gráfica de una recta), por lo tanto, podemos desecharlo de nuestro proceso de selección.

b) Varianza
La varianza está definida por la siguiente ecuación:

Donde:

* V = N - (n+1), siendo “N” el número de datos y “n” el grado del polinomio.
* P(obs) son las presiones experimentales, y P(calc) son las que dan las ecuaciones.

A menor varianza, mayor será la precisión de ajuste de la ecuación. En la tabla 2 se observa que el polinomio de tercer orden es el que tiene la varianza más pequeña (0,0371448), aunque el polinomio de cuarto orden tiene un valor cercano (0,0375345). El polinomio de segundo orden por tener una varianza grande es desechado de nuestro proceso de selección.

c) Intervalos de Confianza
Cuando se analizan ecuaciones con muchos parámetros (a partir de cuarto orden), es importante considerar los intervalos de confianza de los valores de los parámetros porque representan la incertidumbre asociada a dichos valores. Por ejemplo, un parámetro de valor 4 con un intervalo de 0,3 debería ser representado como 4±0,3, es decir, su “verdadero” valor estaría en el rango de 3,7 a 4,3. Un mal ajuste entre el modelo y los datos es a menudo indicado por un intervalo de confianza que incluye al cero dentro del rango.
En la siguiente tabla tenemos los intervalos de confianza para los polinomios de grado 3 y 5:

TABLA 3: Intervalos de confianza obtenidos con POLYMATH

Se observa en la tabla que el intervalo de confianza de los coeficientes a4 (2,30E-08±1,66E-07 da un rango de -1.43E-07 a 1.89E-7) y a5 del polinomio de 5to orden incluyen al cero. Aunque por razones de espacio no está incluida en la tabla, el intervalo de confianza del coeficiente a4 del polinomio de 4to orden también incluye al cero. Sólo el polinomio de 3er orden no presenta ningún intervalo que incluya al cero, por lo tanto es el que mejor se ajusta a los datos.

CONCLUSIONES

· En este ejemplo los polinomios de 3er y 4to orden son los que más se ajustan a los datos, siendo el mejor el de grado 3 por dos características que no tiene el otro: una menor varianza, e intervalos de confianza que no incluyen el valor de cero.
· El polinomio de 5to orden a pesar de tener un índice de correlación cercano a 1 y una buena curva de ajuste, no valdría la pena gastar tiempo y esfuerzo en su utilización.

HENRY FLORES

PD: si desea recibir el artículo completo que incluye el análisis de residuales y las curvas de ajuste, pídala escribiendo al correo : henryabraham@yahoo.com

Condiciones Normales y Condiciones Estándares

En el mundo de la Ingeniería Química hay dudas con respecto a la diferencia entre condiciones normales (N.C.) y condiciones estándares (S.C.); las mismas que han sido sembradas durante años por libros de textos, catálogos y profesionales confundidos.

Todos empleamos ambas condiciones a la presión de 1 atmósfera (1,013 bar; 101,325 kPa; 14,7 psia; etc.). La confusión está en la TEMPERATURA. Muchos asignan e intercambian los siguientes valores para cada condición: 0ºC, 15ºC y 25ºC. En las siguientes líneas aclararemos la disyuntiva.

¿Qué dice la IUPAC?

La Unión Internacional de Química Pura y Aplicada (IUPAC, siglas en inglés) publicó en su página web un glosario revisado en el 2000, en la cual se definen los términos “Condiciones Normales” (Normal Conditions), “Estándar” (Standard) y “Condiciones Estándares para los gases” (Standard Conditions for Gases). Pueden leerlo en el siguiente vínculo:

http://www.iupac.org/reports/1990/6211calvert/glossary.html

He aquí la traducción al castellano:

Condiciones Normales: Es un término cualitativo que depende de la preferencia del investigador; a menudo implica la presión del ambiente y la temperatura del lugar. Es preferible que estas variables de temperatura y presión sean fijadas como valores representativos de las condiciones actuales (o rango de condiciones) empleadas en el estudio.

Estándar: Es un valor exacto o un concepto establecido por autoridad o por acuerdo, que sirve como modelo o regla en la medición de una cantidad o en el establecimiento de una práctica o procedimiento, en el análisis de la contaminación del aire, o el uso de los gases, líquidos y sólidos de referencia estándar para calibrar equipos.

Condiciones Estándares para Gases: A veces se indica con la abreviación STP. Temperatura: 273,15 K (0ºC). Presión: 10e5 pascales. La IUPAC recomienda descontinuar el uso inicial de la presión de 1 atm (equivalente a 1,01325 x 10e5 Pa) como presión estándar.

CONCLUSIÓN

Para fines prácticos no hay una diferencia significativa entre 1,01325 x 10e5 Pa y 10e5 Pa. Podemos seguir empleando la presión de 1 atm para cálculos que no requieran un rigor científico.
De acuerdo a las definiciones anteriores, podemos resumir lo siguiente:

Condiciones Estándares: 1 atm y 0ºC.

Condiciones Normales: Presión y Temperatura del lugar. Depende de las condiciones a las cuales se esté haciendo el experimento, estudio o medición; comúnmente para la presión es 1 atm, y la temperatura: 15ºC, 20ºC, 25ºC ó 27ºC.

HENRY FLORES (henryabraham@yahoo.com)