# @title Google Drive に置いたCSVを読込(ここではMy Drive後の【covid_Oki.csv】を読んでいる)
import pandas as pd
# @title エクセルでいう「A行をindexにしたい」場合。CSVを読み込むときに、『index_col="A行"』とするとOK。『parse_dates=True』でindexは日付データだと認識させる。

df = pd.read_csv('/content/drive/My Drive//covid_Oki.csv',index_col="日付",parse_dates=True, encoding="shift-jis")
df

次に、不要な県コードなどを削除します。👇

# @title 県コード、累計を削除
df = df.drop(columns=df.columns[[0,1,4]])
df

ここで、データの基本統計量を見ていきます👇

# @title df基本統計量
display(df.describe())
print("====================================================")
print(df.info())

2021年10月9日時点の累計感染者数は49,865人。

1日あたりの新規感染者数平均は78人。

といった感じです。

更に、感染状況をグラフにしてみていきます。👇

# @title covidグラフ表示
import matplotlib.pyplot as plt
import japanize_matplotlib
import numpy as np


#グラフの箱サイズ
df.plot(figsize=(10,8),fontsize=16,grid=True)
#凡例の位置
plt.legend(loc = 'upper left', bbox_to_anchor=(0, 1), borderaxespad=0, fontsize=16)

こうして見ると、8月～9月に感染者のピークを向かえ、10月19日現在は落ち着きを取り戻している状況ということが分かります。

しかし、新聞などで専門家も指摘するよう、残念ながら第六波が来ることは避けられないでしょう((+_+))

ここで、第六波が来ると仮定して、時系列分析で有用なfacebookのライブラリで、沖縄県のコロナ感染者数がどのように推移していくかを予測してみます。👇

#@title Prophetで時系列予測(非線形)
periods =  180#@param {type:"number"} 予測期間
from fbprophet import Prophet
import datetime
#トレンドの変化転換点を描画する関数
from fbprophet.plot import add_changepoints_to_plot

#ｘ軸(ds)に時系列日付を。ｙ軸に説明変数を。
covid = pd.DataFrame({"ds": df.index.values, "y": df["沖縄県の感染者数_1日ごとの発表数"].values})
model = Prophet(growth="logistic")
# 非線形の上限を設定
covid['cap'] = 1000 
model.fit(covid)
#空のデータフレームを用意。(periods=に何日先まで空データを作るか？入力。frepは頻度。ここではday)
future = model.make_future_dataframe(periods=periods, freq="d")
# 非線形の上限を設定
future['cap'] = 1000

forecast = model.predict(future)

# データをcsvで保存
forecast.to_csv("dat.csv", encoding="utf-8-sig")

# 予測値と、その上下限を表示
display(forecast[['ds', 'yhat', 'yhat_lower', 'yhat_upper']])
fig = model.plot(forecast)
#垂直に引かれた赤い破線はProphetがトレンドの変化を特定した場所を示し、
#赤い曲線はすべての季節性を取り除いたトレンドを示しています。
a = add_changepoints_to_plot(fig.gca(), model, forecast)
model.plot_components(forecast)
plt.show()

結果は御覧の通り、右肩上がりです。残念です。

予測値は下図のとおり。来年4月頃の1日当たり感染者数は550人程度。

役に立つかはわかりませんが、予測値のCSVファイルも上げておきます。

dat

.csv

Download CSV • 208KB

以上が、沖縄県の新型コロナ感染者数の状況です。

❷琉球銀行の決算状況

ここからは、琉球銀行の決算状況を見ていきます。

まずは同行の企業規模をホームページより抜粋。

https://www.ryugin.co.jp/corporate/about/profile/

2021年3月期　決算説明資料

https://www.ryugin.co.jp/corporate/ir/tanshin/

単体での経常収益は373億円。

連結での経常収益は572億円。

子会社では琉球リースの経常収益が171億円なので、この子会社の業績も同行の業績に大きく影響しています。

因みに、同リース会社の完全子会社化は平成27年。👇

経常収益が約38％跳ね上がった年ですね。

さて、企業規模がある程度把握できたところで、同行の連結四半期損益計算書の推移を見ていきます。👇

8399

.csv

Download CSV • 1KB

# @title 琉球銀行の四半期損益計算書（連結）
PL = pd.read_csv('/content/drive/My Drive//8399.csv',index_col="日付",parse_dates=True, encoding="shift-jis",thousands=',')
display(PL)
PL.info()

次に、沖縄県経済とコロナの関係です。これは同行作成のIRプレゼンテーション2021年3月期より引用します。👇

https://www.ryugin.co.jp/corporate/ir/presentation/

観光立県である沖縄県にとって、コロナの影響は甚大です。

有効求人倍率の低下も著しいです。

業績の要ともいえる貸出金ですが、新型コロナ関連の融資は、「事業者への資金繰り支援を更に徹底する観点から、いわゆる実質※無利子・無担保・据置最大５年の融資」などの貸出金が多く、儲けのための融資というより、政府主導の倒産させない融資が色濃くでたといった印象でしょうか。

コロナによって伸びた貸出金は、同行の収益には結びつかない。

もしくは、貸倒のリスクが高まるため、引当金を多く計上することになり、減益につながる可能性もありそうです。

これを受け同行では、フォワードルッキングな引当を開始したとのことです。一言でいえば、貸倒引当金を念のため、たくさん計上します。ということです。

以上が、琉球銀行の決算状況です。

❸予測結果『琉球銀行（8399）の2022年3月期の経常利益は前期比3倍の11,741百万円！？』

さて、ざっくりと沖縄県のコロナと琉球銀行への影響が見えてきたところで、コロナ感染者数と同行の業績にどのような関係性があるのかを見ていきます。

金融機関は基本的に日銀の金利政策に大きく影響を受けますので、財務省から日本国債の10年利回りを入手します。👇

https://www.mof.go.jp/jgbs/reference/interest_rate/index.htm

jgbcm_all

.csv

Download CSV • 10KB

# @title 財務省より国債10年利回り取得　https://www.mof.go.jp/jgbs/reference/interest_rate/index.htm
IR = pd.read_csv('/content/drive/My Drive//jgbcm_all.csv',index_col="基準日",parse_dates=True, encoding="shift-jis",thousands=',')
display(IR.head().append(IR.tail()))
IR.info()

IR.plot(figsize=(10,8),fontsize=16,grid=True)

10年国債利回りの推移👇9月末時点で0.081％です。

一般的に金利が低いときに融資した場合、その後の金利が上がってくると、収益がひっ迫します。（融資時点の調達金利は低いが、返済時の金利は上昇しているため）

以上を踏まえて、琉球銀行の業績と沖縄県のコロナ感染者数を一つのデータに纏めます。👇

# @title covid_Okiから説明変数追加
df["RB経常利益"] = PL['経常利益']
df["RB当期純利益"] = PL['当期純利益']
df["金利"] = IR["10年"]
df = df.dropna()
df

琉球銀行の経常利益（2022年3月期）を沖縄県の新型コロナ感染者数から予測する

こんにちは。乞食です。

皆さん、ワクチンの接種は終わりましたか？

相変わらず10万人あたり感染者数が全国ワーストの沖縄県の経済状況を見ていくことにしました。

（2021年10月18日までの情報を表示）

沖縄県の経済状況を把握するうえで、欠かせない金融機関の業績推移。

新型コロナでどんな影響が出ているのか？

創業70年を超える「なが～いおつきあい」をモットーとする琉球銀行を例に考察します( ﾟДﾟ)

※本考察はワクチンで熱発中の一個人の見解であり、実際の同行の業績や沖縄県の経済状況・ワクチン接種状況等には、一切関係が無いことを事前に申し上げます。

目次

❶沖縄県の新型コロナ感染状況

❷琉球銀行の決算状況

❸予測結果『琉球銀行（8399）の2022年3月期の経常利益は前期比3倍の11,741百万円！？』

❹おまけ：琉球銀行の株価予測

❶沖縄県の新型コロナ感染状況

まず最初に同県の新型コロナ感染状況のデータをNHKより入手します。

Google ドライブに同ファイルをアップロードして、colabで読み込みます。👇

次に、不要な県コードなどを削除します。👇

ここで、データの基本統計量を見ていきます👇

2021年10月9日時点の累計感染者数は49,865人。

1日あたりの新規感染者数平均は78人。

といった感じです。

更に、感染状況をグラフにしてみていきます。👇

こうして見ると、8月～9月に感染者のピークを向かえ、10月19日現在は落ち着きを取り戻している状況ということが分かります。

しかし、新聞などで専門家も指摘するよう、残念ながら第六波が来ることは避けられないでしょう((+_+))

ここで、第六波が来ると仮定して、時系列分析で有用なfacebookのライブラリで、沖縄県のコロナ感染者数がどのように推移していくかを予測してみます。👇

結果は御覧の通り、右肩上がりです。残念です。

予測値は下図のとおり。来年4月頃の1日当たり感染者数は550人程度。

役に立つかはわかりませんが、予測値のCSVファイルも上げておきます。

以上が、沖縄県の新型コロナ感染者数の状況です。

❷琉球銀行の決算状況

ここからは、琉球銀行の決算状況を見ていきます。

まずは同行の企業規模をホームページより抜粋。

最新の通期連結決算です👇同行IRより

2021年3月期 決算説明資料

単体での経常収益は373億円。

連結での経常収益は572億円。

子会社では琉球リースの経常収益が171億円なので、この子会社の業績も同行の業績に大きく影響しています。

因みに、同リース会社の完全子会社化は平成27年。👇

経常収益が約38％跳ね上がった年ですね。

さて、企業規模がある程度把握できたところで、同行の連結四半期損益計算書の推移を見ていきます。👇

次に、沖縄県経済とコロナの関係です。これは同行作成のIRプレゼンテーション2021年3月期より引用します。👇

観光立県である沖縄県にとって、コロナの影響は甚大です。

有効求人倍率の低下も著しいです。

コロナによって伸びた貸出金は、同行の収益には結びつかない。

もしくは、貸倒のリスクが高まるため、引当金を多く計上することになり、減益につながる可能性もありそうです。

これを受け同行では、フォワードルッキングな引当を開始したとのことです。一言でいえば、貸倒引当金を念のため、たくさん計上します。ということです。

以上が、琉球銀行の決算状況です。

❸予測結果『琉球銀行（8399）の2022年3月期の経常利益は前期比3倍の11,741百万円！？』

さて、ざっくりと沖縄県のコロナと琉球銀行への影響が見えてきたところで、コロナ感染者数と同行の業績にどのような関係性があるのかを見ていきます。

金融機関は基本的に日銀の金利政策に大きく影響を受けますので、財務省から日本国債の10年利回りを入手します。👇

10年国債利回りの推移👇9月末時点で0.081％です。

一般的に金利が低いときに融資した場合、その後の金利が上がってくると、収益がひっ迫します。（融資時点の調達金利は低いが、返済時の金利は上昇しているため）

以上を踏まえて、琉球銀行の業績と沖縄県のコロナ感染者数を一つのデータに纏めます。👇

2021年3月期　決算説明資料