Skip to content

Instantly share code, notes, and snippets.

@genkuroki
Created March 5, 2018 07:28
Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save genkuroki/a6a6901d4594084c79b740dbc0ed967c to your computer and use it in GitHub Desktop.
Save genkuroki/a6a6901d4594084c79b740dbc0ed967c to your computer and use it in GitHub Desktop.
quantile(Chisq(2), 0.95)
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"metadata": {
"trusted": true
},
"cell_type": "code",
"source": "using Distributions\nusing Plots\n@show d = Chisq(2)\n@show q = quantile(d, 0.95)\nplot(x->pdf(d,x), 0.1, 10, label=\"pdf of d = Chisq($(d.ν))\")\nvline!([q], label=\"quantile(d, 0.95)\", ls=:dash)",
"execution_count": 1,
"outputs": [
{
"output_type": "execute_result",
"execution_count": 1,
"data": {
"text/html": "<?xml version=\"1.0\" encoding=\"utf-8\"?>\n<svg xmlns=\"http://www.w3.org/2000/svg\" xmlns:xlink=\"http://www.w3.org/1999/xlink\" width=\"600\" height=\"400\" viewBox=\"0 0 600 400\">\n<defs>\n <clipPath id=\"clip0800\">\n <rect x=\"0\" y=\"0\" width=\"600\" height=\"400\"/>\n </clipPath>\n</defs>\n<polygon clip-path=\"url(#clip0800)\" points=\"\n0,400 600,400 600,0 0,0 \n \" fill=\"#ffffff\" fill-opacity=\"1\"/>\n<defs>\n <clipPath id=\"clip0801\">\n <rect x=\"120\" y=\"0\" width=\"421\" height=\"400\"/>\n </clipPath>\n</defs>\n<polygon clip-path=\"url(#clip0800)\" points=\"\n32.4745,375.869 580.315,375.869 580.315,11.811 32.4745,11.811 \n \" fill=\"#ffffff\" fill-opacity=\"1\"/>\n<defs>\n <clipPath id=\"clip0802\">\n <rect x=\"32\" y=\"11\" width=\"549\" height=\"365\"/>\n </clipPath>\n</defs>\n<polyline clip-path=\"url(#clip0802)\" style=\"stroke:#000000; stroke-width:0.5; stroke-opacity:0.1; fill:none\" points=\"\n 137.875,375.869 137.875,11.811 \n \"/>\n<polyline clip-path=\"url(#clip0802)\" style=\"stroke:#000000; stroke-width:0.5; stroke-opacity:0.1; fill:none\" points=\"\n 250.252,375.869 250.252,11.811 \n \"/>\n<polyline clip-path=\"url(#clip0802)\" style=\"stroke:#000000; stroke-width:0.5; stroke-opacity:0.1; fill:none\" points=\"\n 362.63,375.869 362.63,11.811 \n \"/>\n<polyline clip-path=\"url(#clip0802)\" style=\"stroke:#000000; stroke-width:0.5; stroke-opacity:0.1; fill:none\" points=\"\n 475.007,375.869 475.007,11.811 \n \"/>\n<polyline clip-path=\"url(#clip0802)\" style=\"stroke:#000000; stroke-width:0.5; stroke-opacity:0.1; fill:none\" points=\"\n 32.4745,300.598 580.315,300.598 \n \"/>\n<polyline clip-path=\"url(#clip0802)\" style=\"stroke:#000000; stroke-width:0.5; stroke-opacity:0.1; fill:none\" points=\"\n 32.4745,222.527 580.315,222.527 \n \"/>\n<polyline clip-path=\"url(#clip0802)\" style=\"stroke:#000000; stroke-width:0.5; stroke-opacity:0.1; fill:none\" points=\"\n 32.4745,144.455 580.315,144.455 \n \"/>\n<polyline clip-path=\"url(#clip0802)\" style=\"stroke:#000000; stroke-width:0.5; stroke-opacity:0.1; fill:none\" points=\"\n 32.4745,66.384 580.315,66.384 \n \"/>\n<polyline clip-path=\"url(#clip0800)\" style=\"stroke:#000000; stroke-width:1; stroke-opacity:1; fill:none\" points=\"\n 32.4745,375.869 580.315,375.869 \n \"/>\n<polyline clip-path=\"url(#clip0800)\" style=\"stroke:#000000; stroke-width:1; stroke-opacity:1; fill:none\" points=\"\n 32.4745,375.869 32.4745,11.811 \n \"/>\n<polyline clip-path=\"url(#clip0800)\" style=\"stroke:#000000; stroke-width:1; stroke-opacity:1; fill:none\" points=\"\n 137.875,375.869 137.875,370.408 \n \"/>\n<polyline clip-path=\"url(#clip0800)\" style=\"stroke:#000000; stroke-width:1; stroke-opacity:1; fill:none\" points=\"\n 250.252,375.869 250.252,370.408 \n \"/>\n<polyline clip-path=\"url(#clip0800)\" style=\"stroke:#000000; stroke-width:1; stroke-opacity:1; fill:none\" points=\"\n 362.63,375.869 362.63,370.408 \n \"/>\n<polyline clip-path=\"url(#clip0800)\" style=\"stroke:#000000; stroke-width:1; stroke-opacity:1; fill:none\" points=\"\n 475.007,375.869 475.007,370.408 \n \"/>\n<polyline clip-path=\"url(#clip0800)\" style=\"stroke:#000000; stroke-width:1; stroke-opacity:1; fill:none\" points=\"\n 32.4745,300.598 40.6921,300.598 \n \"/>\n<polyline clip-path=\"url(#clip0800)\" style=\"stroke:#000000; stroke-width:1; stroke-opacity:1; fill:none\" points=\"\n 32.4745,222.527 40.6921,222.527 \n \"/>\n<polyline clip-path=\"url(#clip0800)\" style=\"stroke:#000000; stroke-width:1; stroke-opacity:1; fill:none\" points=\"\n 32.4745,144.455 40.6921,144.455 \n \"/>\n<polyline clip-path=\"url(#clip0800)\" style=\"stroke:#000000; stroke-width:1; stroke-opacity:1; fill:none\" points=\"\n 32.4745,66.384 40.6921,66.384 \n \"/>\n<g clip-path=\"url(#clip0800)\">\n<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:12; text-anchor:middle;\" transform=\"rotate(0, 137.875, 389.669)\" x=\"137.875\" y=\"389.669\">2</text>\n</g>\n<g clip-path=\"url(#clip0800)\">\n<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:12; text-anchor:middle;\" transform=\"rotate(0, 250.252, 389.669)\" x=\"250.252\" y=\"389.669\">4</text>\n</g>\n<g clip-path=\"url(#clip0800)\">\n<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:12; text-anchor:middle;\" transform=\"rotate(0, 362.63, 389.669)\" x=\"362.63\" y=\"389.669\">6</text>\n</g>\n<g clip-path=\"url(#clip0800)\">\n<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:12; text-anchor:middle;\" transform=\"rotate(0, 475.007, 389.669)\" x=\"475.007\" y=\"389.669\">8</text>\n</g>\n<g clip-path=\"url(#clip0800)\">\n<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:12; text-anchor:end;\" transform=\"rotate(0, 26.4745, 305.098)\" x=\"26.4745\" y=\"305.098\">0.1</text>\n</g>\n<g clip-path=\"url(#clip0800)\">\n<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:12; text-anchor:end;\" transform=\"rotate(0, 26.4745, 227.027)\" x=\"26.4745\" y=\"227.027\">0.2</text>\n</g>\n<g clip-path=\"url(#clip0800)\">\n<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:12; text-anchor:end;\" transform=\"rotate(0, 26.4745, 148.955)\" x=\"26.4745\" y=\"148.955\">0.3</text>\n</g>\n<g clip-path=\"url(#clip0800)\">\n<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:12; text-anchor:end;\" transform=\"rotate(0, 26.4745, 70.884)\" x=\"26.4745\" y=\"70.884\">0.4</text>\n</g>\n<polyline clip-path=\"url(#clip0802)\" style=\"stroke:#009af9; stroke-width:1; stroke-opacity:1; fill:none\" points=\"\n 32.4745,11.811 33.8326,16.2179 35.1907,20.5718 36.5488,24.8735 49.1064,62.2796 61.6641,95.7308 74.2217,125.645 86.7793,152.397 101.873,180.835 116.966,205.699 \n 129.727,224.266 142.488,240.841 156.35,256.836 170.212,270.975 182.992,282.551 195.772,292.883 223.174,311.446 254.922,327.991 283.137,339.243 307.878,347.034 \n 334.852,353.785 361.987,359.124 391.989,363.703 420.534,367.06 450.214,369.755 475.199,371.532 505.703,373.229 534.377,374.454 580.315,375.869 \n \"/>\n<polyline clip-path=\"url(#clip0802)\" style=\"stroke:#e26f46; stroke-width:1; stroke-opacity:1; fill:none\" stroke-dasharray=\"8, 5\" points=\"\n 362.15,375.869 362.15,11.811 \n \"/>\n<polygon clip-path=\"url(#clip0800)\" points=\"\n392.955,78.051 562.315,78.051 562.315,32.691 392.955,32.691 \n \" fill=\"#ffffff\" fill-opacity=\"1\"/>\n<polyline clip-path=\"url(#clip0800)\" style=\"stroke:#000000; stroke-width:1; stroke-opacity:1; fill:none\" points=\"\n 392.955,78.051 562.315,78.051 562.315,32.691 392.955,32.691 392.955,78.051 \n \"/>\n<polyline clip-path=\"url(#clip0800)\" style=\"stroke:#009af9; stroke-width:1; stroke-opacity:1; fill:none\" points=\"\n 398.955,47.811 434.955,47.811 \n \"/>\n<g clip-path=\"url(#clip0800)\">\n<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:12; text-anchor:start;\" transform=\"rotate(0, 440.955, 52.311)\" x=\"440.955\" y=\"52.311\">pdf of d = Chisq(2.0)</text>\n</g>\n<polyline clip-path=\"url(#clip0800)\" style=\"stroke:#e26f46; stroke-width:1; stroke-opacity:1; fill:none\" stroke-dasharray=\"8, 5\" points=\"\n 398.955,62.931 434.955,62.931 \n \"/>\n<g clip-path=\"url(#clip0800)\">\n<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:12; text-anchor:start;\" transform=\"rotate(0, 440.955, 67.431)\" x=\"440.955\" y=\"67.431\">quantile(d, 0.95)</text>\n</g>\n</svg>\n"
},
"metadata": {}
}
]
},
{
"metadata": {
"trusted": true
},
"cell_type": "code",
"source": "",
"execution_count": null,
"outputs": []
}
],
"metadata": {
"kernelspec": {
"name": "julia-0.6",
"display_name": "Julia 0.6.2",
"language": "julia"
},
"toc": {
"nav_menu": {},
"number_sections": true,
"sideBar": true,
"skip_h1_title": false,
"title_cell": "Table of Contents",
"title_sidebar": "Contents",
"toc_cell": false,
"toc_position": {},
"toc_section_display": true,
"toc_window_display": false
},
"language_info": {
"file_extension": ".jl",
"name": "julia",
"mimetype": "application/julia",
"version": "0.6.2"
},
"gist": {
"id": "",
"data": {
"description": "quantile(Chisq(2), 0.95)",
"public": true
}
}
},
"nbformat": 4,
"nbformat_minor": 2
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment