课件协方差分析2015.pptxVIP

1;学习内容：协方差分析;控制混杂的方法;问题的提出;5;6;7;质疑结果…;协方差分析的概念;观测值=一般均值+水平影响+协变量影响+随机误差;协方差分析数据结构（单因素随机区组设计试验资料）;基本思想;;回到研究实例;协方差分析的应用条件 ;SAS语句：检验应用条件;一般线性模型General Linear Model (GLM);Proc GLM与其它程序的比较;Proc GLM与其它程序的比较（续）;Proc GLM 选项alpha= (默认0.05); Class 分类变量表； 必选Model句前！ Model 因变量 =自变量表 /选项; 仅允许1Model句！ alpha= (默认0.05) cli 显示预测的个体Y值的 容许区间 clm 显示预测总体均值置信区间 No同时用cli 与clm！ Int或noint 是否包含截距项 P: 对自变量不缺失的观测显示观测值、预测值、残差值 Solution： 输出回归方程中各参数估计值 SS1, SS2, SS3, SS4: 输出每个效应的各型估计方程所对应的平方和（默认SS1）。 ;Model后面的2个可选语句;估计效应参数（回归系数）涉及的ＳＳ;I型SS：与模型中自变量的进入顺序有关，即模型中每一个自变量对因变量的影响效应仅仅校正了在它之前进入模型的自变量（不包括交叉变量） II型SS：与模型中自变量的进入顺序无关，即模型中每一个自变量对因变量的影响效应都校正了模型中其它自变量（不包括交叉变量） ;III型SS和IV型SS：与模型中自变量的进入顺序无关。即模型中每一个自变量（包括交叉变量）对因变量的影响效应都校正了模型中其它自变量（不包括交叉变量）。 这两种类型SS的区别只有当数据不平衡，且有空组的设计时才显示出来。;4种SS算法的选择;SAS分析实例1:两种药物治疗癫疯病菌效果观察;;检验方差齐性： Proc discrim data=aa pool=test; Class drug; Var y; Run; ;检验回归直线斜率是否相等： proc glm data=aa; class drug; model y=drug x drug * x; run; ;SAS结果解读;研究核黄素缺乏对蛋白质利用的影响。将36只大白鼠按某些重要的非处理因素配成12个区组，用随机的方法决定每组中的3只分别进入3个饲料组。 饲料1为缺乏核黄素的饲料；饲料2为含核黄素的饲料，但限制食量使其与第一组的食量相近；饲料3为含核黄素的饲料、但不限制食量。 3组大白鼠的进食量为x，同期内所增体重为y。 分析3种饲料对大白鼠体重的影响有无差别。 ;区组号;cards; 256.9 27.0 260.3 32.0 544.7 160.3 271.6 41.7 271.1 47.1 481.2 96.1 210.2 25.0 214.7 36.7 418.9 114.6 300.1 52.0 300.1 65.0 556.6 134.8 262.2 14.5 269.7 39.0 394.5 76.3 304.4 48.8 307.5 37.9 426.6 72.8 272.4 48.0 278.9 51.5 416.1 99.4 248.2 9.5 256.2 26.7 549.9 133.7 242.8 37.0 240.8 41.0 580.5 147.0 342.9 56.5 340.7 61.3 608.3 165.8 356.9 76.0 356.3 102.1 559.6 169.8 198.2 9.2 199.2 8.1 371.9 54.3 ; Run;;proc sort; by a; run; proc univariate normal; var y; by a; run; proc sort; by b; run; proc univariate normal; var y; by b; run;;proc glm; class a b; model y=x a b x*a; run; proc glm; class a b; model y=x a b; lsmeans b/stderr pdiff; run; ;练习：镉作业工人接触烟尘年数与肺活量关系的研究;练习; 练习 filename aah:\lecture-02\sas\bp962.dbf; proc dbf db3=aa out=bp; Data bb; Set bp; Proc sort data=bb; By group; Run;