SAS上机练习题(习题部分).docVIP

重庆医科大学卫生统计学统计软件包SAS上机实习题重庆医科大学卫生统计学教研室 彭斌 编写2010年12月卫生统计学统计软件包上机实习题SAS上机练习题（一）1、SAS常用的窗口有哪三个？请在三个基本窗口之间切换并记住这些命令或功能键。2、请在PGM窗口中输入如下几行程序，提交系统执行，并查看OUTPUT窗和LOG窗中内容，注意不同颜色的含义；并根据日志窗中的信息修改完善程序。 DATS EX01; INPUTT NAME $ AGE SEX; CARDS; XIAOMIN 19 1 LIDONG 20 1 NANA 18 2 ; PROD PRONT DATS=EX1; RUN; PROC PRINT DATA=EX1; VAR NAME AGE; RUN; 3、将第2题的程序、结果及日志保存到磁盘。4、试根据某班12名学生3门功课成绩表完成后面的问题： 表1 某班12名学生3门功课成绩表学号生化物理病理083687165084746168085737546087757168084858587085787975086807679089677371118705475083706684084627365099827079问题：1） 建立数据集；2） 打印至少有1门功课不及格同学的信息；(提示，使用if语句)5、（选做）根据下列数据建立数据集表2 销售数据开始时间终止时间费用2005/04/2825MAY2009$123,345,0002005 09 1805OCT2009$33,234,5002004050830JUN2009$432,334,500提示：(格式化输入；数据之间以空格分隔，数据对齐；注意格式后面的长度应以前一个位置结束开始计算，如果读入错误，可试着调整格式的宽度；显示日期需要使用输出格式) 开始时间，输入格式yymmdd10. 终止时间，输入格式date10. 费用，输入格式dollar12. 例如：input x1 yymmdd10. X2 date10. X3 dollar12. 6、手机号码一编码规则一般是：YYY-XXXX-ZZZZ，其YYY为号段；XXXX一般为所在地区编码；ZZZZ为对应的个人识别编号。下面有一组电话号码（来源于网络，末位以X替换），请用程序完成下列要求：（1）分别列出属于联通、移动、电信的号码；（2）分别提取地区编号及个人识别编号。（提示：列输入方式或者字符串操作substr()函数） 说明：移动：134-139、150、151、152、157、158、159、188联通：130、131、132、155、156电信：133、153、180、189手机号码手机号码1508320464X1313654836X1510291126X1351425709X1339644156X1303612231X1532464828X1383824855X1365721647X1522322763X1523105754X1357851051X1592624347X1508311759X1331237668X1327313520X1370048578X1556443719X1507244457X1804346016X1321246707X1513441713X1368464734X1308279203X1394913242X1897876466X7、下表是某班学生几门功课的成绩，其中形势课是考查课，其它均为考试课。请完成以下处理并且保存程序、结果和日志：（注意函数在DATA STEP中使用的位置） （1）用Mean()函数求出每位同学的平均分，将其保存在变量中；（2）用sum()函数或者表达式求出每位同学的总分，将其保存在变量中；（3）如果每门考试课权重是0.75，每门考查课权重为0.25，试求每位同学的加权平均分；（4）将数据集保存为永久数据集，保存位置为“D:mydata”；（5）调用上面的永久数据集，并将男、女生的成绩分别保存在数据集Dmale和Dfemale中。表3 某班同学几门功课的成绩学号性别(0=女，1=男)高数生理人解数理统计形势(考查)0083173736474750414190797185780564197878991800774140606165750873168656076750874174685660750875173466566870876179797489870877175685560700878176606471800880181736075820881176818489780883177606066750884085878485900885078756266800886080797773830887085828184800888077757467800889067717169900890075707569800891069767679900892066716060788、下面是3个大类疾病的ICD-10编码及对应的疾病名。请完成以下任务： （1）建立数据集；（2）提取每种疾病的大类编码；（3）分别将3个大类的疾病存入3个数据集。（提示：ICD10编码中小数点前面的三位表示大类；length语句定义字符变量长度；字符串取子串函数）表4 4类疾病的ICD10编码及对应疾病名ICD10DiseaseA01.001伤寒A01.002伤寒杆菌性败血症A01.101甲型副伤寒A01.201乙型副伤寒A01.301丙型副伤寒A01.401副伤寒A02.001B群沙门氏菌肠炎A02.002C群沙门氏菌肠炎A02.004沙门氏菌性肠炎A02.006沙门氏菌胃肠炎A02.007鼠伤寒沙门氏菌性肠炎A02.008婴儿沙门氏菌肠炎A02.101沙门氏菌败血症A02.901沙门氏菌感染A02.902鼠伤寒沙门氏菌感染A02.903沙门氏菌属食物中毒A03.001什密氏志贺菌痢疾A03.101弗氏志贺菌痢疾A03.201鲍氏志贺菌痢疾A03.301宋内氏志贺菌痢疾A03.802菌痢混合感染A03.901细菌性痢疾A03.902慢性细菌性痢疾急性发作A03.903慢性细菌性痢疾A03.906中毒性痢疾9、现有两个文件，内容如下：A文件学号 统计 体育 年龄1903 87 83 231902 56 96 221901 93 75 191904 77 84 241908 88 55 18B文件学号 性别 班级1901 男 11902 女 11903 男 31904 男 2 1905 女 21908 女 3要求：(1)根据A、B文件分别建立数据集A、B；(2)打印至少有一门功课不及格的同学的年龄、性别和班级。 10、某研究员欲分析急性染毒对肝脏功能的影响，将40只小鼠分为两组，雄雌各半，试验组进行急性染毒试验，染毒后2小时测定血液中的ALT，整理的结果见下表。 ALT(丙氨酸转氨酶)NOMaleFemaleControlTreatmentControlTreatment138.457.839.755.7230.662.434.663.8325.557.540.559.9429.949.944.349.7529.055.641.248.8633.647.325.651.3739.450.130.653.8842.548.737.642.6946.757.939.561.71049.360.425.165.3要求：根据上面的数据建立恰当的SAS数据集11、某职业病防治所对30名矿工分别测定血清铜蓝蛋白含量（mol/L），资料如下。问各期血清铜蓝蛋白含量的测定结果有无差别？ 疾病分期测定结果0期8.0 9.0 5.8 6.3 5.4 8.5 5.6 5.4 5.5 7.2 5.6 I期8.5 4.3 11.0 9.0 6.7 9.0 10.5 7.7 7.7II期11.3 7.0 9.5 8.5 9.6 10.8 9.0 12.6 13.9 6.5要求：根据上面的数据建立恰当的SAS数据集 12、24名志愿者随机分成两组，每组12人，接受降胆固醇试验，甲组为特殊饮食组，乙组为药物治疗组。受试者试验前后各测量一次血清胆固醇（mmol/L），数据见下表： 甲组乙组受试者试验前试验后受试者试验前试验后16.116.0016.906.9326.816.8326.406.3536.484.4936.486.4147.597.2847.007.1056.426.3056.536.4166.946.6466.706.6879.178.4279.109.0587.337.0087.316.8396.946.5896.966.91107.677.22106.816.73118.156.57118.167.65126.606.17126.986.52要求：根据上面的数据建立恰当的SAS数据集(彭斌，2010-10-10)SAS上机练习题（二）1、某农村地区1999年14岁女孩的身高资料列于表1。表1 某农村地区1999年14岁女孩身高资料（cm）142.3148.8142.7144.4144.7145.1143.3154.2152.3142.7156.6137.9143.9141.2139.3145.8142.2137.9141.2150.6142.7151.3142.4141.5141.9147.9125.8139.9148.9154.9145.7140.8139.6148.8147.8146.7132.7149.7154.0158.2138.2149.8151.1140.1140.5143.4152.9147.5147.7162.6141.6143.6144.0150.6138.9150.8147.9136.9146.5130.6142.5149.0145.4139.5148.9144.5141.8148.1145.4134.6130.5145.2146.2146.4142.4137.1141.4144.0129.4142.8132.1141.8143.3143.8134.7147.1140.9137.4142.5146.6135.5146.8156.3150.0147.3142.9141.4134.7138.5146.6134.5135.1141.9142.1138.1134.9146.7138.5139.6139.2148.8150.3140.7143.5140.2143.6138.7138.9143.5139.9134.4133.1145.9139.2137.4142.3160.9137.7142.9126.8问题：（1）计算均数、中位数；计算均数的95%可信区间；（2）计算四分位间距、标准差、变异系数；计算标准误；（3）请进行正态性检验；（4）观察频数分布情况；2、（t检验、配对t检验）24名志愿者随机分成两组，每组12人，接受降胆固醇试验，甲组为特殊饮食组，乙组为药物治疗组。受试者试验前后各测量一次血清胆固醇（mmol/L），数据见下表：甲组乙组受试者试验前试验后受试者试验前试验后16.116.0016.906.9326.816.8326.406.3536.484.4936.486.4147.597.2847.007.1056.426.3056.536.4166.946.6466.706.6879.178.4279.109.0587.337.0087.316.8396.946.5896.966.91107.677.22106.816.73118.156.57118.167.65126.606.17126.986.52问题：（1）试建立能够满足以下统计分析的SAS数据集。（2）试判断两组受试者试验前血清胆固醇水平是否相等？（3）分别判断两种降胆固醇措施是否有效？（4）试判断两种降胆固醇措施的效果是否相同？ 3、（秩和检验）某职业病防治所对30名矿工分别测定血清铜蓝蛋白含量（mol/L），资料如下。问各期血清铜蓝蛋白含量的测定结果有无差别？ 疾病分期测定结果0期8.0 9.0 5.8 6.3 5.4 8.5 5.6 5.4 5.5 7.2 5.6 I期8.5 4.3 11.0 9.0 6.7 9.0 10.5 7.7 7.7II期11.3 7.0 9.5 8.5 9.6 10.8 9.0 12.6 13.9 6.5请用秩和检验完成分析。 4、（频数表等）以下是对某年级的问卷调查结果，问题表内容如下：id机器号sex性别：1=male,2=femaledate记忆中影响最大的日子是哪一天？height身高是多少weight体重是多少money昨天吃饭花多少钱（元）？question146级改革之后学习英语的积极性：a增加，b不变，c减弱question2结核病的早期症状有(可以多选)：a=咳嗽,b=痰中带血,c=盗汗,d=贫血,e=胸痛,f=午后潮热,g=胃肠不适请进行以下分析：1） 此次调查中男性的比例为多少；2） 对46级改革之后学习英语的积极性的情况进行描述，三种选择占的比例分别是多少？3） 男生与女性的对46级改革的反应有无差异；4） 计算身体质量指数（Body Mass Index，），比较男女生的BMI有无差异？5） 关于BMI与人体的肥胖程度通常可以根据下面的标准来判断，请列出每一类的人数及比例；偏瘦： BMI指数 30中度肥胖： BMI指数 35 重度肥胖： BMI指数 40原始数据机器号性别身高体重开支英语改革TB症状idsexheightweightmoneyquestion1question2cnw1l011179705.7aABDEFcnw1l021175707.5aABEcnw1l032157474.5aABEcnw1l042163485cABDFcnw1l052161525bABFcnw1l062150445.5bABDcnw1r012158564.4bADFGcnw1r022158566.2bABDEcnw1r031168646.5aABEcnw1r052150445.5bABDcnw1r061162537bABDEcnw2l022160538cABEFcnw2l032164586cABEcnw2l041167607cABCDEFcnw2l052164545cABcnw2r022151446bBcnw2r032150517.8bABDFcnw2r042160555bABEcnw2r052160473.9bABEFcnw2r062156488.5cABCDEFGcnw3l012162488cABDFcnw3l021163607.2aABEcnw3l031161507.3bADFGcnw3r011169575.3aABDFcnw3r022162.546.57aABCDEFcnw3r032155537.5bABCDEFcnw3r052163504.4bABEDcnw3r061163658bABCDEFcnw4l011168605cABDGcnw4l021178688aACcnw4l032158456cADEFcnw4l041173688.5bABDEFcnw4r011175535bABDEFcnw4r021178696bAcnw4r032158508bABDEFcnw4r042165534.5bABDEcnw4r052157465.3bABEcnw4r062155486.7cABDcnw5l012160494.2bABDEFcnw5l022160497aABEcnw5l032170578.5cABcnw5r012160494.2bABDEcnw5r021168659.2aABDEFcnw5r022170578.5cADcnw6l022160504bABEcnw6l031181633.5bADEFcnw6l0411755320aABcnw6l051174606.5bADcnw6l061162524.7bABEcnw6r0121545410aABDFcnw6r022166544.5cABDFcnw6r032155435.4cABEFcnw6r042153455.4cEGcnw6r052160485aABcnw6r062164467.4bABEFcnw7l022160468bAcnw7l051174707bADFcnw7r012155558cABcnw7r0211686010cABEFcnw7r0311816010cABEcnw7r051166607cABDEcnw7r062164467.4bABEFcnw8l011160509cADFcnw8r011174708aADEGcnw8r022172535cABEcnw1r012154494.2CABDEFcnw1r0211696410.4BABDEcnw1r0411696715AABDEFGcnw1r0521535010BABCDcnw1r062160503.6AABEcnw2r022156464AADGcnw2r031173627.5BABFGcnw2r041168578.9BABDFcnw2r052155456BABDcnw2r061170607CABDcnw3r021162566.3CABDcnw3r032161507BABEcnw3r042154485CABEcnw3r052155466BABDFcnw3r062155515AABDEGcnw4r022153503.5CBEcnw4r032150485AABDEcnw4r042152.5465.2AADEcnw4r051185707.5AABDGcnw4r061170586.5BADcnw5r011173626BABDEcnw5r021168506.1AABDEcnw5r032160507CABCEFcnw5r042155426BABDFcnw5r051168557BADFcnw5r061173618BABDcnw6r012158508.1CADFcnw6r0211756110.5BABDEcnw6r0321554411CABDEcnw6r0411626010CADFcnw6r0511736415CADFcnw7r0111758611CABEcnw7r0211726212BABEGcnw7r031172656.3BABEGcnw7r052151487BABDEFGcnw8r021167509BABDEFcnw8r032160505BABDEFGcnw8r0421706510BABDEFGcnw1l011165477.3BABcnw1l022152444.2BABcnw1l03216249.54.5BABEcnw1l042157523.2BABEcnw1l052150453.8AABEcnw2l022170554.9BADFcnw2l032155475.9BABDFcnw2l042163485.8AADFGcnw2l05215047.58AABEcnw3l012153464BABDFGcnw3l0211655435AADFcnw4l0121605212.4AABDFcnw4l021170655ADcnw4l041165515BAFcnw4l052155428BABCDEcnw5l012160457CABDFcnw5l0211796310CABDFcnw5l041170679CADFcnw6l011170589.6BABCDEFcnw6l022160545.5CABCDEFcnw6l032155484.4AABCDEFcnw6l0411706510CABCDEFcnw6l052165477CABDFcnw7l012150485.8CABEcnw7l021160507AABCDEFcnw7l0311756610AABCDEFcnw7l042163518BABDcnw7l052165574.9AABDcnw8l0111606010BABEcnw8l022154504.3ABEcnw8l032160607AABDEF 补充练习题(摘自医学统计学，颜虹主编，人民卫生出版社，第一版)P132：（t检验）第2，4题2. 某医生研究脑缺氧对脑组织中生化指标的影响，将出生状况相近的乳猪按出生体重配成7对共2组，随机接受两种处理，一组设为对照组，一组设为脑缺氧模型组，试验结果如表。试比较两种处理的脑组织钙泵的含量有无差别。乳猪号对照组实验组10.35500.275520.20000.254530.31300.180040.36300.180050.35440.311360.34500.295570.30500.28704. 随机抽取20只小鼠分配到A，B两个不同饲料组，每组10只，在喂养一定时间后，测得鼠肝中铁的含量（g/g），数据如下。试问不同饲料对鼠肝中铁的含量有无影响。A组3.590.963.891.231.612.941.963.681.542.59B组2.231.142.631.001.352.011.641.131.011.70P149：（方差分析）第1，3，4题1. 将12份同一个体的血液标本随机等分成3组，分别用甲，乙，丙3种抗凝剂处理后测红细胞沉降率。试验结果如下，问3种抗凝剂作用的差异是否有统计学意义？抗凝剂红细胞沉降率（mm/h）甲17161615乙10111212丙119893. 比较下列试验结果（如下表）雌激素3个剂量组之间的差异有无统计学意义。种系雌激素注射剂量0.20.40.8甲106116145乙4268115丙70111133丁4263874.分析下列拉丁方设计的试验结果（下表），其中字母A，B，C，D表示胰岛素注射的4种剂量（U/kg），分别为A：0.32，B：0.47，C：0.62，D：0.77。问4种剂量下的测量结果的差异有无统计学意义？实验日期家兔编号12341C32.7A11.2B23.2D48.12B26.2D31.8C28.9A18.73A14.0C14.0D27.5B25.64D33.2B16.5A21.2C40.2P168：（卡方检验）第1，3题1.某医师用两种疗法治疗脑血管梗死，结果如下表，是比较两种疗法的疗效是否相同？疗法有效无效合计有效率（%）甲疗法2563180.65乙疗法2933290.63合计5496385.713.有50份痰液标本，每份分别接种在甲乙两种培养基中，观察结核杆菌的生产情况，结果如下表，是比较两种培养基的效果。甲培养基乙培养基合计+-+23（a）12（b）35-7（c）8（d）15合计302050P184：（秩和检验）第2，3，6题2. 某研究中心观察溶脲脲原体（UU）感染对家兔精子质量的影响，分别测定十只兔子感染UU前后的精子密度（10E9/L），数据见下表，试分析溶脲脲原体是否影响家兔的精子密度。编号12345678910感染前336371386364377292288304333302感染后2582913002852983033122603392903. 用中草药治疗不同类型的小儿肺炎，其疗效分为4个等级，结果见表，试比较该药物治疗不同类型的小儿肺炎的疗效有无差别？疗效病毒性肺炎细菌性肺炎合计控制6542107显效18624有效302353无效131124合计126822086. 4种疾病患者痰液内嗜酸性粒细胞的检查结果如下表，问4种疾病患者痰液内嗜酸性粒细胞的等级分布有无差别？白细胞等级支气管扩张肺水肿肺癌病毒性呼吸道感染-1475+3697+10655+7241合计21182518 (彭斌，2010-10-10)SAS上机练习题（三）1、（方差分析）某研究人员研究3种刺激因素（高葡萄糖、高胰岛素和过氧化氢）刺激大鼠肾小球系膜细胞对单核细胞趋化蛋白表达的影响，在四种条件下对大鼠肾小球系膜细胞进行培养，检测结果如表1。表1 不同刺激因素作用下蛋白表达结果分组表达结果（光密度值）对照组0.35470.36520.38420.35980.37420.3925高葡萄糖组1.22471.23541.21361.22861.24151.2275高胰岛素组1.14531.15261.13751.15471.12941.1659过氧化氢组1.35821.36571.35181.32591.37261.3318问题1：计算每一组的平均光密度，观察各个组间均数的差异大小；问题2：进行方差分析，将结果整理成方差分析表，观察变异、均方及F值之间的关系；问题3：分别采用SNK法和Dunnett法进行多重比较；问题4：对上述结果进行解释。2、（方差分析）研究枸杞多糖对脂肪肝的预防作用，按窝别作为区组标志，每一区组3只大鼠，随机分配到三个组：生理盐水组、酒精组、酒精+枸杞多糖组，三组灌胃5周，检测肝脏中谷胱甘肽（GSH）的含量（mg/gprot），结果见表2。表2 三组小鼠肝脏中谷胱甘肽（GSH）的含量（mg/gprot）区组酒精组LBP预防组生理盐水130.4865.0879.15231.2563.0475.46333.2867.5979.32434.6168.5875.98528.3564.1276.55629.1766.5580.34727.3466.8984.35830.5867.1588.14934.2568.0587.351027.3165.4872.151128.0964.3874.611230.4565.0486.331333.2566.8494.351434.0467.5692.051534.2567.4696.42问题1：按照随机区组设计进行分析，处理组间、区组间效应是否有差异？问题2：忽略区组，按照完全随机设计进行分析，处理组间效应是否有差异？问题3：对上述两个方差分析表进行比较，找出相同项与不同项？观察处理组间比较的F值变化，会得出怎样的结论？问题4：根据分析结果对该研究下一个结论。3、（方差分析）研究人员研究某种物质的毒性，将40只小鼠分为两组，雄雌各半，试验组给予研究物质，2小时后测定血液中碱性磷酸酶的含量，结果如下。表3 40只小鼠给予不同物质后血液中碱性磷酸酶的含量性别分组碱性磷酸酶雄性对照组367.9408.6375.6354.9421.7374.5432.7401.3399.4367.5试验组423.8446.9432.5478.1437.5421.6489.0432.5421.0420.4雌性对照组378.1345.2390.6399.0421.1341.3322.5365.4321.6401.9试验组420.4473.2450.3405.5427.4460.5420.1394.4389.6420.5问题1：研究者考虑了几个影响因素？问题2：性别对碱性磷酸酶有无影响？问题3：该物质对碱性磷酸酶有无影响？4、（回归分析）以下是制作标准曲线时对应的浓度和吸光度值，试求出标准曲线。并求出A=1.15时的浓度。表4 不同浓度对应的吸光度值浓度A值0.970.31.270.561.570.931.881.352.181.515、(相关分析、回归分析)27名糖尿病人的血清总胆固醇、甘油三脂、空腹胰岛素、糖化血红蛋白、空腹血糖的测量值列于表中，试分析血糖与其它几项指标关系的关系。表5 27名糖尿病人的血糖及有关变量的测量结果序号i总胆固醇甘油三脂胰岛素糖化血血 糖（mmol/L）(mmol/L)(U/ml)红蛋白(%)(mmol/L)X1X2X3X4Y15.681.904.538.211.223.791.647.326.98.836.023.566.9510.812.344.851.075.888.311.654.602.324.057.513.466.050.641.4213.618.374.908.5012.608.511.187.083.006.7511.512.193.852.1116.287.99.6104.650.636.597.18.4114.591.973.618.79.3124.291.976.617.810.6137.971.937.579.98.4146.191.181.426.99.6156.132.0610.3510.510.9165.711.788.538.010.1176.402.404.5310.314.8186.063.6712.797.19.1195.091.032.538.910.8206.131.715.289.910.2215.783.362.968