Exchange of Goods and Services in the U.S. for 1947 (in billions of 1947 dollars)
In [58]:
import numpy as np
import pandas as pd
In [2]:
ls
In [16]:
f = open('leontief.csv')
In [17]:
Names = (f.next().strip()).split(',')
data =[]
In [18]:
for row in f:
#np.vstack(data,(row.strip().split(',')))
data.append(row.strip().split(','))
In [35]:
data = np.array(data)
# Extracting External Demand
Ex_Demand= data[:-1][:,-1]
Ex_Demand=Ex_Demand.astype(float)
# Extracting Gross Product
Gross_Product = data[-1][:-1]
Gross_Product=Gross_Product.astype(float)
# Extracting Expenditures Matrix
A = data[:-1][:,:-1]
A=A.astype(float)
print Ex_Demand
print Gross_Product
print A
In [43]:
Demand = []
for i in xrange(len(A[0])):
Demand.append (sum(A[i])+Ex_Demand[i])
Demand= np.array(Demand)
Income = Gross_Product - Demand
Income
Out[43]:
In [68]:
A_tc = A / Gross_Product
print A_tc
In [74]:
L = (np.identity(len(A[0])) - A_tc)
L
Out[74]:
In [73]:
linalg.inv(L)
Out[73]:
In [77]:
P = np.dot(linalg.inv(L), Ex_Demand)
In [78]:
P
Out[78]:
$ Turisti_{t} = \alpha * Turisti_{t-1} + \beta_{0} +\beta_{1}GDP_{cap} + \beta_{2}\frac{Prices_{foreign}}{Prices_{home}} + \beta_{3}*ExRate_{\frac{Foreign}{MDL}}+\beta_{4}CO_{2}$
In [ ]:
In [ ]:
In [ ]: