🔖 キーワード索引

この用語と一緒に検索・参照されやすいタグ。 関連ページに飛ぶときの手がかりにも使えます。

🔖 キーワード索引 — 完全強化版

「ノルム」を理解するうえで必要なキーワードを 10 件以上提示します。 各チップから対応セクションへ移動できます。

30 秒結論 文脈 直感 数式 記号読み解き 実値計算 Python 実装 落とし穴 関連手法 関連用語 グループ教材 概念マップ

🎨 直感で掴む

「ノルム」を最初に学ぶときは、 厳密な定義よりイメージを優先しましょう。 以下は具体例・比喩を用いた直感的理解の入口です。

• L1 = タクシー距離（碁盤目状の街）／ L2 = 直線距離（鳥が飛ぶ）／ L∞ = チェス王様の距離。
• L1 はスパース解を生む。 Lasso 回帰が変数選択になる理由。
• L2 は滑らか。 Ridge 回帰が係数を一律に縮める理由。

🎨 直感で掴む — 完全強化版

ノルム を一言でいえば「47都道府県の特徴ベクトル（人口・所得）の L2 ノルム」。 47 都道府県という小さな母集団でも、 SSDSE-B-2026 の A1101 列に注目すると、 大都市圏と地方の差・人口規模に伴う相対比較など、 様々なパターンが見えてきます。

比喩でいうと、 ノルム はデータ分析の「眼鏡」のようなもの。 同じデータでも眼鏡を変えれば、 平均（中心）・分散（ばらつき）・相関（連動）・因果（影響）と、 異なる情報が浮かび上がります。 SSDSE-B-2026 を題材に、 この眼鏡をかけてみるのが本ページの狙いです。

• 具体例 1：47 都道府県の A1101 の中央値 ± 四分位を見ると、 規模感が掴める。
• 具体例 2：A4101 と組み合わせて散布図にすると、 群構造が見える。
• 具体例 3：政令指定都市の有無で群分けすると、 介入研究の素材になる。

🎨 ノルムの幾何学的解釈

単位球の形（2 次元での視覚化）

Lasso が L1、 Ridge が L2 — 幾何学的に

「損失関数の等高線」と「罰則ボール」の 初接触点が解。 L1 ボールは頂点（軸上の点）を持つため、 接触点が軸上＝係数の一部が 0 になりやすい。 L2 ボールは円なので、 接触点は一般位置＝係数は小さくなるが 0 にはならない。

凸性と最適化

$p \geq 1$ の Lp ノルムは 凸関数。 したがって、 罰則項として使うと最適化問題が凸となり、 局所最適 = 大域最適。 これが「L1, L2 がよく使われ、 L0.5 や L0 はあまり使われない」理由の一つ。

📐 定義・数式

やさしい説明で掴んだ感覚を、ここで ℓ_p ノルムの一般形 の定義式に対応づけます。下の式は左辺 $\|\mathbf{x}\|_p$ が何で決まるかを右辺で書き下したもので、Σ（合計） が現れます。それぞれの記号が何の量を指すのかは、次の「🔬 数式を言葉で読み解く」で 1 つずつ確かめてください。

📐 数式または定義 — 完全強化版

ノルム の代表的な定義式は次のとおりです。

ここで使われる記号や演算の意味は次節で言葉に翻訳します。

📐 数式を言葉で読み解く — Lp ノルムの統一視点

一般的な Lp ノルムは次のように定義されます：

$$ \|\mathbf{x}\|_p \;=\; \left( \sum_{i=1}^{n} |x_i|^p \right)^{1/p} \qquad (1 \le p \le \infty) $$

数式を言葉で読み解く：各成分の絶対値を p 乗してから合計し、 1/p 乗で元のスケールに戻す。 p=1 なら絶対値の和、 p=2 なら二乗和の平方根、 p→∞ なら最大成分。

特殊ケース

🐍 Python 実装 #C1: p を変えてノルムを連続計算

🎯 このコードでやること：東京-沖縄の差ベクトルに対し、 p=1,1.5,2,3,5,10,∞ と変えて Lp 距離を計算。 p が大きくなると最大成分が支配的になる様子を観察。

📥 入力データ：東京-沖縄の差ベクトル（生値）：差成分=[12,618,000, 12,005,000, 1,277,000]。

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16

# ── この抜粋で使うベクトルを用意します ──
import numpy as np
import pandas as pd

df = pd.read_csv('data/raw/SSDSE-B-2026.csv', encoding='cp932', skiprows=[1])
df = df[df['SSDSE-B-2026'] == 2023]        # 2023 年の 47 都道府県
cols = ['A1101', 'A1303', 'A4101']   # 総人口 / 65歳以上人口 / 出生数
vecs = {p: df[df['Prefecture'] == p][cols].values[0].astype(float)
        for p in ['東京都', '沖縄県']}

d = vecs['東京都'] - vecs['沖縄県']
print('差ベクトル:', d.astype(int))
for p in [1, 1.5, 2, 3, 5, 10]:
    norm_p = np.sum(np.abs(d)**p)**(1/p)
    print(f'  L{p:>4}: {norm_p:>14,.0f}')
print(f'  L_inf: {np.max(np.abs(d)):>14,.0f}')

📤 実行結果：

💬 結果の読み方：p を上げると L1=2590 万 → L∞=1262 万へ単調減少。 これは「全成分の合計の方が、 最大成分単独より大きい」ため。 p=2 (ユークリッド) と p=∞ (チェビシェフ) の差は 38%、 用途で選ぶ重要性が分かる。

📐 ノルムの公理 — なぜ「ノルム」と呼べるのか

関数 $\|\cdot\| : \mathbb{R}^n \to \mathbb{R}$ がノルムであるためには、 以下の 3 つの公理を満たす必要があります：

公理 1: 非負性と定値性

$$ \|\mathbf{x}\| \geq 0, \qquad \|\mathbf{x}\| = 0 \iff \mathbf{x} = \mathbf{0} $$

数式を言葉で読み解く：ノルムは常に非負であり、 0 になるのはゼロベクトルのときだけ。 これにより「距離が負」「中身があるのに大きさ 0」という不合理を防ぐ。

公理 2: 同次性（スカラー乗算）

$$ \|c \cdot \mathbf{x}\| = |c| \cdot \|\mathbf{x}\| \qquad \text{for all } c \in \mathbb{R} $$

数式を言葉で読み解く：ベクトルを 2 倍にすると、 ノルムも 2 倍に。 負のスカラー（-3 倍）なら、 ノルムは 3 倍に。 「向きは関係ない、 大きさだけ比例する」。

公理 3: 三角不等式

$$ \|\mathbf{x} + \mathbf{y}\| \leq \|\mathbf{x}\| + \|\mathbf{y}\| $$

数式を言葉で読み解く：直接行く距離は、 経由する距離以下。 これは「三角形の二辺の和は他辺より大きい」の一般化で、 距離空間の根幹。

L0 が「真のノルム」ではない理由

「L0 ノルム」=非ゼロ成分数 は同次性を満たさない：

$$ \|2\mathbf{x}\|_0 = \|\mathbf{x}\|_0 \neq 2 \|\mathbf{x}\|_0 $$

2 倍してもゼロでない成分数は変わらないため。 そのため厳密には「擬ノルム」または「カウント関数」と呼ばれます。 ただし機械学習ではスパース性指標として実用的に L0 と呼ばれます。

Lp が p < 1 でノルムにならない理由

$$ \|(1,0)\|_{0.5} = 1, \quad \|(0,1)\|_{0.5} = 1, \quad \|(1,1)\|_{0.5} = 4 $$

三角不等式 $4 \leq 1 + 1$ は成立しないため、 p < 1 のときは厳密にはノルムでありません（擬ノルム）。 ただし距離尺度として実用的に使われることがあります。

📐 行列ノルムへの拡張

ベクトルのノルムを行列に一般化したものを 行列ノルムと呼びます。 機械学習・最適化で頻出。

代表的な行列ノルム

🐍 Python 実装 #I1: SSDSE データ行列のノルム

🎯 このコードでやること：47 × 6 の SSDSE データ行列について、 各種行列ノルムを計算。

📥 入力データ：z-score 標準化された 47 × 6 の行列 Mz。

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23

import numpy as np
import pandas as pd
from sklearn.preprocessing import StandardScaler

# ── この抜粋で使う 47 × 6 の行列 Mz を用意する ──
_d = pd.read_csv('data/raw/SSDSE-B-2026.csv', encoding='cp932', skiprows=[1])
_d = _d[_d['SSDSE-B-2026'] == 2023].reset_index(drop=True)
COLS6 = ['A1101', 'A1303', 'A4101', 'A4200', 'A9101', 'L3221']
# 総人口 / 65歳以上人口 / 出生数 / 死亡数 / 婚姻件数 / 消費支出
Mz = StandardScaler().fit_transform(_d[COLS6].values.astype(float))

A = Mz   # 47 × 6 行列

print(f'Frobenius  ||A||_F  = {np.linalg.norm(A, "fro"):.3f}')
print(f'スペクトル ||A||_2  = {np.linalg.norm(A, 2):.3f}')
print(f'核ノルム   ||A||_*  = {np.linalg.norm(A, "nuc"):.3f}')
print(f'最大列和   ||A||_1  = {np.linalg.norm(A, 1):.3f}')
print(f'最大行和   ||A||_∞  = {np.linalg.norm(A, np.inf):.3f}')

# 特異値
U, S, Vt = np.linalg.svd(A, full_matrices=False)
print(f'特異値: {S.round(3)}')
print(f'核ノルム = 特異値の和 = {S.sum():.3f}')

📤 実行結果：

💬 結果の読み方：スペクトル = 最大特異値 (15.421)、 核ノルム = 全特異値の和 (24.337)。 第 1 特異値が支配的なので「データはほぼ 1 次元的」、 これは PCA で第 1 主成分の寄与率が 84.3%（= 15.421² / Σs²）という事実と整合する。

📜 ノルムの歴史

📝 「どのノルムを使うか」1 行ガイド

• 線形回帰の罰則（特徴選択したい） → L1 (Lasso)
• 線形回帰の罰則（係数を縮小したい、 多重共線性対策） → L2 (Ridge)
• テキスト類似度 → cosine
• 画像エッジ保存平滑化 → TV ノルム（L1 of gradient）
• 異常検出 → Mahalanobis
• クラスタリングのデフォルト → 標準化後の L2 (Euclidean)
• 外れ値に頑健にしたい → L1 (Manhattan) or Huber
• 最悪ケース誤差 → L∞ (Chebyshev)
• 低ランク行列復元 → 核ノルム
• adversarial robust ML → L∞ ball 内の摂動

📐 Lp ノルム一般化と p の選び方

L1, L2, L∞ は Lp ノルムの特殊例。 $$\|x\|_p = \left(\sum_i |x_i|^p\right)^{1/p}$$ で p=1, 2, ∞ それぞれ。 p を連続的に変えてみる。

このコードでやること: 東京-沖縄ベクトル差について p=0.5, 1, 1.5, 2, 3, 5, 10, ∞ で Lp ノルムを計算し、 p が増えると最大値ノルム L∞ に収束することを確認する。

📥 入力データ: 標準化済み Z（A1101/A1303/A4101） から東京と沖縄の差ベクトル。

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18

# ── この抜粋で使うデータを用意します ──
# Z は 47 都道府県 × 3 指標を標準化した行列。行の並びは SSDSE の順。
import numpy as np
import pandas as pd
from sklearn.preprocessing import StandardScaler

df = pd.read_csv('data/raw/SSDSE-B-2026.csv', encoding='cp932', skiprows=[1])
df = df[df['SSDSE-B-2026'] == 2023].reset_index(drop=True)
cols = ['A1101', 'A1303', 'A4101']   # 総人口 / 65歳以上人口 / 出生数
Z = StandardScaler().fit_transform(df[cols].astype(float).values)

diff = Z[12] - Z[46]  # 東京 - 沖縄
p_list = [0.5, 1, 1.5, 2, 3, 5, 10, 100]
result = []
for p in p_list:
    val = (np.sum(np.abs(diff)**p))**(1/p)
    result.append((p, round(val, 3)))
print(pd.DataFrame(result, columns=['p', 'Lp_norm']))

📤 実行結果:

💬 結果の読み方: p を増やすと L∞≈4.56 に収束。 これは「最大成分が支配する」という直感と一致。 一方 p < 1 では「準ノルム」となり、 三角不等式 $\|x+y\| \leq \|x\| + \|y\|$ が破れる（L0.5=39.19 は L1=13.07 と L0=3 の和を超える可能性）ため、 厳密にはノルムではない。 p < 1 はスパース性をさらに強める用途で稀に使われる。

表 p 値ごとのノルムの性質

🔬 数式を言葉で読み解く — 数式を「言葉」に翻訳

数式を眺めるだけでは身につかないので、 各記号がどんな役割を担っているかを言葉で押さえます。 「数式を音読する習慣」がつくと、 論文や教科書を読むスピードが体感で 2 倍ほど上がります。

x

n 次元ベクトル

‖·‖_p

p ノルム

p=1

L1 = 絶対値和

p=2

L2 = 二乗和の平方根

p=∞

成分の最大絶対値

🔬 数式を言葉で読み解く — 完全強化版

数式の各記号を、日本語の意味に変換します。

• X, Y — 観測されたデータ系列（SSDSE-B-2026 の A1101, A4101 列など）
• n — 標本サイズ（47 都道府県データなら 47、 多年度パネルなら 47×年数）
• β, θ, α — 推定すべき未知パラメータ
• ε, u — 観測されない誤差項・撹乱項
• E[·] — 期待値（母集団平均）
• Cov, Var — 共分散・分散
• arg min — 「最小化するパラメータ」の意味
• L (ラグ演算子) — 時系列で 1 期前の値を返す操作

🔬 SSDSE-B-2026 実データでノルムを徹底比較

ここからは SSDSE-B-2026（47 都道府県 × 112 指標）を素材に、 L1 / L2 / L∞ / cosine の各ノルムを実値で比較します。 「東京・大阪・神奈川・沖縄」の 4 県を 3 つの人口指標（A1101 総人口、 A1102 日本人人口、 A1301 年少人口）でベクトル化し、 県間距離をすべての尺度で計算します。

🧮 シナリオ：4 県の 3 次元ベクトル化

🐍 Python 実装 #A1: 各ベクトルの L1 / L2 / L∞ ノルム

🎯 このコードでやること：SSDSE-B-2026 から 2023 年・4 県を抽出し、 各都道府県ベクトルの L1 / L2 / L∞ ノルム（大きさ）を計算。 ノルムは「原点 0 からの距離」、 つまり「県全体の規模」を 1 つの数値に集約。

📥 入力データ（読み込み直後の先頭 3 行）：

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12

import pandas as pd
import numpy as np

df = pd.read_csv('data/raw/SSDSE-B-2026.csv', encoding='cp932', skiprows=[1])
df = df[df['SSDSE-B-2026']==2023].reset_index(drop=True)

cols = ['A1101','A1102','A1301']
for pref in ['東京都','大阪府','神奈川県','沖縄県']:
    v = df[df['Prefecture']==pref][cols].values[0].astype(float)
    print(f'{pref:>6}: L1={np.sum(np.abs(v)):>14,.0f}, '
          f'L2={np.linalg.norm(v):>14,.0f}, '
          f'Linf={np.max(np.abs(v)):>14,.0f}')

📤 実行結果：

💬 結果の読み方：L1 は全成分の絶対値の合計、 L2 はユークリッド距離（直線距離）、 L∞ は最大成分。 どの尺度でも東京 > 大阪 ≈ 神奈川 > 沖縄 という順序は変わらないが、 値の比は L1 で 9.2 倍、 L∞ で 9.6 倍となるなど、 尺度ごとに「規模差の強調具合」が異なる。

🐍 Python 実装 #A2: 県間距離（L1 / L2 / L∞）

🎯 このコードでやること：東京・大阪・沖縄の 3 県ペアに対し、 ベクトル差のノルムを計算。 これが「県の特徴がどれだけ違うか」の定量化。

📥 入力データ：先ほどの 3 次元ベクトル 4 県分。

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14

prefs = ['東京都','大阪府','神奈川県','沖縄県']
vecs = {p: df[df['Prefecture']==p][cols].values[0].astype(float)
        for p in prefs}

pairs = [('東京都','大阪府'),('東京都','沖縄県'),('東京都','神奈川県'),
         ('大阪府','沖縄県'),('神奈川県','沖縄県')]

print(f'{"ペア":>16} {"L1":>13} {"L2":>13} {"Linf":>13}')
for a,b in pairs:
    d = vecs[a] - vecs[b]
    print(f'{a:>5}-{b:<5}: '
          f'{np.sum(np.abs(d)):>13,.0f} '
          f'{np.linalg.norm(d):>13,.0f} '
          f'{np.max(np.abs(d)):>13,.0f}')

📤 実行結果：

💬 結果の読み方：東京-神奈川が最も近く（L2=662 万）、 東京-沖縄が最も遠い（L2=1746 万）。 大阪と神奈川は東京からほぼ等距離。 ただし、 単位を揃えずに生人口で計算しているため「規模差を距離として認識」している状態。

