miércoles, 26 de junio de 2013

Problema34: Análisis de Regresión

Ej34. Los datos reflejados en la siguiente tabla relacionan los niveles de concentración de una sustancia, presentes en 6 muestras, con la respuesta observada en un colorímetro:


x.4.24505.948012.865026.666035.186050.6510
y0.02800.04400.06800.09200.10950.1295

El objetivo es encontrar una curva de calibración del tipo: y = α·eβx, que permita determinar la concentración a partir de la respuesta observada en el colorímetro.

Determinar el valor de α y de β.

NOTA: La variable x representa la respuesta del colorímetro y la variable y representa la concentración de la sustancia de la muestra. 


Este problema comparte los mismos datos que el Problema14 pero con distinta curva de calibración.




El modelo que nos ofrece el problema no es un modelo lineal simple, por lo que tenemos que adecuarlo, en este caso, actuamos mediante logaritmo neperiano:


Ln(y) = Ln(α·eβx) = Ln(α) + β·x

Cambio de Variable:
¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨
· y* = Ln(y)
· b0 = Ln(α)

El modelo ajustado, con el cambio de variable, es:

y* = b0 + b1·x

Y la tabla quedará tal y como sigue:

x*.4.24505.948012.865026.666035.186050.6510
y*Ln(0.0280).Ln(0.0440).Ln(0.0680).Ln(0.0920).Ln(0.1095)Ln(0.1295)

Ahora, podemos hacer una recopilación de datos que se extraen de la tabla una vez realizado el cambio de variable.

· n = 6

·

·

·

·

·


Para calcular la pendiente, la expresión matemática es:


Para obtener su valor, necesitamos saber los valores de Sxy y Sxx:

·

·


Por lo tanto, la pendiente es:


Una vez obtenida la pendiente, podemos tener el valor del estimador para la ordenada. Sustituimos valores:


Por lo tanto, la ecuación de regresión ajustada es:

y*(x) = -3.331880 + 0.029227·x

Siendo:

· y* = Ln(y)
· b0 = Ln(α) = -3.331880
· b1 = β = 0.029227

El modelo no lineal quedará:

 y = 0.035726·e0.029227·x

Siendo:

· α = e-3.331880 ≈ 0.035726

2 comentarios:

Anónimo dijo...

¿El modelo no lineal seria
y(x)= 0.035726·e(elevado a x.0.029227)

Siendo:

· α = e-3.331880 ≈ 0.035726. ?

AqueronteBlog dijo...

Buenas:

Tienes razón, ha sido una errata, además, incluso el enunciado del problema te indica como debe ser el modelo no lineal, y evidentemente, debe ser del tipo exponencial.

Ya está subsanado, muchas gracias por el apunte.


Un saludo.