1．随机取组有无重复试验的两种本题是无重复DATAPGM15G;DOA=1TO4;/*A为窝别*/DOB=1TO3;/*B为雌激素剂量*/INPUTX@@;/*X为子宫重量*/OUTPUT;END;END;CARDS;106116145426811570111133426387;RUN;odshtml;/*将结果输出成网页格式，SAS9.0以后版本可用*/PROCGLMDATA=PGM15G;CLASSAB;MODELX=AB/SS3;MEANSAB;/*给出因素A、B各水平下的均值和标准差*/MEANSB/SNK;/*对因素B（即剂量）各水平下的均值进行两两比较*/RUN;ODSHTMLCLOSE;2．2*3析因设计两因素完全随机统计方法2*3析因设计tiff=f的开方DATAaaa;DOzs=125,200;DOrepeat=1TO2;/*每种试验条件下有2次独立重复试验*/dojs=0.015,0.030,0.045;INPUTcl@@;OUTPUT;END;END;END;CARDS;2.702.452.602.782.492.722.832.852.862.862.802.87;run;PROCGLM;CLASSzsjs;MODELcl=zsjszs*js/SS3;MEANSzs*js;LSMEANSzs*js/TDIFFPDIFF;/*对zs和js各水平组合而成的试验条件进行均数进行两两比较*/RUN;ODSHTMLCLOSE;练习一：2*2横断面研究列链表方法：卡方矫正卡方FISHERDATAPGM19A;DOA=1TO2;DOB=1TO2;INPUTF@@;OUTPUT;END;END;CARDS;226821;run;PROCFREQ;WEIGHTF;TABLESA*B/CHISQ;RUN;样本大小=57练习二：对裂列连表结果变量换和不换三部曲1横断面研究P《0.05RDATAPGM19B;DOA=1TO2;DOB=1TO2;INPUTF@@;OUTPUT;END;END;CARDS;403414119252;run;odshtml;PROCFREQ;WEIGHTF;TABLESA*B/CHISQcmh;RUN;odshtmlclose;样本大小=57练习三：病例对照2*2病例组中有何没有那个基因是正常的3.8倍，则有可能导致痴呆要做前瞻性研究用对裂DATAPGM20;DOA=1TO2;DOB=1TO2;INPUTF@@;OUTPUT;END;END;CARDS;24060360340;run;odshtml;PROCFREQ;WEIGHTF;TABLESA*B/CHISQcmh;RUN;odshtmlclose;总样本大小=1000练习四：配对设计隐含金标准2*2MC卡方检验34和0在总体上（B+C《40用矫正卡方）是否相等则可得甲培养基优于乙培养基一般都用矫正因卡方为近似计算DATAPGM19F;INPUTbc;chi=(ABS(b-c)-1)**2/(b+c);p=1-PROBCHI(chi,1);求概率1减掉从左侧积分到卡方的值chi=ROUND(chi,0.001);IFp>0.0001THENp=ROUND(p,0.0001);FILEPRINT;PUT（打印在输出床口）#2@10'Chisq'@30'Pvalue'（#表示行）#4@10chi@30p;CARDS;340;run;odshtmlclose;练习五：双向有序R*C列连表用KPAdataaaa;doa=1to3;dob=1to3;inputf@@;output;end;end;cards;582314278917;run;odshtml;*简单kappa检验;procfreqdata=aaa;weightf;（频数）tablesa*b;testkappa;run;*加权kappa检验;procfreq;weightf;tablesa*b;testwtkap;run;odshtmlclose;SAS系统FREQ过程频数百分比行百分比列百分比a*b表ab合计12315839.4621.3632.046342.86a*b表ab合计12392.0686.573.173.774.7611.11210.682.001.494228.5784.0079.2574.7614.0025.935034.01385.4423.5311.9496.1226.4716.981711.5650.0062.963423.13合计6745.585336.052718.37147100.00a*b表的统计量对称性检验统计量(S)2.8561自由度3Pr>S0.4144对称性检验指总体上主对角线的上三角数相加是否与下三角三个数相加对称性检验与KPA检验是否一致是否一个可以代替另一个检验Pe理论观察一致率独立假设性基础上计算的相互独立简单Kappa系数Kappa0.6809渐近标准误差0.050095%置信下限0.5830简单Kappa系数95%置信上限0.7788H0检验:Kappa=0总体的H0下的渐近标准误差0.0597Z11.4112单侧Pr>Z<.0001双侧Pr>|Z|<.0001总体的KPA是否为0KPA大于0两种方法的一致性有统计学意义小于0不一致性有统计学意义加权的Kappa系数加权的Kappa0.6614渐近标准误差0.056095%置信下限0.551695%置信上限0.7711置信区间不包括0拒绝H0按此计算结果可以用一种取代另一种方法但要看专业要求达到多少才可以观测一致率达到多少才可以代替样本大小=147FREQ过程频数百分比行百分比列百分比a*b表ab合计12315839.4692.0621.363.1732.044.766342.86a*b表ab合计12386.573.7711.11210.682.001.494228.5784.0079.2574.7614.0025.935034.01385.4423.5311.9496.1226.4716.981711.5650.0062.963423.13合计6745.585336.052718.37147100.00a*b表的统计量对称性检验统计量(S)2.8561自由度3Pr>S0.4144简单Kappa系数Kappa0.6809渐近标准误差0.05009...