🐍 Python 実装 #A3: cosine 類似度（方向だけを見る）

🎯 このコードでやること：規模を無視して「成分比のパターン」だけが似ているかを cosine 類似度で見る。 値域は [-1, 1]、 1 なら完全に同じ方向、 0 なら直交、 -1 なら正反対。

📥 入力データ：4 県の 3 次元ベクトル。

1
2
3
4
5
6

def cos_sim(a, b):
    return np.dot(a, b) / (np.linalg.norm(a) * np.linalg.norm(b))

for a, b in pairs:
    s = cos_sim(vecs[a], vecs[b])
    print(f'cos({a:>5}, {b:<5}) = {s:.5f}')

📤 実行結果：

💬 結果の読み方：すべてのペアで cos ≈ 1.0。 つまり「規模を除けば、 県の人口指標は 3 次元空間でほぼ同じ方向を向いている」＝3 つの指標の比率が県によらず似通っている、 ということ。 cosine は「規模差を完全に無視」してしまうので、 用途を選ぶ必要があります。

🔬 数式を言葉で読み解く — L0/L1/L2/L∞ を SSDSE-B で全比較

「どのノルムを使うべきか」は分析目的で決まる。 SSDSE-B-2026 の 47 都道府県を「総人口・高齢人口・出生数」（A1101/A1303/A4101）の 3 次元ベクトルで表し、 東京と沖縄の距離を 4 種類のノルムで測ると、 それぞれの「距離観」の違いが鮮明になる。

このコードでやること: SSDSE-B-2026 で「東京」と「沖縄」の 3 次元ベクトル（総人口・高齢人口・出生数）を標準化後に L0/L1/L2/L∞ で測り、 4 種類の距離値を比較する。

📥 入力データ:

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15

import pandas as pd
import numpy as np
from sklearn.preprocessing import StandardScaler

df = pd.read_csv('data/raw/SSDSE-B-2026.csv', encoding='cp932', skiprows=[1])
df = df[df['SSDSE-B-2026']==2023].reset_index(drop=True)
cols = ['A1101', 'A1303', 'A4101']   # 総人口 / 高齢人口 / 出生数
X = StandardScaler().fit_transform(df[cols].values)
tokyo = X[df['Prefecture'].tolist().index('東京都')]
okina = X[df['Prefecture'].tolist().index('沖縄県')]
diff = tokyo - okina
print(f'L0 (非ゼロ要素数): {np.count_nonzero(diff)}')
print(f'L1 (Manhattan)  : {np.linalg.norm(diff, 1):.3f}')
print(f'L2 (Euclidean)  : {np.linalg.norm(diff, 2):.3f}')
print(f'L-inf (Chebyshev): {np.linalg.norm(diff, np.inf):.3f}')

📤 実行例:

💬 結果の読み方: 標準化後の東京-沖縄距離は L1=13.07, L2=7.55, L∞=4.56。 L∞ は「総人口の差」が支配しており（東京の総人口 z-score が突出）、 L1 は 3 軸の絶対差の総和、 L2 は中庸。 ノルム選択がそのままモデル感度を決める典型例。

🎯 まとめ — ノルム完全強化版

本ページでは「ノルム」を 12 必須セクション（🔖 索引 / 💡 結論 / 📍 文脈 / 🎨 直感 / 📐 数式 / 🔬 数式を言葉で読み解く / 🧮 実値計算 / 🐍 Python / ⚠️ 落とし穴 / 🌐 関連手法 / 🔗 関連用語 / 📚 グループ教材）で整理しました。 SSDSE-B-2026 で東京・大阪・神奈川・沖縄の 3 次元ベクトルを L1 / L2 / L∞ / cosine / Mahalanobis のすべての尺度で計算し、 「標準化前は人口が距離を独占（東京-沖縄 L2=1746 万）、 標準化後は意味のある相対距離（7.6 σ）」という具体結論まで到達。 統計データ分析コンペでは「単位と尺度の前処理」「標準化の有無」「次元数」を必ず明示することが重要です。

🧮 実値で計算してみる

数式だけでは「実感」が湧きにくいので、 具体的な数値で 1 度手計算してみると理解が定着します。 以下の例は、 本サイトで扱う SSDSE-B-2026 や公開教材に近い形式で用意しました。

ベクトル x = (3, 4) のノルム：

同じベクトルでもノルムにより値が違う。 用途で選ぶ。

手計算で得た値と、 後述の Python 実装で算出した値が一致することを確認すると、 「数式とコードの対応関係」がクリアに見えるようになります。

🧮 実値で計算してみる — SSDSE-B-2026 で ノルム（完全強化版）

SSDSE-B-2026（公的統計の社会・教育系データセット、 47 都道府県 × 10 年分超 × 100 以上の列）を用いて、 「ノルム」を体感します。 ファイル名は SSDSE-B-2026.csv、 読み込みは下記の Python コードで行います。

1
2
3
4
5
6
7
8

import pandas as pd

# SSDSE-B-2026 を読み込む（cp932 / Shift_JIS）
df = pd.read_csv('data/raw/SSDSE-B-2026.csv', skiprows=[1], encoding='cp932')
print(df.shape)          # (564, 112)
print(df['SSDSE-B-2026'].unique())  # 含まれる年度
latest = df[df['SSDSE-B-2026'] == df['SSDSE-B-2026'].max()].copy()
print(latest[['Prefecture', 'A1101', 'A4101']].head())

ここで使った中心列 A1101 は SSDSE-B-2026 における 47都道府県の特徴ベクトル（人口・所得）の L2 ノルム に関連する指標です。 算出例：

• 47 都道府県・最新年度の A1101 平均と標準偏差を求める
• A1101 と A4101 の相関（線形・順位）を比較する
• 群分け（例：人口上位 10 県 vs それ以外）で平均差を見る

🧮 47 都道府県のクラスタリング

標準化済みの 6 次元データを使い、 47 都道府県を k=5 のクラスターに分けます。 距離（ノルム）が異なれば、 クラスタも変わる、 という重要な性質を確認します。

🐍 Python 実装 #H1: K-means クラスタリング

🎯 このコードでやること：標準化済み 6 次元データで K-means (k=5) を実行し、 各クラスタに属する都道府県を出力。

📥 入力データ：47 × 6 の z-score 化された行列。

1
2
3
4
5
6
7
8

from sklearn.cluster import KMeans

km = KMeans(n_clusters=5, n_init=10, random_state=0).fit(Mz)
df['cluster'] = km.labels_

for c in range(5):
    members = df[df['cluster']==c]['Prefecture'].tolist()
    print(f'クラスタ {c} (n={len(members)}): {", ".join(members)}')

📤 実行結果：

💬 結果の読み方：東京都がぽつんと「自分だけのクラスタ 3」になる（z=4σ 級の外れ値）。 次が埼玉・神奈川・愛知・大阪の大都市圏（クラスタ 1、 n=4）、 続いて北海道・千葉・静岡・兵庫・福岡の準大都市圏（クラスタ 4、 n=5）。 残りは中規模県 23 件（クラスタ 0）と小規模県 14 件（クラスタ 2）に分かれる。 L2 距離を使った結果、 「人口規模」が支配的な軸になっていることが分かる。 なお random_state=0 を指定しているので、 何度実行しても同じ分割になる（クラスタ番号も固定）。

🐍 Python 実装 #H2: 距離を L1 に変えるとどうなる？

🎯 このコードでやること：K-medoids (PAM) で L1 距離を使い、 結果を K-means と比較。

📥 入力データ：同じ 47 × 6 z-score 化済み。

 1
 2
 3
 4
 5
 6
 7
 8
 9
10

from sklearn_extra.cluster import KMedoids  # pip install scikit-learn-extra

km_l1 = KMedoids(n_clusters=5, metric='manhattan', random_state=0).fit(Mz)
df['cluster_l1'] = km_l1.labels_

# L2 と L1 のクラスタ割り当ての一致度
from sklearn.metrics import adjusted_rand_score
ari = adjusted_rand_score(df['cluster'], df['cluster_l1'])
print(f'L2-K-means と L1-K-medoids の ARI = {ari:.3f}')
print('(1.0 = 完全一致、 0.0 = ランダム)')

📤 実行結果（イメージ）：

💬 結果の読み方：ARI=0.82 で「ほぼ似ているが、 一部の県が異なるクラスタに移動」。 これがノルム選択の影響。 外れ値の多いデータでは L1 (Manhattan / median ベース) が L2 (Euclidean / mean ベース) より頑健になる。

🧮 単位球の形と最適化の幾何イメージ

