统计SAS实例及答案-何飞武.docxVIP

统计SAS实例及答案（专业结论待完善）目录第一部分 计量资料3（一）t检验31、单样本——实例3-532、两组配对——实例3-643、两组独立——实例3-75（二）方差分析61、完全随机——实例4-262、随机区组——实例4-473、两阶段交叉——实例4-684、析因设计92×2——补充1个92×2——实例11-1103×3——实例11-2115×2×2——实例11-3125、正交设计14有空白列，无重复——实例11-414无空白列，有重复——补充打印116有空白列，有重复——补充打印217（三）非参181、单样本——类似3-5——实例8-2182、两组配对——类似3-6——实例8-1193、两组独立——类似3-7——实例8-3204、多组——非参——实例8-521第二部分 计数资料22（一）四格表（2×2）221、独立四格表22基本卡方?2——实例7-122连续校正?c2或Fisher精确概率——实例7-223Fisher精确概率——实例7-4242、配对四格表25Kappa值——实例7-325（二）行×列表（r×c）261、教材三个263×2——率——实例7-6262×3——构成比——实例7-7274×3——列联系数c——实例7-8282、补充三个292×3——列有序-CMH?2——补充1个293×3——双向有序-CMH?2——补充1个304×3——分层卡方-CMH?2——补充1个31（三）非参322×5——列有序-非参——实例8-4324×4——列有序-非参——实例8-733第三部分 回归34（一）相关34计量均正态-Pearson r（r）——实例9-534计数秩相关-Spearman r（rs）——实例9-835（二）回归361、直线回归36一元线性回归——实例9-136多元线性回归——实例15-1372、曲线回归39对数y=a+blog10(x)——实例9-1339指数y=exp(a+bx)——实例9-1440抛物线y=c+bx+axx——补充41（三）回归诊断-略41（四）逻辑回归421、y二值42一元无哑变量——补充42一元有哑变量——补充43多元无交互作用——补充44多元有交互作用——补充45多元有交互作用-筛选危险因子——实例16-2462、y多值-略47第四部分 其他481、聚类-略482、判别-略483、主成分分析-略484、因子分析——补充48实例第一部分 计量资料（一）t检验1、单样本——实例3-5原始数据：112 105 172 157 165 140 137 178 116 151 171 163129 130 125 135 128 100 126 128 90 113 128 12988 126 96 175 160 116 90 103 162 129 110 127SAS：data eg; input x@@;cards;112……127;proc univariate normal;/*正态性检验*/ var x;proc ttest h0=140;/*t检验，均值和140比较*/ var x;run;答案：H0：服从正态分布P=0.0988＞0.05不拒绝H0，服从正态分布H0：等于140t=-2.14，P=0.0397P＜0.05，拒绝H0，差异有统计学意义x=130.83＜140因此，所有从事铅工作男工人血红蛋白含量均数低于正常人均数140g/L2、两组配对——实例3-6SAS：data eg; input x1 x2@@; d=x1-x2;cards;0.840.58……0.870.506;proc univariate normal; var d;/*对d值做正态性检验*/proc ttest; paired x1*x2;/*配对t检验*/proc means n mean; var x1 x2;/*算两组的均值*/run;答案：H0：d值服从正态分布P=0.9993＞0.10不拒绝H0，d值服从正态分布H0：两种相同t=7.93，P＜0.0001拒绝H0，差异有统计学意义x1=0.7952＞x2=0.5228因此，哥特里-罗紫法偏高3、两组独立——实例3-7SAS：data eg; input c x@@;cards;1-0.7……2-2;proc univariate normal; var x; class c;/*按c进行分组*/proc ttest; var x; class c;run;答案：两组正态性检验：P=0.7527P=0.7754均＞0.10，均正态H0：方差齐Equality of VariancesP=0.