How does scipy calculate the binomial CDF so fast?












1












$begingroup$


I ran scipy.stats.binom.cdf(500006, 1000000, 0.5) and it took less than a milisecond. This is crazy as binomial CDF involves summing up a bunch of binomial coefficients. What approximation algorithm goes behind the implementation? Thanks!






binomial computational-statistics scipy







share|cite|improve this question









$endgroup$

















    1












    $begingroup$


    I ran scipy.stats.binom.cdf(500006, 1000000, 0.5) and it took less than a milisecond. This is crazy as binomial CDF involves summing up a bunch of binomial coefficients. What approximation algorithm goes behind the implementation? Thanks!






    binomial computational-statistics scipy







    share|cite|improve this question









    $endgroup$















      1












      1








      1





      $begingroup$


      I ran scipy.stats.binom.cdf(500006, 1000000, 0.5) and it took less than a milisecond. This is crazy as binomial CDF involves summing up a bunch of binomial coefficients. What approximation algorithm goes behind the implementation? Thanks!






      binomial computational-statistics scipy







      share|cite|improve this question









      $endgroup$




      I ran scipy.stats.binom.cdf(500006, 1000000, 0.5) and it took less than a milisecond. This is crazy as binomial CDF involves summing up a bunch of binomial coefficients. What approximation algorithm goes behind the implementation? Thanks!






      binomial computational-statistics scipy




      binomial computational-statistics scipy






      share|cite|improve this question













      share|cite|improve this question











      share|cite|improve this question




      share|cite|improve this question










      asked 2 hours ago









      John JiangJohn Jiang

      42628




      42628






















          1 Answer
          1






          active

          oldest

          votes


















          4












          $begingroup$

          There's a connection between the binomial cdf and the beta cdf. As a result, most software implementations of binomial cdf use a call to a regularized incomplete beta function. That's what scipy.stats does also, as it says in the documentation of special.bdtr (which is what gets called when you ask for a binomial cdf)



          https://docs.scipy.org/doc/scipy/reference/generated/scipy.special.bdtr.html






          share|cite|improve this answer











          $endgroup$













            Your Answer





            StackExchange.ifUsing("editor", function () {
            return StackExchange.using("mathjaxEditing", function () {
            StackExchange.MarkdownEditor.creationCallbacks.add(function (editor, postfix) {
            StackExchange.mathjaxEditing.prepareWmdForMathJax(editor, postfix, [["$", "$"], ["\\(","\\)"]]);
            });
            });
            }, "mathjax-editing");

            StackExchange.ready(function() {
            var channelOptions = {
            tags: "".split(" "),
            id: "65"
            };
            initTagRenderer("".split(" "), "".split(" "), channelOptions);

            StackExchange.using("externalEditor", function() {
            // Have to fire editor after snippets, if snippets enabled
            if (StackExchange.settings.snippets.snippetsEnabled) {
            StackExchange.using("snippets", function() {
            createEditor();
            });
            }
            else {
            createEditor();
            }
            });

            function createEditor() {
            StackExchange.prepareEditor({
            heartbeatType: 'answer',
            autoActivateHeartbeat: false,
            convertImagesToLinks: false,
            noModals: true,
            showLowRepImageUploadWarning: true,
            reputationToPostImages: null,
            bindNavPrevention: true,
            postfix: "",
            imageUploader: {
            brandingHtml: "Powered by u003ca class="icon-imgur-white" href="https://imgur.com/"u003eu003c/au003e",
            contentPolicyHtml: "User contributions licensed under u003ca href="https://creativecommons.org/licenses/by-sa/3.0/"u003ecc by-sa 3.0 with attribution requiredu003c/au003e u003ca href="https://stackoverflow.com/legal/content-policy"u003e(content policy)u003c/au003e",
            allowUrls: true
            },
            onDemand: true,
            discardSelector: ".discard-answer"
            ,immediatelyShowMarkdownHelp:true
            });


            }
            });














            draft saved

            draft discarded


















            StackExchange.ready(
            function () {
            StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fstats.stackexchange.com%2fquestions%2f389663%2fhow-does-scipy-calculate-the-binomial-cdf-so-fast%23new-answer', 'question_page');
            }
            );

            Post as a guest















            Required, but never shown

























            1 Answer
            1






            active

            oldest

            votes








            1 Answer
            1






            active

            oldest

            votes









            active

            oldest

            votes






            active

            oldest

            votes









            4












            $begingroup$

            There's a connection between the binomial cdf and the beta cdf. As a result, most software implementations of binomial cdf use a call to a regularized incomplete beta function. That's what scipy.stats does also, as it says in the documentation of special.bdtr (which is what gets called when you ask for a binomial cdf)



            https://docs.scipy.org/doc/scipy/reference/generated/scipy.special.bdtr.html






            share|cite|improve this answer











            $endgroup$


















              4












              $begingroup$

              There's a connection between the binomial cdf and the beta cdf. As a result, most software implementations of binomial cdf use a call to a regularized incomplete beta function. That's what scipy.stats does also, as it says in the documentation of special.bdtr (which is what gets called when you ask for a binomial cdf)



              https://docs.scipy.org/doc/scipy/reference/generated/scipy.special.bdtr.html






              share|cite|improve this answer











              $endgroup$
















                4












                4








                4





                $begingroup$

                There's a connection between the binomial cdf and the beta cdf. As a result, most software implementations of binomial cdf use a call to a regularized incomplete beta function. That's what scipy.stats does also, as it says in the documentation of special.bdtr (which is what gets called when you ask for a binomial cdf)



                https://docs.scipy.org/doc/scipy/reference/generated/scipy.special.bdtr.html






                share|cite|improve this answer











                $endgroup$



                There's a connection between the binomial cdf and the beta cdf. As a result, most software implementations of binomial cdf use a call to a regularized incomplete beta function. That's what scipy.stats does also, as it says in the documentation of special.bdtr (which is what gets called when you ask for a binomial cdf)



                https://docs.scipy.org/doc/scipy/reference/generated/scipy.special.bdtr.html







                share|cite|improve this answer














                share|cite|improve this answer



                share|cite|improve this answer








                edited 6 mins ago

























                answered 2 hours ago









                Glen_b♦Glen_b

                210k22401743




                210k22401743






























                    draft saved

                    draft discarded




















































                    Thanks for contributing an answer to Cross Validated!


                    • Please be sure to answer the question. Provide details and share your research!

                    But avoid …



                    • Asking for help, clarification, or responding to other answers.

                    • Making statements based on opinion; back them up with references or personal experience.


                    Use MathJax to format equations. MathJax reference.


                    To learn more, see our tips on writing great answers.




                    draft saved


                    draft discarded














                    StackExchange.ready(
                    function () {
                    StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fstats.stackexchange.com%2fquestions%2f389663%2fhow-does-scipy-calculate-the-binomial-cdf-so-fast%23new-answer', 'question_page');
                    }
                    );

                    Post as a guest















                    Required, but never shown





















































                    Required, but never shown














                    Required, but never shown












                    Required, but never shown







                    Required, but never shown

































                    Required, but never shown














                    Required, but never shown












                    Required, but never shown







                    Required, but never shown







                    -

                    Popular posts from this blog

                    SQL Server 17 - Attemping to backup to remote NAS but Access is denied

                    Image
                    .everyoneloves__top-leaderboard:empty,.everyoneloves__mid-leaderboard:empty,.everyoneloves__bot-mid-leaderboard:empty{ margin-bottom:0; } 0 I was attemping to use Ola Hallengrens scripts to backup to a remote NAS, but it was failing and not giving any direct errors, so I tried doing it manually using Server Management Studio and got the following error: TITLE: Microsoft SQL Server Management Studio Failed to retrieve data for this request. (Microsoft.SqlServer.Management.Sdk.Sfc) For help, click: http://go.microsoft.com/fwlink?ProdName=Microsoft%20SQL%20Server&LinkId=20476 ADDITIONAL INFORMATION: An exception occurred while executing a Transact-SQL statement or batch. (Microsoft.SqlServer.ConnectionInfo) The operating system returned the error '0x80070005(Access is denied.)&#
                    Read more

                    Always On Availability groups resolving state after failover - Remote harden of transaction...

                    Image
                    0 My explanation will be opinion based as it is based on my findings. If more information is needed, I am always happy to add it. The setup We have 4 SQL Server 2016 Instances, and all 4 underlying servers are in a WSFC (In my eyes a strange setup since Always On AG's are 2 by 2. But that is a different question.) 4 nodes + Fileshare witness Version: 13.0.5026.0 (all nodes) The instances that are having the failover issues are instances 1 & 2. They have 3 AG's between them. AG1 is primary on instance 1 and is having the issue. AG2 and AG3 are primary on instance 2. Instances 3 and 4 have their own separate AG's (2). AG1's settings The week before Last week, after checking for long running queries (none) and the synchronization state (green accross the board) we initiated a
                    Read more

                    Restoring from pg_dump with foreign key constraints

                    Image
                    5 2 In restoring a database from a pg_dump , a number of errors are being generated and the whole table is subsequently being ignored. An example: ERROR: insert or update on table "channelproducts" violates foreign key constraint "fk_rails_dfaae373a5" DETAIL: Key (channel_id)=(1) is not present in table "channels". Interestingly enough, I've noted that all these instances are popping up because of the sequence of loading. channels is after channelproducts both alphabetically and in the file and thus I can understand why postgres complains about having to create a child without a parent. caveat: the foreign key is being generated automatically by rails 4.2: I could remove the problem at the source but that still does not really solve the problem... version: PostgreSQL
                    Read more