28:null 29:["$","title",null,{"children":"Z-distribution | Mattias Villani | Observable"}] 2f:[] 2c:{"id":"cc5c11b1c594861c","slug":null,"title":"Z-distribution","thumbnail":"bcfe365ae2c9d64104598d853f47c72a6c3781536dbfabdaf58e633ed8d8c232","access_level":"public","create_time":"$D2024-08-22T20:28:59.608Z","update_time":"$D2024-08-22T20:29:18.865Z","view_time":null,"first_public_time":"$D2024-08-22T20:28:59.608Z","fork":{"id":"811384ce5401e838","slug":"z-distribution","title":"Fisher Z-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":887,"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":0,"collection":null,"type":"notebook","body":{"title":"Z-distribution","cells":[{"id":38,"value":"# *Z*-distribution","pinned":false,"hidden":false,"mode":"md"},{"id":722,"value":"html`

${tex`X \\sim \\mathcal{Z}(\\alpha, \\beta, \\mu, \\sigma) \\text{ for }X\\in\\mathbb{R}, \\alpha>0, \\beta>0, \\mu \\in \\mathbb{R}, \\sigma >0\\\\[1em]`}`","pinned":false,"hidden":false,"mode":"ojs"},{"id":118,"value":"viewof myinputs = Inputs.form([\n Inputs.range([0.1, 100], {value: 0.5, step: 0.1, label: tex`\\alpha`}),\n Inputs.range([0.1, 100], {value: 0.5, step: 0.1, label: tex`\\beta`}),\n Inputs.range([-10, 10], {value: 0, step: 0.1, label: tex`\\mu`}),\n Inputs.range([0.1, 10], {value: 1, step: 0.1, label: tex`\\sigma`}),\n Inputs.range([-10, 10], {value: -1.96, step: 0.01, label: \"quantile\"})\n ])","pinned":false,"hidden":false,"mode":"ojs"},{"id":450,"value":"viewof options = Inputs.form([\n Inputs.toggle({value: true, label: \"show closest normal\"}),\n Inputs.toggle({label: \"show pdf in logs\", value: false})\n ])","pinned":false,"hidden":false,"mode":"ojs"},{"id":757,"value":"tailprobs = \n{\n var tailprobs;\n if (shownormal){\n tailprobs = tex`\\text{Normal distribution: } P(X\\leq ${myinputs[4]})=${normcdf.toPrecision(4)} \\\\ \\text{Z-distribution: } P(X\\leq ${myinputs[4]})=${zcdf.toPrecision(4)}`\n }\n else{\n tailprobs = tex`\\text{Z-distribution: } P(X\\leq ${myinputs[4]})=${zcdf.toPrecision(4)}`\n }\n return tailprobs;\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":305,"value":"plt = Plot.plot({\n color: {\n legend: shownormal\n },\n x: {\n label: \"x\",\n axis: true\n },\n y: {\n axis: false,\n label: yaxislabel,\n transform: (x) => {if (showlogpdf){return Math.log(x)}else{return x}}\n //domain: [0,1.05*jstat.normal.pdf(mu, mu, sigma)]\n },\n title: yaxislabel,\n marks: [\n Plot.ruleY([0]),\n Plot.line(data, {filter: d => d.dist == \"z\", x: \"x\", y: \"pdf\", stroke : \"dist\", strokeWidth: 2}),\n Plot.line(data, {filter: d => d.dist == \"normal\" && shownormal, x: \"x\", y: \"pdf\", stroke : \"dist\", strokeWidth: 2}),\n Plot.areaY(data, {filter: d => d.x <= myinputs[4] && d.dist == \"normal\" && shownormal, x: \"x\", y: \"pdf\", fill: \"steelblue\", opacity: 0.2}),\n Plot.areaY(data, {filter: d => d.x <= myinputs[4] && d.dist == \"z\", x: \"x\", y: \"pdf\", fill: \"orange\", opacity: 0.2}),\n \n ]\n })","pinned":false,"hidden":false,"mode":"ojs"},{"id":828,"value":"viewof n = Inputs.range([1, 100], {value: 50, label: \"sample size, n\", step: 1})","pinned":false,"hidden":false,"mode":"ojs"},{"id":832,"value":"pltmax = Plot.plot({\n color: {\n legend: true\n },\n x: {\n label: \"z_max\",\n axis: true\n },\n y: {\n axis: false,\n //domain: [0,1.05*jstat.normal.pdf(mu, mu, sigma)]\n },\n title: \"pdf for the maximum of sample of size n\",\n marks: [\n Plot.ruleY([0]),\n Plot.line(datapdfmax, {filter: d => d.dist == \"z\", x: \"x\", y: \"pdf\", stroke : \"dist\", strokeWidth: 2}),\n Plot.line(datapdfmax, {filter: d => d.dist == \"normal\", x: \"x\", y: \"pdf\", stroke : \"dist\", strokeWidth: 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 only code and details that you probably don't 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":791,"value":"digamma = require( 'https://cdn.jsdelivr.net/gh/stdlib-js/math-base-special-digamma@umd/browser.js' )","pinned":true,"hidden":false,"mode":"ojs"},{"id":804,"value":"trigamma = require( 'https://cdn.jsdelivr.net/gh/stdlib-js/math-base-special-trigamma@umd/browser.js' )","pinned":true,"hidden":false,"mode":"ojs"},{"id":223,"value":"#### Data","pinned":false,"hidden":false,"mode":"md"},{"id":874,"value":"yaxislabel = {if (showlogpdf){return \"log pdf\"}else{return \"pdf\"}}","pinned":true,"hidden":false,"mode":"ojs"},{"id":862,"value":"shownormal = options[0]","pinned":false,"hidden":false,"mode":"ojs"},{"id":863,"value":"showlogpdf = options[1]","pinned":false,"hidden":false,"mode":"ojs"},{"id":511,"value":"alpha = myinputs[0]","pinned":false,"hidden":false,"mode":"ojs"},{"id":531,"value":"beta = myinputs[1]","pinned":false,"hidden":false,"mode":"ojs"},{"id":533,"value":"mu = myinputs[2]","pinned":false,"hidden":false,"mode":"ojs"},{"id":535,"value":"sigma = myinputs[3]","pinned":false,"hidden":false,"mode":"ojs"},{"id":819,"value":"quantile = myinputs[4]","pinned":false,"hidden":false,"mode":"ojs"},{"id":17,"value":"data = {\n const x = d3.range(-20, 20, 0.01);\n const normpdf = x.map(x => ({x: x, pdf: jstat.normal.pdf(x, Zmean, Zsd), dist: \"normal\"}));\n const zpdf = x.map(x => ({x: x, pdf: Zpdf(x, alpha, beta, mu, sigma), dist: \"z\"}));\n const data = normpdf.concat(zpdf)\n return data\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":735,"value":"logdata = {\n const x = d3.range(-20, 20, 0.01);\n const normpdf = x.map(x => ({x: x, pdf: Math.log(jstat.normal.pdf(x, Zmean, Zsd)), dist: \"normal\"}));\n const zpdf = x.map(x => ({x: x, pdf: Zpdf(x, alpha, beta, mu, sigma, true), dist: \"z\"}));\n const data = normpdf.concat(zpdf)\n return data\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":19,"value":"normcdf = jstat.normal.cdf(quantile, Zmean, Zsd);","pinned":false,"hidden":false,"mode":"ojs"},{"id":675,"value":"zcdf = jstat.beta.cdf(sigmoid((quantile-mu)/sigma), alpha, beta);","pinned":false,"hidden":false,"mode":"ojs"},{"id":834,"value":"datapdfmax = {\n const x = d3.range(0, 20, 0.01);\n const normpdf = x.map(x => ({x: x, pdf: normalpdfmax(x, Zmean, Zsd), dist: \"normal\"}));\n const zpdf = x.map(x => ({x: x, pdf: Zpdfmax(x, alpha, beta, mu, sigma), dist: \"z\"}));\n const data = normpdf.concat(zpdf)\n return data\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":839,"value":"function normalpdfmax(x, mu, sigma){\n return jstat.normal.pdf(x, mu, sigma)*jstat.normal.cdf(x, mu, sigma)**(n-1)\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":542,"value":"#### Functions","pinned":false,"hidden":false,"mode":"md"},{"id":544,"value":"function sigmoid(x){\n return 1/(1 + Math.exp(-x))\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":558,"value":"function logsigmoid(x){\n return -Math.log(1 + Math.exp(-x))\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":551,"value":"function Zpdf(x, alpha, beta, mu = 0, sigma = 1, log = false){\n const z = (x - mu)/sigma\n if (log){\n return alpha*logsigmoid(z) + beta*logsigmoid(-z) - Math.log(sigma) - jstat.betaln(alpha, beta)\n }else{\n return sigmoid(z)**alpha*sigmoid(-z)**beta/(sigma*jstat.betafn(alpha,beta))\n }\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":816,"value":"Zcdf = function(x, alpha, beta, mu, sigma){\n return jstat.beta.cdf(sigmoid((x-mu)/sigma), alpha, beta)\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":660,"value":"Zmean = mu + sigma*(digamma(alpha) - digamma(beta))","pinned":false,"hidden":false,"mode":"ojs"},{"id":670,"value":"Zsd = Math.sqrt(sigma**2*(trigamma(alpha) + trigamma(beta)))","pinned":false,"hidden":false,"mode":"ojs"},{"id":826,"value":"function Zpdfmax(x, alpha, beta, mu, sigma){\n return Zpdf(x, alpha, beta, mu, sigma)*Zcdf(x, alpha, beta, mu, sigma)**(n-1)\n}","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":"$W2f","stdlib":"1","titling":"auto"}} 2d:"$2c:body"