检验统计量-清华大学.pptVIP

* * 例子：粒子鉴别中常遇到的问题 带电粒子动量 粒子在每单位长度介质中的能量沉积 如果观察到的实验点在 此区域，它们是什么？ * * 检验拟合优度 为了达到此目的… * * 检验统计量与拟合优度 小的 t 数据与 H 更不符合 大的 t 数据与 H 更符合 * * P-值定义 将拟合优度用P-值表示 (也称为观察的显著水平或置信水平) 注意: 这不是 H 为真的概率。 在经典统计学上，我们从不涉及 P(H) 。 ?(H):H 的先验概率 对所有可能性进行归一化积分 而在贝叶斯统计理论中，则把 H 当成了随机变量，并利用贝叶斯定理 得到 * * P-值与假设检验 如果 H 为真，则对于连续的 ，P 在[0,1]范围内均匀分布. 在P-值定义中不涉及别的假设。 根据 P-值的定义，对 H 假设拟合优度的检验可以通过计算P-值的大小 来完成。但是应注意以下两点： P-值是一个随机变量。前面的显著水平在检验时已经被指定为常数。 如果 H 非真，则 P 的概率密度函数通常很接近零。 例如，在闪烁体叠层中根据每 一层测量的电离能损，并与利 用测量轨迹长度与不同粒子假 设估计出的电离能损进行比较， 可以计算 P-值进行 ?/? 粒子 鉴别。 Phys. Rev. D77，052003（2008） * * 例子:拟合优度检验 投 N 次硬币，观察到 nh 次头朝上的概率服从二项式分布： 假设H：硬币是公平的(朝上的 ph = 朝下的 pt = 0.5) 投 N=20 次硬币，观察到17次头朝上，则 在 t-空间中，具有相同或较少符合的区域为 取拟合优度检验统计量 * * 拟合优度检验中的问题 问题:当 P-值等于0.0026,是否意味着 H 假设是错的? P-值并不回答此问题。它只是给出与观察到的结果一样， 与 H 假设不符或者高于 H 假设( ph=pt=0.5 )的概率。 P-值=“偶然”得到如此奇怪结果的概率 一种实用的检验方法是在同样的假设下，产生同样数目的事例足够多次。检查如此奇怪的结果发生的概率是否与P-值相当。 * * 观测到一个信号的显著程度 假设观测 n 个事例，包含了 nb=已知过程(或本底)的事例数 如果?b=0.5 ，而且观测到 nobs=5 可否就此声称该迹象为新的发现？ 假设 H：?s=0，即只有本底过程出现。 ns=新过程(或信号)的事例数 也就是所谓的“无效假设” * * 观测到一个信号的显著程度(续) 对应的P-值 给出了得到这种极端结果的概率：虽然很小但不为零！ * * 潜在的问题之一 一个误导读者但又常常被使用的结果表示… 信号 即与零有两倍的标准偏差 实际想要的是：均值 ?b=0.5的泊松变量给出观测量大于 5 概率是多少？ 概率为1.7×10-4 但上面的结果表示隐含了均值为4.5，? = 2.2的高斯变量给出零或更少 的概率： “观测的信号” 如果?s 1，没有问题，即 n 服从高斯分布。 * * 潜在的问题之二 实际问题中会涉及系统误差，例如?b=0.8，则概率变为 建议给出与 ?b 合理变化相对应的P-值范围 虽然本底只增大了0.3，但却比 ?b=0.5 时小了一个量级。 * * 信号峰的显著性 假设我们不但测量了总的事例数，还测量了每个事例对应的不变质量。 频数 不变质量 信号峰 在显示信号峰的两个区间，有11个 事例，本底估计为?b=3.2 Q1：在哪寻找信号峰？ 计算任何两相连区间的P-值 Q2：信号宽度与分辨率相符吗？ 将区间增大至分辨率的几倍 Q3：信号峰是人为制造出来的吗？ 调整选择条件，分析新数据 … Qn：能发表信号峰结果吗？ 观察到的实验数据与期待本底大小的直方图，每个区间是泊松分布的一个变量。 * * 例子：?(2230)的观测 陈老师也是 作者之一 为什么不计算任何两相连区间的P-值？ 重复实验得不到先前的结果！ * * 皮尔逊的 ?2 检验 在观测的数据 与预言的期待值 之间进行比 较的检验统计量 如果 ni 是相互独立而且服从均值为?i 泊松分布，所有?i 并不太小(5)， 那么?2 将服从 N 个自由度的最小二乘概率密度函数分布。所观察的?2 可给出P-值 这里，f (z; N) 自由度为N 的最小二乘概率密度函数。 * * 皮尔逊的 ?2 检验(续) 自由度为 N 的最小二乘概率密度函数的期待值为 E(z)=N 通常以 ?2/N 来体现符合的程度 服从N-1自由度的 ?2分布(pi ntot1) 最好分别给出 ?2，N ，例如 * * 例子: ?2 检验 频数 不变质量 共20个区间， 自由度为 20 1）由于许多区间只有很少或根本没有计数，它将不服从 ?2 的概率密度函数分