簡單隨機:
PROC PLAN SEED=210002;
FACTORS n=240; OUTPUT OUT=aaa;
DATA bbb; SET aaa; number=_n_;
IF n<=120 THEN group='A'; ELSE group='B';
PROC PRINT NOOBS;
VAR number group;
RUN;
分層區(qū)組隨機:
PROC PLAN SEED=210002;
FACTORS center=3 block=10 length=10;
OUTPUT OUT=fffff;
DATA ggggg; SET fffff; number=_n_;
IF length<=5 THEN group='A'; ELSE group='B';
PROC PRINT NOOBS;
VAR center number group;
RUN;
數(shù)據(jù)庫比較
proc compare base=sasuser.aaa compare=sasuser.bbb;
title ‘report’;
run;
秩和檢驗
Data PGM20C;
DO A=1 TO 4;
DO B=1 TO 2;
INPUT F; OUTPUT;
END; END; CARS;
43 44
15 15
2 1
0 0
;
DATA TWO; SET PGM20C;
DO I=1 TO F;
OUTPUT TWO;
END; RUN;
PROC NPAR1WAY WILCOXON DATA=TWO;
CLASS B; EXACT; VAR A; RUN;
方差分析:
重復(fù)測量方差分析程序
與某一時間點的比較,組間比較,未分組
data a;
set b
proc glm;
class group;
model T1 T2 T3 T4=group/nouni;
repeated time 4 contrast(2)/printm summary;
lsmeans group/stderr;