L1, L2, L∞ それぞれの「ノルム=1 となる点の集合（単位球）」は、 2 次元平面では下記の形になります。 ノルム最小化問題は「目的関数の等高線と単位球の最初の接点」を探す操作と見なせます。

表 単位球の形と最適化への影響

このコードでやること: 2 次元で L1, L2, L∞ の単位球の頂点座標を生成し、 それぞれの「ノルム=1」を満たす点を列挙する。

📥 入力データ: 角度パラメータ θ = 0 〜 2π を 8 点サンプリング。

1
2
3
4
5
6
7
8
9

theta = np.linspace(0, 2*np.pi, 9)[:-1]
for p in [1, 2, 100]:
    pts = []
    for t in theta:
        x, y = np.cos(t), np.sin(t)
        r = (abs(x)**p + abs(y)**p)**(1/p)
        pts.append((round(x/r, 2), round(y/r, 2)))
    print(f'L{p} 単位球サンプル: {pts}')

📤 実行結果:

💬 結果の読み方: L1 は座標の和が 1（diamond の辺上）、 L2 は半径 1 の円、 L∞ は最大成分が 1（正方形の辺上）。 Lasso がスパース解を出す理由は「diamond の頂点が軸上にあり、 そこで目的関数等高線と接しやすい」幾何に由来する。 一方 Ridge は円の表面で接するため、 係数が滑らかに縮む。

🧪 行列ノルム（フロベニウス・スペクトル・核ノルム）

ベクトルノルムを行列に拡張すると、 主に フロベニウスノルム（全要素の二乗和の平方根）、 スペクトルノルム（最大特異値）、 核ノルム（特異値の和）がある。 ニューラルネット重みの正則化や低ランク近似で重要。

表 行列ノルムの比較

このコードでやること: SSDSE-B-2026 を 47×3 の行列とみなし、 3 種類の行列ノルムを比較する。

📥 入力データ: 標準化済み Z（A1101/A1303/A4101） (47×3 行列)。

1
2
3
4
5
6

fro = np.linalg.norm(Z, 'fro')
spec = np.linalg.norm(Z, 2)
nuc = np.linalg.norm(Z, 'nuc')
print(f'Frobenius = {fro:.3f}')
print(f'Spectral  = {spec:.3f}')
print(f'Nuclear   = {nuc:.3f}')

📤 実行結果:

💬 結果の読み方: フロベニウス 11.87 ≈ スペクトル 11.83 という近接値は、 第 1 主成分（最大特異値）が全体のばらつきの大半を吸収することを意味する。 これは PCA で 1 軸目（規模軸＝総人口）が分散の大半を説明する事実と一致。 核ノルム 13.12 は特異値の総和（11.83 + 0.97 + 0.32 程度）。 核ノルム最小化は推薦システムでユーザー × アイテム行列を低ランク補完する標準手法。

🛠 実装チェックリスト（コンペ用）

統計データ分析コンペでノルムを使う前に、 以下の項目を必ず確認することで「再現性」「説明可能性」「数値安定性」を担保できます。

表 ノルム計算前後のチェックリスト

このコードでやること: SSDSE-B-2026 で 47 都道府県の L2 距離行列を計算し、 三角不等式と対称性を実際に検証する。

📥 入力データ: 標準化済み Z（A1101/A1303/A4101） (47×3)。

1
2
3
4
5
6
7
8
9
10

from scipy.spatial.distance import pdist, squareform

D = squareform(pdist(Z, metric='euclidean'))
assert np.allclose(D, D.T), '対称性 NG'
# 三角不等式: D[i,k] <= D[i,j] + D[j,k]
i, j, k = 12, 26, 46  # 東京, 大阪, 沖縄
print(f'd(東京, 沖縄) = {D[i,k]:.3f}')
print(f'd(東京, 大阪) + d(大阪, 沖縄) = {D[i,j]+D[j,k]:.3f}')
print(f'三角不等式: {D[i,k] <= D[i,j]+D[j,k]}')

📤 実行結果:

💬 結果の読み方: 東京→沖縄を直行する距離 7.55 が、 大阪経由 7.62 より短い（実空間の三角不等式）。 これが破れる距離尺度（準ノルム）を使うと k-means などのクラスタリングが暴走する。 検証はランダム 3 点抽出を 100 回程度繰り返せば十分。

🎓 歴史的経緯と命名の由来

「ノルム」という語は ラテン語 norma（直角定規）に由来し、 「基準・規範」を意味する。 数学的には 1900 年代初頭にバナッハ（S. Banach）が「バナッハ空間」を定義した際、 ベクトル空間に「長さ」を導入する公理として整備された。 L1, L2, L∞ という記号は Henri Lebesgue（ルベーグ）の積分論における Lp 空間に由来する。

表 ノルムの歴史年表

SSDSE-B-2026 を題材としたノルム分析は、 これら 130 年の理論蓄積の応用にあたる。 「単位を揃え、 軸間相関を考慮し、 適切な p を選ぶ」という基本は変わらない。

🧮 数式に値を入れて手で計算する: L1/L2/L∞ ノルム

合成ベクトル x=[3, -4, 5] で 3 ノルムを計算する。

Step 1: 各ノルム

Step 2: 大小関係

🐍 Python で再現

1
2
3
4
5

import numpy as np
x = np.array([3, -4, 5])
print(f"L1: {np.linalg.norm(x, 1)}")
print(f"L2: {np.linalg.norm(x, 2):.3f}")
print(f"L∞: {np.linalg.norm(x, np.inf)}")

📤 実行結果

💬 手計算 (Step 1) と Python 出力が完全一致。

🐍 Python 実装

公的統計（SSDSE-B-2026）を題材に、 最小限の Python コードで動作させます。 ファイルパス（data/raw/SSDSE-B-2026.csv）は自分の環境に合わせて変更してください。 まずはこのまま動かすことが理解の最短ルートです。

1
2
3
4
5
6

import numpy as np

x = np.array([3, 4])
print('L1:', np.linalg.norm(x, ord=1))   # 7
print('L2:', np.linalg.norm(x, ord=2))   # 5
print('Linf:', np.linalg.norm(x, ord=np.inf))  # 4

▶ 実行 を押せばこのページの中でそのまま動きます（ライブラリもデータも同梱済みで、 準備は要りません）。 手元の Python に移して動かすときは pip install numpy pandas scikit-learn scipy statsmodels が必要です。 読んでいるデータは data/raw/SSDSE-B-2026.csv。 日本語を含むので encoding='cp932' の指定を落とさないでください。

本サイトの全コードは 論文一覧ページ から実例として確認できます。 自分のデータで試したい場合は、 列名・欠損記号・単位の違いだけ調整すれば、 ほぼそのまま流用できます。

👣 ステップバイステップ実例

「ノルム」を初めて使う方向けに、 ハンズオン的な実行手順を整理します。 上の Python 実装と組み合わせて、 1 度自分の手でなぞってみることを強く推奨します。

1. 環境準備：このページのコードは ▶ 実行 ボタンでそのまま動くので、 まずは何も入れずに試す。 手元で動かしたくなったら Python 3.9 以上に pandas・scipy・matplotlib を入れ、 Jupyter Notebook か Google Colab を使うと試行錯誤しやすい。
2. データ取得：本サイト題材の SSDSE-B-2026 を data/raw/ に配置（または自分のデータを用意）。 列名と単位を確認。
3. 探索的に観察：df.head()、 df.describe()、 df.isna().sum() で全体像を把握。 ここで欠損や外れ値の見当を付ける。
4. 前提検証：ノルム をこのデータに当てはめてよいか（このページの「⚠️ よくある落とし穴」で挙げた 正規化忘れ・L1 と L2 の取り違え など）を確認。 NG なら別手法を検討。
5. 本処理：上のコードブロックを参考に、 関数を呼び出して値を取得。 中間出力をその都度プリントして合っているか確認。
6. 結果可視化：散布図、 棒グラフ、 ヒートマップなど、 解釈しやすい図を 1〜2 枚作る。 タイトルには結論を書く。
7. 解釈・記録：「📝 レポートでの報告」の 5 点セットに沿って Notebook に書き残す。 後の自分のために結論・限界・次の一手を明記。
8. 共有：Notebook を GitHub や Drive に置き、 関係者にレビュー依頼。 ピアレビューで穴が見つかることが多いので大事。

この 8 ステップを 1 度回すと、 「用語を読んで分かった気になる」段階から「実際に使える」段階に進めます。 知識は身体で覚えるのが結局のところ最速です。

🐍 Python 実装 — 完全強化版

scipy / pandas / scikit-learn / statsmodels を中心とした標準的な実装例です。 まず CSV を読み込み、 次に ノルム の解析を行います。

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20

import pandas as pd
import numpy as np
from scipy import stats

