SAS练习题及程序答案.docx
《SAS练习题及程序答案.docx》由会员分享,可在线阅读,更多相关《SAS练习题及程序答案.docx(20页珍藏版)》请在冰豆网上搜索。
![SAS练习题及程序答案.docx](https://file1.bdocx.com/fileroot1/2023-7/3/89c5ffac-94ed-4668-8525-93fe35da7410/89c5ffac-94ed-4668-8525-93fe35da74101.gif)
SAS练习题及程序答案
1.随机取组有无重复试验的两种本题是无重复
DATAPGM15G;
DOA=1TO4;/*A为窝别*/
DOB=1TO3;/*B为雌激素剂量*/
INPUTX@@;/*X为子宫重量*/
OUTPUT;
END;
END;
CARDS;
106116145
4268115
70111133
426387
;
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.60
2.782.492.72
2.832.852.86
2.862.802.87
;
run;
PROCGLM;
CLASSzsjs;
MODELcl=zsjszs*js/SS3;
MEANSzs*js;
LSMEANSzs*js/TDIFFPDIFF;/*对zs和js各水平组合而成的试验条件进行均数进行两两比较*/
RUN;
ODSHTMLCLOSE;
练习一:
2*2横断面研究列链表方法:
卡方矫正卡方FISHER
DATAPGM19A;
DOA=1TO2;
DOB=1TO2;
INPUTF@@;
OUTPUT;
END;
END;
CARDS;
226
821
;
run;
PROCFREQ;
WEIGHTF;
TABLESA*B/CHISQ;
RUN;
样本大小=57
练习二:
对裂列连表结果变量换和不换三部曲1横断面研究P《0.05R
DATAPGM19B;
DOA=1TO2;
DOB=1TO2;
INPUTF@@;
OUTPUT;
END;
END;
CARDS;
403414
119252
;
run;
odshtml;
PROCFREQ;
WEIGHTF;
TABLESA*B/CHISQcmh;
RUN;
odshtmlclose;
样本大小=57
练习三:
病例对照2*2病例组中有何没有那个基因是正常的3.8倍,则有可能导致痴呆要做前瞻性研究用对裂
DATAPGM20;
DOA=1TO2;
DOB=1TO2;
INPUTF@@;
OUTPUT;
END;
END;
CARDS;
24060
360340
;
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列连表用KPA
dataaaa;
doa=1to3;
dob=1to3;
inputf@@;
output;
end;
end;
cards;
5823
1427
8917
;
run;
odshtml;
*简单kappa检验;
procfreqdata=aaa;
weightf;(频数)
tablesa*b;
testkappa;
run;
*加权kappa检验;
procfreq;
weightf;
tablesa*b;
testwtkap;
run;
odshtmlclose;
SAS系统
FREQ过程
频数
百分比
行百分比
列百分比
a*b表
a
b
合计
1
2
3
1
58
39.46
92.06
86.57
2
1.36
3.17
3.77
3
2.04
4.76
11.11
63
42.86
a*b表
a
b
合计
1
2
3
2
1
0.68
2.00
1.49
42
28.57
84.00
79.25
7
4.76
14.00
25.93
50
34.01
3
8
5.44
23.53
11.94
9
6.12
26.47
16.98
17
11.56
50.00
62.96
34
23.13
合计
67
45.58
53
36.05
27
18.37
147
100.00
a*b表的统计量
对称性检验
统计量(S)
2.8561
自由度
3
Pr>S
0.4144
对称性检验指总体上主对角线的上三角数相加是否与下三角三个数相加对称性检验与KPA检验是否一致是否一个可以代替另一个检验Pe理论观察一致率独立假设性基础上计算的相互独立
简单Kappa系数
Kappa
0.6809
渐近标准误差
0.0500
95%置信下限
0.5830
95%置信上限
0.7788
H0检验:
Kappa=0
总体的H0下的渐近标准误差
0.0597
Z
11.4112
单侧Pr>Z
<.0001
H0检验:
Kappa=0
双侧Pr>|Z|
<.0001
总体的KPA是否为0KPA大于0两种方法的一致性有统计学意义小于0不一致性有统计学意义
加权的Kappa系数
加权的Kappa
0.6614
渐近标准误差
0.0560
95%置信下限
0.5516
95%置信上限
0.7711
置信区间不包括0拒绝H0按此计算结果可以用一种取代另一种方法但要看专业要求达到多少才可以观测一致率达到多少才可以代替
样本大小=147
FREQ过程
频数
百分比
行百分比
列百分比
a*b表
a
b
合计
1
2
3
1
58
39.46
92.06
86.57
2
1.36
3.17
3.77
3
2.04
4.76
11.11
63
42.86
2
1
0.68
2.00
1.49
42
28.57
84.00
79.25
7
4.76
14.00
25.93
50
34.01
3
8
5.44
23.53
11.94
9
6.12
26.47
16.98
17
11.56
50.00
62.96
34
23.13
合计
67
45.58
53
36.05
27
18.37
147
100.00
a*b表的统计量
对称性检验
统计量(S)
2.8561
自由度
3
Pr>S
0.4144
简单Kappa系数
Kappa
0.6809
渐近标准误差
0.0500
95%置信下限
0.5830
95%置信上限
0.7788
加权的Kappa系数
加权的Kappa
0.6614
渐近标准误差
0.0560
95%置信下限
0.5516
95%置信上限
0.7711
H0检验:
加权的Kappa=0
H0下的渐近标准误差
0.0646
Z
10.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;
431490902
388410800
495587950
137179325
;
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;
1541
49915
315045
52224
;
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;
215131148
67101128
4463132
;
run;
odshtml;
PROCCORRSPEARMAN;
VARAB;
FREQF;
RUN;
odshtmlclose;
统计分析与SAS实现第1次上机实习题
一、定量资料上机实习题
要求:
(1)先判断定量资料所对应的实验设计类型;
(2)假定资料满足参数检验的前提条件,请选用相应设计的定量资料的方差分析,并用SAS软件实现统计计算;
(3)摘录主要计算结果并合理解释,给出统计学结论和专业结论。
【练习1】取4窝不同种系未成年的大白鼠,每窝3只,随机分配到三个实验组中,分别注射不同剂量雌激素,经过一定时间后处死大白鼠测子宫重量,资料见表1。
问剂量和窝别的各自水平下子宫重量之间的差别有无统计学意义?
若剂量间差别有统计学意义,请作两两比较。
表1未成年大白鼠注射不同剂量雌激素后的子宫重量
窝别
子宫重量(g)
剂量(μg/100g):
0.2
0.4
0.8
合计
1
106
116
145
367
2
42
68
115
225
3
70
111
133
314
4
42
63
87
192
合计
260
358
480
1098
【SAS程序】:
程序1
【练习2】一位工程师研究由钻头压力产生的冲力。
考察了A(钻孔速度)和B(进料速度),两因素分别取2与3水平,各水平组合下均做了两次独立重复实验,资料见表2。
假定资料满足参数检验的前提条件,且两因素对观测结果的影响地位平等,已知冲力越小越好,试作分析,尽可能给出较为明确的统计和专业结论。
表2在钻孔速度和进料速度取不同水平的条件下冲力的测定结果
钻孔冲力(单位)
速度进料速度:
0.0150.0300.045
1252.702.452.60
2.782.492.72
2002.832.852.86
2.862.802.87
【SAS程序】:
程序2.
二、定性资料上机实习题
要求:
(1)若题目中未给出表格,请列出标准的列联表,并对其命名;
(2)若题目中已列出不规范的表格,先修改,然后对其命名;
(3)根据分析目的或自己提出分析目的、资料的前提条件选用相应的统计分析方法,并用SAS软件实现计算;
(4)将主要计算结果摘录出来,给出统计学和专业结论。
【练习1】某卫生防疫站对屠宰场及肉食零售点的猪肉,检查其表层沙门氏菌带菌情况,结果如下表。
试比较屠宰场与肉食零售点猪肉表层沙门氏菌的带菌率之间差别有无统计学意义?
表1屠宰场及肉食零售点猪肉表层沙门氏菌抽检结果
采样地点
检查数
阳性数值
带菌率(%)
屠宰场
28
2
7.14
零售点
29
8
27.59
【SAS程序】:
练习1
【练习2】有人对某部门22707名雇员中,普查了HBsAg,其中3454名阳性,19253名为阴性。
从1975年起,追踪了3年,发现在阳性组有40名患了肝癌,阴性组仅一名患肝癌。
试选用合适的方法对资料进行全面分析。
【SAS程序】:
练习2
【练习3】APOE-4等位基因与老年痴呆性的关联研究:
以600名晚发及散发老年痴呆患者和400名正常对照为研究对象,分析APOE-4等位基因与老年痴呆性的关系。
表4APOE-4等位基因与老年痴呆性病例对照关联研究
AGT等位基因
例数
合计
病例组
对照组
APOE-4
240
60
300
非APOE-4
360
340
700
合计
600
400
1000
【SAS程序】:
练习3
【练习4】请分析下表资料。
已从专业上认定培养的阳性结果就是“真阳性”,而不会出现假阳性。
表3两种培养基对同一批痰液标本同时培养的结果
甲培养基
培养结果
例数
乙培养基结果:
+
-
合计
+
36
34
70
-
0
135
135
合计
36
169
205
配对设计2×2列联表资料总体率差异性检验统计量的计算公式
若b+c≥40时
若b+c<40时
【SAS程序】:
练习4
【练习5】请分析下表资料。
表6两法检查室壁收缩运动的符合情况
━━━━━━━━━━━━━━━━━━━
对比法测冠心病人数
定的结果核素法∶正常减弱异常
───────────────────
正常5823
减弱1427
异常8917
───────────────
合计675327
━━━━━━━━━━━━━━━━━━━
【SAS程序】:
练习5
简单kappa检验和加权kappa检验这两种方法都是用来检验两种评价方法是否具有一致性的方法。
其主要的区别是两种方法计算的公式不一样,更具体地说是对两个变量的打分不一样,简单kappa检验主要是利用对角线上的信息,加权kappa检验除了利用对角线上的数据外,还将对角线外的数据进行加权打分,将对角线外的信息也充分利用。
所以在选择方法时应根据专业知识,如果两个变量取值的界线比较明确,如“+”“++”“+++”等,这时可以选用简单的kappa检验,如果两个变量的取值不十分明确,人为的因素较多时,就可选用加权kappa检验。
【练习6】请分析下表资料。
表5某地6094人按2种血型系统划分的结果
━━━━━━━━━━━━━━━━━━━
ABO人数
血型MN血型:
MNMN
───────────────────
O431490902
A388410800
B495587950
AB137179325
───────────────
合计145116662977
━━━━━━━━━━━━━━━━━━━
【SAS程序】:
练习6
•分析方法选择
不超过1/5的格子理论频数<5时,此类资料应采用
检验处理。
其他情形,可采用Fisher精确检验处理。
【练习7】请分析下表资料。
表73种药物疗效的观察结果
━━━━━━━━━━━━━━━━━━━
疗人数
效药物∶ABC
───────────────────
治愈1541
显效49915
好转315045
无效52224
───────────────
合计1008585
━━━━━━━━━━━━━━━━━━━
【SAS程序】:
练习7
【练习8】请分析下表资料。
表8眼晶状体混浊度与年龄之关系
━━━━━━━━━━━━━━━━━━━
晶状体混眼数
浊程度年龄∶20~30~40~
───────────────────
+215131148
++67101128
+++4463132
───────────────
合计326295408
━━━━━━━━━━━━━━━━━━━
(注:
本资料素材和资料部分来自网络,仅供参考。
请预览后才下载,期待你的好评与关注!
)