声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 5442|回复: 5

[综合讨论] 如何统计出某列数据符合什么分布

[复制链接]
发表于 2007-8-26 09:33 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
问题1 : 给定一列数据,如何统计出这列数据符合什么分布呢?
如:
0 0.41214 0.76525 0.68232 0.73807 1.019 0.60452 1.0834 1.1135 0.94951 0.69712 1.0637 1.3658 0.8992 1.715 1.604 2.5476 2.3084 2.3531 2.6823 2.8665 2.8017 2.6988 3.1254 2.9164 2.9204 3.1865 3.1029 3.1999 3.2573 3.6025 3.8058 4.2056 4.5434 5.7013 5.2751 6.3063 6.8964 7.0063 7.1357 6.9776 6.6779 6.5496 6.5395 7.257 7.5914 7.4754 7.8587 8.1931 8.3006 8.2573 8.0932 8.6925 9.2447 9.61 9.6346 9.959 10.524 11.109 11.381 11.447 12.109 12.437 12.933 13.222 13.256 13.347 13.432 13.539 13.431 13.557 13.722 13.804 13.337 13.924 13.628 14.11 14.188 14.151 14.239 13.839 14.313 14.444 14.276 15.167 16.276 17.003 17.519 17.242

先谢了!!!!
回复
分享到:

使用道具 举报

发表于 2007-8-26 10:21 | 显示全部楼层
先画图,再估计.
 楼主| 发表于 2007-8-26 11:55 | 显示全部楼层
大侠能详细点吗?如果只是plot(x)肯定不行


===================================
可以参考下面bainhome的程序;
也可以先画图、初步判断,再用相应的Matlab函数判断.
By xjzuo
==================================

[ 本帖最后由 xjzuo 于 2007-8-26 22:06 编辑 ]
发表于 2007-8-26 12:14 | 显示全部楼层
你的数据服从正态分布。
以前写过一个假设检验的函数,主要是利用K-S检验的一个MATLAB函数kstest。其他分布自己按照格式照猫画虎填写即可:
  1. function f=p_judge(A,alpha)
  2. % 本程序用于判别所给数据源在置信率为0.05时的概率分布形式。A的形式为n×1。
  3. A=A(:);
  4. [mu,sigma]=normfit(A);
  5. p1=normcdf(A,mu,sigma);
  6. [H1,s1]=kstest(A,[A,p1],alpha)
  7. n=length(A);
  8. if H1==0
  9.     disp('该数据源服从正态分布。')
  10. else
  11.     disp('该数据源不服从正态分布。')
  12. end
  13. phat=gamfit(A,alpha);
  14. p2=gamcdf(A,phat(1),phat(2));
  15. [H2,s2]=kstest(A,[A,p2],alpha)
  16. if H2==0
  17.     disp('该数据源服从γ分布。')
  18. else
  19.     disp('该数据源不服从γ分布。')
  20. end
  21. lamda=poissfit(A,alpha);
  22. p3=poisscdf(A,lamda);
  23. [H3,s3]=kstest(A,[A,p3],alpha)
  24. if H3==0
  25.     disp('该数据源服从泊松分布。')  
  26. else
  27.     disp('该数据源不服从泊松分布。')  
  28. end
  29. mu=expfit(A,alpha);
  30. p4=expcdf(A,mu);
  31. [H4,s4]=kstest(A,[A,p4],alpha)
  32. if H4==0
  33.     disp('该数据源服从指数分布。')  
  34. else
  35.     disp('该数据源不服从指数分布。')  
  36. end
  37. [phat, pci] = raylfit(A, alpha)
  38. p5=raylcdf(A,phat);
  39. [H5,s5]=kstest(A,[A,p5],alpha)
  40. if H5==0
  41.     disp('该数据源服从rayleigh分布。')  
  42. else
  43.     disp('该数据源不服从rayleigh分布。')
  44. end
复制代码

评分

1

查看全部评分

 楼主| 发表于 2007-8-26 14:44 | 显示全部楼层
谢谢bainhome , 我先学习一下你的程序,如果不明白还有多多请教
发表于 2009-12-17 16:42 | 显示全部楼层
呱呱,太感谢了,正是我想要的!!找了好久了!!:lol :lol :lol :lol
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-11 15:39 , Processed in 0.058228 second(s), 20 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表