内容发布更新时间 : 2024/12/24 2:50:51星期一 下面是文章的全部内容请认真阅读。
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;
13．
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
run;
ods html close;
练习五：双向有序R*C列连表 用KPA
data aaa;
do a=1 to 3;
do b=1 to 3;
input f @@; output;
end;
end; cards;
58 2 3 1 42 7 8 9 17 ; run; ods html;
*简单kappa检验;
proc freq data=aaa;
*加权kappa检验; proc freq;
weight f; tables a*b; test wtkap; weight f;（频数） tables a*b; test kappa;
run;
run;
ods html close; SAS 系统
FREQ 过程
频数 百分比 行百分比 列百分比 a * b 表 a 1 1 58 2 b 2 3 3 63 合计 39.46 92.06 86.57 2 1 0.68 2.00 1.49 8 5.44 23.53 11.94 1.36 3.17 3.77 42 28.57 84.00 79.25 9 6.12 26.47 16.98 53 36.05 2.04 4.76 11.11 7 4.76 14.00 25.93 17 11.56 50.00 62.96 27 18.37 42.86 50 34.01 34 23.13 147 100.00 3 合计 67 45.58
a * b 表的统计量
对称性检验 统计量 (S) 自由度 Pr > S 2.8561 3 0.4144 对称性检验指 总体上主对角线的上三角数相加是否与下三角三个数相加 对称性检验与KPA 检验是否一致 是否一个可以代替另一个检验 Pe理论观察一致率 独立假设性基础上计算的 相互独立
简单 Kappa 系数 Kappa 渐近标准误差 95% 置信下限 95% 置信上限
0.6809 0.0500 0.5830 0.7788