Co-evolution of human and technology
総務省の公開しているe-Statには社会疫学的指標が多く含まれる．今回熱中症搬送人員数に様々な指標を加えて解析してみた．
説明変数として日最高気温，日平均水蒸気圧，都道府県人口に加えて過去30日間の平均気温，エアコン保有台数，年間収入のジニ係数，光熱・水道費，実収入，第1次産業就業者比率，第2次産業就業者比率，都市公園数，都市緑化割合，自然公園割合，自然公園数，生活保護被保護人員を加えた．
すべての変数が有意であったが，VIFを見ると多重共線性を疑わせる変数もあり，良いモデルとは言えない結果となった．
以前の記事ではポアソン回帰モデルおよび負の二項分布モデルを用いて熱中症搬送人員数に対する日最高気温と平均水蒸気圧の回帰係数を推定した．
人口10万人あたり何名の罹患者数，というのは割り算値である．総務省消防庁の公開している熱中症搬送人員数は都道府県ごとの搬送数であり，もともと都道府県別人口が異なるのだから搬送人員数を都道府県人口で割った割合のほうが指標として適切なのではないか，という指摘は一理ある．
しかし，割り算値ではなく実数を解析すべきである．変形した観測値を統計モデルの応答変数にするのは不必要であるばかりか，誤った結果を導きかねないからである．割り算値からは確からしさの情報が失われること，変換された値の分布が不明であることから，割り算値は避けるべきである．その代わりに割り算の分母をオフセット項として線形予測子に組み込む手法がある．
熱中症搬送人員数はカウントデータであり，その期待値は集計ゾーンの集計対象人口に依存する．都道府県人口をオフセット項とすることで，都道府県の人口規模の影響を調整した回帰分析ができる．今回は都道府県人口をオフセット項として線形予測子に組み込み，一般化線形回帰分析を行ってみた．
地球温暖化に伴い，台風の大型化や洪水被害の拡大が懸念されている．今回は国土交通省の公表している洪水浸水想定区域データを SQL Server にアップロードしたのでその経緯を報告する．
このデータは2022年9月時点で令和3年のものであり，毎年更新されていくものと思われる．関係各位はデータ更新が大変だろうが，頑張っていただきたい．
eStatには納税義務者数（所得割）と課税対象所得のデータが1985年から2020年まである．今回は都道府県別の自治体ごとの納税義務者数と課税対象所得の散布図を示す．
2020 年の国勢調査の結果がようやくeSTATに反映された．日本の市区町村よりも粒度の細かい小地域（町丁・字等別）の人口構成が公表されたのは2022年6月24日付である．今回はこのデータをSQL Serverに取り込んでみたい．
“eSTATの小地域（町丁・字等別）毎の年齢（５歳階級、４区分）別、男女別人口をSQL ServerにBULK INSERTする” の続きを読む
件の統計は復興庁の全国の避難者の数（所在都道府県別・所在施設別の数にあるが，このページはトップページから辿ることができず，検索からのみ到達できる．時系列でのデータは必須と思われるが，トップページから辿れるのは最新の情報のみであり，これは国民の利益に反する．
熱中症の搬送人員と最高気温との相関関係を可視化し閾値をχ二乗検定するでは最高気温と搬送人員との相関関係を解析した．今回は水蒸気圧と搬送人員との関係を可視化し，閾値を求めた．重症度別の搬送人員についての検討は日平均蒸気圧と熱中症の重症度別搬送人員との関連を調べるに記述した．
河川データを都道府県別ではなく，水域別に抽出したい．そんな動機から QGIS と EXCEL の間を行ったり来たりしている．QGISで都道府県ごとの河川データをマージするではかなり無謀なことをやった．今回はもう少し丁寧にデータを扱ってみたい．
前回の記事ではフォルダーから一括してデータを読み込む方法を紹介した．今回は月次の更新ファイルを読み込んで既存のテーブルにデータを追記する方法を紹介する．
既に読み込んだファイルは拒否したい．監査としてのワークシートが必要だ．リレーショナルデータベースならデータの一意制約から可能だが，EXCEL では自前で作る必要がある．
色々と挑戦しがいのある課題であった．
“EXCEL VBA で既存のテーブルにネットワークドライブ上の Workbook からデータを追記する” の続きを読む
Power Query が使えないと不便である．先日 EXCEL 2010 の素の環境でフォルダ内のブックをすべて開き，データを読み込む必要があったのだが，Power Query が使えなかったため，VBA でブックを開いて読み込まなければならなかった．備忘録としての記事である．
この記事はPower Query でフォルダから複数ファイルを一括インポートすると対応している．やっていることは同じだが，.xls 形式だとクエリの検証に時間がかかるため，VBA で読み込んだほうが動作は早いかもしれない．
フォルダー内のファイル一覧を取得するには FileSystemObject を使う場合と Dir() 関数を使う方法とがある．ここでは FileSystemObject を使うことにする．
データは Range オブジェクトに格納されているため，Range オブジェクトを取得するのが当面の目標となる．