求助，如何确定数据的分布？谢谢。

hitman

想问下一，SPSS或者EXCEL能确定一组数据是属于什么分布的吗？比如有最近10年的每年交通事故量，能看出是属于什么分布的吗？谢谢。

hitman

多谢版主，每次都能学到新东西~~~~~~~~~~~~

$ }; ^- X5 s' W# m; G开始研究下。

olh2008

数据分布的估计SPSS有这样的功能，而Excel这方面功能不怎么强。
$ `9 c& k4 s  C0 e8 \其实你应该多看看SPSS的相关书籍的。因为不同的分布判定不一样。
4 O% ~# G+ v" Y% r6 [5 @这里给你一个能直接用于判定的Matlab函数，好好研究一下吧。
& `8 U" ^7 {" l2 k9 N8 R  R& n! M

1. - n! }0 @$ B* a7 m: l' s1 d
   
2. function f=p_judge(A,alpha)7 H$ s2 S\" j2 _\" K/ k
   
3. %判定数据的概率分布形式9 x$ ?7 n( x9 Y) z\" |9 Q' V( }; R
   
4. A=A(:);
   $ j: B7 g* K7 U& i% G
5. [mu,sigma]=normfit(A);. G( n) |6 h' K5 y
   
6. p1=normcdf(A,mu,sigma);+ X( c9 S2 [% J& v
   
7. [H1,s1]=kstest(A,[A,p1],alpha);0 M# u\" J/ Z/ {, X9 K* D9 w4 ~
   
8. n=length(A);3 `: C6 ^: A  W% [, A0 B
   
9. if H1==0- C4 C( R. a9 \6 M% Y1 f
   
10.     disp('该数据服从正态分布。')% y2 Y7 M0 c& Q$ J3 J' h3 a, X
    
11. end
    # v2 [3 [9 I; K4 q+ d# m& U2 T
12. phat=gamfit(A,alpha);
    $ p5 ~* t- ]% m# d+ k# [
13. p2=gamcdf(A,phat(1),phat(2));
    & c! t* Z* B# f- w
14. [H2,s2]=kstest(A,[A,p2],alpha);
    4 w4 X2 I3 J, h
15. if H2==0  x# K3 D( g0 f6 m4 W9 r' x( |
    
16.     disp('该数据服从γ分布。')
    ) S: n\" N1 Z$ Z$ ~  k1 g
17. end
    4 L, a# k; N( u4 ^8 j7 }3 y# _
18. lamda=poissfit(A,alpha);* h; C/ q3 a* l
    
19. p3=poisscdf(A,lamda);' `# `7 \# L- @8 ^: I& C- |$ p
    
20. [H3,s3]=kstest(A,[A,p3],alpha);
    4 D! S% q1 L- e6 ~' V1 n* i$ A, t$ K
21. if H3==01 @( i3 ?! n+ s6 }4 S
    
22.     disp('该数据服从泊松分布。')  / X9 |6 }- a9 N3 ^9 D
    
23. end\" Y6 X; u# s' m9 H  o
    
24. mu=expfit(A,alpha);% P6 `9 {; Z( V/ j
    
25. p4=expcdf(A,mu);
    4 L- j8 @! u( k* j4 K3 @& L! f
26. [H4,s4]=kstest(A,[A,p4],alpha);
    - ~3 M4 t9 P& E+ n, r; z4 O
27. if H4==0
    4 A9 v/ F9 M+ `1 i+ B
28.     disp('该数据服从指数分布。')  7 D$ T) j6 z$ b7 z+ [
    
29. end
    ( B. u1 [  ?) \# U5 |
30. [phat, pci] = raylfit(A, alpha);
    5 K; b. K; g7 }  k- H
31. p5=raylcdf(A,phat);% H1 I/ C; @, v# E$ W% s; X
    
32. [H5,s5]=kstest(A,[A,p5],alpha);1 Y9 O7 a2 b\" e  a9 f1 q$ }9 ~
    
33. if H5==0
    2 @( a- W3 I/ d+ Y5 C
34.     disp('该数据服从rayleigh分布。')  
    / g. r& m3 d7 w. e; B
35. end' l+ B) z' ^+ |1 w4 h
    

复制代码