医学统计学孙振球(第三版)SAS程序及结果.docx

《医学统计学孙振球(第三版)SAS程序及结果.docx》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、1.两独立样本inputcx@@;procunivariatenormal;varx;classc;时间资料、百分率资料不服从正态分布T检验（正态）procttest;varx;classc;run;TheTTESTProcedure均数、标准差、标准误T-Tests:Pooled(方差齐，t)Satterthwaite（方差不齐，近似t）EqualityofVariances方差齐性检验非参（非正态）procnpar1waywilcoxon;varx;classc;run;WilcoxonTwo-Sa

2、mpleTestZOne-SidedPr>ZTwo-SidedPr>

3、Z

4、2.配对样本t检验inputx1x2@@;d=x1-x2;cards;procunivariatedata=ex3_6;vard;run;非参（非正态）单样本中位数与总体中位数inputx1@@;median=45.30;d=x1-median;cards;procunivariatenormal;vard;run;procmeansp50;varx1;run;MomentsBasicStatisticalMeasuresStdD

5、eviation标准偏差TestsforLocation:Mu0=0ØStudent'stt检验SignØSignedRank非参（非正态时）Quantiles(Definition5)ExtremeObservations3.方差分析1）完全随机procanova;classc;modelx=c;meansc/lsdsnkdunnett;meansc/hovtestwelch;run;方差齐性检验:Levene'sTestforHomogeneityofxVarianceANOVAofSquaredD

6、eviationsfromGroupMeans若p>0.05,齐用F检验若<0.05，不齐用Welch'sANOVAforx非参：Kruskal-WallishH检验H0：总体分布位置相同procnpar1waywilcoxon;varx;classc;run;Kruskal-WallisTest2）随机区组设计H0：μ1=μ2=μ3inputxab@@;cards;procanova;classab;modelx=ab;meansa/snk;run;注意：SS组间=SS处理间+SS区组间非参：Frie

7、dmanM检验procfreq;tablesb*a*x/scores=rankcmh2;run;注：区组*处理*指标3）两阶段交叉H0：μ1=μ2inputrtimetreat$x@@;cards;procanova;classrtimetreat;modelx=rtimetreat;meansrtimetreat/snk;run;4）析因设计方差分析2×2析因dataf22;inputabx@@;cards;000.80000.90000.70101.30101.20101.10010.90011.1

8、0011.00112.10112.20112.00;procprint;procanova;classab;modelx=aba*b;meansaba*b;run;H0：A因素效应=0，B因素效应=0A*B因素效应=0a:F=168.75p<.001拒绝Hb:F=90.75p<0.001拒绝Ha*b:F=36p=0.0003拒绝Ho:.甲有效，乙有效，甲乙有交互作用，甲乙都用时血红细胞增加数最多。3×3析因设计dataex11_2;inputxab@@;cards;procanova;classab;m

9、odelx=aba*b;meansaba*b/snk;run;5*2*2析因设计procanova;classabcx@@;modelx=abca*ba*cb*ca*b*c;run;snk中只有5和4只对一个字母说明5.4差异大，1.2.3介于之问5)正交实验设计有空白列,无重复dataex11_4;inputababcn1n2dx@@;cards;111111186111222295122112291122221194212121291212212196221122183221211288;procp

10、rint;procanova;classabcd;modelx=abca*bd;meansabcda*b;run;结果误差是由空白列造成的，一级分解中p>0.05不用看只有Ca*b的p<0.05，看F值，F越大作用越大,c为主ab为次因x越大越好，看x的mean值选择较大者，所以a2b1c2d2最优无空白列,有重复dataex1_2;inputabcdx@@;cards;11111222133321232231231231323213332