#181 Pythonを用いて2014年（平成26年）大阪府産業連関表最終需要部門を推計（その４）
宮崎 稔也投稿#178から、Pythonを用いて2014年（平成26年）大阪府産業連関表における最終需要部門の推計に取り掛かっています。
今回の投稿は、移輸入額の推計についてです。
不要な列の削除と独自の列の挿入
最終需要部門において、大阪府産業連関表にはないが全国延長表にはある列、もしくは全国延長表にはないが大阪府産業連関表には独自に存在する列がいくつか見受けられます。
前者に該当する列は削除し、後者に該当する列を新たに追加していきます。
コード181-1
# 大阪府延長表で不要な列を削除
h26_fd_df = h26_fd_df.drop(["801300", "820000", "841101", "841102", "841200"], axis=1) # 不要な列である801300を削除
# 大阪府延長表で独自の列を挿入
# 列 851101 輸入（普通貿易）を挿入
h26_fd_df.insert(21, "851101", 0)
# 列 851102 輸入（特殊貿易）を挿入
h26_fd_df.insert(22, "851102", 0)
# 列 851200 輸入（直接購入）を挿入
h26_fd_df.insert(23, "851200", 0)
# 列 870000を871100 輸入品関税品に変更
h26_fd_df = h26_fd_df.rename({"870000": "871100"}, axis=1)
# 列 "891100" 移入を挿入
h26_fd_df.insert(27, "891100", 0)
# 列 "900000" 最終需要部門計を挿入
h26_fd_df.insert(28, "900000", 0)（１５）輸入（普通貿易）
推計方法の確認
全国延長表の国内需要額に対する輸入（普通貿易）率を、府内需要額に乗じて推計した。
大阪府総務部統計課「平成25年（2013年）大阪府産業連関表（延長表）報告書」より引用
実際の推計
コード182-2
# 851100 輸入（普通貿易）の推計
for idx in h26_manu_stocks_df.index:
h26_totalfd = h26_fd_df.loc[idx, "711100":"771100"].sum() # 780000 府内最終需要計を算出
h26_totalintsec = h26_adjusted_mid_df.loc[idx, :].sum() # 7000 府内内生部門計の算出
h26_totaldemand = h26_totalintsec + h26_totalfd # 790000 府内需要合計の算出
h26_fd_df.loc[idx, "780000"] = h26_totalfd # 算出した780000 府内最終需要計を代入
h26_fd_df.loc[idx, "790000"] = h26_totaldemand # 算出した7000 内生部門計を代入
h26_japan_totaldemand = h26_japan_enl_fd_df.loc[idx, "790000"] # 全国延長表の国内需要額
if h26_japan_totaldemand == 0:
h26_fd_df.loc[idx, "851101"] = 0
continue
else:
h26_japan_import = h26_japan_enl_fd_df.loc[idx, "841101"] # 全国延長表の輸入（普通貿易）
h26_japan_imp_ratio = h26_japan_import / h26_japan_totaldemand # 全国延長表の国内需要額に対する輸入（普通貿易）の率
h26_fd_df.loc[idx, "851101"] = h26_totaldemand * h26_japan_imp_ratio（１６）輸入（特殊貿易）
推計方法の確認
全国延長表の輸入（特殊貿易）額を、出入国管理統計から得られる出国日本人の対全国比で按分した。
大阪府総務部統計課「平成25年（2013年）大阪府産業連関表（延長表）報告書」より引用
実際の推計
法務省出入国在留管理庁の「出入国管理統計統計表」より、全国及び大阪府における2014年（平成26年）の出国日本人数の値を得ました。
推計に利用したPythonのコードは、以下のようになります。
コード180-3
# 輸入（特殊貿易）の推計
japan_dep = 16903388 # 出入国管理統計における出国日本人の全国値
osaka_dep = 3224562 + 3959 # 出入国管理統計における出国日本人の大阪の値(関空 + 大阪港)
for idx in h26_manu_stocks_df.index:
japan_imp_spe = h26_japan_enl_fd_df.loc[idx, "841102"] # 全国延長表の輸入（特殊貿易）
# 全国延長表の輸入（特殊貿易）額を出入国管理統計から得られる出国日本人の対全国比で按分
osaka_imp_spe = japan_imp_spe * osaka_dep / japan_dep
h26_fd_df.loc[idx, "851102"] = osaka_imp_spe（１７）輸入（直接購入）
推計方法の確認
全国延長表の輸入（直接購入）額を、出入国管理統計から得られる出国日本人の対全国比で按分した。
大阪府総務部統計課「平成25年（2013年）大阪府産業連関表（延長表）報告書」より引用
実際の推計
（１６）と同様にして求めることができるので、コードは割愛します。
（１８）関税
推計方法の確認
全国延長表の輸入額に対する関税率を、平成25年大阪府延長表の輸入額に乗じて推計した。
大阪府総務部統計課「平成25年（2013年）大阪府産業連関表（延長表）報告書」より引用
実際の推計
域内（国内または府内）輸入額は以下の式から求めます。
$$
域内輸入額 = 域内輸入（普通貿易）+ 域内輸入（特殊貿易） + 域内輸入（直接購入）
$$
推計に利用したPythonのコードは、以下のようになります。
コード180-4
# 関税の推計
for idx in h26_manu_stocks_df.index:
# 全国延長表の輸入額 = 輸入（普通貿易） + 輸入（特殊貿易）+ 輸入（直接購入）
japan_total_imp = h26_japan_enl_fd_df.loc[idx, "841101":"841200"].sum()
if japan_total_imp == 0:
h26_fd_df.loc[idx, "861100"] = 0
continue
else:
japan_customduties = h26_japan_enl_fd_df.loc[idx, "851100"] # 関税額
japan_duty_ratio = japan_customduties / japan_total_imp # 全国延長表の輸入額に対する関税率
osaka_total_imp = h26_fd_df.loc[idx, "851101":"851200"].sum()# H26大阪府延長表の輸入額
osaka_customduties = osaka_total_imp * japan_duty_ratio # 関税率をH26大阪府延長表の輸入額に乗じる
h26_fd_df.loc[idx, "861100"] = osaka_customduties # 列861100に大阪府延長表の関税額を入力（１９）輸入品商品税
推計方法の確認
全国延長表の輸入額に対する輸入品商品税率を、平成25年大阪府延長表の輸入額に乗じて推計した。
大阪府総務部統計課「平成25年（2013年）大阪府産業連関表（延長表）報告書」より引用
実際の推計
（１８）関税と推計方法が同じなので、Pythonのコードは割愛します。
（２０）移入
推計方法の確認
各部門について、平成23年大阪府基本表の移入率（移入額／府内需要合計）を使用し、平成25年大阪府延長表の府内需要合計に乗じることにより求めた。
大阪府総務部統計課「平成25年（2013年）大阪府産業連関表（延長表）報告書」より引用
実際の推計
推計に用いたPythonのコードは、以下のようになります。
コード180-5
# 移入の推計
for idx in h26_manu_stocks_df.index:
h23_total_demand = h23_basic_df.loc[idx, "7900"] # H23大阪府基本表の府内需要合計
if h23_total_demand == 0:
h26_fd_df.loc[idx, "891100"] = 0
continue
else:
h23_immig = h23_basic_df.loc[idx, "8911"] # H23大阪府基本表の移入額
h23_immig_ratio = h23_immig / h23_total_demand # H23大阪府基本表の移入率
h26_total_demand = h26_fd_df.loc[idx, "790000"] # H26大阪府延長表の府内需要合計
# H26大阪府延長表の府内需要合計にH23大阪府基本表の移入率を乗じる
h26_immig = h26_total_demand * h23_immig_ratio
h26_fd_df.loc[idx, "891100"] = h26_immigFollow me!
