[Scipy-tickets] [SciPy] #1425: kurtosistest requires at least 4 observations

SciPy Trac scipy-tickets@scipy....
Sun Sep 11 22:26:41 CDT 2011


#1425: kurtosistest requires at least 4 observations
-------------------------+--------------------------------------------------
 Reporter:  bsouthey     |       Owner:  somebody   
     Type:  defect       |      Status:  needs_work 
 Priority:  normal       |   Milestone:  Unscheduled
Component:  scipy.stats  |     Version:  devel      
 Keywords:               |  
-------------------------+--------------------------------------------------

Comment(by josefpktd):

 just to make sure

 {{{
 >>> data = np.array([ 393.,  353.,  334.,  336.,  327.,  300.,  300.,
 308.,  283.,
         285.,  270.,  270.,  272.,  278.,  278.,  263.,  264.,  267.,
         267.,  267.,  268.,  254.,  254.,  254.,  256.,  256.,  258.,
         240.,  243.,  246.,  247.,  248.,  230.,  230.,  230.,  230.,
         231.,  232.,  232.,  232.,  234.,  234.,  236.,  236.,  238.,
         220.,  225.,  225.,  226.,  210.,  211.,  212.,  215.,  216.,
         217.,  218.,  200.,  202.,  192.,  198.,  184.,  167.])
 >>> np.testing.assert_equal([np.round(z, 2), np.round(p, 4)], (2.21,
 0.0269))
 >>> #z_aba, p_aba = (2.21, 0.0269)
 }}}

 example and results from
 Ralph B. D'Agostino, Albert Belanger and Ralph B. D'Agostino, Jr., A
 Suggestion for Using Powerful and Informative Tests of Normality,
 The American Statistician
 Vol. 44, No. 4 (Nov., 1990), pp. 316-321
 Published by: American Statistical Association
 Article Stable URL: http://www.jstor.org/stable/2684359

 Results with R moments z = 2.2226, p-value = 0.02624 (Package moments
 version 0.12) are different than the published results

-- 
Ticket URL: <http://projects.scipy.org/scipy/ticket/1425#comment:5>
SciPy <http://www.scipy.org>
SciPy is open-source software for mathematics, science, and engineering.


More information about the Scipy-tickets mailing list