#180　Pythonを用いて2014年（平成26年）大阪府産業連関表最終需要部門を推計（その３）

2023年1月23日 2023年1月23日

宮崎稔也

投稿#178から、Pythonを用いて2014年（平成26年）大阪府産業連関表における最終需要部門の推計に取り掛かっています。

今回の投稿も、投稿#179に引き続き、最終需要部門の推計についてです。

（１１）輸出（普通貿易）

推計方法の確認

各部門について、平成23年全国基本表から平成25年全国延長表への輸出率（輸出額／生産額）の増減率を、平成23年大阪府基本表の輸出率に乗じて平成25年輸出率を求め、平成25年大阪府延長表の生産額に乗じることにより求めた。

大阪府総務部統計課「平成25年（2013年）大阪府産業連関表（延長表）報告書」より引用

実際の推計

コード180−1

# (11)輸出（普通貿易）の推計
h23_basic_df = pd.read_csv("H23大阪府基本表.csv") # 平成23年大阪府基本表を読み込む
h23_basic_df = to_str_in_code(h23_basic_df, 3) # indexを文字列に型変換

h23_japan_basic_df = pd.read_csv("平成23年産業連関表_統合小分類.csv") # 平成23年全国基本表を読み込む
h23_japan_basic_df = to_str_in_code(h23_japan_basic_df, 3) # indexを文字列に型変換

for idx in h26_manu_stocks_df.index: # 各部門について
    h23_japan_pro = h23_japan_basic_df.loc[idx, "9700"]
    h26_japan_pro = h26_japan_enl_fd_df.loc[idx, "970000"]
    h23_pro = h23_basic_df.loc[idx, "9700"]
    if (h23_japan_pro == 0 or h26_japan_pro == 0 or h23_pro == 0):
        h26_fd_df.loc[idx, "801101"] = 0
        continue
    else:
        h23_japan_exrate = h23_japan_basic_df.loc[idx, "8011"] / h23_japan_pro
        # 平成26年全国基本表における輸出率（輸出額／生産額）を算出
        h26_japan_exrate = ((h26_japan_enl_fd_df.loc[idx, "801101"] + h26_japan_enl_fd_df.loc[idx, "801102"]) 
                                / h26_japan_pro)
        roc = (h26_japan_exrate / h23_japan_exrate) # 輸出率の増減率rocを算出
        h23_exrate = h23_basic_df.loc[idx, "8011"] / h23_pro
        if h23_exrate == 0:
            h26_fd_df.loc[idx, "801101"] = 0
        else:
            h26_exrate = h23_exrate * roc # 平成23年大阪府基本表の輸出率に増減率をかける
            h26_fd_df.loc[idx, "801101"] = h26_pro_df.loc["production", idx] * h26_exrate

（１２）輸出（特殊貿易）および（１３）輸出（直接購入）

推計方法の確認

（１２）輸出（特殊貿易）
全国延長表の輸出（特殊貿易）額を、宿泊旅行統計から得られる外国人宿泊者数の対全国比で按分した。

１３）輸出（直接購入）
全国延長表の輸出（直接購入）額を、宿泊旅行統計から得られる外国人宿泊者数の対全国比で按分した。

大阪府総務部統計課「平成25年（2013年）大阪府産業連関表（延長表）報告書」より引用

実際の推計

外国人宿泊者数は、観光庁「宿泊旅行統計調査」の「平成26年1月～12月分（年の確定値）」の集計結果の「第2表」における外国人延べ宿泊者数の値を使用しました。

コード180-2

# (12)輸出（特殊貿易）と(13)輸出（直接購入）の推計

foreigner_japan = 44824600 # 宿泊旅行統計から得られる外国人延べ宿泊者（日本全体） 
foreigner_osaka = 6200160 # 宿泊旅行統計から得られる外国人延べ宿泊者（大阪府）
ratio = foreigner_osaka / foreigner_japan # 宿泊旅行統計から得られる外国人延べ宿泊者の対全国比
# H26全国延長表の輸出（特殊貿易）額およびと輸出（直接購入）額を、対全国比で按分
for idx in h26_manu_stocks_df.index:
    h26_fd_df.loc[idx, "801102"] = h26_japan_enl_fd_df.loc[idx, "801102"] * ratio # 輸出（特殊貿易）を対全国比で按分
    h26_fd_df.loc[idx, "801200"] = h26_japan_enl_fd_df.loc[idx, "801200"] * ratio # 輸出（直接購入）を対全国比で按分

（１０）調整項

順番が前後しましたが、（１０）調整項をここで推計します。

推計方法の確認

全国延長表の輸出額に対する調整項の率を、平成25年大阪府延長表の輸出額に乗じて推計した。

大阪府総務部統計課「平成25年（2013年）大阪府産業連関表（延長表）報告書」より引用

実際の推計
輸出額は以下の式で求めます。
$$
輸出 = 輸出（普通貿易）+ 輸出（特殊貿易）
$$

コード180-3

# (10)調整項の推計
# 列771100 調整項を追加
h26_fd_df.insert(11, "771100",0)

for idx in h26_manu_stocks_df.index:
    japan_exp = h26_japan_enl_fd_df.loc[idx, "801101":"801102"].sum() # 全国延長表の輸出額
    if japan_exp == 0:
        h26_fd_df.loc[idx, "771100"] = 0
        continue
    else:
        japan_adj = h26_japan_enl_fd_df.loc[idx, "801300"] # 全国延長表の調整項
        japan_ratio = japan_adj / japan_exp # 全国延長表の輸出額に対する調整項の率
        osaka_exp = h26_fd_df.loc[idx, "801101":"801102"].sum() # H26大阪府延長表の輸出額
        h26_fd_df.loc[idx, "771100"] = osaka_exp * japan_ratio # 調整項＝輸出額＊率

（１４）移出

推計方法の確認

各部門について、平成23年大阪府基本表の移出率（移出額／生産額）を使用し、平成25年大阪府延長表の生産額に乗じることにより求めた。

大阪府総務部統計課「平成25年（2013年）大阪府産業連関表（延長表）報告書」より引用

実際の推計

コード180-4

# (14)移出の推計

# 列821100 移出を追加
h26_fd_df.insert(19, "821100", 0)

for idx in h26_manu_stocks_df.index:
    h23_pro = h23_basic_df.loc[idx, "9700"] # H23大阪府基本表の生産額
    if h23_pro == 0:
        h26_fd_df.loc[idx, "821100"] = 0
        continue
    else:
        h23_ship = h23_basic_df.loc[idx, "8211"] # H23大阪府基本表の移出額
        ship_ratio = h23_ship / h23_pro # H23大阪府移出率
        # H23大阪府移出率にH26大阪府延長表の生産額を乗じる
        h26_fd_df.loc[idx, "821100"] = h26_pro_df.loc["production", idx] * ship_ratio

Follow me!