docs: function(size, p) {
           return {
              d: function(x, logp) { return dbinom(size, p, logp)(x); },
              p: function(q, lowerTail, logp) {
                 return pbinom(size, p, lowerTail, logp)(q);
              },
              q: function(prob, lowerTail, logp) {
                 return qbinom(size, p, lowerTail, logp)(prob);
              },
              r: function() { return rbinom(size, p)(); }
           };
        }