定量指标组间比较的SAS编程操作.ppt

* means语句中可设置众多的选项，均与指定效应的组间多重比较有关。 * 例 某地方病研究所调查了8名正常儿童的尿肌酐含量（mmol/24h）如下表，试分析尿肌酐含量和其年龄的相关关系。 1. 直线相关分析 表 8名正常儿童的年龄（岁）与尿肌酐含量（mmol/24h） 程序： data r19; input x y; cards; 13 3.54 11 3.01 9 3.09 6 2.48 8 2.56 10 3.36 12 3.18 7 2.65 ; proc plot; plot y*x=‘#; run; proc corr; var x y; run; proc plot; plot y*x=‘#; run; proc corr; var x y; run; r19 结果： 直线相关系数 P值 ? 结论： 本例Pearson相关系数为r=0.88177， 所对应的P=0.00380.05，说明两个变量之 间存在正相关关系，即一个变量的值增大 时，另一个变量的值也相应地增大。 例9-1 某地方病研究所调查了8名正常儿童的尿肌酐含量（mmol/24h）如表9-1，估计尿肌酐含量对其年龄的回归方程。 2. 直线回归分析 程序： data r19; input x y; cards; 13 3.54 11 3.01 9 3.09 6 2.48 8 2.56 10 3.36 12 3.18 7 2.65 ; proc reg; model y=x; run; 校正决定系数 误差均方的平方根 应变量变异系数 应变量均数 决定系数 方差分析结果 回归方程为： y? = 1.661667 + 0.139167 x 常用选项: ? stb: 输出标准化偏回归系数 , 语句为“model y=x/stb;” 常用选项: ? p: 输出每个观测的实际值、预测值和残差。语句为： “model y=x/p;” 应变量原始值 原始值的预测值 残 差 程序： data r19; input x y@@; cards; 13 3.54 11 3.01 9 3.09 6 2.48 8 2.56 10 3.36 12 3.18 7 2.65 ; proc reg; model y=x/stb p; run; 练习： 给10只中年大鼠注射内毒素后，测得每只大鼠的红细胞含量x与血红蛋白含量y分别如下： X：654 786 667 605 761 642 652 706 602 539 Y: 130 168 143 130 158 129 151 153 149 109 试对x、y两变量进行线性相关和回归分析。 程序: data r20; input x y@@; cards; 654 130 786 168 667 143 605 130 761 158 642 129 652 151 706 153 602 149 539 109 ; proc plot; plot y*x=‘#; run; Proc corr， Var x y； Run； proc reg; model y=x; run; 绘制散点图 相关分析 回归分析 结果： 截距项与0直接的差别无统计学意义。 程序: data r20; input x y@@; cards; 654 130 786 168 667 143 605 130 761 158 642 129 652 151 706 153 602 149 539 109 ; proc plot; plot y*x=‘#; run; Proc corr， Var x y； Run； proc reg; model y=x; run; proc reg; model y=x /noint r cli clm; run; 绘制散点图 相关分析 回归分析 去掉截距项， 进行残差分析 计算各x值所对应的Y值的容许区间。 计算各x值所对应的Y的总体均值的置信区间。 结果: 秩相关 例9.8 某省调查了1995年到1999年当地居民18类死因的构 成以及每种死因导致的潜在工作损失年数WYPLL的构成， 结果见表77.33。以死因构成为x，WYPLL构成为y，试作秩相 关分析。 表7.3 某省1995年到1999年居民死因构成与WYPLL构成 死因类别 1 2 3 4 5 6 7 8 9 各死因构成(%) 0.03 0.14 0.20 0.43 0.44 0.45 0.47 0.65 0.95 WYPLL构成 0.05 0.34 0.93 0.69 0.38 0.79 1.19 4.74 2.31 死因类别