2a:null 2b:["$","title",null,{"children":"Gumbel distribution | Mattias Villani | Observable"}] 31:[] 2e:{"id":"bb517b53087c21d2","slug":"gumbel-distribution","title":"Gumbel distribution","thumbnail":"d8a1a001c908dbd86fd64cb58237e67d95afafa42be022b4e12b2143b6d907a6","access_level":"public","create_time":"$D2025-01-24T09:22:20.594Z","update_time":"$D2025-01-24T09:51:40.115Z","view_time":null,"first_public_time":"$D2025-01-24T09:22:20.594Z","fork":{"id":"0527b03232426322","slug":"generalized-extreme-value-distribution","title":"Generalized extreme value distribution","owner":{"id":"7592a7ab121f74a0","login":"mattiasvillani","name":"Mattias Villani"}},"stars":0,"starred":false,"owner":{"type":"team","id":"7592a7ab121f74a0","login":"mattiasvillani","name":"Mattias Villani","avatar_url":"https://avatars-next.observableusercontent.com/a440aaf39a3abefb86cbcc758fd7a5ab42555ae68d6825b52ec006f385a94d1d","tier":"pro_2024"},"version":1039,"publish_version":null,"role":null,"authored":false,"authors":[{"type":"individual","id":"56499be8cf43f0ef","login":"mattiasvillani","name":"Mattias Villani","avatar_url":"https://avatars-next.observableusercontent.com/a440aaf39a3abefb86cbcc758fd7a5ab42555ae68d6825b52ec006f385a94d1d","tier":null}],"members":[{"type":"individual","id":"56499be8cf43f0ef","login":"mattiasvillani","name":"Mattias Villani","avatar_url":"https://avatars-next.observableusercontent.com/a440aaf39a3abefb86cbcc758fd7a5ab42555ae68d6825b52ec006f385a94d1d","tier":null,"role":"editor"},{"type":"team","id":"7592a7ab121f74a0","login":"mattiasvillani","name":"Mattias Villani","avatar_url":"https://avatars-next.observableusercontent.com/a440aaf39a3abefb86cbcc758fd7a5ab42555ae68d6825b52ec006f385a94d1d","tier":"pro_2024","role":"editor"}],"comments":[],"custom_thumbnail":false,"trashed":false,"subscribed":null,"forks":1,"collection":{"id":"40c8ef5af1a51aff","slug":"distributions","title":"Statistical distributions","owner":{"id":"7592a7ab121f74a0","login":"mattiasvillani","name":"Mattias Villani"}},"type":"notebook","body":{"title":"Gumbel distribution","cells":[{"id":38,"value":"# Gumbel distribution","pinned":false,"hidden":false,"mode":"md"},{"id":629,"value":"

${tex`X\\sim \\mathrm{Gumbel}(\\mu,\\beta)`}","pinned":false,"hidden":false,"mode":"md"},{"id":496,"value":"The density function of a ${tex`X\\sim \\mathrm{Gumbel}(\\mu,\\beta)`} random variable is\n\n${tex`f(x)=\\frac{1}{\\beta}e^{-(z+e^{-z})}`}\n\nwhere ${tex`z = \\frac{x-\\mu}{\\beta}`}","pinned":false,"hidden":false,"mode":"md"},{"id":118,"value":"viewof inputs = Inputs.form({\n mu: Inputs.range([-1, 5], {value: 2, step: 0.1, label: tex`\\mu`}),\n beta: Inputs.range([0, 5], {value: 1, step: 0.01, label: tex`\\beta`}),\n quantile: Inputs.range([-10,10], {value: 2, step: 0.01, label: \"quantile:\"})\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":461,"value":"dist_quantile = tex`\\text{If } X \\sim \\mathrm{Gumbel}(${inputs.mu},${(inputs.beta).toPrecision(2)}) \\text{ then } `","pinned":false,"hidden":false,"mode":"ojs"},{"id":555,"value":"moments = {\n \n return tex`\n \\begin{aligned}\n &\\mathrm{E}( X) = ${(mean).toPrecision(3)} \\\\[1em]\n &\\mathrm{Stdev( X)} =${(Math.sqrt(variance)).toPrecision(3)}\\\\[1em]\n &P( X\\leq ${inputs.quantile}) =${cdf.toPrecision(4)}\n \\end{aligned}\n `\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":305,"value":"plt_pdf = Plot.plot({\n color: {\n legend: true\n },\n x: {\n label: \"x\",\n axis: true,\n //domain: support(inputs.mu, inputs.sigma)\n },\n y: {\n label: \"f(x)\",\n //axis: false,\n //domain: [0,jstat.chisquare.pdf(mode, nu1, nu2)]\n },\n marks: [\n Plot.ruleY([0]),\n //Plot.ruleX([-5]),\n Plot.line(pdfdata, {x: \"x\", y: \"pdf\", stroke : \"steelblue\", strokeWidth: 2}),\n Plot.areaY(pdfdata, {filter: d => d.x <= inputs.quantile, x: \"x\", y: \"pdf\", fill: \"steelblue\", opacity: 0.2})\n ]\n })","pinned":false,"hidden":false,"mode":"ojs"},{"id":109,"value":"
\n
\n
\n
\n
\n
\n
\n
\n
\n
\n
\n
\n
\n
\n
\n
\n\n---\n## 🙈 Appendix 🙈\n\nHere is some code and data details that you probably do not want to read.","pinned":false,"hidden":false,"mode":"md"},{"id":221,"value":"#### Imports","pinned":false,"hidden":false,"mode":"md"},{"id":26,"value":"jstat = require('jstat')","pinned":false,"hidden":false,"mode":"ojs"},{"id":750,"value":"mathjs = require('mathjs')","pinned":false,"hidden":false,"mode":"ojs"},{"id":899,"value":"digamma = require( 'https://cdn.jsdelivr.net/gh/stdlib-js/math-base-special-digamma@umd/browser.js' )","pinned":false,"hidden":false,"mode":"ojs"},{"id":828,"value":"#### Functions","pinned":false,"hidden":false,"mode":"md"},{"id":909,"value":"function t(x){\n return Math.exp(-(x-inputs.mu)/inputs.beta)\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":903,"value":"function pdfgumbel(x, mu, beta){\n const z = (x-mu)/beta\n return (1/beta)*Math.exp(-(z + Math.exp(-z)))\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":907,"value":"function cdfgumbel(x, mu, beta){\n return Math.exp(-Math.exp(-(x-inputs.mu)/inputs.beta))\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":223,"value":"#### Data","pinned":false,"hidden":false,"mode":"md"},{"id":896,"value":"γ = -digamma(1) // Euler's number","pinned":false,"hidden":false,"mode":"ojs"},{"id":850,"value":"mean = inputs.mu + inputs.beta*γ","pinned":false,"hidden":false,"mode":"ojs"},{"id":766,"value":"variance = ((Math.PI**2)/6)*inputs.beta**2","pinned":false,"hidden":false,"mode":"ojs"},{"id":1025,"value":"extrema = [-3,10]","pinned":true,"hidden":false,"mode":"ojs"},{"id":17,"value":"pdfdata = {\n const x = d3.range(extrema[0], extrema[1], 0.01); \n const pdf = x.map(x => ({x: x, pdf: pdfgumbel(x, inputs.mu, inputs.beta)}));\n return pdf\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":19,"value":"cdf = cdfgumbel(inputs.quantile, inputs.mu, inputs.beta);","pinned":false,"hidden":false,"mode":"ojs"},{"id":229,"value":"#### Style","pinned":false,"hidden":false,"mode":"md"},{"id":202,"value":"","pinned":false,"hidden":false,"mode":"html"},{"id":58,"value":"viewof md_colour = Inputs.color({label: \"Highlight color cells\", value: \"#F5F5DC\"})","pinned":false,"hidden":false,"mode":"ojs"}],"files":"$W31","stdlib":"1","titling":"auto"}} 2f:"$2e:body"