06.7.32
モザイクプロットでカテゴリの交差割合を面積で表す
mosaic(df, ["年代", "利用サービス"]) が DataFrame の 2 列から直接この矩形を描き、properties 引数でセルごとの色を指定できます。読み終えると、20 代の列で SNS アプリの矩形が高く、50 代の列で EC が高いという「年代による嗜好の違い」を区切り位置のずれとして読み取り、カイ二乗検定の結果を図で補足できるようになります。直感 #
アンケートで 500 人の年代と「最もよく使うサービス」を聞き、4 × 4 のクロス集計表を作ったとします。16 個の数字から「若い人ほど SNS、年上ほど EC」と読むには比較が要ります。モザイクプロットは、まず横方向を年代で切り、30 代が最も多ければ 30 代の列を最も広く取り、次に各列の中を、その年代の人がどのサービスを選んだかの割合で縦に切ります。関連がなければ区切りはそろいます。年代とサービスに関連があると区切りの高さが列ごとにずれ、そのずれの大きさがカイ二乗統計量の大きさに対応します。
Python で描く #
次のコードは、年代ごとに異なる確率でサービスを選ぶ 500 人分のデータを作り、mosaic でクロス集計を矩形に変えて描きます。年代が上がるほど同じサービスの色を濃くする color_map を properties で渡しています。
| |

実行結果の読み方 #
列の横幅が年代の人数の割合で、30 代の列が最も広いのは p=[0.25, 0.30, 0.25, 0.20] で 30 代を最多にしたためです。各列の中の矩形の高さが、その年代の中でのサービスの割合です。20 代の列では青(SNS アプリ)が最も高く、50 代の列では赤(EC)が最も高いので、列をまたいで区切りの位置が大きくずれています。このずれが「年代とサービスの選択に関連がある」ことの視覚的な証拠で、scipy.stats.chi2_contingency のカイ二乗検定で p 値を添えれば定量的な裏付けになります。
使いどころと落とし穴 #
モザイクプロットは、アンケートの年代別嗜好や AB テストのセグメント別結果のように、2 つのカテゴリ変数の関連をカイ二乗検定とセットで見せたいときに向いています。カテゴリ数が片方でも 6 を超えるとセルが細くなってラベルが読めません。その場合は少数カテゴリを「その他」にまとめるか、グループ化棒グラフに切り替えます。3 変数以上は mosaic に 3 列渡せば描けますが、入れ子が深くなるので 2 変数ずつ複数枚に分けるのが実用的です。色は 1 つの変数(ここではサービス)に割り当て、もう 1 つの変数は軸の位置で区別するのが基本で、両方に色を使うと 16 色の区別を読者に強いることになります。独立性からの乖離をセルごとに見せたいなら、期待度数との差(残差)を色にした残差ヒートマップが代替です。