第十二章圆形分布资料及统计分析.pptVIP

第十二章圆形分布资料的统计分析 * * 上海第二医科大学 生物统计教研室 第一节 角度资料概论 圆形分布(circular distribution)统计方法用于处理角度资料 。医学中有些观察数据常用角度表示,例如心向量图的方位角,脑血流图的上升角,主峰角等；与环境卫生有关的风向也常用罗盘的方向角度来表示。有些数据以一年中的月,日或一昼夜中的时,分来表示,前者如正常人血压值在一年中各月份的变动; 某病的发病率在一年中是否有好发时间,后者如婴儿的出生时刻,心脏病人的发病时刻, 在一天中的任何时刻均有可能,可以研究是否有集中于某一时刻的倾向,这一类时间性的资料可化成角度资料来处理。 圆形分布中的角度,指的都是圆心角, 其特点是周而复始,没有真正的零点,也没有大小之分。习惯上把正北方向定为0,一昼夜中的正午夜(0点0分)也定为0,一年中的1月1日午夜也定为0,但这完全是人为规定的。 圆形分布中最常见的是Von Mises分布, 这是一个单峰圆形分布,相当于线性资料的正态分布,本章所讨论是的都是这类分布。当角度资料在圆上的分布有集中于一个方向的趋势, 所求得的平均角(mean angle)经检验不是均匀分布, 且为一个集中方向时就称之为单峰圆形分布。反之,当角度资料在圆上的分布均匀(uniform circular distribution),无明显的集中趋势,就认为平均角不存在。 第二节 角的均数及其假设检验 二．角离差S和集中趋势r 样本统计量r在圆形分布中是描述α离散程度的一种统计指标,它与α的标准差s的关系如下: s=(180/π) 度 当一组数据中所有αi都等于同一数值时,则这组数据无变异, s =0,而r=1, 当一组数据中的αi均匀地分布在圆周上,则r=0,而s则因平均角不存在而无法计算,但当r趋向于0时,s趋向于无穷大。r值的范围在0～1之间, s值的范围在0～∞无穷大之间,s可称为圆标准差 三．平均角的假设检验 1. 所有αi都均匀分布在圆周的一个总体,其集中趋势量度值ρ=0,但在此总体中随机抽一个样本,所得ρ的估计值r不一定为0,因此, 当同一个样本资料算得平均角与r后,此是否意义(即是否来自ρ=0的总体)必须进行假设检验,称为均匀性检验(test of uniformity)。此时，H0:ρ=0,即为均匀分布,不存在平均角；H1:ρ≠0,即不是均匀分布,存在平均角。 2. 均匀性检验方法很简单,根据样本大小n和算得的r或rc查附表十五,如r(或rc)大于或等于表中界值,则P≤α,即在相应的α水准上拒绝H0,表示存在集中趋势,平均角有意义。如r小于表中界值,则Pα,即在α水准上不拒绝H0,认为是均匀分布, 不存在集中趋势,故均匀角无意义。 检验表12.2二十名妇女的分娩时间所得的均匀角有无意义? 解:H0:ρ=0 H1:ρ≠0 该例中已求得r=0.71645,n=20,查附表十五,n=20时,r0.05=0.3846,r0.01= 0.4718，r0.001=0.5687,现rr0.001故P0.001,即在α=0.001水准上拒绝H0,认为存在集中趋势,此平均角有意义。 第四节 两个或多个样本平均角的比较 两个或多个样本的平均角各自经均匀性检验,如果都拒绝 H0, 则可用 Watson-William检验,判断它们是否来自总体平均角都为ρ的总体,即比较平均角之间是否有显著差别。 二．两样本平均角比较的U2检验 用U2检验法对均匀性及合并r大小等无特殊要求, 故不必作平均角的均匀性检验及求合并r值。 患A病的6个病人,晚上入睡时间分别为20:30,21:00,21:15,21:20,21:45,22:00; 患B病的7个病人,入睡时间分别为21:30,21:45,22:05,22:15,22:20,22:45,22:50, 问两种病病人入睡时间的迟早有无差别? （1）将两样本中的时间化为角度,再按角度大小从上向下排列,但仍分两组, 各自编序号i与j,排列结果见表12.8第2,4列。 （2）两组分别计算i/n1及j/n2,见表12.8第3,6列,A,B两组各有12个可标i或j的格子,而A组有6格填了i,有6个空格,B组有7格填了j,有5个空格,有空格的各行,也要计算i/n1,和j/n2,其i,j的决定方法为: ①i或j等于1以前各空格其值为0,例如表中第6列1～4行,j/n2的计算方法为0/7即等于0。 ②如i或j等于C之后有空格,这些空格在计算i/n1或j/n2时