viernes, 7 de noviembre de 2008

PRONOSTICOS DEL PIB (PATRICIO TAPIA)

--PRONOSTICO DEL PIB PARA EL 1998
--TOMANDO COMO VALORES REFERENCIALES LOS ANTERIORES AÑOS DESDE EL 1993 HASTA EL 1997
zt=[18573 20195 21268 23635 23255]';
n=length(zt)
zt1=[10056;zt(1:n-1)]
zt2=[0;10056;zt(1:n-2)]
x=[zt1,zt2]
b=inv(x'*x)*(x'*zt)
s=poly(0,'s')
s=s^2-b(1)*s-b(2)
raiz=roots(s)
t=[1:n]'
z1=(raiz(1).^t)
z2=(raiz(2).^t)
x=[z1 z2]
a=inv(x'*x)*(x'*zt)
z=a(1)*z1+a(2)*z2
error1=zt-z
sqrt(error1'*error1)
plot(t,error1)
--EN LA SIGUIENTE LINEA SE DA EL PRONOSTICO DEL VALOR PARA EL 1998
t=[1:n+1]';
z1=(raiz(1).^t)

--EN CONCLUSION : EN VALOR RESULTADO DEL PRONOSTICO ES MUY ELEVADO Y DIFERENTE DEL VALOR –
--QUE SE ANUNCIO POR LA PRENSA.EN AQUELLA EPOCA DEBIDO A LOS FACTORES ADVERSOS POR TODOS
--CONOCIDOS .

===============================================

--PRONOSTICO DEL PIB PARA EL 2008
--TOMANDO COMO VALORES REFERENCIALES LOS ANTERIORES AÑOS DESDE EL 2000 HASTA EL 2007
zt=[21250 24899 28636 32642 37187 41763 45789 48508]';
n=length(zt)
zt1=[0;zt(1:n-1)]
zt2=[0;0;zt(1:n-2 )]
x=[zt1,zt2]
b=inv(x'*x)*(x'*zt)
s=poly(0,'s')
s=s^2-b(1)*s-b(2)
raiz=roots(s)
t=[1:n]'
z1=(raiz(1).^t)
z2=(raiz(2).^t)
x=[z1 z2]
a=inv(x'*x)*(x'*zt)
z=a(1)*z1+a(2)*z2
error1=zt-z
sqrt(error1'*error1)
plot(t,error1)
--EN LA SIGUIENTE LINEA SE DA EL PRONOSTICO DEL VALOR PARA EL 2008
t=[1:n+1]';
z1=(raiz(1).^t)
z2=(raiz(2).^t)
z=a(1)*z1+a(2)*z2
plot(t,z)


==================================================================
==================================================================

No hay comentarios: