医療機関コードは10桁の数値からなる．最初の2桁が都道府県コード，次の1桁が区分点数コード，後半7桁が医療機関番号である．医療機関番号は都道府県内では一意であるが，都道府県をまたぐと一意ではなくなる．今回はその医療機関番号から医療機関コードを取得する方法を投稿する．
全国の医療機関コードを取得する
SQL Server Management Studio 21をインストールする
SQL Server Management Studio 21 がリリースされた．ダウンロード先はこちらになる．今回インストール作業を画像で紹介する．
EXCEL VBAで複数のブックからタブ区切りテキストに一括で別名保存する
最近，Lancersを利用する機会があった．月次のファイルを10年分，計120個のブックをダウンロードし，テーブル形式に整形したブックが納品された．これらをSQL Serverにインポートしたいのだが，手動でタブ区切りテキストに別名保存するのは骨が折れる．
結局VBAで処理することにした．その際の備忘録である．
国連の人口データおよび世界銀行のGDPデータをダウンロードする
世界各国の人口推移およびGDP推移を取得したい．そんな場合は国連や世界銀行のデータを活用する．今回は国連から人口推移，世界銀行からGDP推移のデータをそれぞれ取得したので経緯を紹介する．
“国連の人口データおよび世界銀行のGDPデータをダウンロードする” の続きを読む
オーストラリアのジオコーディング
オーストラリアでの住所からのジオコーディングはGeoscape社が担当しており，そのプロジェクト名をG-NAFという．オーストラリア政府から補助金を受けており，2029年まで無料公開されることが決まっている．
データ数は1500万件以上，空間参照系はGDA94(EPSG: 4283)またはGDA2020(EPSG: 7844)である．
今回はSQL Serverでデータベースからテーブル作成，データのインポート，テーブルへの主キーと外部キーの作成までを行う．
e-Statからの社会疫学的指標を加えて熱中症搬送人員数を分析する
総務省の公開しているe-Statには社会疫学的指標が多く含まれる．今回熱中症搬送人員数に様々な指標を加えて解析してみた．
説明変数として日最高気温，日平均水蒸気圧，都道府県人口に加えて過去30日間の平均気温，エアコン保有台数，年間収入のジニ係数，光熱・水道費，実収入，第1次産業就業者比率，第2次産業就業者比率，都市公園数，都市緑化割合，自然公園割合，自然公園数，生活保護被保護人員を加えた．
すべての変数が有意であったが，VIFを見ると多重共線性を疑わせる変数もあり，良いモデルとは言えない結果となった．
日平均気温の過去30日間の移動平均をSQL Serverのウィンドウ関数を用いて計算する
熱中症の搬送人員数が月平均気温と負の相関があるとの情報を得た．普段涼しい地域ほど日最高気温の上昇に弱いという意味である．普段涼しいということを表現するには過去30日間の日平均気温の平均を取ればよいだろうと判断した．こうなるとSQL Serverのウィンドウ関数の出番である．
熱中症搬送人員数に平均風速や平均雲量は影響するか
熱中症搬送人員数に日最高気温と平均水蒸気圧が強く影響することは疑いの余地がない．他の気象条件として風速や雲量が負の影響をおよぼす可能性はないだろうか．言い換えると，風速が強ければ熱中症を発症する可能性が下がることは考えられないか，晴れよりも曇りや雨の日は熱中症を発症する可能性が下がることは考えられないかということである．
前回の記事で熱中症データベースに平均風速をインポートした．詳細は割愛するが，同様の手順で平均雲量のデータもインポートできる．
今回は説明変数として日最高気温，平均水蒸気圧に平均風速および平均雲量を加えて一般化線形モデルにて解析を行い，tree関数で可視化を試みた．
熱中症搬送人員数に都道府県人口をオフセット項として追加し一般化線形回帰分析を行う
以前の記事ではポアソン回帰モデルおよび負の二項分布モデルを用いて熱中症搬送人員数に対する日最高気温と平均水蒸気圧の回帰係数を推定した．
人口10万人あたり何名の罹患者数，というのは割り算値である．総務省消防庁の公開している熱中症搬送人員数は都道府県ごとの搬送数であり，もともと都道府県別人口が異なるのだから搬送人員数を都道府県人口で割った割合のほうが指標として適切なのではないか，という指摘は一理ある．
しかし，割り算値ではなく実数を解析すべきである．変形した観測値を統計モデルの応答変数にするのは不必要であるばかりか，誤った結果を導きかねないからである．割り算値からは確からしさの情報が失われること，変換された値の分布が不明であることから，割り算値は避けるべきである．その代わりに割り算の分母をオフセット項として線形予測子に組み込む手法がある．
熱中症搬送人員数はカウントデータであり，その期待値は集計ゾーンの集計対象人口に依存する．都道府県人口をオフセット項とすることで，都道府県の人口規模の影響を調整した回帰分析ができる．今回は都道府県人口をオフセット項として線形予測子に組み込み，一般化線形回帰分析を行ってみた．
国土数値情報の避難施設のポイントデータをSQL Serverにアップロードする
全国には避難施設が約 12 万箇所ある．それらの位置データを SQL Server にアップロードしたい．何度か試行錯誤の結果アップロードに成功したので備忘録として公開する．