![](data:image/jpeg;base64,iVBORw0KGgoAAAANSUhEUgAAAMgAAACFCAYAAAAenrcsAAAABmJLR0QAOABGAP8xB1bHAAAACXBIWXMAAASaAAAEmgEucM3eAAAAB3RJTUUH3wwYExYagHvHrwAAC/VJREFUeNrtnXt0VMUdx78XA74QEAQKUcQHFU9TAV+gqBU59kBAFOUlUK1VNLagUku1HsEDVWupVtESHoJaXoLyMIRilUdSjYGgBFFIsEC0BCRASLJ5EEjI/PrHZN29+5y92cQk+/2cs2fPnZ2ZO3P3fmfmzr3fO5ZSSlCLZVkwQUTccQWAqv2u8Ypi1YZZAdJateGW1++WT1obpuUiJNrEOUnkdcJaAFrUiqFFAFH4CSXAye4dR/iXkMaEJSLROiklVM9BSFOkRTTFFmyIRAgFQggFQggFQgihQAihQAihQAihQAihQAihQAihQAihQAihQAihQAghFAghFAghFAghFAghFAghFAghFAghFAghFAghhAIhhAIhhAIhhAIhpKGJay4VSUwKH2f9XB1v/Vyz/EzixRqxdlyajUB8/zSe4IRDrAYUHOEQK2aoqRGs+DfwUQZQVAp0ag/ceRswbIAVsAc6USl4dz2QuQMoLAHizgCuuESn6XtV6JfZf/WNYOFq4NuDwPltgHuHAINu8qR5cb6gcwfgwXvs+SxcJTh8DHg2yRNeWKzLvWkL0KY1cPM1wPg7gDNbWT+UefUswbz3gIxsoIUF/LI/8Ju7gapqYP77wKfbAaWA2/oCD48C4s6wp53/PvBZNnCqGkjoATx0D9A93nJcv0aEo2U5YlIgs9/VJ9iLk4H2bYC9/wNeXQS0bS34xXX+x/CVd4DOHWrjtwUqTwK5ecDaNKDvVcH3k18geGE+kDQauP7nQH4B8NcFQJeOgl5X6P38biww8Xng1usEl3XTYfsOCDZuAZKn2vN7ciYwbijwwHCgtBx4azWw7F96203yu0DCT7UoSiuAl98CVm0ADnyvT/j77wTKTujwdenAXQM9aWct1sIf/5ze/uQL4E+vAX//o6BLR8tR/X5kQUQhl2bK4EeCh/8zRfmFZ32l5PG/qIDph09SUn5CRVyGmQuVrPzIni5zh5Kpr9vDtnypZNILSk6f1p+JzyvZuCX8/krLlfz6GXuZ12y0p9u9V8ldk5Rf+K69Sn77Z3vapan++1ySqmTmwsDHxbR+DYQy+ERMTPYgt9/gH3blpcChI4Hjd+sCLFgJjBsquOB885YxNw8YO8QeltBD92De9OtlIXOH4INNutlr3xYY2M++n1NVgpUfA+nbgKNFQPVpvcSjZdkbyv597HlfHA+cqrLQv489Xvd4oOCYPe6Avv51GHA9MOXlutWvgaiXLismBdK5g3/YeedaqKgMHP+ph4A3VwIPTgU6dxD0vAS4oTfQr1foBUaPFQMTnrMC/I/+vf8Dw/VwpqISeGmyf15vrwEOFuiyxHcCzj5Li2b4Y/b8Lzjfnu7cs62g4ZWn7GHt2/rvt0M7wFVW9/rxIr0pTd21iKyx6dLRwrRHgapqQX4B8M23wPIPga07gcn3B093zllA8lRBh3ZW2MZu+YfAr4bpa4u1acCjY+y/Z2wHXn0a6Njek7agMEAzGkSwJisFF7mALh3tYcdLgLbn1b1+nOaNAVq1tHDZRRYSb7EwYxLwyfbQ8XtdoUUUjuwcPWPVv4+FQTcBefk6zDbEqgZatbSn25AZ3fqlbfMPS98G9O5Zt/pRIM2cKX8TbNoqOFYkOF0jKHIJ1mzU4+1QjB2iZ5k2ZApc5YLKk4IduYJpb3hO/rIKwT+W6Zkgd0s/cRzwxlL9m5trE4B57wFFLkFxqWDlx4Ijx6Nbz/wCIGWzoNilPymbBan/AcYOdV4/DrFigPHDgNQ0YP57wMkq4IJ2QP+rgaceDJ2ue7yFGRMF76QAc1YANTVAz0uBMYO9ppyX6Qvhrp08w5KLu1q49TrB7GXA0xN02KOjgbkrgKTpgBLgxt56ePfZjujV87FxWoRLUoHq00DC5cCLjwPxnSzH9WvqRHOddNKE4aM5HGIRQoEQwiEWIexBCKFACKFACKFAGhATey0h7EGaKU2xAWiujRYFQkgIGvWjJt8dEixYBezaC5zZErjpGmDCCMD3adGtOwXLPwS+O6S3u8cD9yb622HD2UOdWmsTk4DU2eFtvCZldbfE3i1yuDvcR44LFq8FsnO0W/CizsDowfjBHZmYBKx5XbBwlbbcusotW54mttlQdt9wZW5Ctlx/Gqsj8NARJWP+oJ1wx0v0Z9UGJS+9qWyutpz9SkZOVrJ5q5Jil5Iil3bjjXhCyZ48j4nswGElo36vZHOWdgfm5mk33pd7PHFmJCuZt0JJQaGSqmolrjIlW3cqeeY1Fda9OGuxkrfXKDl4RMmJSiU79+j807epiMsazA0ZiCOFSsZNUbIkVcn3R5WcqlKye5+S6cn2/F55x1N3b0yOi4jIfU8r+ShDSUWlksPHlLwwT8lbq8OX2TT/xkqjFcjMhSqgBXRRil0g05OVrE3zj5eyWcmf56iI7KFOrbWmNl7TskYikFfe1sIMV75VG1SdbMEmdt9o5t9YaLTXINm5gS2gA/vZt/cdCPzihH5X6ZcxuMnN0w5AbxJ6AHkHPdtua21hceQPF5jYeE3LGgnbcwLv25cbewcONzkup6oES9cJJkwT3DlRkJgEjH7SwtGi8Ps1yZ/XIA4oLQ9uAfWmpBRoF8Dx1q4NUFzq2Taxhzq11gJmNl7TskaCq8zfThsI3+MWyXExtfs6zZ8CcUCb1sEtoL4nV0mZvii2CadMXxC6MbGHOrXWAmY2XtOyRnqcCouB+M6h47WMC1w+k+Niavd1mj+neR1w9ZWBLaCbs+zbl3cDsr7yj5e1U//mJhJ7aKTWWlNMy9oyTlBdbdbCXvszYMMW52UyOS4mdt9gZW7qttxGK5CxQ4HUdLsF9INNgkNH7fFG3A4sSgHSsgQlpfqTliVYnAqMHOSVn4E91Km11hTTsnbpCHy+S78BMhzj7wA2ZgJL1wkKCgVV1YLcPMGMOWYCMzkuJnbfYGVu6rbcRv24e95BPXe/e59uwfr3AR4ZBdz9uH0eP3OHnqP3vrcwJhG4oZe9G99/QNtDd++z20P7XKnj7fxGkJoGfP1fu7V21CCg9Tmh74MEu1fh+5tJWb/YJZizQg9jRKyw90EKCgWLUvTERsUJoFtXYPQg4JZrrbDlMzkupeWCuSuAL3Z77L5Jo4ERT3jKFqrM4fKnQAjhEIsQCoQQCoQQQoEQQoEQQoEQQoH8ONDeS4EQQoLAG4WEhIBvdzcYYrkfm6jLKrJA+FVqAWB/vn68Jmd/rc34amDCSP/Ha+rLPkwokDpRl1Vkw61Se7BA8Ows/YDflAd02KfZwKwl9jKYrC7rdGVe4oMQ49Vy67KKrIlt9aUFSpav90+zdJ2KeHVZp/ZhwlVu64TTVWRNVqn9cg9w3zD/fQ64Xi9q48ZkdVmnK/OSHzoODrGc4HQVWRPbajCbsW9YfduHm8sJ7qSevnNWFEik034OV5E1sa0GsxkXuezb9W0fbg4tf7TgfRDDg17XA29iW+3dU8+C+ZL+uX07Ehtryzjg0guBwTcD0ycKPtnuqU806tVUGrBI83B/WvCEj+ygu9MF+w6GiW113BBgzSZgXbqOU1wqSE0X5B+2l93bxlpSplB5UpCdozDtDU+93PbhwmKgRum3pnywyUJCD/sJ0ByHW3Wpl2+6mBxiOTl47jThvoNhskrthT+x8PwkfX/jzZXAma30fZCHRwIZ2VJ7UW+hezy8Vpe1am2sFsYM9pTD6cq8xOd/5530xk92jr55OHsqZ6MaGl6DNDKmJwu+3iuoqBS4ygQZ2YLXFtlvOBL2IDFLRrZg9Ub91hOlgO5dgTsGAAP7sfegQAjhEIsQCoQQCoQQCoQQCoQQElAgnMQihD0IIc4EEgveAELYgxBCgRBCgRBCgRBCgRBCgRBCgRBCgRBCKBBCKBBCKBBC6hm+ejQGcD+xHc3n7ryfAm/Oz/NRIDFAtE9gpy+G5hCLUHAUCCG8SCeEUCCEUCCEUCCEUCCEUCCEUCCENBpCvRuOAiExT6gbn3zUhDT5Fr8+7+yzByHNZjhUH2m5whQhIYhryO6KkCYnEAqCEF6DEEKBEEKBEEKBEEKBEEKBEEKBEEKBEEKBEEIoEEIoEEIoEEIoEEJ+JP4PkT45sIvGAHQAAAAASUVORK5CYII=) |
Chi2 comparison of X-Y points and a function (chi2, p-values)
Source code name: "stat_chi2_p1d.py"
Programming language: Python
Topic: Statistics/Tests
DMelt Version 1. Last modified: 12/11/2015. License: Free
https://datamelt.org/code/cache/stat_chi2_p1d_2280.py
To run this script using the DataMelt IDE,
copy the above URL link to the menu [File]→[Read script from URL] of the DMelt IDE.
from jhplot import *
p=P1D("Data")
p.add(1, 9, 3); p.add(2, 21, 5); p.add(3, 36, 6)
f=F1D('10+x*x*x',1,3)
print p.compareChi2(f)
c=HPlot()
c.visible()
c.setAutoRange()
c.draw(p); c.draw(f)