Calculo numerico.

Problema#5. Los siguientes datos son para el estudio del  consumo de carne de res en Colombia

INDI: ingreso personal disponible per cápita ($/Hab)

CRES: Consumo per cápita de carne de res (ton/año)

Se pide:
a)Construya un diagrama de dispersión

b) Haga los cálculos para la estimación de
un del modelo exponencial

 c) Determine el grado de ajuste e interprételo

d)Elabore un breve reporte del
modelo exponencial

e)Haga los cálculos para la estimación de un del modelo
potencial

f)Determine el grado de ajuste e interprételo

g)Elabore un breve reporte del modelo potencial

h) Haga los cálculos para la estimación de un del modelo
cuadrático

i)Determine el grado de ajuste e interprételo

j) Elabore un breve reporte
del modelo cuadrático.

K) De los tres modelos hallados cual es el mejor?

l) Haga la curva de cada modelo (con diferente color c/u) y añádala al diagrama de dispersión.

Solución:

a) Diagrama de dispersión
[pic 1]

     Como puede observarse de la gráfica, los datos no muestran una tendencia lineal. Por lo tanto se intentara analizar los datos con otros modelos exponencial, potencial y cuadrático.

COMANDO EN R

plot(x,y,xlab='INDI',ylab='CRES',main='DIAGRAMA DE DISPERSION', col="gold", pch=15)

 MODELO EXPONENCIAL

b) CÁLCULOS

El modelo lineal transformado se linealiza en: yp=ap+b*x, donde yp=log(y) y ap=log(a)

Cálculos de a y b

b =( n ∑xiyi’ – (∑xi) (∑yi’))/ ( n ∑x^2 – (∑xi)^2 ) =
(12 x 4971971.872 – (339780) (160.332651))/( 12 x 11567982384-(339780)^2 )=4.963223x10-5

       _      _

a’ = y' – bx

_

yp = ∑yip /n = 160.332651/12 =13.361054

_
x = ∑xi/n= 339780/12 = 28315

           ___                 ___

ap = 13.361054 – (28315 x 4.963223x10-5 ) =  11.8225=1,8225x101

a=exp (1,8225x101)=82227470.9

Fórmula para calcular R^2

C) Calculo del grado de ajuste
 R^2 = 1-[∑ (ei) ^2 /∑ (yi’ – mean (y’i))^2
=1-(0.05110572/0.1415759) =  0.6390225

Interpretación:

     El porcentaje de variabilidad explicado por el modelo resulto con una bondad de ajuste de   63.9%  lo que indica que el modelo exponencial es muy bueno para explicar el Consumo per cápita de carne de res en función de ingreso personal disponible per cápita

Verificación con el lenguaje r

reg<-lm(yp~x); summary(reg); coef=reg$coefficients; a=coef[1]; b=coef[2]

Call:

Lm (formula = yp ~ x)

Residuals:

      Min                1Q         Median               3Q            Max

-0.13304     -0.01946    0.01411       0.04252     0.07672

...