分类变量的重复测量

分类变量的重复测量第1页，课件共51页，创作于2023年2月

分类变量（categoricalvariable）又称为定性变量（qualitativevariable），在工作中应用甚广。根据其不同的取值性质，又可分为3种类型：第一种是名义刻度（nominalscale）的分类变量，它是按事物属性分类的变量，如性别、职业等。在统计学上为了计算方便，将这些不同的属性进行数量化处理，如男性赋值为1，女性赋值为2。这种数值只是作为属性的代码，其间并无大小之分。第二种为有序刻度(ordinalscale)的分类变量，它是根据事物呈现出的程度或水平不同进行赋值。如临床化验结果用符号“-、+、++、+++”，文化程度用“文盲、小学、中学、大学、研究生”来划分等级，在进行数量化处理时赋值1、2、3、…。这里需要注意的是，1与2之差不一定等于2与3之差。第三种是区间刻度（intervalscale)，如人口学统计中的年龄分组，“0-，10-，20-，…”就是典型的例子。根据资料的性质，区间跨度有等距的，也有不等距的。第2页，课件共51页，创作于2023年2月把分类变量作为反应变量进行重复观察的情形在工作中应用较广。在本书第九章第五节中介绍了二分类反应变量的重复测量资料分析方法。本章将介绍分类反应变量重复测量资料的一般分析方法。主要介绍加权最小二乘法分析方法。第一节一个总体的二分类反应重复测量资料的分析。第3页，课件共51页，创作于2023年2月第一节二分类反应重复测量

资料的分析一、资料结构设有n例受试者，反应变量y，分为r水平（当r=2时称y为二分类反应变量)，有p个观测时间点，每一受试者可能出现的结果共有rp种。当r=2、p=3时的所有可能结果共有8种。当对y用（0，1）赋值方法时，n例受试者在这8种组合情况下出现的边际合计数为n000至n111。详细分配情况见下面示意图图11.1。第4页，课件共51页，创作于2023年2月边际合计数n000至n111是每种可能取值结果组合情况下的观察例数。总例数n为边际合计数之和。如果将上述示意图资料按独立观察资料整理各时间点的频数分布时（表11.1),每一时间点的合计数都为n=n1+=n2+=n3+，总例数变成了3n而不是n。但实际是观察例数为n，观察次数为3n。一般统计学方法是以表11.1的资料类型为出发点进行分析的。而本章介绍的方法是以图11.1的边际合计数为出发点进行分析的。第5页，课件共51页，创作于2023年2月二、模型结构受试者i在时间点j的反应变量为yij

，其取值概率为pj，（j=1,…,p）。单总体重复测量资料的任务为是要分析这p个概率p1，…，pp，是否相等。令反应频数的某种函数即反应函数（responsefunction）F＝f（m），建立的线性模型为：式中X为d×q维的设计矩阵，d=(r-1)p，为反应函数。B为q维未知参数。第6页，课件共51页，创作于2023年2月由于每一个体在不同时间点的反应变量的取值是相关的，因此，在配合模型时应将这种相关性纳入模型中。对式（11.1)求参数估计值的加权最小二乘解法就是在考虑到这种相关性的基础上建立起来的。模型（11.1)的加权最小二乘解法求参数ß的估计值的计算公式为第7页，课件共51页，创作于2023年2月例11.1一项药物对风湿性关节炎的治疗效果观察。用一种药物对46名患风湿性关节炎的病人进行治疗，共三个疗程。每个疗程结束后评价其治疗效果。有效者记为yij=1，否则yij=0i=1，…,46;j=1,2,3。资料列于表11.2中。表中的T1,T2,T3分别代表三个疗程。本项研究的目的是探讨不同疗程的疗差别。第8页，课件共51页，创作于2023年2月三、配合线性模型的步骤

表11.2为资料的原始记录形式，需要将其整理成边际频数表的格式后再配合模型。计算步骤介绍如下。

