SAS过程步课件

1、SAS过程步  2011-03-20 07:11:28|  分类： 默认分类 |  标签：sas  |举报|字号 订阅编号过程名功用必要的语句means计算基本统计量proc means;var x;run;chart制作次数分布表次数分布图proc chart;hbar x;run;univariate plot正态分布检验proc univariate plot;var x;run;PRINT数 据 输 出PROC PRINT ;RUN;SORT数 据 排 序PROC SORT;B

2、Y age;PROC PRINT;RUN;tabulate制作表格PROC PRINT;Proc tabulate format=8.0;Class repname month;var unitsold;table month,repname*unitsold*sum;RUN;anova;单方面分类的方差分析proc anova;classes a;model x = a;anova双方面分类的方差分析proc anova;classes block a;model x =block a;anova三方面分类的方差分析proc anova;classes A B T;model X1=A B

3、T ;  可选择:LSQduncanduncan多重比较可选择:LSQ测验邓肯Q测验图基Q测验proc anova;classes block a;model x =block a;means block a /LSQ alpha=0.05;means block a /LSQ alpha=0.01;  anova2*2 复因子试验的统计分析方法proc anova;classes block pinxi midu;model x=block pinxi midu pinxi*midu;  anova;2*2*2 复因子试验的统计分析方法proc ano

4、va;classes block k p n x;model x=block k p n k*p k*n p*n k*p*n ;  glm协方差分析proc glm;class t ;model y=t x;  anova多元方差分析proc anova;class block tre ;model x1 x2 x3 =block tre;manova h=block tre /summary;run;  anova正交设计的分差分析proc anova;classes block a b c d ;model x = block a b c d ;means b

5、lock a b c d /lsd duncan tukey;run;  nested嵌套设计的方差分析proc nested;classes plant leaf;var calcium;run;  glm带有交叉项的双向不均衡设计的方差分析proc glm;class drug disease;model y=drug disease drug*disease /ss1 ss2 ss3 ss4;run;  corr简 单 相 关 系 数proc corr;var y;with x1;run;  corr一个变量与多个变量的简单相关系数proc cor

6、r;var y;with a b c d e f ;run;  corr多个变量间的简单相关系数的计算proc corr;run;  plot相关点式图的绘制plot y*x='+'rproc plot;un;  corr,partial一 级 净 相 关proc corr;var y;with a; partial f;run;  corr,partial二 级 净 相 关proc corr;var y;with f; partial a b ;run;  corr,partial三 级 净 相 关proc corr;var

7、y;with f; partial a b d ;run;  cancorr典 型 相 关 分 析proc cancorr;var y1 y2 y3 ;with x1 x2 x3 ;run;  reg直 线 回 归proc reg;model y=x1;run;  reg二 元 回 归proc reg;model y=x1 x2;run;  reg可选择:forwardbackwardstepwise多 元 回 归可选择:向 前 回 归向 后 回 归双 重 逐 步 回 归proc reg;model g=a b c d e f  selecti

8、on=forward;run;  cluster可选择:averagewardsdensitysingletwostage聚 类 分 析可选择:类 平 均 法重 心 法离 差 平 均 和 法密 度 估 计 法最 短 距 离 法两阶段密度估计法proc cluster method=average outtree=cate;var a b c d e f ; id bh;proc tree;run;  fastclus动态聚类分析proc fastclus data=iris maxc=2 maxiter=10 out=clus;var sepallen sepalwid p

9、etallen petalwid;proc freq;tables cluster*species;run;         standard数据标准化proc standard mean=0 std=1 out=new;run;  princomp主 成 分 分 析可选择:由 协 差 阵 计 算由 相 关 矩 阵 计 算proc princomp cov;run;或proc princomp;run;  factor 可选择:因 子 分 析可选择:主 分 量 分 析主 因 子 分 析proc

10、factor data=socecon simple corr;run;或proc factor data=socecon priors=smc msa scree residual preplotrotate=promax reorder plotoutstat=fact all;proc print;run;  ttest成 组 法 T 检 验proc ttest data=pinzhs;class pin;var yield;run;  univariate配 对 法 T 检 验proc univariate data=chromat;var methdiff;run

11、;  npar1way秩 和 检 验proc npar1way data=gastric wilcoxon;class group;var lysolevl;run;  univariate符 号 秩 检 验proc print data=chromat;proc univariate data=chromat;var methdiff;run;  freq卡平方测验可选择:2×2 联列表式2×j 联列表式r×c联列表式data;do a=1 to 2;do b=1 to 2;input x;output;end;end;cards;

12、.;proc freq;weight x;tables a*b / chisq ;run;  reg二 次 抛 物 线proc reg data=uspop;var yearsq;model pop=year/r cli clm;plot r. *p. ;add yearsq;print;plot;plot pop*year='a' predicted.*year='p' u95. *year='u'l95. *year='l'/overlay;run;  NLIN可选择:dudMARQUARDT负 指 数 曲

13、 线可选择:试 位 法麦 夸 特 法PROC NLIN BEST=10 METHOD=dud;PARMS B0=0 TO 2 BY .5 B1=.01 TO .09 BY .01;MODEL Y=B0*(1-EXP(-B1*X);DER.B0=1-EXP(-B1*X);DER.B1=B0*X*EXP(-B1*X);OUTPUT OUT=B P=YHAT R=YRESID;PROC PLOT DATA=B;PLOT Y*X='A' YHAT*X='P'/OVERLAY VPOS=25;PLOT YRESID*X/VREF=0 VPOS=25;RUN; 

14、catmod逻 辑 斯 蒂 回 归（标 准 响 应 函 数）proc catmod;weight count;direct heat soak;model y=heat soak / freq ml nogls covb corrb;quit;  sortRSREG二 次 响 应 面 分 析proc sort;by time temp;proc rsreg;model mbt=time temp;ridge max;run;  lackfit欠适测验model mbt=time temp/lackfit;ridge max;run;  score得 分 分 析pr

15、oc score data=fitness score=factout out=fscore;var age weight runtime runpulse rstpulse;run;  stepdisc逐 步 判 别 分 析proc stepdisc data=iris bsscp tsscp;class species ;var sepallen sepalwid petallen petalwid;run;  candisc典 型 判 别 分 析proc candisc data=iris out=outcan distance anova;class species

16、 ;var sepallen sepalwid petallen petalwid;run;proc plot;plot can2*can1=species;format species specchar.;title'Plot of Canonical Variables'run;  RSREG岭 嵴 分 析PROC SORT;BY X1 X2 X3;PROC RSREG;MODEL Y=X1 X2 X3/LACKFIT;RUN;DATA B;*-GET THE ACTUAL VALUES-;SET A END=EOF;OUTPUT;*-CREATE AN X1*X2 GRID FOR PLOTTING-;IF EOF THEN DO;Y=