cdc课件-sas讲义_协方差分析

cdc课件-sas讲义_协方差分析

ID:1682137

大小:161.50 KB

页数:9页

时间:2017-11-13

cdc课件-sas讲义_协方差分析_第1页
cdc课件-sas讲义_协方差分析_第2页
cdc课件-sas讲义_协方差分析_第3页
cdc课件-sas讲义_协方差分析_第4页
cdc课件-sas讲义_协方差分析_第5页
资源描述:

《cdc课件-sas讲义_协方差分析》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库

1、协方差分析用SAS中的glm过程进行协方差分析时,要注意不同试验设计时class语句和model语句的写法。设分类变量为A、B,协变量为X,观察值为Y,则有:①单因素k水平设计的协方差分析模型classA;modelXA;②随机区组设计的协方差分析模型classAB;modelXAB;③两因素析因设计的协方差分析模型classAB;modelXABA*B;实例分析1.一元协方差分析例1研究牡蛎在不同温度的水中不同位置上的生长情况。有人做了如下试验:分别在通向发电站的入口处(温度较低)不同位置(底部和表层)和出口处(温度较高)不同位置(底部和表层)及电站附近的深水处(底

2、部和表层的中间)总共5个不同位置点上,随机地各放4袋牡蛎(每袋中有10个),共5×4=20袋。在将每袋牡蛎放入位置点之前,先洗干净称出每袋的初始体重,放在5个不同点一个月后再称出最后体重。试验结果数据见表26.1所示。表1.1牡蛎在不同温度和位置上的生长数据位置trt重复数rep(x为初始体重,y为最后体重)1234xyxyxyxy1(入口底部)27.232.632.036.633.037.726.831.02(入口顶部)28.633.826.831.726.530.726.830.43(出口底部)28.635.222.429.123.228.924.430.24(出

3、口顶部)29.335.021.827.030.336.424.330.55(附近中部)20.424.619.623.425.130.318.121.8程序如下:datagrowth;dotrt=1to5;dorep=1to4;9inputxy@@;output;end;end;cards;27.232.632.036.633.037.726.831.028.633.826.831.726.530.726.830.428.635.222.429.123.228.924.430.229.335.021.827.030.336.424.330.520.424.619.623

4、.425.130.318.121.8;procanovadata=growth;classtrt;modely=trt;procglmdata=growth;classtrt;modely=trtx/solution;meanstrt;lsmeanstrt/stderrtdiff;contrast'trt12vstrt34'trt-1-1110;estimate'trt1adjmean'intercept1trt10000x25.76;estimate'trt2adjmean'intercept1trt01000x25.76;estimate'adjtrtdiff't

5、rt1-1000;estimate'trt1unadjmean'intercept1trt10000x29.75;estimate'trt2unadjmean'intercept1trt01000x27.175;estimate'unadjtrtdiff'trt1-1000x2.575;run;程序说明:定性变量trt的5个不同位置点对y可能有较大的影响,因此class语句中分组变量为trt,先选用anova过程进行方差分析。然而,牡蛎的初始体重x对牡蛎的最后体重y可能也有一定的影响,故适合选用glm过程进行协方差分析,在model语句中不仅包括分组变量trt,而且应

6、包括协变量x。选择项solution要求输出回归系数的估计值及其标准误差和假设检验等结果。means和lsmeans语句要求输出分组变量trt各水平下y的未修正均值和修正后的均值,选择项stderr要求输出y的修正均值的标准误差、各修正均值与0比较的假设检验结果;选择项tdiff要求输出y的各修正均值之间两两比较所对应的t值和p值。Contrast语句是用来比较入口处底部和顶部均值之和与出口处底部和顶部均值之和是否相等。前三条estimate语句是用来估计入口处底部和顶部调整后的均值及它们之差,并假设检验是否为0,后三条estimate语句是用来估计入口处底部和顶部未

7、调整的均值及它们之差,并假设检验是否为0。程序输出的主要结果见表1.21.31.4所示。9TheSASSystemAnalysisofVarianceProcedureDependentVariable:YSourceDFSumofSquaresMeanSquareFValuePr>FModel4198.4070000049.601750004.640.0122Error15160.2625000010.68416667CorrectedTotal19358.66950000R-SquareC.V.RootMSEYMean0.55317510.5

当前文档最多预览五页,下载文档查看全文

此文档下载收益归作者所有

当前文档最多预览五页,下载文档查看全文
温馨提示:
1. 部分包含数学公式或PPT动画的文件,查看预览时可能会显示错乱或异常,文件下载后无此问题,请放心下载。
2. 本文档由用户上传,版权归属用户,天天文库负责整理代发布。如果您对本文档版权有争议请及时联系客服。
3. 下载前请仔细阅读文档内容,确认文档内容符合您的需求后进行下载,若出现内容与标题不符可向本站投诉处理。
4. 下载文档时可能由于网络波动等原因无法下载或下载错误,付费完成后未能成功下载的用户请联系客服处理。