SAS的卡方检验

卡方检验基本思想

在分类资料统计分析中我们常会遇到这样的资料，如两组大白鼠在不同致癌剂作用下的发癌率如下表，问两组发癌率有无差别？

处理发癌数未发癌数合计发癌率%

甲组52197173.24

乙组3934292.86

合计912211380.33 52 19 39 3 是表中最基本的数据，因此上表资料又被称之为四格表资料。卡方检验的统计量是卡方值，它是每个格子实际频数A与理论频数T差值平方与理论频数之比的累计和。每个格子中的理论频数T是在假定两组的发癌率相等（均等于两组合计的发癌率）的情况下计算出来的，如第一行第一列的理论频数为71*91/113=57.18，故卡方值越大，说明实际频数与理论频数的差别越明显，两组发癌率不同的可能性越大。 利用统计学软件分析结果如下： data kafang; input row column number @@; cards; 1 1 52 1 2 19 2 1 39 2 2 3 ; run; proc freq; tables row*column/chisq; weight number; run;

统计量自由度值概率

卡方16.47770.0109（有统计学意义）

似然比卡方17.31010.0069

连续校正卡方15.28680.0215

Mantel-Haenszel 卡方16.42030.0113

Phi 系数

-0.2394

列联系数

0.2328

Cramer 的 V

-0.2394

二联表的卡方检验方法

假设有两个分类变量X和Y，它们的值域分另为{x1, x2}和{y1, y2}，其样本频数列联表为：

y1y2总计

x1aba+b

x2cdc+d

总计a+cb+da+b+c+d 若要推断的论述为H1：“X与Y有关系”，可以利用独立性检验来考察两个变量是否有关系，并且能较精确地给出这种判断的可靠程度。具体的做法是，由表中的数据算出随机变量K^2的值（即K的平方） K^2 = n (ad - bc) ^ 2 / [(a+b)(c+d)(a+c)(b+d)] 其中n=a+b+c+d为样本容量 K^2的值越大，说明“X与Y有关系”成立的可能性越大。 当表中数据a，b，c，d都不小于5时，可以查阅下表来确定结论“X与Y有关系”的可信程度：

P(K^2≥k)0.500.400.250.150.10

k0.4550.7081.3232.0722.706

P(K^2≥k)0.050.0250.0100.0050.001

k3.8415.0246.6357.87910.828 例如，当“X与Y有关系”的K^2变量的值为6.109，根据表格，因为5.024≤6.109＜6.635，所以“X与Y有关系”成立的概率为1-0.025=0.975，即97.5%。

四格表资料的卡方检验

四格表资料的卡方检验用于进行两个率或两个构成比的比较。 1. 专用公式：

若四格表资料四个格子的频数分别为a，b，c，d，则四格表资料卡方检验的卡方值=n(ad-bc)^2/(a+b)(c+d)(a+c)(b+d)， 自由度v=（行数-1）（列数-1） 2. 应用条件： 要求样本含量应大于40且每个格子中的理论频数不应小于5。当样本含量大于40但理论频数有小于5的情况时卡方值需要校正，当样本含量小于40时只能用确切概率法计算概率。

行X列表资料的卡方检验

行X列表资料的卡方检验用于多个率或多个构成比的比较。 1. 专用公式： r行c列表资料卡方检验的卡方值=n[(A11/n1n1+A12/n1n2+...+Arc/nrnc)-1] 2. 应用条件： 要求每个格子中的理论频数T均大于5或1列联表资料的卡方检验

同一组对象，观察每一个个体对两种分类方法的表现，结果构成双向交叉排列的统计表就是列联表。 1. R*C 列联表的卡方检验： R*C 列联表的卡方检验用于R*C列联表的相关分析，卡方值的计算和检验过程与行X列表资料的卡方检验相同。 2. 2*2列联表的卡方检验： 2*2列联表的卡方检验又称配对记数资料或配对四格表资料的卡方检验，根据卡方值计算公式的不同，可以达到不同的目的。当用一般四格表的卡方检验计算时，卡方值=n（ad-bc)^2/[(a+b)(c+d)(a+c)(b+d)]，此时用于进行配对四格表的相关分析，如考察两种检验方法的结果有无关系；当卡方值=（|b-c|-1）2/(b+c)时，此时卡方检验用来进行四格表的差异检验，如考察两种检验方法的检出率有无差别。 列

联表卡方检验应用中的注意事项同R*C表的卡方检验相同。 卡方检验就是统计样本的实际观测值与理论推断值之间的偏离程度，实际观测值与理论推断值之间的偏离程度就决定卡方值的大小，卡方值越大，越不符合，偏差越小，卡方值就越小，越趋于符合，若量值完全相等时，卡方值就为0，表明理论值完全符合。 一年以来,我们的货运工作在车务段和车站领导的大力支持和全体货运人员的共同努力下，主要开展了以下几个方面的工作，现将2007年度货运工作总结如下：

一、提高装载质量，保一路平安。