1.首先用下列SAS程序计算边际合计数程序中的subj为受试者号，time1、time2、time3代表3个疗程。第9页，课件共51页，创作于2023年2月第10页，课件共51页，创作于2023年2月计算出的边际频数列于表11.3中。从而得到边际频数向量n及其函数边际频率向量m为：

n=(616242466)’M=(6/460.34780.04350.08700.04350.08700.13040.1304)’

用nj及mj表示向量n及m中的元素。

2.给出转换矩阵A本例的转换矩阵A的结构为

A的作用是将反应频率M转换成反应函数F。第11页，课件共51页，创作于2023年2月

3.求出反应函数F它是每一时间点反应变量各个水平的边际概率。但它不包含该反应变量最后一个水平的反应概率。本例有3个时间点（time1time2time3），每个时间点的反应变量为r=2水平，故每一时间点只有一个反应函数。共有3个反应函数为：如f1=0.6086957=0.1304348+0.3478261+0.0434783+0.0869565第12页，课件共51页，创作于2023年2月第13页，课件共51页，创作于2023年2月第14页，课件共51页，创作于2023年2月第15页，课件共51页，创作于2023年2月第16页，课件共51页，创作于2023年2月以上分析可用SAS的CATMOD过程。程序如下：

PROCCATMODDATA=exmp11_1order=DATA;RESPONSEMARGINALS;MODELtime1*time2*time3=_RESPONSE_/ONEWAYCOV;REPEATEDtime3/_RESPONSE_=time;RUN;

程序中的RESPONSEMARGINALS语句是指定反应函数为在MODEL中规定的反应变量的边际概率。输出结果见表11.5、表11.6、表11.7。第17页，课件共51页，创作于2023年2月

从输出的表11.5反应函数和协方差矩阵中看出，第一和第二疗程的有效率都为0.6087（60.87%），第三疗程的有效率为0.3478（34.78%）。从表11.6方差分析看，标记‘time’的有效率分析，第三疗程的效果显著低于前两个疗程。对各疗程疗效一致性的假设检验可用下列公式计算。令对比矩阵C为：

第18页，课件共51页，创作于2023年2月从表11.7加权最小二乘估计值分析可见，有关TIME的两个参数估计值都为正值，反映药物在疗程1、2期的效果优于第3期的效果。第19页，课件共51页，创作于2023年2月在程序中的REPEATED语句的下面加入CONTRAST语句：CONTRAST‘time1VStime3’_RESPONSE_21;

或CONTRAST‘time1VStime3’ALL_PARMS021;为了检验第一疗程与第三疗程疗效的相等性，所用的检验假设为：Ho：2＝－2－3或Ho:22＋3＝0；H1：22＋30这是因为在参数化时用了关系式4＝－2－3（4代表第三疗程疗效），故有2＝4等价于2＝－2－3第20页，课件共51页，创作于2023年2月

上述SAS程序是用每个病例的记录建立数据集配合模型用的，也可以用频数表资料（表11.3)配合模型。所用SAS程序为：得到输出结果见表11.8。从表11.8对比分析看出，第一期的疗效与第三期疗效间有显著差异（P＝0.0161）。第21页，课件共51页，创作于2023年2月DATAexmpl1_1;INPUTtime1time2time3count;CARDS;111611016101210040112010400160006;PROCCATMODDATA=exmpl1_1ORDER=DATA;RESPONSEMARGINALS;WEIGHTcount;MODELtime1*time2*time3=_RESPONSE_/ONEWAYCOV;REPEATEDtime3/_RESPONSE_=time;RUN;第22页，课件共51页，创作于2023年2月第二节多组资料的分析方法

