CHICOM: Code for comparing weighted or unweighted histograms in Fortran-77, C++, R and Python
| Autor: | Helmut Neukirchen, Helgi Halldorsson, Nikolay D. Gagunashvili |
|---|---|
| Rok vydání: | 2019 |
| Předmět: |
Discrete mathematics
Fortran General Physics and Astronomy Python (programming language) 01 natural sciences 010305 fluids & plasmas Goodness of fit Hardware and Architecture Histogram Conjugate gradient method 0103 physical sciences Test statistic 010306 general physics computer Statistic computer.programming_language Mathematics Statistical hypothesis testing |
| Zdroj: | Computer Physics Communications. 245:106872 |
| ISSN: | 0010-4655 |
| Popis: | Improved a program that calculates test statistics to compare weighted and unweighted histograms. The program is presented in Fortran-77, C++, R and Python. The code calculates test statistics for histograms with either normalized or unnormalized weights of events. New version program summary Program Title: CHICOM Program Files doi: http://dx.doi.org/10.17632/424sd4fhj8.1 Licensing provisions: GPLv3 Programming language: Fortran-77, C++, Python, R Journal reference of previous version: CHICOM: A code of tests for comparing unweighted and weighted histograms and two weighted histograms, N. D. Gagunashvili, Comput. Phys. Commun. 183 (2012) 193-196 Does the new version supersede the previous version?: Yes Reasons for the new version: To use an improved version of the chi-square test with better statistical properties, instead of the median statistic [3] as in the previous version. Summary of revisions: An improved version of the test statistic [2] was used that uses an improved chi-square test [4]. The algorithm has been implemented in four commonly used programming languages (Fortran-77, C++, Python and R). Nature of problem: The program calculates test statistics for comparing weighted or unweighted histograms. Solution method: An improved test statistic for comparing weighted histograms is calculated using the formulas presented in Ref. [2]. In order to find the test statistic, we must find the probability that minimizes the sum of the goodness of fit test statistic of each histogram. To do so, the Polak–Ribiere conjugate gradient method is used to converge to the minimum from an initial guess suggested in the article. References [1] N. D. Gagunashvili, Comput. Phys. Commun. 183(2012)193. [2] N. D. Gagunashvili, Eur. Phys. J. Plus (2017) 132: 196. [3] N. D. Gagunashvili, Nucl. Instrum. Meth. A 596(2008)439. [4] N. D. Gagunashvili, J. Instrum. 10(2015)P05004. |
| Databáze: | OpenAIRE |
| Externí odkaz: |
|
Full Text from ScienceDirect