df = pd.read_csv('data/raw/SSDSE-B-2026.csv', skiprows=[1], encoding='cp932')
df = df[df['SSDSE-B-2026'] == df['SSDSE-B-2026'].max()].copy()

x = df['A1101'].astype(float).values
y = df['A4101'].astype(float).values

# 基本統計量
print('n            =', len(x))
print('mean(x)      =', np.mean(x))
print('std(x)       =', np.std(x, ddof=1))

# ノルム の代表的計算（用途に応じて scipy/statsmodels を切替える）
r, p = stats.pearsonr(x, y)
print(f'Pearson r = {r:.4f}, p = {p:.4g}')
rs, ps = stats.spearmanr(x, y)
print(f'Spearman rho = {rs:.4f}, p = {ps:.4g}')

用途別の追加実装：

1
2
3
4
5
6
7
8
9

# 標準化と簡易クラスタリングの例
from sklearn.preprocessing import StandardScaler
from sklearn.cluster import KMeans

X = df[['A1101', 'A4101']].astype(float).values
Xs = StandardScaler().fit_transform(X)
km = KMeans(n_clusters=4, n_init=10, random_state=0).fit(Xs)
df['cluster'] = km.labels_
print(df[['Prefecture', 'A1101', 'A4101', 'cluster']].head(10))

1
2
3
4
5
6
7

# 時系列（北海道の A1101）— 例として ARIMA 系の前処理
import statsmodels.api as sm

ts = df.sort_values('SSDSE-B-2026').groupby('SSDSE-B-2026')['A1101'].mean()
print(ts.tail())
res = sm.tsa.stattools.adfuller(ts)
print('ADF stat:', res[0], 'p:', res[1])

🐍 ノルム・距離計算の主要ライブラリ比較

🐍 Python 実装 #G1: 47 都道府県の全ペア距離行列

🎯 このコードでやること：47 都道府県の標準化済み 6 次元ベクトル間で、 全ペア L2 距離を計算し、 最も似た 5 ペア・最も異なる 5 ペアを抽出。

📥 入力データ：47 都道府県 × 6 列 (A1101, A1102, A1301, A4101, E1101, B4101)。

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19

from scipy.spatial.distance import pdist, squareform

feat = ['A1101','A1102','A1301','A4101','E1101','B4101']
M = df[feat].values.astype(float)
Mz = (M - M.mean(0)) / (M.std(0) + 1e-9)
D = squareform(pdist(Mz, 'euclidean'))

pairs_d = []
for i in range(len(df)):
    for j in range(i+1, len(df)):
        pairs_d.append((D[i,j], df['Prefecture'][i], df['Prefecture'][j]))
pairs_d.sort()

print('--- 最も似た 5 ペア ---')
for d, a, b in pairs_d[:5]:
    print(f'  {d:.3f}  {a} ⇔ {b}')
print('--- 最も異なる 5 ペア ---')
for d, a, b in pairs_d[-5:]:
    print(f'  {d:.3f}  {a} ⇔ {b}')

📤 実行結果（イメージ）：

💬 結果の読み方：最も似た 5 ペアは「近接県＋同規模県」で直感に合う。 最も異なる 5 ペアはすべて「東京 vs 中小県」。 標準化後の L2 距離が「県の総合的な類似度」を 1 つの数字で表現でき、 クラスタリングの自然な入力となる。

⚠️ よくある落とし穴

ノルム (L1/L2/L∞) を実データに適用するとき、 スケール未調整・L1 と L2 の使い分けの誤り・行列ノルムとベクトルノルムの混同 が典型的な失敗。 SSDSE-B-2026 で「人口 (百万人単位) と高齢化率 (%) の混在ベクトル」 をそのまま L2 ノルムにかけると、 ほぼ人口だけが効くため、 標準化後に計算するのが鉄則。

⚠️ 落とし穴 — 完全強化版

ノルム を実務で扱う際に踏みやすい落とし穴を 5 件挙げます。

• L1 vs L2: L1 は外れ値に頑健、L2 は微分しやすい。用途で使い分け。
• 標準化前のノルム: スケールが異なる変数を素のままノルム計算すると大きな変数が支配。
• L0 は厳密にはノルムでない: 0 ノルム（非ゼロ要素数）は擬似ノルム。
• 無限ノルム: ‖x‖∞ = max|x_i|。最悪値を見たいとき有用。
• 正則化のλ調整: ノルム正則化の強さで結果は大きく変わる。交差検証を使う。

⚠️ スケーリングの落とし穴と対策

SSDSE 47 都道府県の人口（百万単位）と高齢化率（%、 0-30 程度）を同じベクトルで扱うと、 人口が L2 ノルムを 独占します。 これが scale dependence という落とし穴。

🐍 Python 実装 #B1: 標準化（z-score）後の距離

🎯 このコードでやること：各列を平均 0・分散 1 に標準化してから L2 距離を取り直す。 これで「全成分が同じ重み」になる。

📥 入力データ：47 都道府県全体での A1101/A1102/A1301 の平均と標準偏差。

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15

all47 = df[cols].values.astype(float)
mu = all47.mean(axis=0)
sd = all47.std(axis=0, ddof=0)

# z-score 化したベクトル
zvecs = {p: (vecs[p] - mu) / sd for p in prefs}

print('z-score 化したベクトル:')
for p in prefs:
    print(f'  {p:>6}: {zvecs[p].round(2)}')

print('\n標準化後 L2 距離:')
for a, b in pairs:
    d = np.linalg.norm(zvecs[a] - zvecs[b])
    print(f'  {a:>5}-{b:<5}: {d:.3f}')

📤 実行結果：

💬 結果の読み方：標準化前は「東京-沖縄」の絶対距離が 1746 万と巨大だったが、 標準化後は 7.6 標準偏差。 「沖縄が全体的に -0.4σ、 東京が +4.1σ」と直感的に解釈できる。 PCA、 k-means、 k-NN など距離ベースの手法は標準化が ほぼ必須。

⚠️ 次元の呪い — 高次元でノルムが「平らになる」

高次元（数百〜数千次元）では、 ランダムな点同士の距離がすべて似た値になる現象 (concentration of measure)。 k-NN が機能しなくなる本質的な理由。

🐍 Python 実装 #E1: 次元 d を変えて L2 距離の分布を観察

🎯 このコードでやること：SSDSE-B-2026 の数値列を 2 / 5 / 10 / 50 / 100 次元に切り取り、 47 県の全ペア距離の 変動係数（std/mean）を計算。 小さいほど「距離が平らに」なっている。

📥 入力データ：47 都道府県 × 最大 100 指標。

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12

from scipy.spatial.distance import pdist

num_cols = df.select_dtypes('number').columns.tolist()
num_cols = [c for c in num_cols if c != 'SSDSE-B-2026']

for d in [2, 5, 10, 50, 100]:
    sel = num_cols[:d]
    M = df[sel].values.astype(float)
    M = (M - M.mean(0)) / (M.std(0) + 1e-9)
    dists = pdist(M, 'euclidean')
    cv = dists.std() / dists.mean()
    print(f'  d={d:>3}: mean={dists.mean():.3f}, std={dists.std():.3f}, CV={cv:.3f}')

📤 実行結果（イメージ）：

💬 結果の読み方：d が 2 → 100 と上がると CV（変動係数）は 0.78 → 0.28 に減少。 つまり「全ペア距離が平均の近くに集中」してくる。 d=100 ともなると k-NN の判別力が弱くなる。 対策：PCA / 特徴選択 / cosine 類似度 / Lp (p=1 や 0.5) の使用。

⚠️ ノルム使用時の追加落とし穴 5 件

🧠 マハラノビス距離との比較

L2 ノルムは「軸が独立」を前提とするが、 SSDSE-B-2026 では 総人口 と 高齢人口 が高相関（r=0.99）。 Mahalanobis 距離は共分散逆行列で「相関を考慮した距離」を計算する。

このコードでやること: 東京-沖縄について L2 と Mahalanobis 距離を比較し、 相関考慮の有無で値がどう変わるかを示す。

📥 入力データ: 標準化済み Z（A1101/A1303/A4101） と共分散行列。

1
2
3
4
5
6
7

from scipy.spatial.distance import mahalanobis

cov = np.cov(Z.T)
cov_inv = np.linalg.pinv(cov)
m_dist = mahalanobis(Z[12], Z[46], cov_inv)
print(f'L2 (Euclidean)   = {np.linalg.norm(Z[12]-Z[46]):.3f}')
print(f'Mahalanobis      = {m_dist:.3f}')

📤 実行結果:

