SAS练习题及程序答案

1．随机取组有无重复试验的两种本题是无重复

DATA PGM15G;

DO A=1TO4; /*A为窝别*/

DO B=1TO3; /*B为雌激素剂量*/

INPUT X @@; /*X为子宫重量*/

OUTPUT;

END;

END;

CARDS;

106 116 145

42 68 115

70 111 133

42 63 87

;

RUN;

ods html; /*将结果输出成网页格式，SAS9.0以后版本可用*/ PROC GLM DATA=PGM15G;

CLASS A B;

MODEL X=A B / SS3;

MEANS A B; /*给出因素A、B各水平下的均值和标准差*/

MEANS B / SNK; /*对因素B（即剂量）各水平下的均值进行两两比较*/ RUN;

ODS HTML CLOSE;

2．2*3析因设计两因素完全随机统计方法 2*3析因设计 tiff =f的开方DATA aaa;

DO zs=125,200;

DO repeat=1TO2; /*每种试验条件下有2次独立重复试验*/

do js=0.015,0.030,0.045;

INPUT cl @@;

OUTPUT;

END;

END;

END;

CARDS;

2.70 2.45 2.60

2.78 2.49 2.72

2.83 2.85 2.86

2.86 2.80 2.87

;

run;

PROC GLM;

CLASS zs js;

MODEL cl=zs js zs*js / SS3;

MEANS zs*js;

LSMEANS zs*js / TDIFF PDIFF; /*对 zs和js各水平组合而成的试验条件进行均数进行两两比较*/

RUN;

ODS HTML CLOSE;

练习一：2*2横断面研究列链表方法：卡方矫正卡方 FISHER

DATA PGM19A;

DO A=1TO2;

DO B=1TO2;

INPUT F @@;

OUTPUT;

END;

END;

CARDS;

2 26

8 21

;

run;

PROC FREQ;

WEIGHT F;

TABLES A*B / CHISQ;

RUN;

样本大小 = 57

练习二：对裂列连表结果变量换和不换三部曲 1横断面研究 P《0.05 R

DATA PGM19B;

DO A=1TO2;

DO B=1TO2;

INPUT F @@;

OUTPUT;

END;

END;

CARDS;

40 3414

1 19252

;

run;

ods html;

PROC FREQ;

WEIGHT F;

TABLES A*B / CHISQ cmh;

RUN;

ods html close;

样本大小 = 57

练习三：病例对照2*2 病例组中有何没有那个基因是正常的3.8倍，则有可能导致痴呆要做前瞻性研究用对裂

DATA PGM20;

DO A=1TO2;

DO B=1TO2;

INPUT F @@;

OUTPUT;

END;

END;

CARDS;

240 60

360 340

;

run;

ods html;

PROC FREQ;

WEIGHT F;

TABLES A*B / CHISQ cmh;

RUN;

ods html close;

总样本大小 = 1000

练习四：配对设计隐含金标准2*2 MC卡方检验 34和0在总体上（B+C《40 用矫正卡方）是否相等则可得甲培养基优于乙培养基一般都用矫正因卡方为近似计算

DATA PGM19F;

INPUT b c;

chi=(ABS(b-c)-1)**2/(b+c);

p=1-PROBCHI(chi,1);求概率 1减掉从左侧积分到卡方的值

chi=ROUND(chi, 0.001);

IF p>0.0001THEN p=ROUND(p,0.0001);

FILE PRINT;

PUT（打印在输出床口） #2 @10'Chisq' @30'P value'（#表示行）

#4 @10 chi @30 p;

CARDS;

34 0

;

run;

ods html close;

练习五：双向有序R*C列连表用KPA data aaa;

do a=1to3;

do b=1to3;

input f @@;

output;

end;

end;

cards;

58 2 3

1 4

2 7

8 9 17

;

run;

ods html;

*简单kappa检验;

proc freq data=aaa;

weight f;（频数）

tables a*b;

test kappa;

run;

FREQ 过程