安全是铁路运输永恒的主题，是铁路其它各项工作的基础，只有保住了安全，铁路运输才会有生存和持续发展的可能。为了能让全体货运人员牢固树

立安全第一的思想，我们利用点名会、规章提问、业务学习以及他站的一些事故教训，向大家宣讲货运安全的重要性，使得每个人都明白无论是谁只要砸了安全生产这个大锅，也就同时砸掉了自己的饭碗，搞好安全生产，首先是对自己负责。

今年4月18日铁路实施第六次大提速，为了确保大提速以后的货运安全，车站首先选货运派骨干人员参加上级组织的培训班，提前组织全体人员学习大提速的有关规定，认真落实装载加固的要求，由货运主任、货运值班员带队深入专用线现场指导，确保了大提速期间的货装安全。

车站办理的货运业务中，机械厂发运的雷达车，每个月虽然量不大，但因为涉及到军事运输及国防安全，为了确保其运输安全优质的完成军交任务，货运室每次在装车时，均严格按照定型图及有关要求装车，不但外勤货运员要检查，货运主任、值班员及安全员都要到现场复查，合格以后才能放行。一年以来，累计装军用车辆150余车，未发生一起整

倒装，装载质量受到好评。

对于普通车辆的装运，对苫盖蓬布的货车，严格按照铁路局“装载加固补充规定”的要求，每车均加盖规定的蓬布绳网。对葛洲坝水泥厂自备蓬布的腰边绳要求企业全部予以更换符合规定的新绳索，同时对无押运人的篷布苫盖货物均使用绳卡，保证运输安全。

在日常工作中加大检查和考核力度，狠抓作业标准化，完善了《车站货管细则》，从现场作业到台帐报表均制定了“标准化作业考核办法”。在日常检查中，对于违章违纪现象，发现一起考核一起。不管是整车、集装箱只要发生超载，每一吨罚款200元，严重超载责任者将给与待岗甚至下岗的处理，从而使得货运室以前懒散的工作作风有了彻底的好转。

二、积极服从大局，保增收创收。

今年以来，根据上级指示，积极开展货车增载工作，结合货运室现状，制定了增载考核的细化办法，确保一车不漏。

今年铁路局加大了对停时的考核力度。为此，货运室会同车站制定了“车站压缩停时考核办法”，并且按照上级要求，与专用线签订停时压缩协议，加大考核力度，从车辆到达起，按照入线前、出线后及货物线作业时间分段计算停留时间，按部门考核，谁造成的停时过大，由谁负责。通过考核，停留时间大幅度压缩，保证了货车的利用率。

今年五月份以来,集装箱需求非常紧张，缺口达300余箱。积压的100余个集装箱不到20天时间全部完毕。货运室积极与多方联系，调运空集装箱，以缓解集装箱运输能力的不足，同时为了缓解篷布的不足，积极与上级联系，先后调运了900床篷布。

三、搞好优质服务，保路风稳定。

货运工作是车站对外服务的窗口之一，我们只有用优质的服务、良好的业务素质，才能使这个窗口让货主感到信任、亲切，才能真正的实现“人民铁路为人民”的宗旨。为了让大家明白优质服务的重要性，彻底摆脱多年形成的“铁老大”意识，我们经常利用点名会等多种场合，结合电视和报纸的新闻向大家宣传服务的重要性和路风问题危害的严重性。使大家明白我们的收入是与各货主单位息息相关、密不可分的，货主就是上帝，是我们大家的摇钱树，也是我们赖以生存的命根子。只有尽可能多的留住货主单位，尽可能多的发运货物，才能确保我们主、副业收入的相对稳定。不管是谁，得罪了货主，受影响的不光是他本人，大家也都会受到影响。特别是在去年的竞争上岗过程中，货运室由20人减至14人，减员率达到30％，针对特殊情况，车站及货运室及时加强宣传教育，耐心做好货运职工的思想工作，防止了由于货运员岗位调整导致的情绪波动，而极有可能发生的服务质量问题。同时多次组织召开货主座谈会，向货主单位征求意见和建议，积极改进、完善我们工作中的不足。

任何一项工作要确保其良性发展，制度是保证，落实是关键。货运室在重新修订《货管细则》，同时在“标准化考核作业办法”中进一步细化。

货运室由于货场内办理煤炭和石膏装卸车作业，灰尘很大，工作环境也比较脏，为了从根本上解决环境脏、乱、差的现象，我们对于大的改造方面积极向上级汇报以求得上级的支持和理解，给予资金上的支持。小的方面，我们不等不靠，积极自筹资金逐步予以解决。在这方面我们货运室树立了漂亮、醒目的导示牌，重新标划了货物线以及货位，绿化了货运室的小花坛，采用租花的方式绿化了营业厅和办公室，调节了营业厅和办公室的枯燥气氛。特别是在上级的大力支持下，将货运大厅整体改造，各种揭挂统一上墙，改造了货运安全室，同时增加了货运微机自助查询系统，从根本上改变美化了环境方便了货主，

给货