统计学上机练习题.docx

上传人:b****1 文档编号:10881395 上传时间:2023-05-28 格式:DOCX 页数:20 大小:25.59KB
下载 相关 举报
统计学上机练习题.docx_第1页
第1页 / 共20页
统计学上机练习题.docx_第2页
第2页 / 共20页
统计学上机练习题.docx_第3页
第3页 / 共20页
统计学上机练习题.docx_第4页
第4页 / 共20页
统计学上机练习题.docx_第5页
第5页 / 共20页
统计学上机练习题.docx_第6页
第6页 / 共20页
统计学上机练习题.docx_第7页
第7页 / 共20页
统计学上机练习题.docx_第8页
第8页 / 共20页
统计学上机练习题.docx_第9页
第9页 / 共20页
统计学上机练习题.docx_第10页
第10页 / 共20页
统计学上机练习题.docx_第11页
第11页 / 共20页
统计学上机练习题.docx_第12页
第12页 / 共20页
统计学上机练习题.docx_第13页
第13页 / 共20页
统计学上机练习题.docx_第14页
第14页 / 共20页
统计学上机练习题.docx_第15页
第15页 / 共20页
统计学上机练习题.docx_第16页
第16页 / 共20页
统计学上机练习题.docx_第17页
第17页 / 共20页
统计学上机练习题.docx_第18页
第18页 / 共20页
统计学上机练习题.docx_第19页
第19页 / 共20页
统计学上机练习题.docx_第20页
第20页 / 共20页
亲,该文档总共20页,全部预览完了,如果喜欢就下载吧!
下载资源
资源描述

统计学上机练习题.docx

《统计学上机练习题.docx》由会员分享,可在线阅读,更多相关《统计学上机练习题.docx(20页珍藏版)》请在冰点文库上搜索。

统计学上机练习题.docx

统计学上机练习题

一、某医院收集了一批糖尿病患者,随机分成2组,对2组慢性肾炎病人分别用A药(group=1),和B药(group=0),评价指标为连续性变量,用Y表示(Y值越低表明疗效越好),结果见Stata文件(ex1X.dta,X=1,2,…,100)。

问这2种药物治疗慢性肾炎的疗效是否相同?

依题意,本题属完全随机设计,所得指标数据为计量资料,且相互独立,在满足各自应具备的条件的前提下,可以对其进行两组计量资料的t检验或者完全随机设计资料的方差分析以推断两个样本所代表的两总体均数μ1,μ2是否相等。

首先,分别对两组计量资料进行正态性检验:

(文件ex14)

H0:

资料正态分布

H1:

资料非正态分布

α=0.05

.sktestyifgroup==1

Skewness/KurtosistestsforNormality

-------joint------

Variable|Pr(Skewness)Pr(Kurtosis)adjchi2

(2)Prob>chi2

-------------+-------------------------------------------------------

y|0.3000.7871.250.5357

.sktestyifgroup==0

Skewness/KurtosistestsforNormality

-------joint------

Variable|Pr(Skewness)Pr(Kurtosis)adjchi2

(2)Prob>chi2

-------------+-------------------------------------------------------

y|0.5560.7790.440.8028

两组资料的正态性检验的P值均大于0.05,故认为两组资料分别服从正态分布。

其次,对两样被资料作方差齐性检验

H0:

σ12=σ22A药组和B药组总体方差相同

H1:

σ12≠σ22A药组和B药组总体方差不同

α=0.10

.sdtesty,by(group)

------------------------------------------------------------------------------

Group|ObsMeanStd.Err.Std.Dev.[95%Conf.Interval]

---------+--------------------------------------------------------------------

0|17-.0260403.40781321.681457-.8905656.838485

1|22.00305.43157762.024279-.8944648.9005648

---------+--------------------------------------------------------------------

combined|39-.0096304.29764781.85881-.6121869.5929261

------------------------------------------------------------------------------

Ho:

sd(0)=sd

(1)

F(16,21)observed=F_obs=0.690

F(16,21)lowertail=F_L=F_obs=0.690

F(16,21)uppertail=F_U=1/F_obs=1.449

Ha:

sd(0)

(1)Ha:

sd(0)~=sd

(1)Ha:

