総務省消防庁の公開している熱中症搬送人員数は都道府県ごとに毎日データを反復抽出しているとも言える．複数の都道府県から繰り返しデータを取るのは独立した反復ではなく，疑似反復と考えられる．このような場合，都道府県単位で差が生じると考えられ，一般化線形混合モデルを用いて回帰係数を推定する必要がある．
今回はRのglmmML()関数を用いて一般化線形混合モデルを用いた回帰係数の推定を行った．
熱中症搬送人員数に都道府県人口をオフセット項として追加し一般化線形回帰分析を行う
以前の記事ではポアソン回帰モデルおよび負の二項分布モデルを用いて熱中症搬送人員数に対する日最高気温と平均水蒸気圧の回帰係数を推定した．
人口10万人あたり何名の罹患者数，というのは割り算値である．総務省消防庁の公開している熱中症搬送人員数は都道府県ごとの搬送数であり，もともと都道府県別人口が異なるのだから搬送人員数を都道府県人口で割った割合のほうが指標として適切なのではないか，という指摘は一理ある．
しかし，割り算値ではなく実数を解析すべきである．変形した観測値を統計モデルの応答変数にするのは不必要であるばかりか，誤った結果を導きかねないからである．割り算値からは確からしさの情報が失われること，変換された値の分布が不明であることから，割り算値は避けるべきである．その代わりに割り算の分母をオフセット項として線形予測子に組み込む手法がある．
熱中症搬送人員数はカウントデータであり，その期待値は集計ゾーンの集計対象人口に依存する．都道府県人口をオフセット項とすることで，都道府県の人口規模の影響を調整した回帰分析ができる．今回は都道府県人口をオフセット項として線形予測子に組み込み，一般化線形回帰分析を行ってみた．
回帰曲線の95%信頼区間をRで求める
前回の記事では熱中症搬送人員数に対する日最高気温の回帰曲線を描いた．今回はポアソン分布に基づく搬送人員数の95%信頼区間を求める．
ポアソン回帰モデルの回帰曲線をRで描く
前回の記事では熱中症搬送人員数に対する日最高気温と平均水蒸気圧の影響をポアソン回帰モデルまたは負の二項分布モデルを用いて回帰係数を求めた．今回はポアソン回帰モデルの回帰曲線をRで描く．
ポアソン回帰モデルおよび負の二項分布モデルを用いて熱中症搬送人員数に対する日最高気温と平均水蒸気圧の回帰係数を推定する
以前，熱中症搬送人員数は日最高気温と相関関係があり，片対数グラフで直線になると述べた．今回はポアソン回帰モデルおよび負の二項分布モデルで熱中症搬送人員数に対する日最高気温と平均水蒸気圧の回帰係数を推定する．
“ポアソン回帰モデルおよび負の二項分布モデルを用いて熱中症搬送人員数に対する日最高気温と平均水蒸気圧の回帰係数を推定する” の続きを読む
RおよびR Studioをインストールする
統計解析のため，SPSS version21.0をかつて使用していた．IBMによると，2018年に同バージョンのサポートは終了している．インストールはできるがアクチベーションができないため問い合わせたところ，上記のような次第であった．そのため，やむなくRを使用することにした．
血圧制御における住宅の役割：日本のスマートウェルネス住宅調査からのエビデンスのレビュー
日本におけるスマートウェルネス住宅調査のレビュー論文であり，介入試験の総括である．短期的には断熱改修により血圧が低下したことを示している．長期コホートについては今後データを収集する予定だそうであるが，エンドポイントに死亡を設定していないのが悔やまれる．英国の統計では死亡統計も検討しており，まさに全省庁を跨いだ介入プログラムの効果を検証する体制が整っている．日本も見習ってほしい．
前回は屋内温度変化と血圧変動との相関関係を見た．今回は介入試験の結果を示しており，非無作為化比較試験であるためエビデンスレベルとしては低いものの，因果関係に迫る内容である．
Role of housing in blood pressure control: a review of evidence from the Smart Wellness Housing survey in Japan Hypertens Res. 2023 Jan;46(1):9-18. doi: 10.1038/s41440-022-01060-6.から本文は読める．読みやすい英語であり，原文のまま挑戦しても問題ない．
室内温度の不安定性が冬季の家庭血圧の日内変動および日間変動に与える影響：日本全国のスマートウェルネス住宅調査から
スマートウェルネス住宅調査の結果から一つの論文を紹介する．前回の投稿では血圧値と屋内温度変化との関係を見たが，今回はベースライン調査から住宅の屋内温度と血圧変動との相関関係を見たものである．表題の論文はここから読める．
“室内温度の不安定性が冬季の家庭血圧の日内変動および日間変動に与える影響：日本全国のスマートウェルネス住宅調査から” の続きを読む
冬季における家庭血圧と屋内温度との関連の横断解析 日本全国のスマートウェルネス住宅調査
表題の原文はCross-Sectional Analysis of the Relationship Between Home Blood Pressure and Indoor Temperature in Winter: A Nationwide Smart Wellness Housing Survey in Japan DOI: 10.1161/HYPERTENSIONAHA.119.12914で読める．家が寒いと血圧が上がるという経験則をデータで示した論文である．血圧が上がれば心血管疾患リスクが上がり，死亡率も上がる．したがって冬には死者が増える，という冬季超過死亡率の上昇まで見てあればよいのだが，残念ながら死亡はエンドポイントとして見ていない．それがこの研究の限界である．
世界保健機関の住宅と健康ガイドラインの死亡率に関する参考文献
これまでの投稿においても参照されている世界保健機関の住宅と健康ガイドラインでは476もの文献が参照されている．この中から死亡率に言及した参考文献を抽出し，要訳とともに掲載する．なお，翻訳は機械翻訳によるものに手動で訂正を加えたものであり，日本語として難解な箇所もある．