SAS练习题及程序答案.doc

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;240 60360 340;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.4692.0686.5721.363.173.7732.044.7611.116342.86  210.682.001.494228.5784.0079.2574.7614.0025.935034.01  385.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.583095%置信上限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.0686.5721.363.173.7732.044.7611.116342.86  210.682.001.494228.5784.0079.2574.7614.0025.935034.01  385.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.050095%置信下限0.583095%置信上限0.7788加权的Kappa系数加权的Kappa0.6614渐近标准误差0.056095%置信下限0.551695%置信上限0.7711H0检验:加权的Kappa=0H0下的渐近标准误差0.0646Z10.2406单侧Pr>Z<.0001双侧Pr>|Z|<.0001对加权的KPA检验与简单的（利用对角线上的数据分析）加权还要利用对角线以外的数据分析样本大小=147练习六：双向无序R*C列连表用卡方理论频数小于5没有超过五分之一，一般用卡方实在不行用FISHER检验超过用KPA两种血型都是按小中大排列相互不影响独立的接受H0不一致行与列变量相互不影响DATAPGM20A; DOA=1TO4; DOB=1TO3; INPUTF@@; OUTPUT; END; END; CARDS;431490902388410800495587950137179325 ;run;odshtml;PROCFREQ; WEIGHTF; TABLESA*B/CHISQ;*exact;RUN;odshtmlclose;样本大小=6094练习七：单向有序R*C秩和检验*方法1;（单因素非参数HO三个药物疗效相同H1不完全相等）DATAPGM20C; DOA=1TO4; DOB=1TO3; INPUTF@@; OUTPUT; END; END; CARDS;15414991531504552224 ;run;odshtml;PROCNPAR1WAYWILCOXON;FREQF; CLASSB; VARA;RUN;*方法2;（FIQCHIM）procfreqdata=PGM20C; weightf; tablesb*a/cmhscores=rank;run;odshtmlclose;总样本大小=270练习八：双向有序属性不同R*C4种目的4种方法SPEARMAN秩相关分析DATAPGM20E; DOA=1TO3; DOB=1TO3; INPUTF@@; OUTPUT; END; END; CARDS;215131148671011284463132 ;run;odshtml;PROCCORRSPEARMAN; VARAB; FREQF;RUN;odshtmlclose;统计分析与SAS实现第1次上机实习题一、定量资料上机实习题要求：（1） 先判断定量资料所对应的实验设计类型；（2） 假定资料满足参数检验的前提条件，请选用相应设计的定量资料的方差分析，并用SAS软件实现统计计算；摘录主要计算结果并合理解释，给出统计学结论和专业结论。

练习1】取4窝不同种系未成年的大白鼠，每窝3只，随机分配到三个实验组中，分别注射不同剂量雌激素，经过一定时间后处死大白鼠测子宫重量，资料见表1问剂量和窝别的各自水平下子宫重量之间的差别有无统计学意义？若剂量间差别有统计学意义，请作两两比较表1未成年大白鼠注射不同剂量雌激素后的子宫重量窝别子宫重量(g)剂量(μg/100g):0.20.40.8合计1106116145367242681152253701111333144426387192合计2603584。