sd(0)>sd

(1)

PF_U=0.4368P>F_obs=0.7735

两样本方差比F=0.690,P=0.4368>0.10,按α=0.10水准不拒绝H0,差别无统计学意义,尚不能认为两组总体方差不同。

由此,认为两组资料满足两组计量资料的t检验和方差分析的要求。

建立假设:

H0:

2种药物治疗慢性肾炎的疗效均数相同

H1:

2种药物治疗慢性肾炎的疗效均数不相同

α=0.05

t检验:

.ttesty,by(group)

Two-samplettestwithequalvariances

------------------------------------------------------------------------------

Group|ObsMeanStd.Err.Std.Dev.[95%Conf.Interval]

---------+--------------------------------------------------------------------

0|17-.0260403.40781321.681457-.8905656.838485

1|22.00305.43157762.024279-.8944648.9005648

---------+--------------------------------------------------------------------

combined|39-.0096304.29764781.85881-.6121869.5929261

---------+--------------------------------------------------------------------

diff|-.0290903.6082877-1.2615981.203418

------------------------------------------------------------------------------

Degreesoffreedom:

37

Ho:

mean(0)-mean

(1)=diff=0

Ha:

diff<0Ha:

diff~=0Ha:

diff>0

t=-0.0478t=-0.0478t=-0.0478

P|t|=0.9621P>t=0.5189

完全随机设计的方差分析:

.onewayygroup

AnalysisofVariance

SourceSSdfMSFProb>F

------------------------------------------------------------------------

Betweengroups.0081152811.0081152810.000.9621

Withingroups131.288523373.54833847

------------------------------------------------------------------------

Total131.296639383.4551747

Bartlett'stestforequalvariances:

chi2

(1)=0.5952Prob>chi2=0.440

两种检验方法最终的P=0.9621〉0.05所以按照α=0.05不拒绝H0,可以认为两种药物治疗慢性肾炎的疗效均数相同

二、某医院收集了一批慢性肾炎患者,随机分成3组,对3组慢性肾炎病人分别用A药(group=1),B药(group=2)和C药(group=3),评价指标为连续性变量,用Y表示(Y值越高表明疗效越好),结果见Stata文件(ex2X.dta,X=1,2,…,100)。

问这3种药物治疗慢性肾炎的有效率是否相同?

依题意,本题属完全随机设计,所得指标数据为计量资料,且相互独立,在满足各自应具备的条件的前提下,可以对其进行完全随机设计资料的方差分析以推断三个样本所代表的总体均数μ1,μ2μ3是否相等。

首先,分别对各组计量资料进行正态性检验:

(文件ex24)

H0:

资料正态分布

H1:

资料非正态分布

α=0.05

.sktestyifgroup==1

Skewness/KurtosistestsforNormality

-------joint------

Variable|Pr(Skewness)Pr(Kurtosis)adjchi2

(2)Prob>chi2

-------------+-------------------------------------------------------

y|0.9170.7930.080.9609

.sktestyifgroup==2

Skewness/KurtosistestsforNormality

-------joint------

Variable|Pr(Skewness)Pr(Kurtosis)adjchi2

(2)Prob>chi2

-------------+-------------------------------------------------------

y|0.5110.3071.640.4411

.sktestyifgroup==3

Skewness/KurtosistestsforNormality

-------joint------

Variable|Pr(Skewness)Pr(Kurtosis)adjchi2

(2)Prob>chi2

-------------+-------------------------------------------------------

y|0.6350.6920.390.8225

三组资料正态性检验的P值都大于0.05,不拒绝H0,认为三组资料分别服从正态分布。

然后对三组数据做Bartkett方差齐性检验:

H0:

σ12=σ22=σ32三个组总体方差相同

H1:

三个组总体方差不全相同

α=0.10

经计算检验统计量χ2=3.6108,P=0.164〉0.10在0.10水准处,不能拒绝H0,故不能认为三组数据的方差不齐。

建立假设进行方差分析:

H0:

3种药物治疗慢性肾炎的疗效均数相同

H1:

3种药物治疗慢性肾炎的疗效均数不全相同

α=0.05

.onewayygroup

AnalysisofVariance

SourceSSdfMSFProb>F

