完全自習用テキスト / 審査・査読対策

論文の統計処理・数的根拠を
自分でゼロから説明できるようになる本

本研究(就業者2,035名の完全決定論的マイクロデータ解析)で実施した全ての数式、検定、効果量、多重比較の背後にある数学的カラクリを、統計検定3級の知識から順を追ってマスターします。

Chapter 00
  1. 第0章: 本書の目的と「統計的ディフェンス」の基本思想
  2. 第1章: 生データから分析変数への変換(データ前処理と尺度構成)
    1. 1.1 リッカート尺度を平均して良いのか?(順序尺度 vs 間隔尺度)
    2. 1.2 クロンバックの $\alpha$ 係数(内的一貫性信頼性)の直感理解
  3. 第2章: 記述統計と95%信頼区間(平均・SD・SE・t分布)
    1. 2.1 標準偏差(SD)と標準誤差(SE)の混同を解消する
    2. 2.2 なぜ正規分布の1.96ではなく「Student-t分布の厳密解」を使うのか?
  4. 第3章: なぜ $t$ 検定ではなく「二元配置分散分析(Two-way ANOVA)」なのか?
    1. 3.1 第一種の過誤の爆発(Family-wise Error Rate)
    2. 3.2 主効果(Main Effect)と交互作用効果(Interaction Effect)
  5. 第4章: 分散分析の計算ロジック(平方和・F値・p値)を紐解く
    1. 4.1 不均衡データ(Unbalanced Design)と Type II 平方和
    2. 4.2 $F$ 値とは「要因シグナル $\div$ 個人差ノイズ」の倍率である
  6. 第5章: 効果量(偏イータ二乗): 職制の説明力が規模の「9.3倍」の根拠
    1. 5.1 サンプルサイズ $N=2,035$ の罠(ビッグサンプル・バイアス)
  7. 第6章: 分散分析の前提条件と監査・ロバスト性検証
  8. 第7章: 事後検定と「現場収束現象」の数学的証明
  9. 第8章: 相関分析と「2次元直交性」の証明
  10. 日常の事例で直感的にわかる4大コラム(同じ統計分析の適用例)
    1. ☕ 新商品カフェラテの満足度(甘さ × 温度)
    2. 💻 在宅勤務の生産性(職種 × 住宅環境)
    3. 🏃 ダイエットサプリの効果(サプリ摂取 × 運動習慣)
    4. 🏥 仕事のやりがい vs 業務のストレス(1本のシーソーではない)
  11. 第9章: 実戦!口頭試問・査読ディフェンス想定問答集(Q&Aスクリプト)
  12. 学習に有用な書籍・Webサイト一覧(推薦リソース)
    1. 📚 推薦書籍(初級〜中級)
    2. 🌐 推薦Webサイト・オンライン教材
  13. 学習達成度チェックリスト

第0章: 本書の目的と「統計的ディフェンス」の基本思想

統計検定3級では「平均値」「中央値」「標準偏差」「相関係数の初歩」を学びます。これらはデータを要約する上で素晴らしい基礎ですが、学術論文や学会口頭発表の審査では、以下のような厳しい質問が飛んできます。

⚠️ 審査員・査読者が投げかけてくる典型的な疑問
  • 「大企業と中小企業を単純比較して大企業が高いと言いますが、役職の比率が違うだけの見かけの差(交絡)ではありませんか?」
  • 「3つの役職を比較するのに、$t$ 検定を3回繰り返してはいけない数学的な理由を説明できますか?」
  • 「$N=2,035$ もあれば、どんなゴミのような小さな差でも $p < 0.05$ で有意になりますよね? その差の実質的な重要性は何で証明したのですか?」

本書の目標は、数式の丸暗記ではなく、「なぜその処理をしないと誤魔化しになるのか」という必然性を理解し、上記の質問に自分の言葉で澱みなく即答できるようになることです。

Chapter 01

第1章: 生データから分析変数への変換(データ前処理と尺度構成)

1.1 リッカート尺度を平均して良いのか?(順序尺度 vs 間隔尺度)

アンケートの「1. 全くそう思わない 〜 5. 非常にそう思う」という選択肢は、厳密な数学上は「順序尺度(順序関係しかない尺度)」です。「1点から2点の距離」と「4点から5点の距離」が等しい保証はありません。

💡 学術論文における解決策: 多項目合成スコア(Composite Score)
単一の1問だけで平均を計算するのは危険ですが、同一の概念を測定する複数の項目(本研究では9項目)を加算平均して「合成スコア」を作成することで、実質的に連続的な「間隔尺度」とみなして分散分析に供することが社会科学・経営学の確立された国際標準です(Norman, 2010)。

1.2 クロンバックの $\alpha$ 係数(内的一貫性信頼性)の直感理解

9問の質問を足し合わせるとき、「質問同士がみんな同じ方向を向いて連動しているか?」を測定する係数がクロンバックの $\alpha$ 係数(Cronbach's Alpha)です。

$$\alpha = \frac{K}{K - 1} \left( 1 - \frac{\sum_{k=1}^K \sigma_{X_k}^2}{\sigma_Y^2} \right)$$

  • $K$: 項目数($K=9$)
  • $\sum \sigma_{X_k}^2$: 各項目ごとの個別の散らばりの合計
  • $\sigma_Y^2$: 9項目を合算した合成スコア全体の散らばり
✅ 本研究の実測値と判定基準($\alpha > 0.80$ で高信頼性)
  • SC18(DX推進意欲・心理的所有感): $\alpha = 0.859$(極めて良好)
  • SC15(DX組織摩擦認知・推進リスク): $\alpha = 0.818$(極めて良好)

項目全体相関(CITC)も全問で $0.43 \sim 0.68 > 0.30$ をクリアしており、足並みを乱す不良項目は1つもありません。

Chapter 02

第2章: 記述統計と95%信頼区間(平均・SD・SE・t分布)

2.1 標準偏差(SD)と標準誤差(SE)の混同を解消する

統計検定3級で最も混同されやすいのが、標準偏差(SD: Standard Deviation)と標準誤差(SE: Standard Error)です。

指標 意味 サンプル数 $n$ との関係 本研究の実測例
標準偏差 (SD) 個人のデータの散らばり $n$ を増やしてもほぼ一定 SC18全体 $SD = 0.628$
標準誤差 (SE) 標本平均のブレやすさ(精度) $SE = \frac{SD}{\sqrt{n}}$($n$ が大きいほどゼロへ収束) 大企業経営層: $SE = 0.1115$
中小部課長: $SE = 0.0265$
📊 視覚的理解 1: サンプルサイズ $n$ と SD vs SE の挙動シミュレーション
💡 グラフの解説: サンプルサイズ $n$ が 10 から 2,035 へ増えても、青線の「標準偏差(SD)」は個人の個性なので約0.63で平坦です。一方、赤線の「標準誤差(SE)」は $\sqrt{n}$ で割るため急激に縮小し、標本平均のブレが極めて小さくなる(推計が精密になる)ことが一目で分かります。

2.2 なぜ正規分布の1.96ではなく「Student-t分布の厳密解」を使うのか?

統計検定3級では「95%信頼区間は 平均 $\pm 1.96 \times SE$」と習いますが、これは母集団の分散が完全に分かっているか、標本が無限大のときの近似です。 大企業経営層($n=31$)のような現実的な標本サイズでは、分散の推定誤差を考慮して自由度 $df = n - 1$ の Student-t 分布臨界値を用いなければなりません。

$$95\% CI = \left[ \bar{x} - t_{0.025, df} \times SE, \quad \bar{x} + t_{0.025, df} \times SE \right]$$

$n=31$ の大企業経営層では $t_{crit} = 2.042$ と $1.96$ より約4.2%大きく、信頼区間が $[3.539, 3.995]$ と慎重に広めに描かれます。これが学術査読で求められる「厳密解」です。

Chapter 03

第3章: なぜ $t$ 検定ではなく「二元配置分散分析(Two-way ANOVA)」なのか?

3.1 第一種の過誤の爆発(Family-wise Error Rate)

本研究には「大企業 vs 中小企業(2水準)」×「経営層 vs 部課長 vs 一般社員(3水準)」= 6つのグループがあります。 初心者は「2グループずつ選んで $t$ 検定を15回やればいいのでは?」と考えます。しかし、15回検定を繰り返すと、誤って「有意差あり」と誤判定してしまう全体確率は以下のように跳ね上がります。

$$\alpha_{family} = 1 - (1 - 0.05)^{15} = 1 - (0.95)^{15} \approx 0.537 \quad (53.7\%!)$$

📈 視覚的理解 2: 検定回数の増加に伴う「第一種の過誤(偽陽性)」の爆発
⚠️ グラフの解説: 1回の検定なら5%だったエラー率が、検定を繰り返すたびに赤色の危険領域へ突入し、15回(6グループ総当たり)では 53.7% に達します。コイントス以下になってしまうため、全体を1回の検定で裁く「分散分析(ANOVA)」が必須になるのです。

3.2 主効果(Main Effect)と交互作用効果(Interaction Effect)

「企業規模(大 vs 中小)によるスコアの差が、組織階層(役員 vs 社員)によって違っている現象」を交互作用効果と呼びます。

🔄 視覚的理解 3: 交互作用効果の有無(平行線 vs 交差・収束)
💡 グラフの解説: 上のボタンを押して切り替えてみてください。「もし平行なら」では役職が変わっても大企業が常に+0.3点リードしますが、「本研究の実態」では一般社員層で大企業と中小企業が 3.18 でピタリと一致します。この「線の傾きの違い」を立証したのが交互作用検定($F=4.32, p=0.0134$)です。
Chapter 04

第4章: 分散分析の計算ロジック(平方和・F値・p値)を紐解く

4.1 不均衡データ(Unbalanced Design)と Type II 平方和

本調査データは、大企業経営層($n=31$)から中小部課長($n=511$)まで人数がバラバラです。 このような不均衡データでは、モデルへの投入順序で結果が変わる Type I 平方和は使えません。主効果の検出力が高く、周辺性の原理に合致した Type II 平方和 が統計学上の最適標準です。

4.2 $F$ 値とは「要因シグナル $\div$ 個人差ノイズ」の倍率である

$$F = \frac{MS_{要因}}{MS_{残差}} = \frac{\text{要因によって生じた散らばり(平均平方: MS)}}{\text{同じグループ内の個人差ノイズ(平均平方: MS_Residual)}}$$

要因 (Source) 平方和 ($SS$) 自由度 ($df$) 平均平方 ($MS$) $F$ 値 $p$ 値 判定
企業規模 (Size) 3.047 1 3.047 8.04 0.0046 有意 ($p < 0.01^{**}$)
職制階層 (Position) 29.252 2 14.626 38.58 < 0.001 極めて有意 ($p < 0.001^{***}$)
交互作用 (Size×Pos) 3.275 2 1.638 4.32 0.0134 有意 ($p < 0.05^*$)
残差 (Residual) 769.231 2,029 0.379 — — —
Chapter 05

第5章: 効果量(偏イータ二乗): 職制の説明力が規模の「9.3倍」の根拠

5.1 サンプルサイズ $N=2,035$ の罠(ビッグサンプル・バイアス)

サンプル数が2,000人を超えると、差がわずか0.01点しかなくても $p < 0.05$ で有意になってしまいます。 そこで、「その要因が結果の散らばりを何%説明できているか」という実質的な重要度を測る指標が、分散分析の効果量「偏イータ二乗(Partial $\eta^2$: $\eta_p^2$)」です。

$$\eta_p^2 = \frac{SS_{要因}}{SS_{要因} + SS_{残差}}$$

📊 視覚的理解 4: 偏イータ二乗(説明力)の圧倒的格差(職制は規模の9.3倍)
💡 グラフの解説: 企業規模主効果は $\eta_p^2 = 0.0039$(わずか0.39%の説明力)に過ぎません。一方、職制階層主効果は $\eta_p^2 = 0.0366$(3.66%)に達し、9.28倍(約9.3倍)の説明力を誇ります。「大企業か中小企業かという議論は、社内の階層差に比べれば微小である」と断言できる数学的根拠です。
Chapter 06

第6章: 分散分析の前提条件と監査・ロバスト性検証

分散分析を査読に耐えうるものにするため、3大前提条件を厳密に監査しています。

  • 正規性の仮定: 全数 $N=2,035$ 名。最小セルでも $n=31 \ge 30$ であり、中心極限定理(CLT)により標本平均の正規性は保証される。
  • 等分散性の仮定: SC18は Brown-Forsythe検定で $p = 0.672$(非有意 $\rightarrow$ 等分散成立)。
  • HC3不均一分散頑健標準誤差: 等分散性が崩れた場合を想定した感度分析でも、ANOVAの主効果・交互作用の有意性は完全に不変。
Chapter 07

第7章: 事後検定と「現場収束現象」の数学的証明

交互作用効果($F=4.32, p=0.0134$)が有意となったため、組織階層ごとにデータをスライスして企業規模の差を検定する層別単純主効果検定(Holm法補正)を実施しました。

組織階層スライス 大企業平均 中小企業平均 平均差 (大 - 中小) $t$ 値 補正後 $p$ 値 効果量 Hedges' $g$ 判定
1. 経営層 3.767 3.481 +0.286 2.45 0.0294* 0.443 (中) 大企業優位
2. 部課長 3.371 3.224 +0.147 3.37 0.0024** 0.243 (小〜中) 大企業優位
3. 一般社員 3.181 3.176 +0.005 0.12 0.9032 (ns) 0.008 (ゼロ) 完全一致・完全消失

経営層では大企業が勝っていますが、一般社員層になると差はわずか0.005点、有意確率 $p = 0.9032$(有意差なし)となります。大企業がどれほど巨額のIT投資を行っても、現場社員の意識は中小企業と全く同一水準に収束しています。これが【現場収束現象】です。

Chapter 08

第8章: 相関分析と「2次元直交性」の証明

「DXを推進したい意欲(SC18)」と「業務負荷や組織摩擦(SC15)」の関係性をピアソンの積率相関係数で検証しました。

📐 相関係数と共有分散(決定係数)の実測値
  • 相関係数: $r = -0.0901$
  • 決定係数(共有分散): $r^2 = (-0.0901)^2 = 0.008118 \quad (\mathbf{0.81\%})$

共有分散が1%未満(99%以上が無関係に独立して変動)であることは、意欲と摩擦が「表裏一体の1本のシーソー」ではなく、「完全に直交する独立した2次元」であることを数学的に証明しています。

Special Columns

日常の事例で直感的にわかる4大コラム(同じ統計分析の適用例)

本論文で使った「二元配置分散分析」「交互作用」「効果量」「2次元直交性」は、DXだけでなく、日常のビジネスや生活のあらゆる現象にそのまま応用できる強力な思考フレームワークです。

コラム 1: 交互作用の直感

☕ 新商品カフェラテの満足度(甘さ × 温度)

あるカフェが新作ラテの味覚テストを行いました。要因Aは「甘さ(微糖 vs 通常)」、要因Bは「温度(ホット vs アイス)」です。

  • ホットで飲んだ場合: 香りが立つため「微糖(甘さ控えめ)」の方が満足度が高い(M = 4.5点 vs 3.2点)。
  • アイスで飲んだ場合: 冷たいと甘みを感じにくいため、「通常(甘め)」の方が満足度が高い(M = 4.4点 vs 3.1点)。

💡 二元配置ANOVAの視点:
もし主効果だけを見て「全体平均すると甘さに差はない(微糖3.8点 ≒ 通常3.8点)」と判断してしまったら、カフェは大失敗します。「ホットかアイスかによって、甘さのベストが正反対に逆転する」という【交互作用効果】を見抜くことで、初めて「ホットは微糖、アイスは甘めで提供する」という正しい商品設計ができるのです。本論文の「経営層では大企業優位だが、一般社員では差が消える」という交互作用と全く同じ構造です。

コラム 2: 単純主効果のスライス

💻 在宅勤務の生産性(職種 × 住宅環境)

リモートワークの生産性を調査しました。要因Aは「住宅環境(個室・書斎あり vs リビング)」、要因Bは「職種(法人営業 vs エンジニア)」です。

  • 法人営業: リビングでも個室でも、生産性にあまり大きな差はない(M = 3.6点 vs 3.8点, 差はわずか)。
  • エンジニア: リビングだと家族の物音で集中が途切れ、個室(4.6点)に比べてリビング(2.8点)では生産性が大崩壊する!

💡 層別単純主効果の視点:
全体平均すると「個室の方が良い(主効果)」ですが、職種ごとにスライス(層別単純主効果検定)すると、「営業にとっては書斎有無の差は非有意だが、エンジニアにとっては絶大な差(g = 1.2)になる」という非対称な交互作用が浮き彫りになります。本論文の「一般社員層だけをスライスすると規模差が完全に消える(p=0.903)」という分析と同じ手法です。

コラム 3: 効果量 vs 有意差

🏃 ダイエットサプリの効果(サプリ摂取 × 運動習慣)

被験者10,000人を集めてダイエットサプリの臨床試験を行いました。

  • サプリを飲んだ群は、飲まなかった群に比べて体重が平均「0.05kg(50グラム)」減りました。
  • サンプル数が10,000人もいるため、$p < 0.001$ で統計的には「超高度に有意差あり」と出ました!

💡 偏イータ二乗(効果量)の視点:
50グラムの差など、コップ1杯の水を飲めば吹き飛ぶ誤差です。偏イータ二乗を計算すると $\eta_p^2 = 0.0001$(説明力0.01%)です。$p$ 値だけを見せられた消費者は「すごい効果だ!」と騙されますが、効果量を見れば「実質的には何の意味もないゴミのような差」だと即座に見破れます。本論文で「企業規模は有意だが $\eta_p^2 = 0.0039$ だから本質的ではない」と喝破したロジックそのものです。

コラム 4: 2次元直交性の直感

🏥 仕事のやりがい vs 業務のストレス(1本のシーソーではない)

「やりがいがある仕事なら、ストレスや不満はゼロになるはずだ」と多くの人が思い込んでいます(1次元のシーソー観)。しかし実際はどうでしょうか?

  • 救急救命医・起業家: 「人命を救うやりがい(超高い)」と同時に「当直や重責のストレス(超高い)」が同居する(第1象限: 高意欲・高ストレス)。
  • 窓際で放置された社員: 「やりがい(超低い)」だが「ノルマも責任もないのでストレス(超低い)」(第3象限: 低意欲・低ストレス)。

💡 概念的直交性($r \approx 0$)の視点:
やりがいとストレスは表裏一体ではなく、完全に独立した2つの軸(X軸とY軸)です。本論文で「DX推進意欲(SC18)」と「組織摩擦・抵抗(SC15)」の共有分散がわずか0.81%($r = -0.090$)であり、両利きのチェンジマネジメントが必要であると主張した根拠と全く一致します。

Chapter 09

第9章: 実戦!口頭試問・査読ディフェンス想定問答集(Q&Aスクリプト)

Question 1 なぜ $t$ 検定の総当たりではなく、二元配置ANOVAを用いたのですか?
【回答】 本研究は企業規模2水準×職制3水準の計6セルから構成されます。$t$ 検定を繰り返すと15回の検定が必要となり、第一種の過誤(偽陽性)の全体確率が $1 - (0.95)^{15} \approx 53.7\%$ まで跳ね上がります。これを防ぎ、規模主効果、職制主効果、および両者の交差パターンである交互作用効果を分離して包括的に評価するために二元配置分散分析を採用しました。
Question 2 企業規模も $p = 0.0046$ で有意なのに、なぜ「規模より階層が重要」と言い切れるのですか?
【回答】 $N=2,035$ という大標本では僅かな差でも $p$ 値が有意になる「ビッグサンプルバイアス」が生じます。実質的な説明力を表す偏イータ二乗($\eta_p^2$)を算出したところ、企業規模主効果は $\eta_p^2 = 0.0039$(説明力約0.4%)に過ぎないのに対し、職制主効果は $\eta_p^2 = 0.0366$(説明力約3.7%)と約9.3倍の説明力を示しました。さらに組織摩擦(SC15)では企業規模は非有意($p=0.160$)であり、階層こそが現象の支配的要因であることが効果量から論証されます。
Question 3 「現場収束現象」はどのような統計処理で証明されたのですか?
【回答】 二元配置ANOVAで交互作用が有意($F=4.32, p=0.0134$)であったことを受け、職制階層ごとの層別単純主効果検定(Holm法補正)を実施しました。その結果、経営層(差 $+0.286, g=0.443, p=0.029$)と部課長(差 $+0.147, g=0.243, p=0.002$)では大企業優位が確認された一方、一般社員層では大企業3.181 vs 中小3.176、差はわずか0.005、Hedges' $g = 0.008$、$p = 0.9032$ と完全に有意差が消滅しました。大企業の推進優位が現場層に移転していない事実が実証されています。
Question 4 大企業経営層が $n=31$ 名しかいませんが、信頼区間や検出力は大丈夫ですか?
【回答】 $n=31$ は中心極限定理の経験則($n \ge 30$)を満たしています。さらに95%信頼区間の算出では漸近的な正規近似(1.96)を用いず、自由度30のStudent-t分布厳密解($t_{crit} = 2.042$)を適用して慎重に推定幅を広げています。また不均一分散に対してもHC3頑健標準誤差による感度分析を行い、統計的結論の不変性を確認済みです。
Learning Resources

学習に有用な書籍・Webサイト一覧(推薦リソース)

統計検定3級から本研究レベル(分散分析・多重比較・効果量)へステップアップするために、本当に役立つ定評ある教材を厳選しました。

📚 推薦書籍(初級〜中級)

1. 『心理学のためのデータ解析テクニック』
著者: 森敏昭・吉田寿夫(福村出版)
★本書で最もおすすめの1冊。数式を最小限に抑え、「分散分析とは何をしているのか」「主効果と交互作用の違い」「なぜ単純主効果検定が必要なのか」を言葉と図で世界一わかりやすく解説しています。第4章〜第6章の分散分析パートを読むだけで、本論文のロジックが手に取るように分かります。
2. 『統計学入門』(赤本)
著者: 東京大学教養学部統計学教室 編(東京大学出版会)
★日本の統計学教科書の金字塔。第11章「分散分析法」にて、平方和の分解の数学的証明、F分布の性質が厳密に書かれています。数式レベルでANOVAの証明を理解したいときのバイブルです。
3. 『社会科学のための統計分析(新版)』
著者: 清水裕士 ほか
★アンケート調査や組織行動データの解析に特化。リッカート尺度の扱い方、クロンバックα係数による尺度信頼性、不均衡データでのType II/III平方和の選び方が実践的に学べます。
4. 『効果量とその活用―検定から効果量へ』
著者: 水本篤・竹内理(外国語教育メディア学会)
★なぜビッグサンプルで $p$ 値だけを見てはいけないのか、偏イータ二乗($\eta_p^2$)や Cohen's $d$ / Hedges' $g$ の正しい使い分けと計算方法がコンパクトにまとまっています。

🌐 推薦Webサイト・オンライン教材

1. BellCurve 統計WEB「統計学の時間」
★初級(統計検定3級相当)から中級(2級相当)まで網羅した日本最高峰の無料Webサイト。特に「中級編: 第11章〜第13章 分散分析」「二元配置分散分析」「事後検定(テューキー法)」のページは図解が秀逸で、スマホでいつでも復習できます。
2. ハンバーガー統計学にようこそ(立命館大学・向後千春先生)
★ハンバーガーショップのポテトの味や売上を題材にした超有名入門教材。分散分析の直感的なイメージを小学生でもわかるストーリー仕立てで掴むことができます。
3. アイスクリーム統計学(立命館大学・向後千春先生)
★上記ハンバーガー統計学の続編。相関分析、二元配置分散分析、交互作用効果の考え方を楽しく学べます。
Checklist

学習達成度チェックリスト

以下の項目にすべてチェックがついたら、あなたはこの論文の統計処理を完全にマスターしています!

  • クロンバック $\alpha$ 係数が「質問同士の足並みの揃い方」であり、0.80以上で合格であることを説明できる。
  • 標準偏差(個人差)と標準誤差(標本平均のブレ)の違いを式($SE = SD/\sqrt{n}$)とともに説明できる。
  • 6グループの比較で $t$ 検定を繰り返すと、エラー率が53.7%に跳ね上がる理由を説明できる。
  • 「主効果」と「交互作用効果」の違いをグラフの傾き(平行 vs 交差)を用いて説明できる。
  • $F$ 値が「要因シグナル $\div$ 個人差ノイズ」の倍率であることを説明できる。
  • 職制主効果の説明力が規模主効果の「9.3倍」である計算ロジック($\eta_p^2$ の比)を示せる。
  • 一般社員層で $p = 0.9032$(非有意)となり、大企業と中小企業が3.18で一致する「現場収束現象」を説明できる。
  • 相関係数 $r = -0.090$(共有分散0.81%)から「2次元直交性」を論理展開できる。
IDM デジタル変革調査研究グループ / 統計解析・研究方法論教育リファレンス
就業者 2,035名 完全決定論的マイクロデータ解析 (Two-way ANOVA & Post-hoc Tests)