原创力文档- 1、本文档共10页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
查看更多
SAS统计课
SAS统计课
SAS
一个语句是一行。
一个语句是以关键词开始，分号结束。如果语句不对，则颜色不对。
语句是分组，组叫块，一块叫步，
数据是以data，过程是以proc开始的。提交运行。
每一步是分步执行的，日志在log上面。
主要的方式是程序方式，但也有非程序的方式：利用分析源建立数据并进行分析。
与一般的统计软件相似。
数据比较少时，可以直接放在编辑的程序里。数据比较多时，则不太好。
1、Data 建立数据集，可以放在work库里，可以放到永久库里，在最上面加上libname my test ‘C:\’ 库名是八个字符，是以字母开头。Data mytest.bodyfat;后面是input规定输入的变量以及变量按怎样的格式读数据。Input 变量：以字母开始，最多32个（V8.2）字符。
Dataline or cards下面是数据，结束时加上一个分号，和run;，永久库的好处是：系统退出后，文件依然存在，下一次可以直接调用。
2、利用已有的数据集进行分析 。
如：
libname myt ‘c:\practice’;
data bodyf;
set myt.bodyfat;
run;
运行就可以了。
利用文本文件。Data bodyf;
infile ‘c:\practice\body.txt’;
input age bodyfat;
run;
这个可以用于datalog的数据管理。
第一行输入变量名，下面的是数据。当然可以不用变量名，但不太好。输入数据时会少一个数据，但要注意一个选择框。
单向分组资料的方差分析：
固定模型效应：
5个处理，每个处理有4个观察值。
＠行跟踪符，多余的数据供后面的input语句做。
过程语句：分类变量class，模型语句model, （误差），多重比较，来筛选最佳数据。
data fertilization;
input trt$ @ ;
do r=1 to 4;
input y@;
output;
end;
cards;
NH4 24 30 28 26
NH42 27 24 21 29
NH4Co3 31 28 25 30
CONH2 32 33 33 28
No 21 22 16 21
;
proc anova;
class trt;
model y=trt;
means trt/ t snk duncan ;
means trt/ t /snk dunca alpha=0.01;
run;
两向分组资料的方差分析：组合内只有单个观察值的两向分组资料
data pea;
length trt $ 10;
input trt $ @;
do block=’ⅰ’, ‘ⅱ’, ‘ⅲ’,’ⅳ’;
input y@;
output;
end;
cards;
ck 60 62 61 60
Chi 65 65 67 65
Dong 63 61 61 60
Dou 64 67 63 61
Lou 62 65 62 64
Ma 61 62 62 65
;
proc anova;
class block trt;
model y=block trt;
means trt /duncan dunnett(‘ck’);
run;
有缺少使用glm而不用anova
Coeff Var10%正常（15%）
原因：模型没有用对（其他因素在里面），实验数据差。
、
二级各级系统分组资料的方差分析
先一个植株，从植株上取叶片，这是一个嵌套。
Proc glm;
Class cultur pot;
Model y=culture pot(culture);
Test h=culture e=pot(culture);
Means culture / duncan e=pot*(culture);
Means culture / duncan e=pot(culture) alpha=0.01;
Random pot(culture) /test;
Run;
Bartlett测验
data fertilization;
input trt$ @ ;
do r=1 to 4;
input y@;
output;
end;
cards;
NH4 24 30 28 26
NH42 27 24 21 29
NH4Co3 31 28 25 30
CONH2 32 33 33 28
No 21 22 16 21
;
proc glm;
class trt;
model y=trt;
means trt/ hovtest=bartlett ;
output out=temp resi
文档评论（0）