------------------------------------------------------------------------

Betweengroups58.9638576229.48192883.250.0450

Withingroups598.58731669.0695047

------------------------------------------------------------------------

Total657.551168689.66987012

Bartlett'stestforequalvariances:

chi2

(2)=3.6108Prob>chi2=0.164

方差分析的结果P=0.045〈0.05按照α=0.05拒绝H0,可以认为三种药物治疗慢性肾炎的疗效不全相同。

比较三种药物处理后治疗慢性肾炎的疗效的均数两两之间的差异:

H0:

第i组与第j组的疗效总体均数相同(i≠j)

H1:

第i组与第j组的疗效总体均数不同

α=0.05

采用LSD-t检验:

探索性研究,灵敏度高,两两之间的差异较灵敏的被发现

.lsdygroup

LSD方法进行两两均数比较

输入效应变量分组变量

mean1-mean2=0.421402P-value=0.631813

mean1-mean3=-1.828885P-value=0.040535

mean2-mean3=-2.250286P-value=0.021088

两两比较结果说明,A药组和B药组的疗效的差异没有统计学意义。

A药组和B药组的疗效都比C药组疗效差,差异具有统计学意义。

采用比较保守的灵敏度比较低的Bonferroni检验,差异较难被发现:

Comparisonofybygroup

(Bonferroni)

RowMean-|

ColMean|12

---------+----------------------

2|-.421402

|1.000

|

3|1.828882.25029

|0.1220.063

Bonferroni检验两两比较的结果是三个药组两两比较互相之间的差异不具有统计学意义。

可见此方法检验结果较保守,灵敏度不高。

三、某医院收集了一批高血脂患者,为了控制年龄与性别对疗效的影响,进行年龄和性别配对,同一对进行随机分为2组,A组用A药,B组用B药,用某个高血脂呈正相关的的指标描述治疗高血脂的疗效,疗效观察值越低表示疗效越好。

用Y1表示A组的疗效观察值,用Y2表示B组的疗效观察值,同一配对为一行。

结果见Stata文件(ex3X.dta,X=1,2,…,100)。

问这2种药物治疗高血脂的疗效是否相同?

本题属异体配对设计,所得数据为配对计量资料,在满足应具备的条件的前提下可以采用配对t检验进行处理:

首先计算配对差值并进行正态性检验(文件ex34)

H0:

配对差值服从正态分布

H1:

配对差值非正态分布

α=0.05

gend=x1-x2

sktestd

Skewness/KurtosistestsforNormality

-------joint------

Variable|Pr(Skewness)Pr(Kurtosis)adjchi2

(2)Prob>chi2

-------------+-------------------------------------------------------

d|0.5690.8810.350.8380

P=0.8380〉0.05,故认为配对的差值近似服从正态分布。

对配对设计资料进行t检验:

H0:

2种药物治疗高血脂的疗效差异的均数为0

H1:

2种药物治疗高血脂的疗效差异的均数不为0

α=0.05

.ttestx1=x2

Pairedttest

------------------------------------------------------------------------------

Variable|ObsMeanStd.Err.Std.Dev.[95%Conf.Interval]

---------+--------------------------------------------------------------------

x1|23.9293222.53317832.557034-.1764222.035066

x2|23.885844.50607642.427057-.16369441.935382

---------+--------------------------------------------------------------------

diff|23.0434782.89839174.308535-1.8196721.906629

------------------------------------------------------------------------------

Ho:

mean(x1-x2)=mean(diff)=0

Ha:

mean(diff)<0Ha:

mean(diff)~=0Ha:

mean(diff)>0

t=0.0484t=0.0484t=0.0484

P|t|=0.9618P>t=0.4809

P值为0.9618,在α=0.05水平不拒绝H0,认为两种药物治疗高血脂的疗效差别没有统计学意义。

同时也可以认为数据为随机区组设计资料,由上文可知,满足方差分析的条件。

源假设不变,采用随机区组设计资料的方差分析,数据格式做出改变,改为3个变量,x表示疗效的数值,变量b为两种不同的药物处理,变量a.为分组情况,变量值1、2…23分别表示23组患者

.anovaxab

Numberofobs=46R-squared=0.2533

