31:T489,viewof ϕvect = { switch (input.p) { case 0: return Inputs.form([Inputs.text({value: "No AR component", label: "", disabled: true})]) case 1: return Inputs.form([ Inputs.range([-1, 2], {value: 0.5, step: 0.01, label: tex`\phi_1`}) ]) break; case 2: return Inputs.form([ Inputs.range([-2, 2], {value: 0.5, step: 0.01, label: tex`\phi_1`}), Inputs.range([-1, 1], {value: 0.0, step: 0.01, label: tex`\phi_2`}) ]) break; case 3: return Inputs.form([ Inputs.range([-1, 1], {value: 0.5, step: 0.01, label: tex`\phi_1`}), Inputs.range([-1, 1], {value: 0.0, step: 0.01, label: tex`\phi_2`}), Inputs.range([-1, 1], {value: 0.0, step: 0.01, label: tex`\phi_3`}) ]) break; case 4: return Inputs.form([ Inputs.range([-1, 1], {value: 0.5, step: 0.01, label: tex`\phi_1`}), Inputs.range([-1, 1], {value: 0.0, step: 0.01, label: tex`\phi_2`}), Inputs.range([-1, 1], {value: 0.0, step: 0.01, label: tex`\phi_3`}), Inputs.range([-1, 1], {value: 0.0, step: 0.01, label: tex`\phi_4`}) ]) break; } }32:T49d,viewof θvect = { switch (input.q) { case 0: return Inputs.form([Inputs.text({value: "No MA component", label: "", disabled: true})]) case 1: return Inputs.form([ Inputs.range([-1, 1], {value: 0.5, step: 0.01, label: tex`\theta_1`}) ]) break; case 2: return Inputs.form([ Inputs.range([-2, 2], {value: 0.5, step: 0.01, label: tex`\theta_1`}), Inputs.range([-1, 1], {value: 0.0, step: 0.01, label: tex`\theta_2`}) ]) break; case 3: return Inputs.form([ Inputs.range([-1, 1], {value: 0.5, step: 0.01, label: tex`\theta_1`}), Inputs.range([-1, 1], {value: 0.0, step: 0.01, label: tex`\theta_2`}), Inputs.range([-1, 1], {value: 0.0, step: 0.01, label: tex`\theta_3`}) ]) break; case 4: return Inputs.form([ Inputs.range([-1, 1], {value: 0.5, step: 0.01, label: tex`\theta_1`}), Inputs.range([-1, 1], {value: 0.0, step: 0.01, label: tex`\theta_2`}), Inputs.range([-1, 1], {value: 0.0, step: 0.01, label: tex`\theta_3`}), Inputs.range([-1, 1], {value: 0.0, step: 0.01, label: tex`\theta_4`}) ]) break; } }33:[] 2f:{"title":"Spectral density - ARMA(p,q) process","cells":[{"id":0,"value":"# Spectral density - ARMA(p,q) process","pinned":false,"hidden":false,"mode":"md"},{"id":194,"value":"The ARMA(p,q) process is of the form\n\n${tex`\n\\phi_p(L)(y_t-\\mu) = \\theta_q(L)\\varepsilon_t,\n`}\n\nwhere\n\n${tex`\\phi_p(L) = 1 - \\phi_1L - \\cdots - \\phi_pL^p`} is the AR lag polynomial,
\n${tex`\\theta_q(L) = 1 + \\theta_1L + \\cdots + \\theta_qL^q`}is the MA \nlag polynomial and
\n${tex`\\{ \\varepsilon_t \\}_{t=1}`} is an iid sequence from ${tex`N(0,\\sigma_\\varepsilon^2)`}.
\n
\n\nThe spectral density of an ARMA(p,q) process is
\n\n${tex`\nf(\\omega)=\\frac{\\sigma^{2}}{2\\pi}\\left|\\frac{\\theta_q(e^{-i\\omega})}{\\phi_{p}(e^{-i\\omega})}\\right|^{2}\n\\text{ for }-\\pi <\\omega\\leq \\pi`}","pinned":false,"hidden":false,"mode":"md"},{"id":247,"value":"#### Model settings","pinned":false,"hidden":false,"mode":"md"},{"id":162,"value":"viewof input = Inputs.form({\n p: Inputs.range([0, 4], {value: 1, step: 1, label: \"number of AR lags\"}),\n q: Inputs.range([0, 4], {value: 1, step: 1, label: \"number of MA lags\"})\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":249,"value":"#### Parameter settings","pinned":false,"hidden":false,"mode":"md"},{"id":244,"value":"viewof σ = Inputs.range([0.01, 5], {label: tex`\\text{error stdev, }\\sigma`, step: 0.01})","pinned":false,"hidden":false,"mode":"ojs"},{"id":116,"value":"$31","pinned":false,"hidden":false,"mode":"ojs"},{"id":121,"value":"$32","pinned":false,"hidden":false,"mode":"ojs"},{"id":152,"value":"Plot.plot({\n y: {label: \"Spectral density f(ω)\"},\n x: {label: \"Angular frequency ω\"},\n title: md`Spectral density of an ARMA(${input.p},${input.q}) process`,\n marks: [\n Plot.ruleY([0]),\n Plot.ruleX([0]),\n Plot.lineY(specdensdata, {x: \"ω\", y: \"specdens\", stroke: \"steelblue\", strokeWidth: 2})\n ]\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":13,"value":"
\n
\n
\n
\n
\n
\n
\n
\n
\n
\n
\n
\n
\n
\n
\n
\n\n---\n## 🙈 Appendix 🙈\nHere follow code and details that you probably do not want to read.","pinned":false,"hidden":false,"mode":"md"},{"id":15,"value":"#### Imports","pinned":false,"hidden":false,"mode":"md"},{"id":11,"value":"math = require(\"mathjs\");","pinned":false,"hidden":false,"mode":"ojs"},{"id":96,"value":"#### Data","pinned":false,"hidden":false,"mode":"md"},{"id":254,"value":"ϕvect_value = {\n if (ϕvect[0] == \"No AR component\"){\n return [0]\n }else{\n return ϕvect\n }\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":256,"value":"θvect_value = {\n if (θvect[0] == \"No MA component\"){\n return [0]\n }else{\n return θvect\n }\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":99,"value":"specdensdata = {\n const ωgrid = d3.range(0,Math.PI, 0.01);\n return ωgrid.map(d => ({ω: d, specdens: SpecDens(d, ϕvect_value, θvect_value, σ)}))\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":94,"value":"#### Functions","pinned":false,"hidden":false,"mode":"md"},{"id":5,"value":"function SpecDens(ω, ϕvect, θvect, σ){\n const eiw = math.exp(math.complex(0,-ω))\n var ARpoly = math.complex(1,0);\n var lag = 0;\n for (let ϕ of ϕvect){\n lag = lag + 1;\n ARpoly = math.add(ARpoly, math.multiply(math.complex(-ϕ,0),math.pow(eiw,lag)))\n } \n var MApoly = math.complex(1,0);\n var lag = 0;\n for (let θ of θvect){\n lag = lag + 1;\n MApoly = math.add(MApoly, math.multiply(math.complex(θ,0),math.pow(eiw,lag)))\n }\n return (σ**2/(Math.PI))*(math.abs(MApoly)**2 / math.abs(ARpoly)**2)\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":229,"value":"function ACF(k, ϕvect, θvect, σ){\n const p = ϕvect.length;\n const q = θvect.length;\n if (p == 1 && q ==0){\n return ϕvect[1]**k\n }\n else{\n return 0\n }\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":232,"value":"d3.range(1,5).map(k => ACF(k, ϕvect, θvect, σ))","pinned":true,"hidden":false,"mode":"ojs"}],"files":"$W33","stdlib":"1","titling":"auto"}