2f:T6df,function controls(id, params) { let form function getParamOptions (type) { switch (type) { case 'unit': return { min: 0.01, max: 0.99, step: 0.01, value: 0.5 } case 'scale': return { min: 0.01, max: 3, step: 0.01, value: 1 } case 'location': return { min: -2, max: 2, step: 0.01, value: 0 } case 'degree': return { min: 1, max: 10, step: 1, value: 1 } default: return { min: 0, max: 3, step: 0.01, value: 1 } } } function addControl(param) { const options = getParamOptions(param.type) const value = typeof param.value === 'number' ? param.value : options.value const min = typeof param.min === 'number' ? param.min : options.min const max = typeof param.max === 'number' ? param.max : options.max return `

${param.name} = ${value}

` } var formHtml = "" params.forEach(function(pi) { formHtml += addControl(pi) }) form = html``; form.value = {} params.forEach(function(pi) { const options = getParamOptions(pi.type) form.value[pi.name] = typeof pi.value === 'number' ? pi.value : options.value form[pi.name].oninput = () => { form.value[pi.name] = form[pi.name].valueAsNumber document.getElementById(`${id}-${pi.name}`).textContent = form.value[pi.name] }; }) return form; }30:T90f,plot = function plot(options) { options = Object.assign({ min: -1, max: 100, size: 1e5, binSize: 10 }, options) const height = 200 const margin = ({top: 20, right: 20, bottom: 30, left: 40}) const binWidth = (options.max - options.min) / options.bins // Generate samples. const data = options.dist.sample(options.size) .filter(d => d >= options.min && d < options.max) // X scale and boundaries var x = d3.scaleLinear() .domain([options.min, options.max]) .range([margin.left, width - margin.right]) // Binned data var bins = d3.bin() .domain([options.min, options.max]) .thresholds(options.max - options.min)(data) // Build PDF/PMF. const pdf = d3.scaleLinear() .domain([options.min, options.max]) .ticks(options.max - options.min) .map(x => ({ x, y: options.size * options.dist.pdf(x) })) .filter(d => isFinite(d.y)) .slice(0, -1) // Y scale var y = d3.scaleLinear() .domain([0, d3.max(bins, d => d.length)]) .range([height - margin.bottom, margin.top]) const svg = d3.select(DOM.svg(width, height)); const bar = svg.append("g") .attr("fill", "#ccc") .attr('stroke', 'none') .selectAll("rect") .data(bins) .enter().append("rect") .attr("x", d => x(d.x0 - 0.5)) .attr("width", d => x(d.x1) - x(d.x0) - 1) .attr("y", d => y(d.length)) .attr("height", d => y(0) - y(d.length)) svg.append("g") .selectAll('.lollipop') .data(pdf).enter() .append('circle') .attr('cx', d => x(d.x)) .attr('cy', d => y(d.y)) .attr('r', 4) .attr('fill', 'red') svg.append("path") .datum(pdf.filter(d => d.y > 0)) .attr("fill", "none") .attr("stroke", "red") .attr("stroke-width", "2px") .attr("d", d3.line() .x(d => x(d.x)) .y(d => y(d.y))) // Add axes var axes = { x: g => g .attr('transform', `translate(0,${height - margin.bottom})`) .call(d3.axisBottom(x).tickSizeOuter(0)), y: g => g .attr('transform', `translate(${margin.left},0)`) .call(d3.axisLeft(y)) .call(g => g.select('.domain').remove()) } svg.append('g') .call(axes.x) svg.append('g') .call(axes.y) return svg.node() }31:[] 2d:{"title":"ranjs.dist discrete distributions","cells":[{"id":0,"value":"# ranjs.dist discrete distributions\n\nThe catalogue for continuous probability distributions implemented in `ranjs.dist`.","pinned":false,"hidden":false,"mode":"md"},{"id":448,"value":"toc()","pinned":false,"hidden":false,"mode":"ojs"},{"id":68,"value":"## Bernoulli","pinned":false,"hidden":false,"mode":"md"},{"id":75,"value":"viewof pBernoulli = controls('bernoulli', [\n {name: 'p', type: 'scale', max: 1, value: 0.5}\n])","pinned":false,"hidden":false,"mode":"ojs"},{"id":79,"value":"plot({\n min: -1,\n max: 2,\n dist: new ranjs.dist.Bernoulli(pBernoulli.p)\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":18,"value":"## beta-binomial","pinned":false,"hidden":false,"mode":"md"},{"id":20,"value":"viewof pBetaBinomial = controls('betabinomial', [\n {name: 'n', type: 'degree', max: 30, value: 20},\n {name: 'α', type: 'scale', value: 1.5},\n {name: 'β', type: 'scale', value: 2}\n])","pinned":false,"hidden":false,"mode":"ojs"},{"id":24,"value":"plot({\n min: -1,\n max: 40,\n dist: new ranjs.dist.BetaBinomial(pBetaBinomial.n, pBetaBinomial.α, pBetaBinomial.β)\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":100,"value":"## binomial","pinned":false,"hidden":false,"mode":"md"},{"id":103,"value":"viewof pBinomial = controls('binomial', [\n {name: 'n', type: 'degree', max: 30},\n {name: 'p', type: 'scale', max: 1, value: 0.5}\n])","pinned":false,"hidden":false,"mode":"ojs"},{"id":107,"value":"plot({\n min: -1,\n max: 40,\n dist: new ranjs.dist.Binomial(pBinomial.n, pBinomial.p)\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":111,"value":"## Borel","pinned":false,"hidden":false,"mode":"md"},{"id":123,"value":"viewof pBorel = controls('borel', [\n {name: 'μ', type: 'scale', max: 1, value: 0.5}\n])","pinned":false,"hidden":false,"mode":"ojs"},{"id":114,"value":"plot({\n min: 0,\n max: 20,\n dist: new ranjs.dist.Borel(pBorel.μ)\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":132,"value":"## Borel-Tanner","pinned":false,"hidden":false,"mode":"md"},{"id":135,"value":"viewof pBorelTanner = controls('boreltanner', [\n {name: 'μ', type: 'scale', max: 1, value: 0.5},\n {name: 'n', type: 'degree'}\n])","pinned":false,"hidden":false,"mode":"ojs"},{"id":138,"value":"plot({\n min: 0,\n max: 20,\n dist: new ranjs.dist.BorelTanner(pBorelTanner.μ, pBorelTanner.n)\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":166,"value":"## categorical","pinned":false,"hidden":false,"mode":"md"},{"id":169,"value":"viewof pCategorical = controls('categorical', [\n {name: 'w1', type: 'scale', max: 1, value: 0.5},\n {name: 'w2', type: 'scale', max: 1, value: 0.5},\n {name: 'w3', type: 'scale', max: 1, value: 0.5},\n])","pinned":false,"hidden":false,"mode":"ojs"},{"id":173,"value":"plot({\n min: -1,\n max: 3,\n dist: new ranjs.dist.Categorical([\n pCategorical.w1,\n pCategorical.w2,\n pCategorical.w3,\n ])\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":177,"value":"## Delaporte","pinned":false,"hidden":false,"mode":"md"},{"id":180,"value":"viewof pDelaporte = controls('delaporte', [\n {name: 'α', type: 'scale', max: 5, value: 0.5},\n {name: 'β', type: 'scale', max: 5, value: 0.5},\n {name: 'λ', type: 'scale', max: 5, value: 0.5},\n])","pinned":false,"hidden":false,"mode":"ojs"},{"id":184,"value":"plot({\n min: -1,\n max: 20,\n dist: new ranjs.dist.Delaporte(pDelaporte.α, pDelaporte.β, pDelaporte.λ)\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":188,"value":"## discrete uniform","pinned":false,"hidden":false,"mode":"md"},{"id":190,"value":"viewof pDiscreteUniform = controls('discreteuniform', [\n {name: 'xmin', type: 'degree', max: 10, value: 5},\n {name: 'xmax', type: 'degree', min: 11, max: 20, value: 15},\n])","pinned":false,"hidden":false,"mode":"ojs"},{"id":194,"value":"plot({\n min: 0,\n max: 22,\n dist: new ranjs.dist.DiscreteUniform(pDiscreteUniform.xmin, pDiscreteUniform.xmax)\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":205,"value":"## discrete Weibull","pinned":false,"hidden":false,"mode":"md"},{"id":208,"value":"viewof pDiscreteWeibull = controls('discreteweibull', [\n {name: 'q', type: 'scale', max: 0.99, value: 0.5},\n {name: 'β', type: 'scale', max: 5, value: 1},\n])","pinned":false,"hidden":false,"mode":"ojs"},{"id":211,"value":"plot({\n min: -1,\n max: 22,\n dist: new ranjs.dist.DiscreteWeibull(pDiscreteWeibull.q, pDiscreteWeibull.β)\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":217,"value":"## Flory Schulz","pinned":false,"hidden":false,"mode":"md"},{"id":220,"value":"viewof pFlorySchultz = controls('floryschultz', [\n {name: 'a', type: 'scale', min: 0.01, max: 0.99, value: 0.5},\n])","pinned":false,"hidden":false,"mode":"ojs"},{"id":223,"value":"plot({\n min: 0,\n max: 40,\n dist: new ranjs.dist.FlorySchulz(pFlorySchultz.a)\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":229,"value":"## generalized Hermite","pinned":false,"hidden":false,"mode":"md"},{"id":232,"value":"viewof pGeneralizedHermie = controls('generalizedhermite', [\n {name: 'a1', type: 'scale', max: 5, value: 0.5},\n {name: 'am', type: 'scale', max: 5, value: 0.5},\n {name: 'm', type: 'degree', min: 2, max: 10, value: 2},\n])","pinned":false,"hidden":false,"mode":"ojs"},{"id":235,"value":"plot({\n min: -1,\n max: 40,\n dist: new ranjs.dist.GeneralizedHermite(pGeneralizedHermie.a1, pGeneralizedHermie.am, pGeneralizedHermie.m)\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":240,"value":"## geometric","pinned":false,"hidden":false,"mode":"md"},{"id":247,"value":"viewof pGeometric = controls('geometric', [\n {name: 'p', type: 'scale', max: 1, value: 0.5},\n])","pinned":false,"hidden":false,"mode":"ojs"},{"id":250,"value":"plot({\n min: -1,\n max: 40,\n dist: new ranjs.dist.Geometric(pGeometric.p)\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":253,"value":"## heads-minus-tails","pinned":false,"hidden":false,"mode":"md"},{"id":256,"value":"viewof pHeadsMinusTails = controls('headsminustails', [\n {name: 'n', type: 'degree', max: 50, value: 1},\n])","pinned":false,"hidden":false,"mode":"ojs"},{"id":259,"value":"plot({\n min: -1,\n max: 20,\n dist: new ranjs.dist.HeadsMinusTails(pHeadsMinusTails.n)\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":267,"value":"## hypergeometric","pinned":false,"hidden":false,"mode":"md"},{"id":270,"value":"viewof pHypergeometric = controls('hypergeometric', [\n {name: 'N', type: 'degree', min: 30, max: 50, value: 10},\n {name: 'K', type: 'degree', max: 30, value: 5},\n {name: 'n', type: 'degree', max: 30, value: 5},\n])","pinned":false,"hidden":false,"mode":"ojs"},{"id":277,"value":"plot({\n min: -1,\n max: 20,\n dist: new ranjs.dist.Hypergeometric(pHypergeometric.N, pHypergeometric.K, pHypergeometric.n)\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":283,"value":"## log-series","pinned":false,"hidden":false,"mode":"md"},{"id":286,"value":"viewof pLogSeries = controls('logseries', [\n {name: 'p', type: 'scale', max: 0.99, value: 0.5},\n])","pinned":false,"hidden":false,"mode":"ojs"},{"id":290,"value":"plot({\n min: 0,\n max: 20,\n dist: new ranjs.dist.LogSeries(pLogSeries.p)\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":296,"value":"## negative-binomial","pinned":false,"hidden":false,"mode":"md"},{"id":299,"value":"viewof pNegativeBinomial = controls('negativebinomial', [\n {name: 'r', type: 'degree', max: 10, value: 0.5},\n {name: 'p', type: 'scale', max: 1, value: 0.5},\n])","pinned":false,"hidden":false,"mode":"ojs"},{"id":303,"value":"plot({\n min: -1,\n max: 50,\n dist: new ranjs.dist.NegativeBinomial(pNegativeBinomial.r, pNegativeBinomial.p)\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":313,"value":"## negative hypergeometric","pinned":false,"hidden":false,"mode":"md"},{"id":316,"value":"viewof pNegativeHypergeometric = controls('negativehypergeometric', [\n {name: 'N', type: 'degree', min: 40, max: 50, value: 40},\n {name: 'K', type: 'degree', min: 20, max: 30, value: 20},\n {name: 'r', type: 'degree', max: 10, value: 5},\n])","pinned":false,"hidden":false,"mode":"ojs"},{"id":320,"value":"plot({\n min: -1,\n max: 50,\n dist: new ranjs.dist.NegativeHypergeometric(pNegativeHypergeometric.N, pNegativeHypergeometric.K, pNegativeHypergeometric.r)\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":326,"value":"## Neyman type A","pinned":false,"hidden":false,"mode":"md"},{"id":329,"value":"viewof pNeymanA = controls('neymana', [\n {name: 'λ', type: 'scale', max: 5, value: 0.5},\n {name: 'ϕ', type: 'scale', max: 5, value: 0.5},\n])","pinned":false,"hidden":false,"mode":"ojs"},{"id":333,"value":"plot({\n min: -1,\n max: 50,\n dist: new ranjs.dist.NeymanA(pNeymanA.λ, pNeymanA.ϕ)\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":340,"value":"## Poisson","pinned":false,"hidden":false,"mode":"md"},{"id":343,"value":"viewof pPoisson = controls('poisson', [\n {name: 'λ', type: 'scale', max: 10, value: 0.5},\n])","pinned":false,"hidden":false,"mode":"ojs"},{"id":348,"value":"plot({\n min: -1,\n max: 20,\n dist: new ranjs.dist.Poisson(pPoisson.λ)\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":352,"value":"## Pólya-Aeppli","pinned":false,"hidden":false,"mode":"md"},{"id":355,"value":"viewof pPolyaAeppli = controls('polyaaeppli', [\n {name: 'λ', type: 'scale', max: 10, value: 0.5},\n {name: 'θ', type: 'scale', max: 0.99, value: 0.5},\n])","pinned":false,"hidden":false,"mode":"ojs"},{"id":358,"value":"plot({\n min: -1,\n max: 50,\n dist: new ranjs.dist.PolyaAeppli(pPolyaAeppli.λ, pPolyaAeppli.θ)\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":367,"value":"## Rademacher","pinned":false,"hidden":false,"mode":"md"},{"id":370,"value":"plot({\n min: -2,\n max: 3,\n dist: new ranjs.dist.Rademacher()\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":375,"value":"## Skellam","pinned":false,"hidden":false,"mode":"md"},{"id":378,"value":"viewof pSkellam = controls('skellam', [\n {name: 'μ1', type: 'scale', max: 10, value: 1},\n {name: 'μ2', type: 'scale', max: 10, value: 1},\n])","pinned":false,"hidden":false,"mode":"ojs"},{"id":382,"value":"plot({\n min: -20,\n max: 20,\n dist: new ranjs.dist.Skellam(pSkellam.μ1, pSkellam.μ2)\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":394,"value":"## soliton","pinned":false,"hidden":false,"mode":"md"},{"id":397,"value":"viewof pSoliton = controls('soliton', [\n {name: 'N', type: 'degree', max: 30, value: 1},\n])","pinned":false,"hidden":false,"mode":"ojs"},{"id":401,"value":"plot({\n min: 0,\n max: 20,\n dist: new ranjs.dist.Soliton(pSoliton.N)\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":407,"value":"## Yule-Simon","pinned":false,"hidden":false,"mode":"md"},{"id":410,"value":"viewof pYuleSimon = controls('yulesimon', [\n {name: 'ρ', type: 'scale', max: 10, value: 1},\n])","pinned":false,"hidden":false,"mode":"ojs"},{"id":413,"value":"plot({\n min: 0,\n max: 20,\n dist: new ranjs.dist.YuleSimon(pYuleSimon.ρ)\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":416,"value":"## zeta","pinned":false,"hidden":false,"mode":"md"},{"id":419,"value":"viewof pZeta = controls('zeta', [\n {name: 's', type: 'scale', min: 1.1, max: 10, value: 1.5},\n])","pinned":false,"hidden":false,"mode":"ojs"},{"id":422,"value":"plot({\n min: 0,\n max: 20,\n dist: new ranjs.dist.Zeta(pZeta.s)\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":431,"value":"## Zipf","pinned":false,"hidden":false,"mode":"md"},{"id":434,"value":"viewof pZipf = controls('zipf', [\n {name: 's', type: 'scale', max: 3, value: 0.5},\n {name: 'N', type: 'degree', max: 50, value: 10},\n])","pinned":false,"hidden":false,"mode":"ojs"},{"id":437,"value":"plot({\n min: 0,\n max: 52,\n dist: new ranjs.dist.Zipf(pZipf.s, pZipf.N)\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":5,"value":"ranjs = require('ranjs')","pinned":false,"hidden":false,"mode":"ojs"},{"id":8,"value":"$2f","pinned":false,"hidden":false,"mode":"ojs"},{"id":14,"value":"$30","pinned":false,"hidden":false,"mode":"ojs"},{"id":11,"value":"function toc(selector = \"h2,h3\", heading = \"Table of Contents\") {\n return Generators.observe(notify => {\n let headings = [];\n\n function observed() {\n const h = Array.from(document.querySelectorAll(selector));\n if (h.length !== headings.length || h.some((h, i) => headings[i] !== h)) {\n notify(html`${heading}

${DOM.text(h.textContent)}`,\n {onclick: e => (e.preventDefault(), h.scrollIntoView())}\n );\n })}`);\n }\n }\n\n const observer = new MutationObserver(observed);\n observer.observe(document.body, {childList: true, subtree: true});\n observed();\n return () => observer.disconnect();\n });\n}","pinned":false,"hidden":false,"mode":"ojs"}],"files":"$W31","stdlib":"1","titling":"auto"}