RootMSE=3.04659AdjR-squared=-0.5274

Source|PartialSSdfMSFProb>F

-----------+----------------------------------------------------

Model|69.2620935233.011395370.320.9952

a|69.2403544223.147288830.340.9929

b|.0217390941.0217390940.000.9618

Residual|204.198249229.28173861

-----------+----------------------------------------------------

Total|273.460343456.07689651

区组处理效应的F检验P值为0.9618与t检验结果一致。

四、某研究者收集一批牛皮癣患者,随机分为2组,A组用A药(group=1),B组用B药(group=0)评价指标为有效(Y=1)和无效(Y=0),结果见Stata文件(ex4X.dta,X=1,2,…,100)。

问;两种药物的疗效如何?

本题属于成组设计两样本率比较的资料,首先对本题的数据资料进行整理(文件ex44)绘制四格表如下:

group|有效无效|Total

-----------+----------------------+----------

B药|1341|54

-----------+----------------------+----------

A药|1329|42

-----------+----------------------+----------

Total|2670|96

经验算,所有理论数T>5,采用成组设计两样本率比较的正态近似法μ检验或者Pearsonχ2检验。

首先建立假设:

H0:

2种药物治疗牛皮癣的疗效有效率相同。

H1:

2种药物治疗牛皮癣的疗效有效率不同。

α=0.05

μ检验的统计量:

Sp1-p2=0.09143,μ=-0.7525查μ界值表,得P>0.05,不拒绝H0,尚不能认为两种药物的疗效不同。

Pearsonχ2检验。

.tabi1341\1329,chir

|col

row|12|Total

-----------+----------------------+----------

1|1341|54

|24.0775.93|100.00

-----------+----------------------+----------

2|1329|42

|30.9569.05|100.00

-----------+----------------------+----------

Total|2670|96

|27.0872.92|100.00

Pearsonchi2

(1)=0.5660Pr=0.452

χ2=0.5660P=0.452>0.05不拒绝H0,尚不能认为两种药物的疗效不同。

五、某研究者收集一批老慢支患者,为了控制病程对疗效的影响,进行病程配对,同一对进行随机分为2组,A组用A药,评价指标为有效(Y1=1)和无效(Y1=0),B组用B药,评价指标为有效(Y2=1)和无效(Y2=0),同一配对为一行,结果见Stata文件(ex5X.dta,X=1,2,…,100)。

问;两种药物的疗效如何?

本题属于配对设计资料两样本率的比较,首先将文件数据整理成配对4格表(文件ex54)

|B组

A组|10|Total

-----------+----------------------+----------

1|368|44

0|3335|68

-----------+----------------------+----------

Total|6943|112

进行假设检验:

H0:

两种药物治疗老慢支的有效率无差别

H1:

两种药物治疗老慢支的有效率有差别

α=0.05

由于b+c>40采用McNemar配对检验

.mcci3683335

|Controls|

Cases|ExposedUnexposed|Total

-----------------+------------------------+----------

Exposed|368|44

Unexposed|3335|68

-----------------+------------------------+----------

Total|6943|112

McNemar'schi2

(1)=15.24Prob>chi2=0.0001

ExactMcNemarsignificanceprobability=0.0001

P值小于0.05,按照α=0.05水准,拒绝H0,认为两种药物治疗老慢支的有效率有差别。

六、某医院收集了一批消化道疾病患者,为了控制病程和性别疗效的影响,进行病程和性别配对,同一对进行随机分为3组,A组用A药,B组用B药,C组用C药,用某个定量的观察指标描述疗效,该指标观察值越低,表示疗效越好。

用Y1表示A组的疗效观察值,用Y2表示B组的疗效观察值,用Y3表示C组的疗效观察值,同一配对为一行。

结果见Stata文件(ex6X.dta,X=1,2,…,100)。

问这3种药物治疗该疾病疗效是否相同?

本题属于随机区组设计资料,所得数据为计量资料。

在满足所需条件的

展开阅读全文
相关资源
猜你喜欢
相关搜索
资源标签

当前位置:首页 > 自然科学 > 物理

copyright@ 2008-2023 冰点文库 网站版权所有

经营许可证编号:鄂ICP备19020893号-2