気温が上がるとビールの出荷量は増えるのか?ツイートシェアはてブLINEPocket2017年7月27日に日刊工業新聞に掲載された記事によると「夏場の平均気温が1°Cあがるごとに、ビールの販売数量は1日80万本(大瓶換算)増える」とあります。確かに、夏の暑い日に飲むビールには格別な旨さがありますし、本数などの数値的な部分はおいておいて、感覚的には違和感のない内容のように感じます。気温が上がるとビールの出荷量が増えるのかが気になる?ビールの出荷量を予想してみたいという人に向けて、アサヒスーパードライの実際の出荷量データと気象庁の気温データとを用いた単回帰分析をしてみた結果を紹介していきます。平均気温とビールの出荷量は相関があるのか?次のグラフは、東京都の平均気温とスーパードライの出荷量のデータを散布図で表わしたものです。データは2011年8月から2019年7月までの8年間(96か月分)のデータを用います。出荷量のひと箱は大瓶(633ml)×20本となります。このグラフ(回帰直線)より、気温が高い方がスーパードライの出荷量が多くなる傾向があることが分かります。相関係数は0.358と非常に弱い相関となります。非常に強い相関ここで、改めてグラフを見ると、平均気温が5°Cから10°Cの区間に出荷量が大きい月がいくつか存在すること分かります。これら月が相関係数を低くしている要因となっていそうです。気温が低いのに出荷量が多い月は何月か?次のグラフは、スーパードライの96か月分の出荷量データを時系列に折れ線グラフで表わしたものです。このグラフより、出荷量の推移が毎年同じように繰り返す傾向があることが見て取れます。その中に、毎年、飛びぬけて出荷量が増えている月があることが分かります。その月は何月なのでしょうか?12月です。12月は忘年会シーズンですね。今回の事例で言えば、アサヒグループホールディングス株式会社の決算月が12月ということも、12月の出荷量が大きく伸びている要因なのかもしれません。12月の出荷量が多い反動で、翌月1月の出荷量の落ち込みが大きいのも特徴的だと感じました。いずれにせよ、さきほどの散布図の中で、単回帰直線から離れている点は、12月である可能性が高そうです。次は、12月が散布図の中でどこにプロットされるかを確認していきます。12月は散布図のどこに位置するか?次のグラフは、東京都の平均気温とスーパードライの出荷量を12月とそれ以外の月とで、色分けをしたうえで散布図で表わしたものです。回帰直線から外れ、気温が低いのに出荷量が大きくなっていた月は、12月であることを確認できました。12月を除外した場合の相関係数は?次のグラフは、平均気温とスーパードライ出荷量のデータを12月のみ除外した形で、散布図とその回帰直線を引いたものです。12月のデータを除くことにより、相関係数が0.358から0.789に変わり、気温とビール出荷量の間には、正の強い相関があることを示す結果に変わりました。回帰式の当てはまり度合いを表す決定係数は0.128から0.623に変わりました。回帰係数(単回帰式の傾き)は19.1となりました。増えることが、過去8年のデータより推定できる、回帰式を得ることができました。この月あたり382万本は、一日あたりの本数に換算すると、12.7万本となります。冒頭に紹介した「夏場の平均気温が1°Cあがるごとに、ビールの販売数量は1日80万本(大瓶換算)増える」という記事と比較しても、12.7万本という今回の分析結果は、ひとつの銘柄の結果であることを考慮すると、大きく外れていないのかもしれませんが、今回の分析はここまでとします。まとめ今回の記事では「気温が上がるとビールの出荷量は増えるのか?」という疑問に対し、東京都の平均気温とアサヒスーパードライの出荷量データを単回帰分析した結果を紹介していきました。結果は、平均気温とビールの出荷量の関係は正の弱い相関があるという結果でした。12月は平均気温が低いにも関わらずビールの出荷量が大きく伸びている特殊な月であったことから、12月のみ除外した場合は、正の強い相関がある結果となりました。今回の結果より、12月を考慮に入れなければ、ビールの出荷量を平均気温からある程度推定できるようになりました。12月も考慮に入れて、ビールの出荷量を予測する場合はどのようにすればいいのでしょうか?ひとつの方法は、説明変数の数を増やすことです。説明変数を増やすことで、平均気温だけでは説明できない要因を予想に含めることができます。このように、複数の説明変数で目的変数を予測する手法のことを「重回帰分析」と呼びますが、「重回帰分析」の内容については、また別の記事で紹介したいと思います。求められるビールの出荷量の予測精度は?Pythonで散布図の特異なデータを色分け表示させる方法データサイエンスデータサイエンスビール単回帰分析統計学ツイートシェアはてブLINEPocketfeedly[article-banner-2]2023年スタート!インボイス制度の影響を受ける事業者とは?NISAやiDeCoだけじゃない!個人事業主やひとり社長におすすめの節税対策RECOMMENDこちらの記事も人気です。データサイエンス2021.7.25私の統計検定2級合格まで軌跡。メーカーに入社後、エンジニアとして半導体製品の企画・開発に従事。30代後半に軸ずらし転職でキャリアをシフト。本業の傍ら独学でPython&統計学を学習中。