Skip to content

Instantly share code, notes, and snippets.

@vassilyl vassilyl/Filzbach.ipynb
Last active May 10, 2016

Embed
What would you like to do?
Sample Angara notebook
{
"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
You can’t perform that action at this time.