matlab 单侧检验,假设检验(单侧).ppt

假设检验(单侧).ppt

总体废品率的 假设检验(双侧) H0成立时 取N(0,1)的1-?/2分位数u1-?/2 假设检验(单侧) 问题 H0成立，乙方应接受那批产品. 问：如果将?提高，会有什么结果？ 总体分布的正态性检验 正态分布的偏度g1=0，峰度g2=3，对于一个样本计算其g1和g2，若样本来自正态总体，则g1和g2应分别在 0 和 3 附近。基于这个思想，构造一个包含g1，g2的 统计量。 2)Kolmogorov-Smirnov检验 通过样本的经验分布函数与给定分布函数的比较，推断该样本是否来自给定分布函数的总体。该正态性检验只能做标准正态检验 。 3)Lilliefors检验 它将Kolmogorov-Smirnov检验改进用于一般的正态性检验。 1)Jarque-Bera检验 假设检验的MATLAB实现 MATLAB命令使用说明 输出参数h=0表示接受H0, h=1表示拒绝H0, sig标示对假设的接受和拒绝程度。 ci给出的置信区间，zval是样本统计量z的值。 用N(5,1)随机数产生n=100的样本，在总体方差未知的情况下分别取?=0.05和?=0.01检验总体均值? ?5.2。 x = normrnd(5,1,100,1); m = mean(x), [h1,sig1,ci1] = ttest(x,5.2,0.05,-1) [h2,sig2,ci2] = ttest(x,5.2,0.01,-1) m = 5.0111, [h1,sig1,ci1] = 1 0.0343 -Inf 5.1815 [h2,sig2,ci2] = 0 0.0343 -Inf 5.2537 可知在?=0.05下拒绝H0(此时sig1)，? 的区间估计(-? 5.1815]不包含5.2；而在?=0.01下接受H0(此时sig2>?)，? 的区间估计(-? 5.2537]包含5.2。 计算结果 demo1202.m 四. 示例求解与分析 示例：吸烟对血压的影响 测量了吸烟组(66人)和不吸烟组(62人)两组人群的6项血压指标 吸烟组均值 吸烟组标准差 不吸烟组均值 不吸烟组标准差 24hSBP(mmHg) 119.35 10.77 114.79 8.28 24hDBP(mmHg) 76.83 8.45 72.87 6.20 　dSBP(mmHg) 122.70 11.36 117.60 8.71 　dDBP(mmHg) 79.52 8.75 75.44 6.80 nSBP(mmHg) 109.95 10.78 107.10 10.11 nDBP(mmHg) 69.35 8.60 65.84 7.03 两个总体均值的假设检验 function [h,sig]=pttest2(xbar,ybar,s1,s2,m,n,alpha,tail) spower=((m-1)*s1^2+(n-1)*s2^2)/(m+n-2); t=(xbar-ybar)/sqrt(spower/m+ spower/n); if tail==0 a=tinv(1-alpha/2,m+n-2); sig =2*(1-tcdf(abs(t),m+n-2)); if abs(t)<=a h=0; else h=1; end end if tail==1 a=tinv(1-alpha,m+n-2); sig =1-tcdf(t,m+n-2); if t<=a h=0; else h=1; end end if tail==-1 a=tinv(alpha,m+n-2); sig =tcdf(t,m+n-2); if t>=a h=0; else h=1; end end h Sig 接受或拒绝H0 24hSBP(mmHg) 1 8.5097e-003 拒绝 24hDBP(mmHg) 1 3.1860e-003 拒绝 　dSBP(mmHg) 1 5.3064e-003 拒绝 　dDBP(mmHg) 1 3.9950e-003 拒绝 nSBP(mmHg) 0 1.2597e-001 接受 nDBP(mmHg) 1 1.3027e-002 拒绝 计算结果(?=0.05) 除夜间(10Pm-6Am)平均收缩压(nSBP)外，其余5项指标都拒绝了H0，于是综合起来可以认为，吸烟对血压的影响显著。 Demo1203.m 加油站(乙方)以含硫量不超过0.08%的标准决定是否接受炼油厂(甲方)