総務省消防庁の公開している熱中症搬送人員数は都道府県ごとに毎日データを反復抽出しているとも言える．複数の都道府県から繰り返しデータを取るのは独立した反復ではなく，疑似反復と考えられる．このような場合，都道府県単位で差が生じると考えられ，一般化線形混合モデルを用いて回帰係数を推定する必要がある．
今回はRのglmmML()関数を用いて一般化線形混合モデルを用いた回帰係数の推定を行った．
日最高気温と熱中症の重症度別の搬送人員数を検討する
熱中症の搬送人員と最高気温との相関関係を可視化し閾値をχ二乗検定するでは日最高気温と熱中症搬送人員との相関関係を可視化した．今回は熱中症搬送人員を重症度別に日最高気温との相関関係を可視化する．
熱中症の重症度別搬送人員数を最高気温と平均湿度別にプロットする
熱中症の搬送人員と最高気温との相関関係を可視化し閾値をχ二乗検定するでは最高気温と搬送数をプロットした．今回は気象庁から湿度のデータをダウンロードし，重症度別にプロットして可視化する．
熱中症の搬送人員と最高気温との相関関係を可視化し閾値をχ二乗検定する
最高気温と熱中症の搬送人数との間に相関関係はあるだろうか．熱中症で救急搬送された人数は総務省の消防庁のサイトにある．これと気象庁のデータを結合してみた．