💬 結果の読み方: Mahalanobis 距離は L2 より小さい（7.55 → 6.40）。 これは「総人口 と 高齢人口 の相関を考慮すると、 東京の異常さは見かけより穏当」と修正された結果。 つまり「人口が大きい県は高齢人口も多いのが普通」というベースラインを差し引いたうえでの距離。 異常検知では Mahalanobis が標準。

🌐 関連手法・派生

ノルムは「距離」「正則化」「異常検知」の 3 用途で姿を変える。 L1 (Manhattan) はスパース性誘導、 L2 (Euclidean) は滑らかな縮小、 L∞ (max-norm) はロバスト最適化、 Mahalanobis は共分散考慮で異常検知に使う — 同じ「距離」概念でも背後にあるノルムが違えば挙動も結論も変わる。

「どのノルムを選ぶか」は分析の前提と一蓮托生。 外れ値の影響を抑えたいなら L1、 滑らかな勾配が欲しいなら L2、 最悪値を抑えたいなら L∞、 共分散を考慮したいなら Mahalanobis。

⚖️ 似た用語との使い分け

「ノルム」と隣接する手法を、 ざっと俯瞰できる比較表として再整理します。 場面に応じてどれを採用するか、 まずは「適用条件」「仮定」「強み・弱み」の 3 軸で見比べてください。

「とりあえずデフォルト」で進めてしまうと、 適用条件外でも気付かず使い続ける事故になりがちです。 1 度「なぜこれを選んだか」を 1 文で書く習慣をつけると、 後の説明・査読でも強力な武器になります。

🛠 現場でのワークフロー例

「ノルム」を実際の分析プロジェクトに組み込むときの典型的な作業順序を示します。 教科書の例題と違って、 実データ・実業務では準備と検証に多くの時間を使うことに注意。

数学基礎 カテゴリのほかの用語と組合せて使う場面が多いため、 上記④までで終わらせず、 ⑤⑥まで丁寧に進めることが「結果が伝わる分析」の鍵です。

🔭 立場で変わる「ノルム」の見方

同じ ノルム でも、 どの立場から読むかで「まず気にすること」が入れ替わります。 下の表は、 立場ごとにこのページのどこから読むと近道かを整理したものです。

本ページはすべての立場を意識して構成されていますが、 自分の関心に応じてセクションを取捨選択して読むのが現実的です。

📜 歴史と背景

ノルム は、 統計学と計算機科学の流れの中から生まれました。 下の年表はこの分野全体の流れで、 ノルム 固有の年表ではありません。 この用語がどの時代の産物かを掴むために置いています。

歴史を知っておくと、 「なぜこの用語がこの定義になっているのか」「なぜ似た用語が複数あるのか」が腑に落ちやすくなります。 用語が生まれた動機を理解することが、 応用する力を養う近道です。

📔 ミニ用語集

「ノルム」を読み解く上で出てきた周辺の小用語を、 すぐに引けるよう 1 か所に集めました。 各説明は本ページの記述と整合しています。

x

n 次元ベクトル

‖·‖_p

p ノルム

p=1

L1 = 絶対値和

p=2

L2 = 二乗和の平方根

p=∞

成分の最大絶対値

🌐 関連手法・派生 — 完全強化版

• 内積 — $L_2$ ノルムは自分自身との内積の平方根
• 標準化 — 尺度の違う列を混ぜると、距離が大きい列だけで決まる
• 正則化 — $L_1$ ならスパース、$L_2$ なら縮小。ノルムの選択が効果を決める
• LASSO — $L_1$ ノルムを罰則に使い、係数をちょうど 0 にする
• Ridge 回帰 — $L_2$ ノルムを罰則に使い、0 に近づけるが 0 にはしない

🌐 機械学習でのノルム利用

1. 正則化（Regularization）

過学習を抑える「罰則項」として使われる：

• Lasso (L1)：$\min \|y - X\beta\|_2^2 + \lambda \|\beta\|_1$　→　係数を完全にゼロにする（スパース）
• Ridge (L2)：$\min \|y - X\beta\|_2^2 + \lambda \|\beta\|_2^2$　→　係数を小さくする（縮小）
• Elastic Net (L1+L2)：$\min \|y - X\beta\|_2^2 + \lambda_1 \|\beta\|_1 + \lambda_2 \|\beta\|_2^2$

2. 距離計算（k-NN、 k-means、 階層クラスタ）

• デフォルトはユークリッド距離（L2）
• 離散カテゴリではマンハッタン（L1）が有利
• 高次元テキストでは cosine 類似度がほぼ標準

3. ニューラルネットワーク

• weight decay = L2 正則化
• Batch Normalization は L2 ノルムを利用した正規化
• Adversarial robustness では L∞ 球内の摂動を仮定

🐍 Python 実装 #D1: Lasso vs Ridge の係数比較

🎯 このコードでやること：SSDSE-B-2026 で「人口 (A1101) を他の指標から予測」する Lasso と Ridge を学習させ、 係数のスパース性を比較。

📥 入力データ：47 都道府県 × 6 指標（A1102, A1301, A4101, E1101, C3301, B4101）。

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17

from sklearn.linear_model import Lasso, Ridge
from sklearn.preprocessing import StandardScaler

feat_cols = ['A1102','A1301','A4101','E1101','C3301','B4101']
X = df[feat_cols].values.astype(float)
y = df['A1101'].values.astype(float)

scaler = StandardScaler()
Xz = scaler.fit_transform(X)
yz = (y - y.mean()) / y.std()

lasso = Lasso(alpha=0.1).fit(Xz, yz)
ridge = Ridge(alpha=1.0).fit(Xz, yz)

print(f'{"特徴量":>6} {"Lasso":>14} {"Ridge":>14}')
for name, l, r in zip(feat_cols, lasso.coef_, ridge.coef_):
    print(f'{name:>6}: {l:>14.4f} {r:>14.4f}')

📤 実行結果（SSDSE-B-2026・2023年・47都道府県で実測）：

💬 結果の読み方：Lasso は A1102（総世帯数＝最も人口と相関する変数）にほぼ全係数（0.90）を集中させ、 他は 0 に落とした（A4101 のみ 0.0001 の無視できる残り）。 Ridge は全変数に係数を残しつつ、 値を小さく揃えた。 用途：解釈性重視なら Lasso、 予測精度重視なら Ridge。 係数は説明変数・目的変数とも標準化してあるため、 単位のない相対的な寄与度として比較できる。

🌐 Mahalanobis 距離 — 相関を考慮したノルム

共分散行列 $\Sigma$ を使い、 変数間の相関を考慮した距離：

$$ d_M(\mathbf{x}, \mathbf{y}) \;=\; \sqrt{(\mathbf{x}-\mathbf{y})^\top \Sigma^{-1} (\mathbf{x}-\mathbf{y})} $$

数式を言葉で読み解く：差ベクトルを共分散の逆行列で「相関と分散の影響を取り除いて」から L2 ノルムを取る。 ユークリッドの「相関無視・等方扱い」を補正。

🐍 Python 実装 #F1: SSDSE で Mahalanobis 距離

🎯 このコードでやること：47 都道府県の 3 次元データから共分散行列を推定し、 東京-沖縄 / 東京-神奈川の Mahalanobis 距離を計算。

📥 入力データ：4 県の 3 次元ベクトル + 47 県の共分散。

 1
 2
 3
 4
 5
 6
 7
 8
 9
10

from scipy.spatial.distance import mahalanobis

M = df[cols].values.astype(float)
cov = np.cov(M.T)
cov_inv = np.linalg.pinv(cov)

for a, b in pairs:
    dm = mahalanobis(vecs[a], vecs[b], cov_inv)
    de = np.linalg.norm(vecs[a] - vecs[b])
    print(f'{a:>5}-{b:<5}: Mahalanobis={dm:.3f}, Euclid={de:>14,.0f}')

📤 実行結果（SSDSE-B-2026・2023年・cols=A1101/A1102/A1301 で実測）：

💬 結果の読み方：Mahalanobis 距離は共分散逆行列で相関とスケールを取り除いた「単位なしの相対距離」で、 標準化の有無に依らず一定。 生の Euclid では 東京-神奈川（662 万）が最も近いのに、 Mahalanobis では 6.214 と 2 番目に大きい — 「人口規模という共通の変動」を差し引くと、 東京は隣接する神奈川からもなお際立って離れていることが分かる。 一方 大阪-沖縄 は Mahalanobis=3.785 と最小で、 相関構造に沿った平凡な差にすぎない。 異常検出（外れ値判定）で広く使われ、 χ² 分布の上限値で閾値設定可能。

🌐 信号処理・最適化でのノルム

スパース信号復元

圧縮センシング (Compressed Sensing) では、 観測 $y = A\mathbf{x}$ から元信号 $\mathbf{x}$ を復元する際、 L0 ノルムを最小化するのが理想：

$$ \min_{\mathbf{x}} \|\mathbf{x}\|_0 \quad \text{subject to} \quad A\mathbf{x} = y $$

しかし L0 は NP 困難なので、 L1 で近似（基底追跡 Basis Pursuit）：

$$ \min_{\mathbf{x}} \|\mathbf{x}\|_1 \quad \text{subject to} \quad A\mathbf{x} = y $$

Restricted Isometry Property (RIP) という条件下では、 L1 解と L0 解が一致することが証明されている (Candès-Tao 2005)。

画像復元（Total Variation）

画像 $u$ のノイズ除去で、 隣接画素差の L1 ノルム（=TV ノルム）を罰則とする：

$$ \min_u \frac{1}{2} \|u - u_{\text{noisy}}\|_2^2 + \lambda \|\nabla u\|_1 $$

これによりエッジを保ちつつ平坦領域だけ平滑化される。 ROF モデル (Rudin-Osher-Fatemi 1992) として知られる。

最適化での近接作用素

L1 罰則の近接作用素は閾値関数（soft-thresholding）：

$$ \text{prox}_{\lambda \|\cdot\|_1}(x_i) = \text{sign}(x_i) \cdot \max(|x_i| - \lambda, 0) $$

数式を言葉で読み解く：絶対値が λ 未満は 0 に押しつぶし、 それ以外は λ だけ縮める。 これが L1 のスパース性を生む技術的核心。

🔗 関連用語（前提・並列・発展）

この用語と直接つながる前後関係。 「前提」=先に知っておくと理解が早い／「並列」=同じ問題に別の角度／「発展」=応用・拡張版。 学習計画の参考にどうぞ。

• 🔗 距離 — 関連概念
• 🔗 正則化 — 主応用
• 🔗 ベクトル — 対象
• 🔗 内積 — 関連演算

✅ 実務チェックリスト

分析を提出する前に、 以下を順に確認すると見落としが大きく減ります。 教材として身につけたい「思考の型」でもあります。

• □ 「ノルム」を使う場面かを再確認したか（適用範囲外で無理に使っていないか）
• □ データの尺度・分布・サンプル数を確認したか
• □ 前提条件を満たしているか（独立性・正規性・線形性など）
• □ 欠損値・外れ値の扱い方針が明確か
• □ 計算した値だけでなく不確実性（標準誤差・信頼区間）も把握したか
• □ 結果の解釈と限界を区別したか
• □ 再現性のためにパッケージバージョン・乱数シードを記録したか
• □ 関連グループ教材で全体像を確認したか

❓ よくある質問（FAQ）

📝 レポートでの報告

「ノルム」を用いた分析を文書化する際、 以下の項目を順序立てて記述すると、 読み手が結果を追体験しやすくなります。 学術論文でも実務レポートでも基本構造は共通です。

• 使ったデータ：出典（例: SSDSE-B-2026）、 期間、 サンプル数 n、 取得日
• 前処理の方針：欠損補完、 外れ値処理、 単位統一、 変数変換（対数、 標準化など）
• 適用条件の確認：本用語の前提が満たされているかを明示的に検証した結果
• 推定値：点推定だけでなく、 標準誤差・95% 信頼区間・p 値などの不確実性も併記
• 結果の可視化：図のキャプションに n・期間・変数の単位を含める
• 解釈：「何を意味するか」を、 ドメイン知識と結びつけて記述
• 限界：「何を意味しないか」を率直に書く（相関は因果ではない、 標本の偏り、 時期の特殊性など）
• 再現性：使用パッケージのバージョン、 乱数シード、 解析コードへのリンク

この型に沿うことで、 査読・上司・将来の自分の誰が読んでも追跡できる記述になります。

🔗 関連用語（前提・並列・発展） — 完全強化版

本ページ「ノルム」に関連する用語を 前提（基礎となる概念）／並列（同カテゴリ）／発展（応用先） の 3 区分で整理します。

📥 前提となる用語

• 平均 — 中心傾向の最も基本的な要約値
• 分散 — ばらつきの定量化
• 標準偏差 — 分散の平方根、原データと同じ単位
• 確率 — 不確実性の数学的基礎
• 正規分布 — 連続確率分布の基盤

↔️ 並列・関連の用語

• 相関 — 2 変数の線形関係
• 共分散 — 相関の規格化前の量
• 線形回帰 — 関係性のモデル化
• 仮説検定 — 統計的判断の枠組み
• 信頼区間 — 推定の不確実性
• p 値 — 帰無仮説下での観測の極端さ
• 効果量 — 実質的な大きさ
• 標準誤差 — 推定量のばらつき

🚀 発展・応用の用語

• t 検定 — 平均差の検定
• ANOVA — 多群比較
• カイ二乗検定 — カテゴリ変数の独立性
• ロジスティック回帰 — 二値結果のモデリング
• クラスタリング — 教師なし分類
• 時系列解析 — 時間構造の解析
• パネルデータ — 時間と個体の二次元構造

🔗 関連用語（前提・並列・発展）— 拡張版

前提となる用語

• ベクトル — ノルムは「ベクトルの大きさ」の一般化
• 内積 — L2 ノルムは内積の平方根
• 行列 — 行列ノルムへの拡張の基盤
• 標準化 — スケール依存問題への基本対策
• 距離 — ノルムから距離が定義される

並列・同類の概念

• cosine 類似度 — 方向だけを見る尺度
• Mahalanobis 距離 — 相関を考慮したノルム
• KL ダイバージェンス — 確率分布間の「距離」（対称でない）
• ハミング距離 — 離散ベクトルの不一致数
• Jaccard 指数 — 集合の類似度
• 編集距離 — 文字列の類似度

発展・上位の概念

• 距離空間 — ノルムが定義する公理的構造
• バナッハ空間 — 完備ノルム空間
• ヒルベルト空間 — 内積から導出される L2 空間
• 正則化 — L1/L2 罰則の応用
• Lasso — L1 罰則回帰
• Ridge — L2 罰則回帰
• Elastic Net — L1+L2 の混合罰則

📋 ノルム チートシート

❓ よくある質問

Q1. なぜ L1 はスパース解を生む？

A. L1 ボール（菱形）の頂点が座標軸上にあるため、 損失関数の等高線と最初に接する点が「軸上＝係数 0」になりやすい。 L2 ボール（円）は頂点を持たないため、 係数は小さくなるがゼロにはならない。

Q2. cosine 類似度と相関係数の違いは？

A. Pearson 相関 = 中心化（平均を引いた）ベクトルの cosine 類似度。 つまり「平均を引いてから cosine を取る」かどうかの違い。

Q3. 高次元 (d > 100) では何を使えば良い？

A. 距離崩壊が起きるため、 (a) 次元削減（PCA、 UMAP）後に距離計算、 (b) cosine 類似度、 (c) Lp ノルムで p < 1 を使う（Aggarwal et al., 2001 が推奨）。

Q4. ベクトルのノルムが 0 になるのはどんな時？

A. すべての成分が 0 のときだけ。 これは「ノルムの公理：定値性」の一部。

Q5. 行列のノルムはどう定義する？

A. Frobenius ノルム（全成分の二乗和の平方根）、 スペクトルノルム（最大特異値）、 核ノルム（特異値の和）など複数あり、 用途で使い分ける。

🔗 関連用語・派生概念リンク

前提として知っておく用語

• ベクトル空間 — ノルムが定義される土台
• 内積 — L2 ノルムの基礎、 $\|x\|_2 = \sqrt{\langle x, x \rangle}$
• 標準化 — ノルム使用前の必須前処理
• 分散・共分散 — Mahalanobis 距離で活躍
• 距離公理 — 三角不等式・非負性・対称性

並列概念（同カテゴリの距離・尺度）

• マンハッタン距離 — L1 距離の別名
• ユークリッド距離 — L2 距離の別名
• チェビシェフ距離 — L∞ 距離の別名
• コサイン類似度 — 角度ベースの類似度
• マハラノビス距離 — 相関考慮の距離
• ハミング距離 — 離散ベクトル用

発展・応用先

• Lasso 回帰 — L1 ノルム正則化
• Ridge 回帰 — L2 ノルム正則化
• Elastic Net — L1 と L2 の混合
• k-means クラスタリング — L2 距離が標準
• k-NN 分類 — ノルム選択がモデル感度を左右
• SVM とカーネル — RBF カーネルは L2 距離ベース
• 勾配ノルム — 学習の収束判定に使用
• スペクトラルノルム — 行列の最大特異値

🎯 理解度チェック

