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