当有多组资料时，统计分析所要回答的问题有三点：①资料是否存在组间差异。②反应是否随时间变化。③是否存在处理组与时间之间的交互效应。分析所用公式以及计算步骤与单总体的情形基本一致。下面用一个具有二分类反应变量的2组重复测量资料为例加以介绍。第23页，课件共51页，创作于2023年2月例11.2对65岁及以上老年人参与某项社会活动的6年追踪调查。样本量为男性662人，女性1311人，分别在当年、第3年和第6年调查他们参与社会活动的情况。反应变量y的赋值规则为：

y=1经常参加该项社会活动

0不经常参加该项社会活动资料列于表11.9中。本例的分组变量为性别，分组数a＝2。反应变量y的分类数：r＝2。重复测量时间点数p＝3。计算步骤如下。第24页，课件共51页，创作于2023年2月第25页，课件共51页，创作于2023年2月

1．列出边际频数矩阵N及相应的边际频率矩阵M每一组的边际频数占一列，即N＝(n1’，n2)。与之对应，每一组的边际频率也占频率矩阵中的一列，即M＝(ml’，m2’)。第26页，课件共51页，创作于2023年2月同时可以建立一个16×1维的边际频率向量L为：

L＝(0.690.070.020.040.030.020.020.120.590.050.020.040.020.030.030.22）’

2.建立转换矩阵A由于每一组都是3个反应变量，每一反应变量具有2个水平，所以所建立的转换矩阵A与例11．1的A完全相同。即第27页，课件共51页，创作于2023年2月3．计算反应函数F第28页，课件共51页，创作于2023年2月

4．分组计算边际频率mg（g＝1，2）的协方差矩阵V（ml）和V（m2）5．两组反应频率的协方差矩阵V（M）本例的V（M）为16×16维矩阵，对角元素为V（mg），非主对角元素为0。第29页，课件共51页，创作于2023年2月

6．求F的协方差矩阵V（F）本例的V（F）为6×6维矩阵：7．定设计矩阵X例11.2资料的设计矩阵（含交互效应）见表11.11．第30页，课件共51页，创作于2023年2月故本例含有交互作用项的设计矩阵为：第31页，课件共51页，创作于2023年2月

8．用式（11.2）计算参数估计值（表11.12）第32页，课件共51页，创作于2023年2月9．用对比矩阵C作两组比较的假设检验C的结构为：利用式（11.6）得到2＝31.38，自由度＝3，查附表3：2界值表，有20.05（3）＝7.81，反映两组间的差异有统计学意义。以上计算过程所用SAS程序如下：第33页，课件共51页，创作于2023年2月DATAexmpl1_2;inputgender$year0year3year6count@@;CARDS;f111904f11088f10125f10051f01133f01022f00130f000158m111391m11036m10112m100143m01115m01021m00118m000143;PROCCATMODORDER=DATA;WEIGHTcount;RESPONSEmarginals;MODELyear0*year3*year6=gender|_RESPONSE_/PRED=FREQCOV;REPEATEDyear;RUN;第34页，课件共51页，创作于2023年2月输出结果见表11.13－表11.17。第35页，课件共51页，创作于2023年2月第36页，课件共51页，创作于2023年2月第37页，课件共51页，创作于2023年2月第38页，课件共51页，创作于2023年2月第39页，课件共51页，创作于2023年2月

从输出表11.17方差分析结果看，性别（gender）和年份（year）都有显著的统计学意义（P＜0.0001）；但性别*年份（gender*year）间交互作用无显著的统计学意义（P＝0.6476）。因此，在模型中剔除交互作用项后从新配合下列模型：

MODELyear0*year3*year6＝gender＿RESPONSE＿/FRED＝FREQCOV；

得到输出结果见表11.18，表11.19。第40页，课件共51页，创作于2023年2月第41页，课件共51页，创作于2023年2月第42页，课件共51页，创作于2023年2月

输出结果表11.19表明，性别间的差异和年份间的差异都具有统计学意义。由于无性别关年份交互作用，所以不同性别受试者对参与某项社会活动的时间趋势是类似的。