医学STATA软件及其应用.ppt

主要内容 定量资料的假设检验 t检验单样本t检验配对t检验成组t检验及方差齐性检验方差分析单因素方差分析及两两比较两因素以上的方差分析 样本均数与总体均数比较的t检验 有原始数据时ttest变量名 val仅有数据之统计量时ttesti obs mean sd val 例1 Ex 1 样本 某医生随机抽查10名某病患者的血红蛋白 如下 问题 该病患者的平均Hb含量是否与正常人的平均Hb含量相同 正常人的平均Hb含量为14 02 g dl 例 样本 某医生随机抽查10名某病患者的血红蛋白 求得其均数为12 59 g dl 标准差为1 632619 g dl 问题 该病患者的平均Hb含量是否与正常人的平均Hb含量相同 正常人的平均Hb含量为14 02 g dl 操作 ttest变量名 valttestx 14 02ttesti obs mean sd valttesti1012 591 63261914 02 One samplettest ObsMeanStd Err Std Dev 95 Conf Interval x 1012 59 51627951 63261911 4220913 75791 Degreesoffreedom 9Ho mean x 14 02Ha mean14 02t 2 7698t 2 7698t 2 7698P t 0 0218P t 0 9891 P t 2 7698 0 0218 假设检验与可信区间的分析具有相同的量化判定 其结论是等价的 ttesti1012 51 63261914 02 配对设计 检验 ttest变量1 变量2注意 变量1和变量2必须是成对输入的配对样本昆明治疗羊角风 例3 配对设计 Ex 3 例3现用两种测量肺活量的仪器对12名妇女测得最大呼气率 PEER L min 资料如下 问两种方法的检测结果有无差别 P60 ttestx2 x1 Pairedttest Variable ObsMeanStd Err Std Dev 95 Conf Interval x2 12419 333326 0085590 09624362 0889476 5778x1 12402 166727 4159494 97161341 8246462 5087 diff 1217 1666711 6416440 32782 8 4564142 78974 Ho mean x2 x1 mean diff 0Ha mean diff 0t 1 4746t 1 4746t 1 4746P t 0 1684P t 0 0842 ttestx1 x2 是否还有其他方法 独立性正态性方差齐性方差齐 t检验方差不齐 t 检验 t检验应用条件 两组资料间的方差齐性检验 sdtest变量名1 变量名2 两组数据各用一个变量sdtest变量名 by 分组变量 两组数据用同一个变量输入 用分组变量加以区分sdtesti obs1 mean1 sd1 obs2 mean2 sd2 mean1 mean2可被替换为点号 成组设计 检验 ttest变量1 变量2 unpaired unequal welch ttest变量 by 分组变量 unequal welch ttesti obs1 mean1 sd1 obs2 mean2 sd2 unequalwelch unpaired表示非配对的 如不选就作配对 检验unequal表示假设两组方差不齐 如不选表示假设两组方差达到齐性welch当方差不齐时 用welch法比较注 welch必须和unequal同时选才有效 例 成组设计 Ex 5 某医生研究转铁蛋白对病毒性肝炎诊断的临床意义 测得12名正常人和15名病毒性肝炎患者血清转铁蛋白含量 g dl 结果如下 问患者和正常人转铁蛋白含量是否有差异 同例 Varianceratiotest Group ObsMeanStd Err Std Dev 95 Conf Interval 1 12271 89172 99558910 37702265 2984278 48492 15235 20673 71566214 3907227 2374243 176 combined 27251 51114 31252922 40856242 6466260 3756 Ho sd 1 sd 2 F 11 14 observed F obs 0 520F 11 14 lowertail F L F obs 0 520F 11 14 uppertail F U 1 F obs 1 923Ha sd 1 sd 2 PF U 0 2649P F obs 0 8596 Sdtestx by group ttestx by group Two samplettestwithequalvariances Group ObsMeanStd Err Std Dev 95 Conf Interval 1 12271 89172 99558910 37702265 2984278 48492 15235 20673 71566214 3907227 2374243 176 combined 27251 51114 31252922 40856242 6466260 3756 diff 36 6854 9500326 4902246 87978 Degreesoffreedom 25Ho mean 1 mean 2 diff 0Ha diff0t 7 4111t 7 4111t 7 4111P t 0 0000P t 0 0000 P t 7 4111 0 0000 思考 t检验两个均数比较t 检验那么 三个 四个 均数如何比较 方差齐 方差不齐 方差分析 单因素方差分析及方差齐性检验 oneway响应变量分组变量 选择项 noanova 不打印方差分析表missing 将缺省值作为单独的一组tabulate 打印各组的基本统计量表scheffe Scheffe法bonferroni Bonferroni法sidak Sidak法注意 方差分析时数据必须按特定形式输入 即新定义一分组变量g 而将不同处理组的观察值以某变量x表示 例6 4三组小鼠的FDP酶活力 资料输入形式 XgABC2 7912 793 835 412 6912 693 153 47 3 8323 152 Stackabc into x onewayxg tsch Summaryofxg MeanStd Dev Freq 1 2 7025 5001356982 3 66125 9850806983 3 90125011 01644258 Total 3 4216667 9827320724AnalysisofVarianceSourceSSdfMSFProb F Betweengroups6 4368089723 218404494 280 0275Withingroups15 775724621 751224983 Total22 212533623 965762331Bartlett stestforequalvariances chi2 2 3 4559Prob chi2 0 178 两两比较结果 Comparisonofxbyg Scheffe RowMean ColMean 12 2 95875 0 111 3 1 198