《非参数统计》与matlab编程第四章符号秩和检验法

《非参数统计》与matlab编程第四章符号秩和检验法

ID:32324312

大小:48.00 KB

页数:10页

时间:2019-02-03

《非参数统计》与matlab编程第四章符号秩和检验法_第1页
《非参数统计》与matlab编程第四章符号秩和检验法_第2页
《非参数统计》与matlab编程第四章符号秩和检验法_第3页
《非参数统计》与matlab编程第四章符号秩和检验法_第4页
《非参数统计》与matlab编程第四章符号秩和检验法_第5页
资源描述:

《《非参数统计》与matlab编程第四章符号秩和检验法》由会员上传分享,免费在线阅读,更多相关内容在应用文档-天天文库

1、第四章符号和检验法函数signrank格式p=signrank(x)原假设为x的中位数为0,显著性水平为0.05的双侧检验。p=signrank(x,m)原假设为x的中位数为m,显著性水平为0.05的双侧检验。p=signrank(x,m,alpha)原假设为x的中位数为m,显著性水平为alpha的双侧检验。[p,h]=signrank(...,'alpha',alpha)例:[p,h]=signrank(...,'alpha',0.01)[p,h,stats]=signrank(...,'method',‘exact’)用精

2、确的方法[p,h]=signrank(...,'method',‘approximate’)用正态近似的方法[p,h,stats]=signrank(x,y,'alpha',0.01,'method','exact')[p,h,stats]=signrank(y1,y2,0.01,'method','approximate')所P值除以2,得到相应单侧检验的P值。§4.2x=[20.323.52219.12124.716.118.521.924.223.425];y=[1821.722.51721.224.817.214.92

3、021.122.723.7];[p,h,stats]=signtest(x,y)p=0.3877h=0stats=sign:4length(find((x-y)>0))ans=82*(1-binocdf(7,12,0.5))ans=0.3877p=0.3877与书上算的不一样,书上算错了。符号检验接受原假设。W+个数取的值01不取任何数111212323、(1,2)424、(1,3)535、(1+4)、(2+3))646、(1,5)、(2,4)、(1,2,3)757、(1,6)、(2,5)、(3,4)、(1,2,4)868、(

4、1,7)、(2,6)、(3,5)、(1,2,5)、(1,3,4)989、(1,8)、(2,7)、(3,6)、(4,5)、(1,2,6)、(1,3,5)、(2,3,4)101010、(1,9)、(2,8)、(3,7)、(4,6)、(1,2,7)、(1,3,6)、(2,3,5)、(4,5,1)、(1、2、3、4)111211、(1,10)、(2,9)、(3,8)、(4,7)、(5,6)、(1,2,8)、(1,3,7)、(1,4,6)、(2,3,6)、(2,4,5)、(1,2,3,5)符号秩和检验:[p,h,stats]=signr

5、ank(x,y)p=0.02685546875000h=1stats=signedrank:11P43,表4.5中:3+2+1+5=11,12*13/2-11ans=67a=11122345681012p=2*sum(a)/2^12p=0.02685546875000在显著性水平0.05下,拒绝原假设。符号秩和检验应用的条件:假设总体服从对称分布,而符号检验不需要。习题四1.x1=[22.3225.7624.2321.3523.4326.9718.3620.7524.0726.4325.4127.22];x2=[21.2523

6、.9724.7719.2623.1226.0019.4017.1822.2323.3524.9825.90]符号秩和[p,h,stats]=signrank(x1,x2)p=0.01220703125000h=1stats=signedrank:8在0.05显著性下,拒绝原假设,有差异。2.符号检验法:[p,h,stats]=signtest(x1,x2)p=0.03857421875000h=1stats=sign:2[h,p,ci,tstat]=ttest(x1,x2)h=1p=0.00888227075133ci=0.3

7、79924103111732.10174256355494tstat=tstat:3.17229233575613df:11sd:1.354973722209693.x1=[390390450380400390350400370430]x2=[270280350300300340290320280320]x3=x2+100[phstats]=signrank(x1,x3)p=0.26562500000000h=0stats=signedrank:9或者:[p,h,stats]=signrank(x1-x2,100)p=0.26

8、562500000000h=0stats=signedrank:9接受原假设,差价是100L-H估计:a=x1-x2a=1201101008010050608090110b=sort([mean(nchoosek(a,2),2);a']);中位数:median(b)ans=9

当前文档最多预览五页,下载文档查看全文

此文档下载收益归作者所有

当前文档最多预览五页,下载文档查看全文
温馨提示:
1. 部分包含数学公式或PPT动画的文件,查看预览时可能会显示错乱或异常,文件下载后无此问题,请放心下载。
2. 本文档由用户上传,版权归属用户,天天文库负责整理代发布。如果您对本文档版权有争议请及时联系客服。
3. 下载前请仔细阅读文档内容,确认文档内容符合您的需求后进行下载,若出现内容与标题不符可向本站投诉处理。
4. 下载文档时可能由于网络波动等原因无法下载或下载错误,付费完成后未能成功下载的用户请联系客服处理。