2019精选医学SAS练习题及程序答案.doc 下载本文

内容发布更新时间 : 2024/12/24 9:34:21星期一 下面是文章的全部内容请认真阅读。

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

DATA PGM15G;

DO A=1 TO 4; /*A为窝别*/

DO B=1 TO 3; /*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=1 TO 2; /*每种试验条件下有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=1 TO 2; DO B=1 TO 2; 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 DATA PGM19B; DO A=1 TO 2; DO B=1 TO 2; INPUT F @@;

OUTPUT;

END;

END;

CARDS;

40 3414 1 19252 ; run; ods html;

R

PROC FREQ;

WEIGHT F;

TABLES A*B / CHISQ cmh;

RUN;

ods html close;

样本大小 = 57

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

DATA PGM20;

DO A=1 TO 2;

DO B=1 TO 2;

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.0001 THEN p=ROUND(p,0.0001); FILE PRINT;

PUT(打印在输出床口) #2 @10 'Chisq' @30 'P value'(#表示行) #4 @10 chi @30 p; ;

CARDS; 34 0