| 计算说明(本节列出脚本中用于计算与检验的公式与步骤,便于复现): 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, |
| 1 | Down>Left | 0.002010143 | 10464 | 10051 | 0.012783 |
| 1a | Down>Left | 0.002709737 | 10441 | 10042 | 0.01235 |
| 1b | Down>Left | 0.002547857 | 10450 | 10048 | 0.012442 |
| 2c | Left>Up | 1.03203E-27 | 10256 | 8760 | 0.046304 |
| 2c | Down>Left | 2.36828E-21 | 11650 | 10256 | 0.043147 |
| 2c | Right>Up | 7.39349E-27 | 10230 | 8760 | 0.0455 |
| 2c | Down>Up | 1.57198E-91 | 11650 | 8760 | 0.089452 |
| 2c | Down>Right | 4.13824E-22 | 11650 | 10230 | 0.043952 |
| 2d | Left>Up | 8.92438E-17 | 10629 | 9461 | 0.036153 |
| 2d | Right>Up | 3.24029E-11 | 10382 | 9461 | 0.028508 |
| 2d | Down>Up | 2.81399E-17 | 10649 | 9461 | 0.036772 |
| 3c | Left>Down | 0.000647571 | 10680 | 10214 | 0.014424 |
| 3c | Up>Down | 0.000106943 | 10751 | 10214 | 0.016622 |
| 3c | Right>Down | 9.6121E-06 | 10835 | 10214 | 0.019222 |
| 3d | Up>Down | 0.004746246 | 10809 | 10430 | 0.011731 |
| 3d | Right>Down | 0.001535243 | 10863 | 10430 | 0.013403 |
| 4c | Left>Down | 1.59219E-05 | 10739 | 10137 | 0.018634 |
| 4c | Up>Down | 0.000119097 | 10668 | 10137 | 0.016436 |
| 4c | Right>Down | 1.2844E-07 | 10885 | 10137 | 0.023153 |
| 4d | Left>Down | 0.001414743 | 10783 | 10348 | 0.013465 |
| 4d | Right>Down | 0.000126972 | 10882 | 10348 | 0.016529 |
| 6 | Left>Up | 0 | 13341 | 0 | 0.412945 |
| 6 | Left>Down | 4.62544E-29 | 13341 | 11584 | 0.054384 |
| 6 | Right>Up | 0 | 13140 | 0 | 0.406723 |
| 6 | Down>Up | 0 | 11584 | 0 | 0.35856 |
| 6 | Right>Down | 2.23991E-23 | 13140 | 11584 | 0.048163 |
| 6a | Left>Up | 0 | 13344 | 0 | 0.413037 |
| 6a | Left>Down | 3.76427E-29 | 13344 | 11584 | 0.054477 |
| 6a | Right>Up | 0 | 13140 | 0 | 0.406723 |
| 6a | Down>Up | 0 | 11584 | 0 | 0.35856 |
| 6a | Right>Down | 2.23991E-23 | 13140 | 11584 | 0.048163 |
| 6b | Left>Up | 0 | 13343 | 0 | 0.413006 |
| 6b | Left>Down | 3.23137E-29 | 13343 | 11581 | 0.054539 |
| 6b | Right>Up | 0 | 13140 | 0 | 0.406723 |
| 6b | Down>Up | 0 | 11581 | 0 | 0.358467 |
| 6b | Right>Down | 1.83933E-23 | 13140 | 11581 | 0.048256 |
| 6c | Right>Left | 0.003707046 | 10803 | 10412 | 0.012103 |
| 6c | Right>Down | 0.000837115 | 10803 | 10345 | 0.014176 |
| 6d | Right>Down | 0.005978479 | 10728 | 10362 | 0.011329 |
| 7 | Up>Down | 0.001805944 | 10776 | 10352 | 0.013124 |
| 7 | Right>Down | 0.002241036 | 10766 | 10352 | 0.012815 |
| 7a | Right>Down | 0.000748757 | 10818 | 10355 | 0.014331 |
| 7b | Up>Down | 0.002541597 | 10756 | 10348 | 0.012629 |
| 7b | Right>Down | 0.00147904 | 10781 | 10348 | 0.013403 |
| 7c | Up>Down | 0.002454204 | 10828 | 10417 | 0.012722 |
| 8 | Up>Down | 0.002675187 | 10829 | 10422 | 0.012598 |
| 8a | Up>Down | 0.00256412 | 10830 | 10421 | 0.01266 |
| 8b | Up>Down | 0.002560983 | 10827 | 10418 | 0.01266 |
| 8c | Left>Down | 0.005383004 | 10711 | 10340 | 0.011484 |
| 8c | Up>Down | 0.002934374 | 10741 | 10340 | 0.012412 |
| 8c | Right>Down | 0.004975466 | 10715 | 10340 | 0.011607 |
| 10 | Up>Left | 0.000922779 | 11046 | 10587 | 0.014207 |
| 10 | Left>Down | 0.002640375 | 10587 | 10184 | 0.012474 |
| 10 | Up>Right | 0.005685401 | 11046 | 10672 | 0.011576 |
| 10 | Up>Down | 1.71075E-09 | 11046 | 10184 | 0.026682 |
| 10 | Right>Down | 0.000372661 | 10672 | 10184 | 0.015105 |
| 11 | Left>Up | 0 | 13165 | 0 | 0.407497 |
| 11 | Left>Down | 5.20987E-20 | 13165 | 11731 | 0.044387 |
| 11 | Right>Up | 0 | 13462 | 0 | 0.41669 |
| 11 | Down>Up | 0 | 11731 | 0 | 0.36311 |
| 11 | Right>Down | 5.53052E-28 | 13462 | 11731 | 0.05358 |
| 11a | Left>Up | 0 | 13391 | 0 | 0.414492 |
| 11a | Left>Down | 2.76886E-37 | 13391 | 11391 | 0.061906 |
| 11a | Right>Up | 0 | 13325 | 0 | 0.412449 |
| 11a | Down>Up | 0 | 11391 | 0 | 0.352586 |
| 11a | Right>Down | 4.43756E-35 | 13325 | 11391 | 0.059863 |