- 1、本文档共9页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
数学实验报告利用MALTAB进行统计推断
PAGE \* MERGEFORMAT 10
实验九 统计推断
一、汽油
㈠问题描述
据说某地汽油的价格是115美分/gal,为了验证这种说法,一位司机开车随机选择了一些加油站,得到某年1月和2月的数据如下:
(1)分别用两个月的数据验证这种说法的可靠性;
(2)分别给出1月和2月汽油价格的置信区间;
(3)如何给出1月和2月汽油价格差的置信区间
㈡简要分析
本题要求在总体方差未知的情况下对总体均值进行假设检验,并求置信区间。
1、假设检验
H0:u=u0,H1:u!=u0.
由于样本方差未知,用样本方差代替总体方差,由
知:
当Px-usn≤t1-a2=1-a时,H0成立。使用ttest求解。
2、置信区间
需要根据具体问题具体分析,具体内容见下一节。
㈢结果与分析
1、假设检验
N = 20
计算得到标准正态分布0.975的分位数为1.96。
由两者均值:
Mean1 = 115.1500000000000
Mean2 = 120.7500000000000
方差:
Var1 = 14.976315789473684
Var2 = 13.776315789473685
计算得到1月份接受假设的区间为
[113.3039619640030 116.6960380359971]
可以看出Mean1在此范围中,故接受假设;
同理可得2月份接受假设的区间为
[113.3733293233395 116.6266706766605]
可以看出Mean2不在此范围内,故不接受假设。
下面使用MATLAB函数直接解决此问题:
[h01,sig01,ci01]=ttest(data1,115);
[h02,sig02,ci02]=ttest(data2,115);
计算得到h01 = 0,h02 = 1,与前面得到的结果一致。
2、置信区间
沿用第一的计算方法,把正态分布的分位数改为t分布的分位数,计算得到:一月份平均值的置信区间为[113.3449294088214 116.9550705911786];
二月份平均值的置信区间为[119.0187561613285 122.4812438386715];
3、价格差的置信区间
此时一月份与二月份的总体方差都未知,有两种处理办法
(1)设两月的总体方差相同,于是由统计部分知识得到:
x1-x2 ~ N(u1-u2,1m+1ns2)
mn(m+n-2)m+nx-y-(u1-u2)m-1sx2+(n-1)sy2 ~ t(m+n-2)
最终计算得到u1-u2的置信区间为
[-8.027274417173238 -3.172725582826751]
(2)设两月的总体方差不同,于是由统计部分知识得:
s02=sx2m+sy2n
T=x-y-u1-u2/s0
T近似服从自由度为l的t分布。
l=s04sx4m2(m-1)+sy4n2(n-1)
于是计算得到u1-u2的置信区间为
[-8.928993085768209 -2.271006914231780]
(3)对比
可以发现后者计算得到的区???比前者要大,这是假设带来的影响。
㈣程序清单
clear;clc;
data1 = [119 117 115 116 112 121 115 122 116 118 109 112 119 112 117 113 114 109 109 118];
data2 = [118 119 115 122 118 121 120 122 128 116 120 123 121 119 117 119 128 126 118 125];
mean0 = 115;
norm975 = norminv(0.975,0,1);
t975 = tinv(0.975,20);
t97540 = tinv(0.975,38);
mean1 = mean(data1);
mean2 = mean(data2);
var1 = var(data1);
var2 = var(data2);
[h01,sig01,ci01]=ttest(data1,115);
[h02,sig02,ci02]=ttest(data2,115);
x11 = mean0+norm975*sqrt(var1)/sqrt(20);
x12 = mean0-norm975*sqrt(var1)/sqrt(20);
x21 = mean0+norm975*sqrt(var2)/sqrt(20);
x22 = mean0-norm975*sqrt(var2)/sqrt(20);
y11 = mean1+t975*sqrt(var1/20);
y12 = mean1-t975*sqrt
您可能关注的文档
- 数字图像处理复习基本内容2012年度最终版.doc
- 数字化城市管理系统建设方案25页,9个标准子系统.doc
- 数字图像处理实验4图象处理中的正交变换.doc
- 数字图像处理技术在大蒜切脐根机中的应用研究.doc
- 数字图像处理技术应用除草机器人.doc
- 数字图像处理教学大纲2014新版.doc
- 数字图像处理实验三数字图像分割处理的工程实现.doc
- 数字图像处理抖动加网伽马校正直方图均匀化.docx
- 数字图像处理均值滤波.doc
- 数字图像处理的MATLAB实现新实验指导书.doc
- 上课《归园田居(其一)》课件28张.pptx
- 色彩的魅力课件55.ppt
- 山水田园诗导学练习鉴赏.ppt
- 部编版《白鹭》课件.pptx
- 《Excel在财务中的应用》项目七 全面预算编制方法与应用 教学课件.pptx
- 《Excel在财务中的应用》项目六 利润管理决策模型与应用 教学课件.pptx
- 《Excel在财务中的应用》项目五 运营资金的管理 教学课件.pptx
- 《Excel在财务中的应用》项目四 项目投资管理决策与应用 教学课件.pptx
- 《Excel在财务中的应用》项目三 筹资管理决策模型与应用 教学课件.pptx
- 《Excel在财务中的应用》项目二 Excel 常用函数基础 教学课件.pptx
文档评论(0)