Last active
May 10, 2016 14:01
-
-
Save vassilyl/6d1607d01aebd3a5a62e to your computer and use it in GitHub Desktop.
Sample Angara notebook
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
{ | |
"cells": [ | |
{ | |
"cell_type": "markdown", | |
"metadata": {}, | |
"source": [ | |
"MH vs Filzbach\n", | |
"------\n", | |
"Compares sampling distributions of Filzbach and raw Metropolis-Hastings implementation for Gamma(2,1) distribution." | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": null, | |
"metadata": { | |
"collapsed": false | |
}, | |
"outputs": [], | |
"source": [ | |
"#N \"Angara.Base\"\n", | |
"#N \"Angara.Chart\"\n", | |
"#N \"Angara.Html\"\n", | |
"#N \"Angara.Table\"\n", | |
"#N \"Angara.Reinstate\"\n", | |
"#N \"Angara.Serialization\"\n", | |
"#N \"Angara.Serialization.Json\"\n", | |
"#N \"Angara.Statistics\"\n", | |
"#N \"DotNetZip\"\n", | |
"#N \"Newtonsoft.Json\"\n", | |
"open Angara.Statistics\n", | |
"open Angara.Charting\n", | |
"Angara.Base.Init()\n", | |
"let angaraDisplay (artefact:obj) =\n", | |
" Angara.Html.MakeEmbeddable \"600px\" artefact\n", | |
" |> Util.Html |> Display" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 2, | |
"metadata": { | |
"collapsed": false | |
}, | |
"outputs": [], | |
"source": [ | |
"let hist count samples = \n", | |
" if count < 1 then invalidArg \"count\" \"must be 1 or greater\"\n", | |
" let s = summary samples\n", | |
" let xmin = s.min\n", | |
" let xmax = s.max\n", | |
" let dx = (xmax - xmin) / (float count)\n", | |
" let h = Array.zeroCreate count\n", | |
" if dx > 0. then \n", | |
" let mutable scount = 0\n", | |
" for x in samples do\n", | |
" let idx = min 100 (max 1 (int (ceil ((x - xmin) / dx)))) - 1\n", | |
" h.[idx] <- h.[idx] + 1.\n", | |
" scount <- scount + 1\n", | |
" if scount > 0 then \n", | |
" for i in 0..count - 1 do\n", | |
" h.[i] <- h.[i] / dx / (float scount)\n", | |
" xmin, xmax, h" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 3, | |
"metadata": { | |
"collapsed": true | |
}, | |
"outputs": [], | |
"source": [ | |
"let hlines (xmin:float, xmax:float, h:float[]) = \n", | |
" let c = Array.length h\n", | |
" let dx = (xmax - xmin) / (float c)\n", | |
" let x = Array.init (c * 2 + 1) (fun i -> xmin + dx * (float (i / 2)))\n", | |
"\n", | |
" let y = \n", | |
" Array.init (c * 2 + 1) (fun i -> \n", | |
" if i = 0 then 0.\n", | |
" else float h.[(i - 1) / 2])\n", | |
" x, y" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 4, | |
"metadata": { | |
"collapsed": false | |
}, | |
"outputs": [], | |
"source": [ | |
"let r = MT19937()\n", | |
"let gamma x = log_pdf (Gamma(2., 1.)) x\n", | |
"let s = \n", | |
" Seq.unfold (fun (i, x) -> \n", | |
" if i < 1000000 then \n", | |
" let step = r.normal()\n", | |
" let xs = x + step\n", | |
" let A = gamma xs - gamma x\n", | |
"\n", | |
" let x = \n", | |
" if log (1. - r.uniform_float64()) < A then xs\n", | |
" else x\n", | |
" Some(x,(i + 1, x))\n", | |
" else None) (1, 1.)\n", | |
" |> Seq.toArray" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 5, | |
"metadata": { | |
"collapsed": false, | |
"scrolled": false | |
}, | |
"outputs": [], | |
"source": [ | |
"let hx, hy = \n", | |
" s\n", | |
" |> hist 100\n", | |
" |> hlines\n", | |
"\n", | |
"let chart = \n", | |
" Chart.ofList [ Plot.line (hx, hx |> Array.map (gamma>>exp), stroke = \"red\", displayName = \"Gamma(2,1)\")\n", | |
" Plot.line (hx, hy, thickness = 1.5, displayName = \"Metropolis\") ]" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 6, | |
"metadata": { | |
"collapsed": false | |
}, | |
"outputs": [ | |
{ | |
"data": { | |
"text/html": [ | |
"<h2 style='text-align:center'>Simple MH over analytic Г(2,1)</h2><hr />" | |
] | |
}, | |
"metadata": {}, | |
"output_type": "display_data" | |
}, | |
{ | |
"data": { | |
"text/html": [ | |
"<section>\r\n", | |
" <link href=\"https://cdn.rawgit.com/predictionmachines/Angara/0.2.2/dist/Angara.Html.css\" rel=\"stylesheet\" />\r\n", | |
" <script src=\"https://cdnjs.cloudflare.com/ajax/libs/require.js/2.2.0/require.min.js\"></script>\r\n", | |
" <script>\r\n", | |
" var content_a84758a236e240e28e9f1100c96b62f6 = {\r\n", | |
" \":Chart\": [\r\n", | |
" {\r\n", | |
" \"displayName\": \"Gamma(2,1)\",\r\n", | |
" \"kind\": \"line\",\r\n", | |
" \"properties\": {\r\n", | |
" \"fill_68\": \"#1F497D\",\r\n", | |
" \"fill_95\": \"#1F497D\",\r\n", | |
" \"stroke\": \"red\",\r\n", | |
" \"thickness\": 1.0,\r\n", | |
" \"treatAs\": \"0\",\r\n", | |
" \"x:double array\": \"AJ38QccOaD8AnfxBxw5oPxk3YGXQOcI/GTdgZdA5wj/fPdzWsgnSP9893NayCdI/MmAIe3322j8yYAh7ffbaP0JBmg+k8eE/QkGaD6Tx4T9rUrBhCWjmP2tSsGEJaOY/lWPGs27e6j+VY8azbt7qP7503AXUVO8/vnTcBdRU7z/0QvmrnOXxP/RC+auc5fE/iEsEVc8g9D+ISwRVzyD0PxxUD/4BXPY/HFQP/gFc9j+yXBqnNJf4P7JcGqc0l/g/RmUlUGfS+j9GZSVQZ9L6P9ptMPmZDf0/2m0w+ZkN/T9wdjuizEj/P3B2O6LMSP8/gj+jpf/BAECCP6Ol/8EAQMzDKPqY3wFAzMMo+pjfAUAWSK5OMv0CQBZIrk4y/QJAYcwzo8saBEBhzDOjyxoEQKtQufdkOAVAq1C592Q4BUD11D5M/lUGQPXUPkz+VQZAQFnEoJdzB0BAWcSgl3MHQIrdSfUwkQhAit1J9TCRCEDUYc9Jyq4JQNRhz0nKrglAH+ZUnmPMCkAf5lSeY8wKQGlq2vL86QtAaWra8vzpC0Cz7l9HlgcNQLPuX0eWBw1A/XLlmy8lDkD9cuWbLyUOQEj3avDIQg9ASPdq8MhCD0DKPXgiMTAQQMo9eCIxMBBA7/+6zP2+EEDv/7rM/b4QQBTC/XbKTRFAFML9dspNEUA5hEAhl9wRQDmEQCGX3BFAXkaDy2NrEkBeRoPLY2sSQIMIxnUw+hJAgwjGdTD6EkCoyggg/YgTQKjKCCD9iBNAzoxLyskXFEDOjEvKyRcUQPNOjnSWphRA806OdJamFEAYEdEeYzUVQBgR0R5jNRVAPdMTyS/EFUA90xPJL8QVQGKVVnP8UhZAYpVWc/xSFkCHV5kdyeEWQIdXmR3J4RZArRncx5VwF0CtGdzHlXAXQNLbHnJi/xdA0tsecmL/F0D3nWEcL44YQPedYRwvjhhAHGCkxvscGUAcYKTG+xwZQEEi53DIqxlAQSLncMirGUBm5CkblToaQGbkKRuVOhpAjKZsxWHJGkCMpmzFYckaQLFor28uWBtAsWivby5YG0DWKvIZ++YbQNYq8hn75htA++w0xMd1HED77DTEx3UcQCCvd26UBB1AIK93bpQEHUBFcboYYZMdQEVxuhhhkx1AajP9wi0iHkBqM/3CLSIeQJD1P236sB5AkPU/bfqwHkC1t4IXxz8fQLW3ghfHPx9A2nnFwZPOH0DaecXBk84fQAAeBDawLiBAAB4ENrAuIEASfyWLFnYgQBJ/JYsWdiBAJeBG4Hy9IEAl4EbgfL0gQDdBaDXjBCFAN0FoNeMEIUBKoomKSUwhQEqiiYpJTCFAXAOr36+TIUBcA6vfr5MhQG9kzDQW2yFAb2TMNBbbIUCCxe2JfCIiQILF7Yl8IiJAlCYP3+JpIkCUJg/f4mkiQKeHMDRJsSJAp4cwNEmxIkC56FGJr/giQLnoUYmv+CJAzElz3hVAI0DMSXPeFUAjQN6qlDN8hyNA3qqUM3yHI0DxC7aI4s4jQPELtojiziNABG3X3UgWJEAEbdfdSBYkQBbO+DKvXSRAFs74Mq9dJEApLxqIFaUkQCkvGogVpSRAO5A73XvsJEA7kDvde+wkQE7xXDLiMyVATvFcMuIzJUBhUn6HSHslQGFSfodIeyVAc7Of3K7CJUBzs5/crsIlQIYUwTEVCiZAhhTBMRUKJkCYdeKGe1EmQJh14oZ7USZAq9YD3OGYJkCr1gPc4ZgmQL03JTFI4CZAvTclMUjgJkDQmEaGricnQNCYRoauJydA4/ln2xRvJ0Dj+WfbFG8nQPVaiTB7tidA9VqJMHu2J0AIvKqF4f0nQAi8qoXh/SdAGh3M2kdFKEAaHczaR0UoQC1+7S+ujChALX7tL66MKEA/3w6FFNQoQD/fDoUU1ChAUkAw2nobKUBSQDDaehspQGWhUS/hYilAZaFRL+FiKUB3AnOER6opQHcCc4RHqilAimOU2a3xKUCKY5TZrfEpQJzEtS4UOSpAnMS1LhQ5KkCvJdeDeoAqQK8l14N6gCpAwob42ODHKkDChvjY4McqQNTnGS5HDytA1OcZLkcPK0DnSDuDrVYrQOdIO4OtVitA+alc2BOeK0D5qVzYE54rQAwLfi165StA\",\r\n", | |
" \"y:double array\": \"f2QS27f8Zz9/ZBLbt/xnP/Fsq9Aqnb8/8Wyr0Cqdvz8wL0E7HjfLPzAvQTseN8s/UOIaTmSx0T9Q4hpOZLHRPw/yaGHwe9Q/D/JoYfB71D/l5RqXuT/WP+XlGpe5P9Y/dDVeUAw11z90NV5QDDXXP1OoOqUAitc/U6g6pQCK1z9p0C1DLmTXP2nQLUMuZNc/BEPwFRvi1j8EQ/AVG+LWP96dWFxvHNY/3p1YXG8c1j9AOuYy+CbVP0A65jL4JtU/dbVHQYAR1D91tUdBgBHUP2ZLbw+F6NI/ZktvD4Xo0j9exfmJzrXRP17F+YnOtdE/BSBqZu2A0D8FIGpm7YDQP9A1DcFIn84/0DUNwUifzj/LepJcf0zMP8t6klx/TMw/SzmjBrsPyj9LOaMGuw/KP1IHURhU7cc/UgdRGFTtxz+t0ADUIejFP63QANQh6MU/kR3azMsBxD+RHdrMywHEP0RnWaMMO8I/RGdZoww7wj+2CXKH6JPAP7YJcofok8A/62x4G7MXvj/rbHgbsxe+P6vo0CLnQ7s/q+jQIudDuz85tKwrAqq4Pzm0rCsCqrg/i4WNvTRHtj+LhY29NEe2P+3nFodvGLQ/7ecWh28YtD9H78wtgBqyP0fvzC2AGrI/50GesidKsD/nQZ6yJ0qwP4+yETRXSK0/j7IRNFdIrT9x2FanxUqqP3HYVqfFSqo/raV5ZH6Vpz+tpXlkfpWnP3mZ2kGoIqU/eZnaQagipT/dayRKtuyiP91rJEq27KI/KeyGQ23uoD8p7IZDbe6gPyAO+PfMRZ4/IA7498xFnj+xMSXxIAubP7ExJfEgC5s/9kTjEFwkmD/2ROMQXCSYP4PH30eoiZU/g8ffR6iJlT/6KdrXzTOTP/op2tfNM5M/vT7PGSwckT+9Ps8ZLByRP0HsC+5ieY4/QewL7mJ5jj+DO3bTpR+LP4M7dtOlH4s/bP0njwYhiD9s/SePBiGIPzrQoRFVdIU/OtChEVV0hT+loL1lOhGDP6WgvWU6EYM/Y77OBSjwgD9jvs4FKPCAP+H9xrOPFH4/4f3Gs48Ufj8HDdj22bJ6PwcN2PbZsno/jx+hoAqwdz+PH6GgCrB3P5wZ7g4YAnU/nBnuDhgCdT/XRe70+Z9yP9dF7vT5n3I/fBTg5pGBcD98FODmkYFwP4m7XkkrP20/ibteSSs/bT9eps0p9uZpP16mzSn25mk/A2Qiu8vuZj8DZCK7y+5mP2gwOgwtTGQ/aDA6DC1MZD+QIYsXtvVhP5Ahixe29WE/6u3ORAPGXz/q7c5EA8ZfP2l5XUcfGVw/aXldRx8ZXD/JZcMGXNdYP8llwwZc11g/fzNnC8f0VT9/M2cLx/RVP9A4rWu5ZlM/0Dita7lmUz8w93dNtSNRPzD3d021I1E/XFrViI1GTj9cWtWIjUZOPw5/7YTPu0o/Dn/thM+7Sj/FyROqzZlHP8XJE6rNmUc/bIsG3KDURD9siwbcoNREP7fOpvi0YUI/t86m+LRhQj9cZos5pDdAP1xmizmkN0A/OMl8vyycPD84yXy/LJw8P9rZ0pdGOzk/2tnSl0Y7OT8nLFhacT82PycsWFpxPzY/4Eak8ASdMz/gRqTwBJ0zP04MN8CrSTE/Tgw3wKtJMT9WRhk3enguP1ZGGTd6eC4/Qk82fDfZKj9CTzZ8N9kqPwxHbdEvpyc/DEdt0S+nJz99BhN8qdUkP30GE3yp1SQ/HR5QEGFZIj8dHlAQYVkiP8iAoTFfKCA/yIChMV8oID8Tptn2pXMcPxOm2falcxw/GE0kJeELGT8YTSQl4QsZP4MLkNumCxY/gwuQ26YLFj+eZXzcx2YTP55lfNzHZhM/vg8v2X8SET++Dy/ZfxIRPzRpbJ+XCg4/NGlsn5cKDj+PFHgyim0KP48UeDKKbQo/6pq6SQE/Bz/qmrpJAT8HP5LeK9LjcQQ/kt4r0uNxBD8RcN9EovoBPxFw30Si+gE/jKO8BRKe/z6Mo7wFEp7/PrztiM8vzPs+vO2Izy/M+z6MJjz7vG/4PowmPPu8b/g+3JSfmbh69T7clJ+ZuHr1PrLI2qLJ4PI+ssjaosng8j7roeBYDZfwPuuh4FgNl/A+Mayly9Yn7T4xrKXL1iftPswo0B7enek+\"\r\n", | |
" },\r\n", | |
" \"titles\": {}\r\n", | |
" },\r\n", | |
" {\r\n", | |
" \"displayName\": \"Metropolis\",\r\n", | |
" \"kind\": \"line\",\r\n", | |
" \"properties\": {\r\n", | |
" \"fill_68\": \"#1F497D\",\r\n", | |
" \"fill_95\": \"#1F497D\",\r\n", | |
" \"stroke\": \"#1F497D\",\r\n", | |
" \"thickness\": 1.5,\r\n", | |
" \"treatAs\": \"0\",\r\n", | |
" \"x:double array\": \"AJ38QccOaD8AnfxBxw5oPxk3YGXQOcI/GTdgZdA5wj/fPdzWsgnSP9893NayCdI/MmAIe3322j8yYAh7ffbaP0JBmg+k8eE/QkGaD6Tx4T9rUrBhCWjmP2tSsGEJaOY/lWPGs27e6j+VY8azbt7qP7503AXUVO8/vnTcBdRU7z/0QvmrnOXxP/RC+auc5fE/iEsEVc8g9D+ISwRVzyD0PxxUD/4BXPY/HFQP/gFc9j+yXBqnNJf4P7JcGqc0l/g/RmUlUGfS+j9GZSVQZ9L6P9ptMPmZDf0/2m0w+ZkN/T9wdjuizEj/P3B2O6LMSP8/gj+jpf/BAECCP6Ol/8EAQMzDKPqY3wFAzMMo+pjfAUAWSK5OMv0CQBZIrk4y/QJAYcwzo8saBEBhzDOjyxoEQKtQufdkOAVAq1C592Q4BUD11D5M/lUGQPXUPkz+VQZAQFnEoJdzB0BAWcSgl3MHQIrdSfUwkQhAit1J9TCRCEDUYc9Jyq4JQNRhz0nKrglAH+ZUnmPMCkAf5lSeY8wKQGlq2vL86QtAaWra8vzpC0Cz7l9HlgcNQLPuX0eWBw1A/XLlmy8lDkD9cuWbLyUOQEj3avDIQg9ASPdq8MhCD0DKPXgiMTAQQMo9eCIxMBBA7/+6zP2+EEDv/7rM/b4QQBTC/XbKTRFAFML9dspNEUA5hEAhl9wRQDmEQCGX3BFAXkaDy2NrEkBeRoPLY2sSQIMIxnUw+hJAgwjGdTD6EkCoyggg/YgTQKjKCCD9iBNAzoxLyskXFEDOjEvKyRcUQPNOjnSWphRA806OdJamFEAYEdEeYzUVQBgR0R5jNRVAPdMTyS/EFUA90xPJL8QVQGKVVnP8UhZAYpVWc/xSFkCHV5kdyeEWQIdXmR3J4RZArRncx5VwF0CtGdzHlXAXQNLbHnJi/xdA0tsecmL/F0D3nWEcL44YQPedYRwvjhhAHGCkxvscGUAcYKTG+xwZQEEi53DIqxlAQSLncMirGUBm5CkblToaQGbkKRuVOhpAjKZsxWHJGkCMpmzFYckaQLFor28uWBtAsWivby5YG0DWKvIZ++YbQNYq8hn75htA++w0xMd1HED77DTEx3UcQCCvd26UBB1AIK93bpQEHUBFcboYYZMdQEVxuhhhkx1AajP9wi0iHkBqM/3CLSIeQJD1P236sB5AkPU/bfqwHkC1t4IXxz8fQLW3ghfHPx9A2nnFwZPOH0DaecXBk84fQAAeBDawLiBAAB4ENrAuIEASfyWLFnYgQBJ/JYsWdiBAJeBG4Hy9IEAl4EbgfL0gQDdBaDXjBCFAN0FoNeMEIUBKoomKSUwhQEqiiYpJTCFAXAOr36+TIUBcA6vfr5MhQG9kzDQW2yFAb2TMNBbbIUCCxe2JfCIiQILF7Yl8IiJAlCYP3+JpIkCUJg/f4mkiQKeHMDRJsSJAp4cwNEmxIkC56FGJr/giQLnoUYmv+CJAzElz3hVAI0DMSXPeFUAjQN6qlDN8hyNA3qqUM3yHI0DxC7aI4s4jQPELtojiziNABG3X3UgWJEAEbdfdSBYkQBbO+DKvXSRAFs74Mq9dJEApLxqIFaUkQCkvGogVpSRAO5A73XvsJEA7kDvde+wkQE7xXDLiMyVATvFcMuIzJUBhUn6HSHslQGFSfodIeyVAc7Of3K7CJUBzs5/crsIlQIYUwTEVCiZAhhTBMRUKJkCYdeKGe1EmQJh14oZ7USZAq9YD3OGYJkCr1gPc4ZgmQL03JTFI4CZAvTclMUjgJkDQmEaGricnQNCYRoauJydA4/ln2xRvJ0Dj+WfbFG8nQPVaiTB7tidA9VqJMHu2J0AIvKqF4f0nQAi8qoXh/SdAGh3M2kdFKEAaHczaR0UoQC1+7S+ujChALX7tL66MKEA/3w6FFNQoQD/fDoUU1ChAUkAw2nobKUBSQDDaehspQGWhUS/hYilAZaFRL+FiKUB3AnOER6opQHcCc4RHqilAimOU2a3xKUCKY5TZrfEpQJzEtS4UOSpAnMS1LhQ5KkCvJdeDeoAqQK8l14N6gCpAwob42ODHKkDChvjY4McqQNTnGS5HDytA1OcZLkcPK0DnSDuDrVYrQOdIO4OtVitA+alc2BOeK0D5qVzYE54rQAwLfi165StA\",\r\n", | |
" \"y:double array\": \"AAAAAAAAAADgHhsJzd+wP+AeGwnN37A/DwLEGj8Kxj8PAsQaPwrGP40p9n8/NM8/jSn2fz80zz8Q4iFgv2fTPxDiIWC/Z9M/GGtWsxuX1T8Ya1azG5fVP47ssg820tY/juyyDzbS1j9AVucU30rXP0BW5xTfStc/aAwDF7971z9oDAMXv3vXP7qOiYwtbdc/uo6JjC1t1z/REX9uRH3WP9ERf25EfdY/u29jnvh51T+7b2Oe+HnVP+kt/vZ0ldQ/6S3+9nSV1D8uo4fEoJTTPy6jh8SglNM/7uq4fH9t0j/u6rh8f23SPzM0trp+6NA/MzS2un7o0D8wICH92dbPPzAgIf3Z1s8/TfRhbYmFzT9N9GFtiYXNP1ytw8ezAcs/XK3Dx7MByz9J5gyfXBrJP0nmDJ9cGsk/WSuS/qgNxz9ZK5L+qA3HP5lxCYSX38Q/mXEJhJffxD8UzWAOjCLDPxTNYA6MIsM/P7WPXZ9jwT8/tY9dn2PBP7K6ZNbj5b8/srpk1uPlvz9dLIEr3ue8P10sgSve57w/bwOf+iJDuj9vA5/6IkO6P2SdDWc8mbc/ZJ0NZzyZtz8BcRBxRvC0PwFxEHFG8LQ/eaY+g7QEsz95pj6DtASzP5G6UAh5zbA/kbpQCHnNsD8VB6N5n6GuPxUHo3mfoa4/cQyIytlBqz9xDIjK2UGrP/TAl+bScqg/9MCX5tJyqD8tR+CA1l6mPy1H4IDWXqY/AX/y3ogwpD8Bf/LeiDCkP+RMsYS/XaE/5EyxhL9doT+QR7L0WROfP5BHsvRZE58/zgKCQXEdnT/OAoJBcR2dP5vsSRaBM5k/m+xJFoEzmT9hhcZc+KqWP2GFxlz4qpY/Y6h7b57LlD9jqHtvnsuUP5PLptubG5I/k8um25sbkj9iTi3SI4SPP2JOLdIjhI8/MTzl2F5Xjj8xPOXYXleOP3hrtca0IIk/eGu1xrQgiT9hftUl1wqHP2F+1SXXCoc/7aHsRIO5gj/toexEg7mCP5PLptubG4I/k8um25sbgj86A0Pg9r2APzoDQ+D2vYA/3rbq/Q9YfD/etur9D1h8PwV+6XIorHg/BX7pciiseD8BZSacmQ53PwFlJpyZDnc/sRe0eFQQcj+xF7R4VBByPzVVqK/t9G8/NVWor+30bz84jTBn9/VsPziNMGf39Ww/+dhRUQGMbj/52FFRAYxuPyX0nFmoYWY/JfScWahhZj8N6/jcRctlPw3r+NxFy2U//E5kygBUYj/8TmTKAFRiP9MwGL1jm10/0zAYvWObXT/N+5AKUH1dP837kApQfV0/xcwLYiglVz/FzAtiKCVXP63DZ+XFjlY/rcNn5cWOVj8VWAhHY+pSPxVYCEdj6lI/pJTiPJ42UD+klOI8njZQP/BuQ+zZT04/8G5D7NlPTj8EQ2C2KMhOPwRDYLYoyE4/hOCklCi6Sz+E4KSUKLpLP6NZWYCeUkY/o1lZgJ5SRj8ZjY/5dghDPxmNj/l2CEM/GY2P+XYIQz8ZjY/5dghDP6yI3ijGqjw/rIjeKMaqPD8NrW4bUAQ/Pw2tbhtQBD8/N8u6KO28Mz83y7oo7bwzP9RroXljfzc/1GuheWN/Nz/B0g1sFOswP8HSDWwU6zA/v1z78hQjLT+/XPvyFCMtP4UbLlEoniU/hRsuUSieJT9ec/S8iq0UP15z9LyKrRQ/ECOBlE/MIj8QI4GUT8wiPzVVqK/t9B8/NVWor+30Hz8QI4GUT8wSPxAjgZRPzBI/XnP0vIqtFD9ec/S8iq0UP8HSDWwU6xA/wdINbBTrED/B0g1sFOsQP8HSDWwU6xA/rcNn5cWOBj+tw2flxY4GP8HSDWwU6xA/wdINbBTrED/B0g1sFOsQP8HSDWwU6xA/5gQ1h7IT7j7mBDWHshPuPuYENYeyE94+5gQ1h7IT3j7mBDWHshPuPuYENYeyE+4+5gQ1h7IT7j7mBDWHshPuPuYENYeyE/4+5gQ1h7IT/j7mBDWHshPuPuYENYeyE+4+rcNn5cWOBj+tw2flxY4GPwAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAOYENYeyE94+5gQ1h7IT3j7mBDWHshPePuYENYeyE94+\"\r\n", | |
" },\r\n", | |
" \"titles\": {}\r\n", | |
" }\r\n", | |
" ]\r\n", | |
"};\r\n", | |
" require.config({\r\n", | |
" baseUrl: \"https://cdn.rawgit.com/predictionmachines/Angara/0.2.2/dist\",\r\n", | |
" paths: {\r\n", | |
" \"angara.serializationjs\": \"https://cdn.rawgit.com/predictionmachines/Angara.Serialization/v0.2.0/dist/Angara.Serialization.umd\",\r\n", | |
" \"Chart\": \"https://cdn.rawgit.com/predictionmachines/Angara.Chart/0.2.0/dist/Chart\",\r\n", | |
" \"Table\": \"https://cdn.rawgit.com/predictionmachines/Angara.Table/0.2.1/dist/Table\",\r\n", | |
" \"angara.tablejs\": \"https://cdn.rawgit.com/predictionmachines/Angara.Table/0.2.1/dist/angara.tablejs\",\r\n", | |
" \"jquery.dataTables\": \"https://cdnjs.cloudflare.com/ajax/libs/datatables/1.10.11/js/jquery.dataTables.min\",\r\n", | |
" \"jquery.dataTables-css\": \"https://cdn.rawgit.com/predictionmachines/Angara.Table/0.2.1/dist/angara.tablejs\",\r\n", | |
" \"chartViewer.umd\": \"https://cdn.rawgit.com/predictionmachines/InteractiveDataDisplay/v1.2.1/dist/chartViewer.umd\",\r\n", | |
" \"chartViewer-css\": \"https://cdn.rawgit.com/predictionmachines/InteractiveDataDisplay/v1.2.1/dist/chartViewer\",\r\n", | |
" \"idd\": \"https://cdn.rawgit.com/predictionmachines/InteractiveDataDisplay/v1.2.1/dist/idd\",\r\n", | |
" \"idd-css\": \"https://cdn.rawgit.com/predictionmachines/InteractiveDataDisplay/v1.2.1/dist/idd\",\r\n", | |
" \"domReady\": \"https://cdnjs.cloudflare.com/ajax/libs/require-domReady/2.0.1/domReady.min\",\r\n", | |
" \"jquery\": \"https://cdnjs.cloudflare.com/ajax/libs/jquery/2.2.2/jquery.min\",\r\n", | |
" \"jquery-ui\": \"https://cdnjs.cloudflare.com/ajax/libs/jqueryui/1.11.4/jquery-ui\",\r\n", | |
" \"css\": \"https://cdnjs.cloudflare.com/ajax/libs/require-css/0.1.8/css.min\",\r\n", | |
" \"rx\": \"https://cdnjs.cloudflare.com/ajax/libs/rxjs/4.1.0/rx.lite.min\"\r\n", | |
" }\r\n", | |
" });\r\n", | |
" require(['Angara.Show', 'angara.serializationjs', 'domReady!'], function (Angara, Marshalling, document) {\r\n", | |
" var viewer = document.getElementById('container_a84758a236e240e28e9f1100c96b62f6');\r\n", | |
" var content_is = Marshalling.InfoSet.Unmarshal(content_a84758a236e240e28e9f1100c96b62f6);\r\n", | |
" var content = Marshalling.InfoSet.Deserialize(content_is);\r\n", | |
" return Angara.Show(content, viewer);\r\n", | |
" });\r\n", | |
" </script>\r\n", | |
" <div style=\"height:600px\" id=\"container_a84758a236e240e28e9f1100c96b62f6\" class=\"viewer-container\">\r\n", | |
" </div>\r\n", | |
"</section>\r\n" | |
] | |
}, | |
"metadata": {}, | |
"output_type": "display_data" | |
} | |
], | |
"source": [ | |
"Util.Html \"<h2 style='text-align:center'>Simple MH over analytic Г(2,1)</h2><hr />\" |> Display\n", | |
"angaraDisplay chart" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 7, | |
"metadata": { | |
"collapsed": false | |
}, | |
"outputs": [], | |
"source": [ | |
"open Angara.Filzbach\n", | |
"let s1 = \n", | |
" let p = Parameters.Empty.Add(\"x\", Uniform(0.001, 100.))\n", | |
"\n", | |
" let logl (p : Parameters) = \n", | |
" let x = p.GetValue 0\n", | |
" gamma x\n", | |
"\n", | |
" let out = Sampler.runmcmc (p, logl, 1000, 1000000, 1)\n", | |
" [| for { values = v } in out.samples -> v.[0] |]\n", | |
"\n", | |
"let hx1, hy1 = \n", | |
" s1\n", | |
" |> hist 100\n", | |
" |> hlines\n", | |
"\n", | |
"let c1 = \n", | |
" Chart.ofList [ Plot.line (hx, hx |> Array.map (gamma>>exp), stroke = \"red\", displayName = \"Gamma(2,1)\")\n", | |
" Plot.line (hx1, hy1, thickness = 1.5, displayName = \"Filzbach\") ]\n", | |
"\n" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 8, | |
"metadata": { | |
"collapsed": false | |
}, | |
"outputs": [ | |
{ | |
"data": { | |
"text/html": [ | |
"<h2 style='text-align:center'>Filzbach over analytic Г(2,1)</h2>" | |
] | |
}, | |
"metadata": {}, | |
"output_type": "display_data" | |
}, | |
{ | |
"data": { | |
"text/html": [ | |
"<section>\r\n", | |
" <link href=\"https://cdn.rawgit.com/predictionmachines/Angara/0.2.2/dist/Angara.Html.css\" rel=\"stylesheet\" />\r\n", | |
" <script src=\"https://cdnjs.cloudflare.com/ajax/libs/require.js/2.2.0/require.min.js\"></script>\r\n", | |
" <script>\r\n", | |
" var content_11100160ce744343a7e1af6423d50a06 = {\r\n", | |
" \":Chart\": [\r\n", | |
" {\r\n", | |
" \"displayName\": \"Gamma(2,1)\",\r\n", | |
" \"kind\": \"line\",\r\n", | |
" \"properties\": {\r\n", | |
" \"fill_68\": \"#1F497D\",\r\n", | |
" \"fill_95\": \"#1F497D\",\r\n", | |
" \"stroke\": \"red\",\r\n", | |
" \"thickness\": 1.0,\r\n", | |
" \"treatAs\": \"0\",\r\n", | |
" \"x:double array\": \"AJ38QccOaD8AnfxBxw5oPxk3YGXQOcI/GTdgZdA5wj/fPdzWsgnSP9893NayCdI/MmAIe3322j8yYAh7ffbaP0JBmg+k8eE/QkGaD6Tx4T9rUrBhCWjmP2tSsGEJaOY/lWPGs27e6j+VY8azbt7qP7503AXUVO8/vnTcBdRU7z/0QvmrnOXxP/RC+auc5fE/iEsEVc8g9D+ISwRVzyD0PxxUD/4BXPY/HFQP/gFc9j+yXBqnNJf4P7JcGqc0l/g/RmUlUGfS+j9GZSVQZ9L6P9ptMPmZDf0/2m0w+ZkN/T9wdjuizEj/P3B2O6LMSP8/gj+jpf/BAECCP6Ol/8EAQMzDKPqY3wFAzMMo+pjfAUAWSK5OMv0CQBZIrk4y/QJAYcwzo8saBEBhzDOjyxoEQKtQufdkOAVAq1C592Q4BUD11D5M/lUGQPXUPkz+VQZAQFnEoJdzB0BAWcSgl3MHQIrdSfUwkQhAit1J9TCRCEDUYc9Jyq4JQNRhz0nKrglAH+ZUnmPMCkAf5lSeY8wKQGlq2vL86QtAaWra8vzpC0Cz7l9HlgcNQLPuX0eWBw1A/XLlmy8lDkD9cuWbLyUOQEj3avDIQg9ASPdq8MhCD0DKPXgiMTAQQMo9eCIxMBBA7/+6zP2+EEDv/7rM/b4QQBTC/XbKTRFAFML9dspNEUA5hEAhl9wRQDmEQCGX3BFAXkaDy2NrEkBeRoPLY2sSQIMIxnUw+hJAgwjGdTD6EkCoyggg/YgTQKjKCCD9iBNAzoxLyskXFEDOjEvKyRcUQPNOjnSWphRA806OdJamFEAYEdEeYzUVQBgR0R5jNRVAPdMTyS/EFUA90xPJL8QVQGKVVnP8UhZAYpVWc/xSFkCHV5kdyeEWQIdXmR3J4RZArRncx5VwF0CtGdzHlXAXQNLbHnJi/xdA0tsecmL/F0D3nWEcL44YQPedYRwvjhhAHGCkxvscGUAcYKTG+xwZQEEi53DIqxlAQSLncMirGUBm5CkblToaQGbkKRuVOhpAjKZsxWHJGkCMpmzFYckaQLFor28uWBtAsWivby5YG0DWKvIZ++YbQNYq8hn75htA++w0xMd1HED77DTEx3UcQCCvd26UBB1AIK93bpQEHUBFcboYYZMdQEVxuhhhkx1AajP9wi0iHkBqM/3CLSIeQJD1P236sB5AkPU/bfqwHkC1t4IXxz8fQLW3ghfHPx9A2nnFwZPOH0DaecXBk84fQAAeBDawLiBAAB4ENrAuIEASfyWLFnYgQBJ/JYsWdiBAJeBG4Hy9IEAl4EbgfL0gQDdBaDXjBCFAN0FoNeMEIUBKoomKSUwhQEqiiYpJTCFAXAOr36+TIUBcA6vfr5MhQG9kzDQW2yFAb2TMNBbbIUCCxe2JfCIiQILF7Yl8IiJAlCYP3+JpIkCUJg/f4mkiQKeHMDRJsSJAp4cwNEmxIkC56FGJr/giQLnoUYmv+CJAzElz3hVAI0DMSXPeFUAjQN6qlDN8hyNA3qqUM3yHI0DxC7aI4s4jQPELtojiziNABG3X3UgWJEAEbdfdSBYkQBbO+DKvXSRAFs74Mq9dJEApLxqIFaUkQCkvGogVpSRAO5A73XvsJEA7kDvde+wkQE7xXDLiMyVATvFcMuIzJUBhUn6HSHslQGFSfodIeyVAc7Of3K7CJUBzs5/crsIlQIYUwTEVCiZAhhTBMRUKJkCYdeKGe1EmQJh14oZ7USZAq9YD3OGYJkCr1gPc4ZgmQL03JTFI4CZAvTclMUjgJkDQmEaGricnQNCYRoauJydA4/ln2xRvJ0Dj+WfbFG8nQPVaiTB7tidA9VqJMHu2J0AIvKqF4f0nQAi8qoXh/SdAGh3M2kdFKEAaHczaR0UoQC1+7S+ujChALX7tL66MKEA/3w6FFNQoQD/fDoUU1ChAUkAw2nobKUBSQDDaehspQGWhUS/hYilAZaFRL+FiKUB3AnOER6opQHcCc4RHqilAimOU2a3xKUCKY5TZrfEpQJzEtS4UOSpAnMS1LhQ5KkCvJdeDeoAqQK8l14N6gCpAwob42ODHKkDChvjY4McqQNTnGS5HDytA1OcZLkcPK0DnSDuDrVYrQOdIO4OtVitA+alc2BOeK0D5qVzYE54rQAwLfi165StA\",\r\n", | |
" \"y:double array\": \"f2QS27f8Zz9/ZBLbt/xnP/Fsq9Aqnb8/8Wyr0Cqdvz8wL0E7HjfLPzAvQTseN8s/UOIaTmSx0T9Q4hpOZLHRPw/yaGHwe9Q/D/JoYfB71D/l5RqXuT/WP+XlGpe5P9Y/dDVeUAw11z90NV5QDDXXP1OoOqUAitc/U6g6pQCK1z9p0C1DLmTXP2nQLUMuZNc/BEPwFRvi1j8EQ/AVG+LWP96dWFxvHNY/3p1YXG8c1j9AOuYy+CbVP0A65jL4JtU/dbVHQYAR1D91tUdBgBHUP2ZLbw+F6NI/ZktvD4Xo0j9exfmJzrXRP17F+YnOtdE/BSBqZu2A0D8FIGpm7YDQP9A1DcFIn84/0DUNwUifzj/LepJcf0zMP8t6klx/TMw/SzmjBrsPyj9LOaMGuw/KP1IHURhU7cc/UgdRGFTtxz+t0ADUIejFP63QANQh6MU/kR3azMsBxD+RHdrMywHEP0RnWaMMO8I/RGdZoww7wj+2CXKH6JPAP7YJcofok8A/62x4G7MXvj/rbHgbsxe+P6vo0CLnQ7s/q+jQIudDuz85tKwrAqq4Pzm0rCsCqrg/i4WNvTRHtj+LhY29NEe2P+3nFodvGLQ/7ecWh28YtD9H78wtgBqyP0fvzC2AGrI/50GesidKsD/nQZ6yJ0qwP4+yETRXSK0/j7IRNFdIrT9x2FanxUqqP3HYVqfFSqo/raV5ZH6Vpz+tpXlkfpWnP3mZ2kGoIqU/eZnaQagipT/dayRKtuyiP91rJEq27KI/KeyGQ23uoD8p7IZDbe6gPyAO+PfMRZ4/IA7498xFnj+xMSXxIAubP7ExJfEgC5s/9kTjEFwkmD/2ROMQXCSYP4PH30eoiZU/g8ffR6iJlT/6KdrXzTOTP/op2tfNM5M/vT7PGSwckT+9Ps8ZLByRP0HsC+5ieY4/QewL7mJ5jj+DO3bTpR+LP4M7dtOlH4s/bP0njwYhiD9s/SePBiGIPzrQoRFVdIU/OtChEVV0hT+loL1lOhGDP6WgvWU6EYM/Y77OBSjwgD9jvs4FKPCAP+H9xrOPFH4/4f3Gs48Ufj8HDdj22bJ6PwcN2PbZsno/jx+hoAqwdz+PH6GgCrB3P5wZ7g4YAnU/nBnuDhgCdT/XRe70+Z9yP9dF7vT5n3I/fBTg5pGBcD98FODmkYFwP4m7XkkrP20/ibteSSs/bT9eps0p9uZpP16mzSn25mk/A2Qiu8vuZj8DZCK7y+5mP2gwOgwtTGQ/aDA6DC1MZD+QIYsXtvVhP5Ahixe29WE/6u3ORAPGXz/q7c5EA8ZfP2l5XUcfGVw/aXldRx8ZXD/JZcMGXNdYP8llwwZc11g/fzNnC8f0VT9/M2cLx/RVP9A4rWu5ZlM/0Dita7lmUz8w93dNtSNRPzD3d021I1E/XFrViI1GTj9cWtWIjUZOPw5/7YTPu0o/Dn/thM+7Sj/FyROqzZlHP8XJE6rNmUc/bIsG3KDURD9siwbcoNREP7fOpvi0YUI/t86m+LRhQj9cZos5pDdAP1xmizmkN0A/OMl8vyycPD84yXy/LJw8P9rZ0pdGOzk/2tnSl0Y7OT8nLFhacT82PycsWFpxPzY/4Eak8ASdMz/gRqTwBJ0zP04MN8CrSTE/Tgw3wKtJMT9WRhk3enguP1ZGGTd6eC4/Qk82fDfZKj9CTzZ8N9kqPwxHbdEvpyc/DEdt0S+nJz99BhN8qdUkP30GE3yp1SQ/HR5QEGFZIj8dHlAQYVkiP8iAoTFfKCA/yIChMV8oID8Tptn2pXMcPxOm2falcxw/GE0kJeELGT8YTSQl4QsZP4MLkNumCxY/gwuQ26YLFj+eZXzcx2YTP55lfNzHZhM/vg8v2X8SET++Dy/ZfxIRPzRpbJ+XCg4/NGlsn5cKDj+PFHgyim0KP48UeDKKbQo/6pq6SQE/Bz/qmrpJAT8HP5LeK9LjcQQ/kt4r0uNxBD8RcN9EovoBPxFw30Si+gE/jKO8BRKe/z6Mo7wFEp7/PrztiM8vzPs+vO2Izy/M+z6MJjz7vG/4PowmPPu8b/g+3JSfmbh69T7clJ+ZuHr1PrLI2qLJ4PI+ssjaosng8j7roeBYDZfwPuuh4FgNl/A+Mayly9Yn7T4xrKXL1iftPswo0B7enek+\"\r\n", | |
" },\r\n", | |
" \"titles\": {}\r\n", | |
" },\r\n", | |
" {\r\n", | |
" \"displayName\": \"Metropolis\",\r\n", | |
" \"kind\": \"line\",\r\n", | |
" \"properties\": {\r\n", | |
" \"fill_68\": \"#1F497D\",\r\n", | |
" \"fill_95\": \"#1F497D\",\r\n", | |
" \"stroke\": \"#1F497D\",\r\n", | |
" \"thickness\": 1.5,\r\n", | |
" \"treatAs\": \"0\",\r\n", | |
" \"x:double array\": \"AJ38QccOaD8AnfxBxw5oPxk3YGXQOcI/GTdgZdA5wj/fPdzWsgnSP9893NayCdI/MmAIe3322j8yYAh7ffbaP0JBmg+k8eE/QkGaD6Tx4T9rUrBhCWjmP2tSsGEJaOY/lWPGs27e6j+VY8azbt7qP7503AXUVO8/vnTcBdRU7z/0QvmrnOXxP/RC+auc5fE/iEsEVc8g9D+ISwRVzyD0PxxUD/4BXPY/HFQP/gFc9j+yXBqnNJf4P7JcGqc0l/g/RmUlUGfS+j9GZSVQZ9L6P9ptMPmZDf0/2m0w+ZkN/T9wdjuizEj/P3B2O6LMSP8/gj+jpf/BAECCP6Ol/8EAQMzDKPqY3wFAzMMo+pjfAUAWSK5OMv0CQBZIrk4y/QJAYcwzo8saBEBhzDOjyxoEQKtQufdkOAVAq1C592Q4BUD11D5M/lUGQPXUPkz+VQZAQFnEoJdzB0BAWcSgl3MHQIrdSfUwkQhAit1J9TCRCEDUYc9Jyq4JQNRhz0nKrglAH+ZUnmPMCkAf5lSeY8wKQGlq2vL86QtAaWra8vzpC0Cz7l9HlgcNQLPuX0eWBw1A/XLlmy8lDkD9cuWbLyUOQEj3avDIQg9ASPdq8MhCD0DKPXgiMTAQQMo9eCIxMBBA7/+6zP2+EEDv/7rM/b4QQBTC/XbKTRFAFML9dspNEUA5hEAhl9wRQDmEQCGX3BFAXkaDy2NrEkBeRoPLY2sSQIMIxnUw+hJAgwjGdTD6EkCoyggg/YgTQKjKCCD9iBNAzoxLyskXFEDOjEvKyRcUQPNOjnSWphRA806OdJamFEAYEdEeYzUVQBgR0R5jNRVAPdMTyS/EFUA90xPJL8QVQGKVVnP8UhZAYpVWc/xSFkCHV5kdyeEWQIdXmR3J4RZArRncx5VwF0CtGdzHlXAXQNLbHnJi/xdA0tsecmL/F0D3nWEcL44YQPedYRwvjhhAHGCkxvscGUAcYKTG+xwZQEEi53DIqxlAQSLncMirGUBm5CkblToaQGbkKRuVOhpAjKZsxWHJGkCMpmzFYckaQLFor28uWBtAsWivby5YG0DWKvIZ++YbQNYq8hn75htA++w0xMd1HED77DTEx3UcQCCvd26UBB1AIK93bpQEHUBFcboYYZMdQEVxuhhhkx1AajP9wi0iHkBqM/3CLSIeQJD1P236sB5AkPU/bfqwHkC1t4IXxz8fQLW3ghfHPx9A2nnFwZPOH0DaecXBk84fQAAeBDawLiBAAB4ENrAuIEASfyWLFnYgQBJ/JYsWdiBAJeBG4Hy9IEAl4EbgfL0gQDdBaDXjBCFAN0FoNeMEIUBKoomKSUwhQEqiiYpJTCFAXAOr36+TIUBcA6vfr5MhQG9kzDQW2yFAb2TMNBbbIUCCxe2JfCIiQILF7Yl8IiJAlCYP3+JpIkCUJg/f4mkiQKeHMDRJsSJAp4cwNEmxIkC56FGJr/giQLnoUYmv+CJAzElz3hVAI0DMSXPeFUAjQN6qlDN8hyNA3qqUM3yHI0DxC7aI4s4jQPELtojiziNABG3X3UgWJEAEbdfdSBYkQBbO+DKvXSRAFs74Mq9dJEApLxqIFaUkQCkvGogVpSRAO5A73XvsJEA7kDvde+wkQE7xXDLiMyVATvFcMuIzJUBhUn6HSHslQGFSfodIeyVAc7Of3K7CJUBzs5/crsIlQIYUwTEVCiZAhhTBMRUKJkCYdeKGe1EmQJh14oZ7USZAq9YD3OGYJkCr1gPc4ZgmQL03JTFI4CZAvTclMUjgJkDQmEaGricnQNCYRoauJydA4/ln2xRvJ0Dj+WfbFG8nQPVaiTB7tidA9VqJMHu2J0AIvKqF4f0nQAi8qoXh/SdAGh3M2kdFKEAaHczaR0UoQC1+7S+ujChALX7tL66MKEA/3w6FFNQoQD/fDoUU1ChAUkAw2nobKUBSQDDaehspQGWhUS/hYilAZaFRL+FiKUB3AnOER6opQHcCc4RHqilAimOU2a3xKUCKY5TZrfEpQJzEtS4UOSpAnMS1LhQ5KkCvJdeDeoAqQK8l14N6gCpAwob42ODHKkDChvjY4McqQNTnGS5HDytA1OcZLkcPK0DnSDuDrVYrQOdIO4OtVitA+alc2BOeK0D5qVzYE54rQAwLfi165StA\",\r\n", | |
" \"y:double array\": \"AAAAAAAAAADgHhsJzd+wP+AeGwnN37A/DwLEGj8Kxj8PAsQaPwrGP40p9n8/NM8/jSn2fz80zz8Q4iFgv2fTPxDiIWC/Z9M/GGtWsxuX1T8Ya1azG5fVP47ssg820tY/juyyDzbS1j9AVucU30rXP0BW5xTfStc/aAwDF7971z9oDAMXv3vXP7qOiYwtbdc/uo6JjC1t1z/REX9uRH3WP9ERf25EfdY/u29jnvh51T+7b2Oe+HnVP+kt/vZ0ldQ/6S3+9nSV1D8uo4fEoJTTPy6jh8SglNM/7uq4fH9t0j/u6rh8f23SPzM0trp+6NA/MzS2un7o0D8wICH92dbPPzAgIf3Z1s8/TfRhbYmFzT9N9GFtiYXNP1ytw8ezAcs/XK3Dx7MByz9J5gyfXBrJP0nmDJ9cGsk/WSuS/qgNxz9ZK5L+qA3HP5lxCYSX38Q/mXEJhJffxD8UzWAOjCLDPxTNYA6MIsM/P7WPXZ9jwT8/tY9dn2PBP7K6ZNbj5b8/srpk1uPlvz9dLIEr3ue8P10sgSve57w/bwOf+iJDuj9vA5/6IkO6P2SdDWc8mbc/ZJ0NZzyZtz8BcRBxRvC0PwFxEHFG8LQ/eaY+g7QEsz95pj6DtASzP5G6UAh5zbA/kbpQCHnNsD8VB6N5n6GuPxUHo3mfoa4/cQyIytlBqz9xDIjK2UGrP/TAl+bScqg/9MCX5tJyqD8tR+CA1l6mPy1H4IDWXqY/AX/y3ogwpD8Bf/LeiDCkP+RMsYS/XaE/5EyxhL9doT+QR7L0WROfP5BHsvRZE58/zgKCQXEdnT/OAoJBcR2dP5vsSRaBM5k/m+xJFoEzmT9hhcZc+KqWP2GFxlz4qpY/Y6h7b57LlD9jqHtvnsuUP5PLptubG5I/k8um25sbkj9iTi3SI4SPP2JOLdIjhI8/MTzl2F5Xjj8xPOXYXleOP3hrtca0IIk/eGu1xrQgiT9hftUl1wqHP2F+1SXXCoc/7aHsRIO5gj/toexEg7mCP5PLptubG4I/k8um25sbgj86A0Pg9r2APzoDQ+D2vYA/3rbq/Q9YfD/etur9D1h8PwV+6XIorHg/BX7pciiseD8BZSacmQ53PwFlJpyZDnc/sRe0eFQQcj+xF7R4VBByPzVVqK/t9G8/NVWor+30bz84jTBn9/VsPziNMGf39Ww/+dhRUQGMbj/52FFRAYxuPyX0nFmoYWY/JfScWahhZj8N6/jcRctlPw3r+NxFy2U//E5kygBUYj/8TmTKAFRiP9MwGL1jm10/0zAYvWObXT/N+5AKUH1dP837kApQfV0/xcwLYiglVz/FzAtiKCVXP63DZ+XFjlY/rcNn5cWOVj8VWAhHY+pSPxVYCEdj6lI/pJTiPJ42UD+klOI8njZQP/BuQ+zZT04/8G5D7NlPTj8EQ2C2KMhOPwRDYLYoyE4/hOCklCi6Sz+E4KSUKLpLP6NZWYCeUkY/o1lZgJ5SRj8ZjY/5dghDPxmNj/l2CEM/GY2P+XYIQz8ZjY/5dghDP6yI3ijGqjw/rIjeKMaqPD8NrW4bUAQ/Pw2tbhtQBD8/N8u6KO28Mz83y7oo7bwzP9RroXljfzc/1GuheWN/Nz/B0g1sFOswP8HSDWwU6zA/v1z78hQjLT+/XPvyFCMtP4UbLlEoniU/hRsuUSieJT9ec/S8iq0UP15z9LyKrRQ/ECOBlE/MIj8QI4GUT8wiPzVVqK/t9B8/NVWor+30Hz8QI4GUT8wSPxAjgZRPzBI/XnP0vIqtFD9ec/S8iq0UP8HSDWwU6xA/wdINbBTrED/B0g1sFOsQP8HSDWwU6xA/rcNn5cWOBj+tw2flxY4GP8HSDWwU6xA/wdINbBTrED/B0g1sFOsQP8HSDWwU6xA/5gQ1h7IT7j7mBDWHshPuPuYENYeyE94+5gQ1h7IT3j7mBDWHshPuPuYENYeyE+4+5gQ1h7IT7j7mBDWHshPuPuYENYeyE/4+5gQ1h7IT/j7mBDWHshPuPuYENYeyE+4+rcNn5cWOBj+tw2flxY4GPwAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAOYENYeyE94+5gQ1h7IT3j7mBDWHshPePuYENYeyE94+\"\r\n", | |
" },\r\n", | |
" \"titles\": {}\r\n", | |
" }\r\n", | |
" ]\r\n", | |
"};\r\n", | |
" require.config({\r\n", | |
" baseUrl: \"https://cdn.rawgit.com/predictionmachines/Angara/0.2.2/dist\",\r\n", | |
" paths: {\r\n", | |
" \"angara.serializationjs\": \"https://cdn.rawgit.com/predictionmachines/Angara.Serialization/v0.2.0/dist/Angara.Serialization.umd\",\r\n", | |
" \"Chart\": \"https://cdn.rawgit.com/predictionmachines/Angara.Chart/0.2.0/dist/Chart\",\r\n", | |
" \"Table\": \"https://cdn.rawgit.com/predictionmachines/Angara.Table/0.2.1/dist/Table\",\r\n", | |
" \"angara.tablejs\": \"https://cdn.rawgit.com/predictionmachines/Angara.Table/0.2.1/dist/angara.tablejs\",\r\n", | |
" \"jquery.dataTables\": \"https://cdnjs.cloudflare.com/ajax/libs/datatables/1.10.11/js/jquery.dataTables.min\",\r\n", | |
" \"jquery.dataTables-css\": \"https://cdn.rawgit.com/predictionmachines/Angara.Table/0.2.1/dist/angara.tablejs\",\r\n", | |
" \"chartViewer.umd\": \"https://cdn.rawgit.com/predictionmachines/InteractiveDataDisplay/v1.2.1/dist/chartViewer.umd\",\r\n", | |
" \"chartViewer-css\": \"https://cdn.rawgit.com/predictionmachines/InteractiveDataDisplay/v1.2.1/dist/chartViewer\",\r\n", | |
" \"idd\": \"https://cdn.rawgit.com/predictionmachines/InteractiveDataDisplay/v1.2.1/dist/idd\",\r\n", | |
" \"idd-css\": \"https://cdn.rawgit.com/predictionmachines/InteractiveDataDisplay/v1.2.1/dist/idd\",\r\n", | |
" \"domReady\": \"https://cdnjs.cloudflare.com/ajax/libs/require-domReady/2.0.1/domReady.min\",\r\n", | |
" \"jquery\": \"https://cdnjs.cloudflare.com/ajax/libs/jquery/2.2.2/jquery.min\",\r\n", | |
" \"jquery-ui\": \"https://cdnjs.cloudflare.com/ajax/libs/jqueryui/1.11.4/jquery-ui\",\r\n", | |
" \"css\": \"https://cdnjs.cloudflare.com/ajax/libs/require-css/0.1.8/css.min\",\r\n", | |
" \"rx\": \"https://cdnjs.cloudflare.com/ajax/libs/rxjs/4.1.0/rx.lite.min\"\r\n", | |
" }\r\n", | |
" });\r\n", | |
" require(['Angara.Show', 'angara.serializationjs', 'domReady!'], function (Angara, Marshalling, document) {\r\n", | |
" var viewer = document.getElementById('container_11100160ce744343a7e1af6423d50a06');\r\n", | |
" var content_is = Marshalling.InfoSet.Unmarshal(content_11100160ce744343a7e1af6423d50a06);\r\n", | |
" var content = Marshalling.InfoSet.Deserialize(content_is);\r\n", | |
" return Angara.Show(content, viewer);\r\n", | |
" });\r\n", | |
" </script>\r\n", | |
" <div style=\"height:600px\" id=\"container_11100160ce744343a7e1af6423d50a06\" class=\"viewer-container\">\r\n", | |
" </div>\r\n", | |
"</section>\r\n" | |
] | |
}, | |
"metadata": {}, | |
"output_type": "display_data" | |
} | |
], | |
"source": [ | |
"Util.Html \"<h2 style='text-align:center'>Filzbach over analytic Г(2,1)</h2>\" |> Display\n", | |
"angaraDisplay chart" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": null, | |
"metadata": { | |
"collapsed": true | |
}, | |
"outputs": [], | |
"source": [] | |
} | |
], | |
"metadata": { | |
"kernelspec": { | |
"display_name": "IFSharp", | |
"language": "fsharp", | |
"name": "ifsharp" | |
}, | |
"language": "fsharp", | |
"language_info": { | |
"name": "fsharp", | |
"version": "1.0.0" | |
} | |
}, | |
"nbformat": 4, | |
"nbformat_minor": 0 | |
} |
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment