计算说明(本节列出脚本中用于计算与检验的公式与步骤,便于复现):

1) 概率估计
   - 观测计数 A_i 为某类别的观测次数,总试验次数为 N(或在两类比较中为 m = A + B)。
   - 类别的估计概率: p_hat = A / N (若只比较两类,则条件概率 p_hat = A / m)。

2) Wilson 置信区间(用于单个类别概率的置信区间,双侧)
   - 给定显著性水平 alpha(双侧),先计算 z = Phi^{-1}(1 - alpha/2)。
   - 设 p_hat = A / N,
     denom = 1 + z^2 / N
     center = p_hat + z^2 / (2N)
     rad = z * sqrt( p_hat*(1-p_hat)/N + z^2/(4N^2) )
     下界 lower = (center - rad) / denom
     上界 upper = (center + rad) / denom
   - 区间被截断到 [0,1]。

3) 两类显著性比较(用于判断 A > B)
   - 只看落在 A 或 B 的样本,令 m = A + B。
   - 在原假设 H0: p_A = p_B(条件下 p = 0.5)下,A ~ Binomial(m, 0.5)。
   - 使用精确二项检验(binomtest)做单侧检验:H1: p_A > 0.5(即 A 的条件概率大于 B)。
   - 得到单侧 p_value(脚本中称为原始 p_value)。

4) 多重比较校正(组内 Bonferroni)
   - 在同一组(例如门组或某方向子组)内做所有两两比较,若组内共有 T 个两两比较,
     则每次检验的显著性阈值设为 alpha_per_test = overall_alpha / T(overall_alpha 默认为 0.05)。
   - 仅当单侧 p_value <= alpha_per_test 时,记录显著结论 A>B(或 B>A)。

5) 差值与效应量(使用总样本 N 作为分母)
   - 在两类比较中,差值定义为:
       diff_total = (A/N) - (B/N),其中 N 为总样本量。
   - 这样差值直接反映在总样本中的概率差异,便于与“误差 ≤1%(按 N 分母)”规则对齐。
   - 在结论页中,差值单元格按以下规则着色:
       <1% → 灰色;1%~2% → 黄色;≥2% → 绿色。

6) 方向比较限制(脚本实现细节)
   - 方向分为两组:直角组 = {Left, Up, Right, Down};斜角组 = {Up-Left, Up-Right, Down-Right, Down-Left}。
   - 仅在组内做两两比较,不跨组比较。

7) 概率统计表(概率统计页)
   - 每个类别的概率由原始计数除以 Total 得到:Door_X = Door_X_count / Total;Dir_Y = Dir_Y_count / Total。

8) 均值/方差/标准差(均值偏差页)
   - 组内均值(例如门组): mean = mean(p_i)(忽略为 0 的项以避免 Total=0 的影响)。
   - 以百分比形式计算偏差列: (p_i - mean) * 100。
   - 方差(%^2)使用样本方差(ddof=0),标准差为方差的平方根(%)。

本页(结论-门)说明:
 - 对每个实验,门组内共有 4 个类别,共 T = C(4,2) = 6 个两两比较。
 - 使用 overall_alpha = 0.05,组内 Bonferroni 校正后 alpha_per_test = 0.05 / 6 ≈ 0.008333。
 - 每条结论为单侧精确二项检验显著的比较,格式: 实验编号 | 结论 | p_value | A_count | B_count | 差值(pA-pB, N分母)。
实验编号 结论 p_value A_count B_count 差值(pA-pB, N分母)
1 Down>Right 0.008184426 10380 10036 0.010783
2 Left>Up 0.003421673 10327 9941 0.0121
2 Right>Up 0.007070462 10291 9941 0.010971
2 Down>Up 0.000943795 10385 9941 0.013918
2a Left>Up 0.002717373 10344 9947 0.012444
2a Right>Up 0.007801094 10292 9947 0.010814
2a Down>Up 0.001534633 10370 9947 0.013259
2b Left>Up 0.002435975 10345 9943 0.012601
2b Right>Up 0.007497275 10290 9943 0.010877
2b Down>Up 0.001138306 10379 9943 0.013667
2c Up>Left 0.002486355 10444 10041 0.012632
2c Down>Left 0.00175321 10460 10041 0.013134
2c Up>Right 1.83504E-05 10444 9855 0.018463
2c Down>Right 1.12771E-05 10460 9855 0.018964
2d Left>Down 0.005924034 10241 9883 0.011222
3c Left>Down 0.000130371 10720 10191 0.016582
3c Up>Down 0.000269956 10692 10191 0.015704
3c Right>Down 0.005314144 10560 10191 0.011567
3d Left>Down 0.002468901 10738 10329 0.012821
3d Up>Down 0.002468901 10738 10329 0.012821
4c Left>Down 0.000159213 10654 10134 0.0163
4c Up>Down 0.000430465 10615 10134 0.015077
4c Right>Down 0.000536512 10606 10134 0.014795
4d Up>Down 0.003603662 10674 10284 0.012225
4d Right>Down 0.001313033 10721 10284 0.013698
6 Left>Up 0 13118 0 0.411197
6 Left>Down 1.76452E-31 13118 11303 0.056893
6 Right>Up 0 13153 0 0.412294
6 Down>Up 0 11303 0 0.354304
6 Right>Down 1.38121E-32 13153 11303 0.05799
6a Left>Up 0 13119 0 0.411228
6a Left>Down 1.51873E-31 13119 11302 0.056956
6a Right>Up 0 13153 0 0.412294
6a Down>Up 0 11302 0 0.354272
6a Right>Down 1.27605E-32 13153 11302 0.058021
6b Left>Up 0 13118 0 0.411197
6b Left>Down 1.76452E-31 13118 11303 0.056893
6b Right>Up 0 13153 0 0.412294
6b Down>Up 0 11303 0 0.354304
6b Right>Down 1.38121E-32 13153 11303 0.05799
8c Right>Left 0.006628156 10744 10383 0.011316
10 Left>Down 0.000160121 10622 10103 0.016269
10 Up>Down 1.63949E-07 10843 10103 0.023196
10 Right>Down 0.001150069 10542 10103 0.013761
11 Left>Up 0 13050 0 0.409065
11 Left>Down 3.52935E-22 13050 11542 0.04727
11 Right>Up 0 13224 0 0.414519
11 Down>Up 0 11542 0 0.361795
11 Right>Down 5.93384E-27 13224 11542 0.052724
11a Left>Up 0 13217 0 0.4143
11a Left>Down 1.21613E-36 13217 11248 0.06172
11a Right>Up 0 13161 0 0.412545
11a Down>Up 0 11248 0 0.35258
11a Right>Down 9.01144E-35 13161 11248 0.059965