1. Q1: 東京-沖縄の L2 距離が標準化前は 1294 万、 標準化後は 7.55 になる理由を 50 字以内で説明せよ。
   A: 単位の桁差で人口軸が距離を独占していたが、 標準化により 3 軸が同尺度になり真の差が見える。
2. Q2: Lasso と Ridge の幾何形状（単位球の形）はそれぞれ何か？
   A: Lasso=正八面体（角あり、 スパース解を誘発）、 Ridge=球（滑らか）。
3. Q3: L∞ ノルム = 4.559 となるとき、 3 成分のうちどれが値を決定するか？
   A: 絶対値が最大の 1 成分のみ（東京の場合も差ベクトルの場合も総人口軸）。
4. Q4: Mahalanobis 距離が L2 より小さくなる場面はどういうとき？
   A: 軸間に相関がある時。 共分散逆行列で相関方向の距離を縮める。
5. Q5: Lp ノルムで p < 1 にすると何が問題か？
   A: 三角不等式が成り立たず厳密にはノルムでなくなる（準ノルム）。 ただしスパース性は L1 より強い。

📚 さらに学ぶための入口

本ページは初学者向けの導入に重きを置いています。 もう一段深く学びたい方向けの参考方向性を以下にまとめました。 具体的な書誌情報は出典を確認の上で各自で取得してください。

• 大学教科書レベル：基礎統計・線形代数・確率論の教科書から該当章を確認すると、 本用語の理論的裏付けが押さえられます。
• 専門書・モノグラフ：本用語の名前で和書・英書を検索すると、 数百ページの体系的解説に出会えます。 1 度通読する価値あり。
• 論文・サーベイ：Google Scholar や arXiv で本用語を検索し、 引用数の多いサーベイ論文を読むと、 最新の派生・発展が見渡せます。
• 公的統計：本サイトの題材である SSDSE（教育用標準データセット）や e-Stat を使うと、 実データで手を動かしながら学べます。
• OSS ドキュメント：scikit-learn・statsmodels・PyTorch などの公式ドキュメントは、 アルゴリズム解説と実装例が揃った優良教材です。
• 本サイトの再現論文：用語がどう実問題に使われるかは、 論文一覧 から該当ジャンルを選ぶと具体例が確認できます。

🎯 このページの要点（最終確認）

「ノルム」を 1 行で言える ように整理：

• カテゴリ：数学基礎
• 何をする道具か：ノルムは、 ベクトルの「大きさ」を測る関数。 機械学習・正則化・距離計算の基礎。
• 使う前に必ず確認：適用条件、 サンプル数、 前提仮定
• 結果と一緒に必ず示す：不確実性（標準誤差・信頼区間）、 解釈、 限界
• 関連グループ教材：このページ末尾のリンクから全体像へ

🧭 次に読むなら：土台が怪しいと感じたら 平均・分散 へ戻り、 もう一段進みたければ t 検定・ANOVA へ進んでください。 この用語集は必要になった時に開く前提で作っているので、 今すぐ全部読む必要はありません。

📚 関連グループ教材

この用語の全体像を学ぶには、 横断的な教材から入るのが効率的：

• 📚 数学基礎 — 領域
• 📚 線形代数 — 関連

🗺 概念マップ — 完全強化版

• 基礎：平均／分散／標準偏差／確率
• 分布：正規分布／標本／中心極限定理
• 関係性：相関／共分散／線形回帰
• 推測：仮説検定／信頼区間／p 値／効果量
• 因果：内生性／IV／DID／平行トレンド／RDD
• 分類：クラスタリング／階層クラスタリング／Ward 法
• 時系列：ARIMA／VAR／グレンジャー因果／自己相関

🎮 触って理解する

ノルムは「ベクトルの大きさの測り方」で、 測り方は 1 通りではありません。 下の平面でベクトルの先端を原点からドラッグ（スマホはタッチ）してみてください。 同じ 1 本のベクトルでも、 L1・L2・L∞ で「大きさ」の数値が変わります。 さらに p スライダーを動かすと「ノルム＝1 の単位球」の形が菱形→円→正方形へ連続的に変わり、 なぜ L1 がスパース性を生むのか（＝角が軸上にある）が目で分かります。

🧭 何が起きているか

• L1（赤・菱形）：$|x|+|y|$。 碁盤目の街を歩く距離。 単位球の角が座標軸の上に来るのがポイント。
• L2（青・円）：$\sqrt{x^2+y^2}$。 鳥が飛ぶ直線距離。 どの方向にも角がなく滑らか。
• L∞（緑・正方形）：$\max(|x|,|y|)$。 一番大きい成分だけで決まる。
• ベクトルを長さそのままで回転させても L2 は不変ですが、 L1・L∞ は向きで値が変わります（対角方向で L1 最大・L∞ 最小）。 これを先端のドラッグで確かめてみてください。

🎯 L1 がスパース性を生む「角」の理由

Lasso（L1 正則化）は「損失の等高線」と「L1 単位球（菱形）を膨らませた罰則ボール」が最初に触れる点を解として選びます。 菱形の頂点は座標軸の上にあるため、 接触点が頂点になりやすく＝一方の係数がちょうど 0 になりやすい（スパース化）。 上の図で p を 1 に近づけると尖った角が軸上に現れる様子が、 まさにこの現象の幾何学的な正体です。 一方 L2（円）は角がないので、 係数は小さくなっても 0 にはなりにくいのです。

⚠️ よくある落とし穴

• スケール未調整のままノルムを取る：単位の大きい成分（例：人口）が数値を独占し、 小さい成分（例：失業率）が無視される。 距離・正則化の前に標準化で各成分のスケールを揃えるのが定石。
• 「L2 が常に正解」と思い込む：外れ値に強くしたい・スパース解が欲しいなら L1、 最悪成分を抑えたいなら L∞ が適切なこともある。
• p<1 を安易にノルム扱い：$p<1$ は三角不等式を満たさず厳密にはノルムではない（擬ノルム）。 上の図で p を 1 未満にすると単位球が非凸の星形になり、 最適化が難しくなる理由が視覚化されます。

🚀 発展

ノルムは正則化（Ridge=L2／Lasso=L1）、 距離計算（KNN など）、 損失関数の収束判定と、 分析の至る所で「大きさ・近さ」の物差しとして働きます。 ベクトル・コサイン類似度（向きの物差し）と対にして覚えると、 「大きさ（ノルム）」と「向き（角度）」でベクトルを二面から捉えられます。

🔗 隣接手法への橋渡し

「ノルム」は単独で完結せず、 前後の手法と組み合わさって価値が発揮される。 入力データの準備 (上流)・同目的の代替手法との比較 (並列)・結果の活用 (下流) という 3 軸で隣接領域を整理する。

• 上流 (入力の準備): ベクトル・行列 の準備、 正規化・標準化 によるスケール調整
• 並列 (同目的の代替): L1 ノルム (マンハッタン)、 L2 ノルム (ユークリッド)、 無限ノルム、 マハラノビス距離 を距離計算で交換
• 下流 (結果の活用): 正則化 (Lasso・Ridge) の罰則項、 距離 計算、 損失関数 の収束判定で活用

この上流・並列・下流の対応を地図化することで、 「ノルム」を中核に据えた分析パイプライン (データ準備 → 手法選択 → 結果の検証と展開) の全体像が見えてくる。

🌳 手法選択フロー

「ノルム」を実際に使うとき、 何をどう選ぶかを順に判断する。 上から順に答えていくと、 使うべき手法と評価の仕方が決まる。

1. 外れ値の影響をどうしたいか
   $L_2$（ユークリッド）は差を二乗するので、 大きく外れた 1 点に引きずられる。 外れ値に頑健にしたいなら $L_1$（マンハッタン）。 東京のような外れ値がある都道府県データでは、 この選択が結果を変える。
2. スパースな解が欲しいか
   $L_1$ を正則化に使うと係数がちょうど 0 になり、 変数選択の効果が出る（LASSO）。 $L_2$ は 0 に近づけるが 0 にはしない（Ridge）。 「使う変数を絞りたい」なら $L_1$。
3. 単位の違う列を混ぜていないか
   ノルムは各成分をそのまま足すので、 総人口（百万単位）と出生率（1 前後）を混ぜると 距離が人口だけで決まる。 標準化してから測る。
4. 次元はいくつか
   高次元では、 どの 2 点の距離もほぼ同じ値に集中する（次元の呪い）。 低次元の感覚で閾値を決めると機能しないので、 そのデータの距離分布を見てから決める。

ノルムは「大きさ・距離の測り方」の選択で、 $p$ を変えると何を重く見るかが変わる。 $p$ を決めずに「距離」と言っている場合、 たいてい $L_2$ が暗黙に使われている。