{"id":4101,"date":"2023-12-01T11:02:59","date_gmt":"2023-12-01T02:02:59","guid":{"rendered":"https:\/\/blog.since2020.jp\/?p=4101"},"modified":"2023-12-01T11:42:42","modified_gmt":"2023-12-01T02:42:42","slug":"mcmc3","status":"publish","type":"post","link":"https:\/\/since2020.jp\/media\/mcmc3\/","title":{"rendered":"\u3010MCMC\u2462\u3011MCMC\u3092\u7528\u3044\u3066\u5b9f\u969b\u306b\u30b5\u30f3\u30d7\u30ea\u30f3\u30b0\u3092\u3057\u3066\u307f\u3088\u3046\uff01"},"content":{"rendered":"\n<p>\u3053\u306e\u8a18\u4e8b\u3067\u306f\u3001MCMC\uff08\u30de\u30eb\u30b3\u30d5\u9023\u9396\u30e2\u30f3\u30c6\u30ab\u30eb\u30ed\u6cd5\uff09\u3067\u30b5\u30f3\u30d7\u30ea\u30f3\u30b0\u3092\u884c\u3046\u624b\u9806\u3068\u3001Python\u3092\u7528\u3044\u305f\u5b9f\u969b\u306e\u9069\u7528\u4f8b\u3092\u3054\u7d39\u4ecb\u3057\u307e\u3059\uff01<\/p>\n\n\n<h2>\u30b5\u30f3\u30d7\u30ea\u30f3\u30b0\u624b\u9806\u306b\u3064\u3044\u3066<\/h2>\n<p>MCMC\u3092\u4f7f\u7528\u3059\u308b\u969b\u306e\u57fa\u672c\u7684\u306a\u624b\u9806\u306f\u4ee5\u4e0b\u306e\u3088\u3046\u306b\u306a\u308a\u307e\u3059\u3002<\/p>\r\n<b><strong>1. \u554f\u984c\u306e\u5b9a\u5f0f\u5316<\/strong><\/b>\r\n<p>MCMC\u3092\u4f7f\u7528\u3059\u308b\u76ee\u7684\u3092\u660e\u78ba\u306b\u3057\u307e\u3059\u3002\u4f8b\u3048\u3070\u3001\u300c\u78ba\u7387\u5206\u5e03\u304b\u3089\u306e\u30b5\u30f3\u30d7\u30ea\u30f3\u30b0\u300d\u3084\u300c\u30d9\u30a4\u30ba\u7d71\u8a08\u30e2\u30c7\u30eb\u306e\u4e8b\u5f8c\u5206\u5e03\u306e\u30d1\u30e9\u30e1\u30fc\u30bf\u63a8\u5b9a\u300d\u306a\u3069\u3067\u3059\u3002<\/p>\r\n<b><strong>2. MCMC\u30a2\u30eb\u30b4\u30ea\u30ba\u30e0\u306e\u9078\u629e<\/strong><\/b>\r\n<p>\u300c\u30e1\u30c8\u30ed\u30dd\u30ea\u30b9\u30fb\u30d8\u30a4\u30b9\u30c6\u30a3\u30f3\u30b0\u30b9\u6cd5\u300d\u3001\u300c\u30ae\u30d6\u30b9\u30b5\u30f3\u30d7\u30ea\u30f3\u30b0\u300d\u3001\u300c\u30cf\u30df\u30eb\u30c8\u30cb\u30a2\u30f3\u30fb\u30e2\u30f3\u30c6\u30ab\u30eb\u30ed\u300d\u306a\u3069\u3001\u9069\u5207\u306aMCMC\u30a2\u30eb\u30b4\u30ea\u30ba\u30e0\u3092\u9078\u629e\u3057\u307e\u3059\u3002\uff08\u2192\u30a2\u30eb\u30b4\u30ea\u30ba\u30e0\u306e\u8a73\u3057\u3044\u8aac\u660e\u306f<a href=\"https:\/\/blog.since2020.jp\/data_analysis\/mcmc2\/\">\u3053\u3061\u3089<\/a>\uff09<\/p>\r\n<p><span style=\"text-decoration: underline\">\u30d1\u30e9\u30e1\u30fc\u30bf\u8a2d\u5b9a<\/span><\/p>\r\n<p>\u30a2\u30eb\u30b4\u30ea\u30ba\u30e0\u306b\u3088\u3063\u3066\u7570\u306a\u308a\u307e\u3059\u304c\u3001MCMC\u306e\u4e3b\u8981\u306a\u30d1\u30e9\u30e1\u30fc\u30bf\u306b\u306f\u3001\u300c\u521d\u671f\u5024\u300d\u3001\u300c\u30d0\u30fc\u30f3\u30a4\u30f3\u671f\u9593\u300d\u3001\u300c\u30b9\u30c6\u30c3\u30d7\u30b5\u30a4\u30ba\u300d \u3001\u300c\u30b5\u30f3\u30d7\u30eb\u6570\u300d\u306a\u3069\u304c\u3042\u308a\u307e\u3059\u3002<\/p>\r\n<p>\u203b\u30d0\u30fc\u30f3\u30a4\u30f3\u671f\u9593\u3068\u306f\u3001\u30a2\u30eb\u30b4\u30ea\u30ba\u30e0\u306e\u521d\u671f\u6bb5\u968e\u3067\u751f\u6210\u3055\u308c\u308b\u30b5\u30f3\u30d7\u30eb\u3092\u3001\u5206\u6790\u304b\u3089\u9664\u5916\u3059\u308b\u671f\u9593\u306e\u3053\u3068\u3067\u3059\u3002\u30d0\u30fc\u30f3\u30a4\u30f3\u671f\u9593\u304c\u77ed\u3059\u304e\u308b\u3068\u3001\u672a\u53ce\u675f\u306e\u30b5\u30f3\u30d7\u30eb\u304c\u5206\u6790\u306b\u542b\u307e\u308c\u308b\u30ea\u30b9\u30af\u304c\u3042\u308a\u307e\u3059\u3002\u4e00\u65b9\u3067\u3001\u9577\u3059\u304e\u308b\u3068\u6709\u7528\u306a\u30b5\u30f3\u30d7\u30eb\u3092\u5931\u3046\u3053\u3068\u306b\u306a\u308a\u307e\u3059\u3002<\/p>\r\n<b><strong>3. \u30b5\u30f3\u30d7\u30ea\u30f3\u30b0\u306e\u5b9f\u884c<\/strong><\/b>\r\n<p>\u521d\u671f\u5024\u3092\u8a2d\u5b9a\u3057\u305f\u5f8c\u3001\u9078\u629e\u3057\u305f\u30a2\u30eb\u30b4\u30ea\u30ba\u30e0\u306b\u5f93\u3063\u3066\u3001\u76ee\u7684\u306e\u5206\u5e03\u304b\u3089\u30b5\u30f3\u30d7\u30eb\u3092\u751f\u6210\u3057\u307e\u3059\u3002<\/p>\r\n<b><strong>4. \u53ce\u675f\u306e\u78ba\u8a8d<\/strong><\/b>\r\n<p>\u751f\u6210\u3055\u308c\u305f\u30b5\u30f3\u30d7\u30eb\u304c\u76ee\u7684\u306e\u5206\u5e03\u306b\u53ce\u675f\u3057\u3066\u3044\u308b\u304b\u3092\u8a55\u4fa1\u3057\u307e\u3059\u3002<\/p>\r\n<b><strong>5. \u7d50\u679c\u306e\u5206\u6790<\/strong><\/b>\r\n<ul>\r\n\t<li>\u30d0\u30fc\u30f3\u30a4\u30f3\u671f\u9593\u30b5\u30f3\u30d7\u30eb\u306e\u9664\u5916\uff1a\u521d\u671f\u306e\u30b5\u30f3\u30d7\u30eb\u3092\u6392\u9664\u3057\u307e\u3059\u3002<\/li>\r\n\t<li>\u4e8b\u5f8c\u5206\u5e03\u306e\u5206\u6790\uff08\u5fc5\u8981\u306a\u5834\u5408\uff09\uff1a\u53ce\u675f\u3057\u305f\u30b5\u30f3\u30d7\u30eb\u3092\u4f7f\u7528\u3057\u3066\u3001\u63a8\u5b9a\u3057\u305f\u3044\u30d1\u30e9\u30e1\u30fc\u30bf\u306e\u4e8b\u5f8c\u5206\u5e03\u3092\u5206\u6790\u3057\u307e\u3059\u3002<\/li>\r\n<\/ul>\r\n<p><!-- notionvc: a9bb0d03-224b-4cad-ba2f-d5b5e23de3a3 --><\/p>\n\n<h2>\u5b9f\u969b\u306bPython\u3067\u5b9f\u88c5\u3057\u3066\u307f\u3088\u3046\uff01<\/h2>\n<p>\u5148\u307b\u3069\u306e\uff15\u3064\u306e\u30b9\u30c6\u30c3\u30d7\u306b\u5247\u3063\u3066\u3001\u5b9f\u969b\u306b\u30b5\u30f3\u30d7\u30ea\u30f3\u30b0\u3092\u884c\u3044\u305f\u3044\u3068\u601d\u3044\u307e\u3059\u3002<\/p>\r\n<p>\u3053\u3053\u3067\u306f\u3001\u7c21\u5358\u306a\u7d71\u8a08\u30e2\u30c7\u30eb\uff08\u30d9\u30a4\u30b8\u30a2\u30f3\u7dda\u5f62\u56de\u5e30\u30e2\u30c7\u30eb\uff09\u3092\u53d6\u308a\u4e0a\u3052\u3001MCMC\u3092\u7528\u3044\u305f\u4e8b\u5f8c\u5206\u5e03\u306e\u30d1\u30e9\u30e1\u30fc\u30bf\u63a8\u5b9a\u3092\u884c\u3046\u3053\u3068\u3092\u76ee\u6a19\u3068\u3057\u307e\u3059\u3002<\/p>\n\n<h2>1. \u554f\u984c\u306e\u5b9a\u5f0f\u5316<\/h2>\n<p><!-- notionvc: fb463d31-752e-4bcf-840b-f1b5d412609b --><\/p>\r\n<p>\u4eca\u56de\u306f\u3001\u7dda\u5f62\u56de\u5e30\u30e2\u30c7\u30eb\u306e\u30d1\u30e9\u30e1\u30fc\u30bf\uff08\u50be\u304d\u3068\u5207\u7247\uff09\u3092\u63a8\u5b9a\u3059\u308b\u3053\u3068\u3092\u76ee\u7684\u3068\u3057\u307e\u3059\u3002<\/p>\r\n<p>\u30e2\u30c7\u30eb\u5f0f\u306f\u6b21\u306e\u3068\u304a\u308a\u3067\u3059\u3002<\/p>\r\n<p><img decoding=\"async\" src=\"https:\/\/since2020.jp\/media\/wp-content\/uploads\/2023\/11\/%E3%82%B9%E3%82%AF%E3%83%AA%E3%83%BC%E3%83%B3%E3%82%B7%E3%83%A7%E3%83%83%E3%83%88-2023-12-01-10.53.52-300x99.png\" alt=\"\" data-name=\"image\" class=\"aligncenter\" width=\"176\" height=\"58\" \/><\/p>\r\n<p>\uff08<em>\u03f5<\/em> \u306f\u8aa4\u5dee\u9805\u3067\u3001\u6b63\u898f\u5206\u5e03\u306b\u5f93\u3046\u3068\u4eee\u5b9a\uff09<\/p>\r\n<p>\u4e0a\u5f0f\u306b\u304a\u3044\u3066\u3001\u4eca\u56de\u306f\u771f\u306e\u30d1\u30e9\u30e1\u30fc\u30bf\u3092<\/p>\r\n<p><img decoding=\"async\" src=\"https:\/\/since2020.jp\/media\/wp-content\/uploads\/2023\/11\/%E3%82%B9%E3%82%AF%E3%83%AA%E3%83%BC%E3%83%B3%E3%82%B7%E3%83%A7%E3%83%83%E3%83%88-2023-12-01-10.55.07-300x104.png\" alt=\"\" data-name=\"image\" class=\"aligncenter\" \/><\/p>\r\n<p>\u3068\u8a2d\u5b9a\u3059\u308b\u3053\u3068\u306b\u3057\u307e\u3059\u3002<\/p>\r\n<p>\u306f\u3058\u3081\u306b\u30e9\u30a4\u30d6\u30e9\u30ea\u306e\u30a4\u30f3\u30dd\u30fc\u30c8\u3068\u3001\u4eee\u60f3\u30c7\u30fc\u30bf\u306e\u751f\u6210\u3092\u884c\u3044\u307e\u3059\u3002<\/p>\r\n<div class=\"hcb_wrap\">\r\n<pre class=\"prism line-numbers lang-plain\" data-lang=\"Plain Text\">import numpy as np \r\nimport scipy.stats as stats \r\nimport matplotlib.pyplot as plt \r\n\r\n#random.seed\u3092\u8a2d\u5b9a \r\nnp.random.seed(0) \r\n# \u4eee\u60f3\u30c7\u30fc\u30bf\u3092\u4f5c\u6210 \r\nx = np.linspace(0, 10, 100) y = 2.5 + 0.5 * x + np.random.normal(scale=1.0, size=100)<!-- notionvc: 614512f7-9e66-4e3b-a53f-fef1718a6771 --><\/pre>\r\n<\/div>\r\n<p><!-- notionvc: 61cc5578-ce39-4daa-9933-7167a9a765ed --><\/p>\r\n<p><!-- notionvc: 3b2d5b75-db61-403f-8491-ec73ec9ad7db --><\/p>\n\n<h2>2. MCMC\u30a2\u30eb\u30b4\u30ea\u30ba\u30e0\u306e\u9078\u629e \u304a\u3088\u3073 3. \u30b5\u30f3\u30d7\u30ea\u30f3\u30b0\u306e\u5b9f\u884c<\/h2>\n<p>\u4eca\u56de\u306f\u3001\u30e1\u30c8\u30ed\u30dd\u30ea\u30b9\u30fb\u30d8\u30a4\u30b9\u30c6\u30a3\u30f3\u30b0\u30b9\u6cd5\u3092\u9078\u629e\u3057\u307e\u3059\u3002<\/p>\r\n<p>\u4e0b\u8a18\u306e\u30b3\u30fc\u30c9\u3067\u306f\u3001<\/p>\r\n<ol>\r\n\t<li>\u73fe\u5728\u306e\u30d1\u30e9\u30e1\u30fc\u30bf\u5024\u304b\u3089\u65b0\u3057\u3044\u5024\u3092\u63d0\u6848\uff08\u6b63\u898f\u5206\u5e03\u3092\u4f7f\u7528\uff09<\/li>\r\n\t<li>\u63d0\u6848\u3055\u308c\u305f\u5024\u306e\u53d7\u7406\u78ba\u7387\u3092\u8a08\u7b97\u3057\u3001\u30e9\u30f3\u30c0\u30e0\u306b\u53d7\u7406\/\u62d2\u5426\u3092\u6c7a\u5b9a\u3002<\/li>\r\n\t<li>\u3053\u306e\u30d7\u30ed\u30bb\u30b9\u30925000\u56de\u7e70\u308a\u8fd4\u3059<\/li>\r\n<\/ol>\r\n<p>\u3068\u3044\u3046\u3053\u3068\u304c\u884c\u308f\u308c\u3066\u3044\u307e\u3059\u3002<\/p>\r\n<p><!-- notionvc: 8c5b20f6-b1e2-4ae0-b192-366f23f8f297 --><\/p>\r\n<div>\r\n<div>\r\n<div class=\"hcb_wrap\">\r\n<pre class=\"prism line-numbers lang-python\" data-lang=\"Python\"><code># \u30e1\u30c8\u30ed\u30dd\u30ea\u30b9\u30fb\u30d8\u30a4\u30b9\u30c6\u30a3\u30f3\u30b0\u30b9\u6cd5\r\ndef metropolis_hastings(lr, iterations, x, y):\r\nbeta_0 = 0\r\nbeta_1 = 0\r\nsigma = 1\r\nsamples = []\r\nfor i in range(iterations):\r\nbeta_0_proposal = np.random.normal(beta_0, lr)\r\nbeta_1_proposal = np.random.normal(beta_1, lr)\r\nsigma_proposal = np.random.normal(sigma, lr)\r\n\r\n\r\nlikelihood_current = np.sum(stats.norm.logpdf(y, beta_0 + beta_1 * x, sigma))\r\nlikelihood_proposal = np.sum(stats.norm.logpdf(y, beta_0_proposal + beta_1_proposal * x, sigma_proposal))\r\n\r\n\r\nprior_current = stats.norm.logpdf(beta_0, 0, 10) + stats.norm.logpdf(beta_1, 0, 10) + stats.halfnorm.logpdf(sigma, 0, 1)\r\nprior_proposal = stats.norm.logpdf(beta_0_proposal, 0, 10) + stats.norm.logpdf(beta_1_proposal, 0, 10) + stats.halfnorm.logpdf(sigma_proposal, 0, 1)\r\n\r\n\r\n#\u53d7\u7406\u78ba\u7387\r\np_accept = np.exp(likelihood_proposal + prior_proposal - likelihood_current - prior_current)\r\n\r\n\r\nif np.random.rand() &lt; p_accept:\r\nbeta_0, beta_1, sigma = beta_0_proposal, beta_1_proposal, sigma_proposal\r\n\r\n\r\nsamples.append([beta_0, beta_1, sigma])\r\n\r\n\r\nreturn np.array(samples)\r\n\r\n\r\n# MCMC\u3092\u5b9f\u884c\r\nsamples = metropolis_hastings(0.05, 5000, x, y)<\/code><\/pre>\r\n<\/div>\r\n<\/div>\r\n<\/div>\n\n<h2>4. \u53ce\u675f\u306e\u78ba\u8a8d<\/h2>\n<p>\u672c\u6765\u306a\u3089\u3001\u6b21\u306b\u53ce\u675f\u5224\u5b9a\u3092\u884c\u3044\u307e\u3059\u304c\u3001\u3053\u3053\u3067\u306f\u30b3\u30fc\u30c9\u304c\u9577\u304f\u8907\u96d1\u306a\u3082\u306e\u306b\u306a\u3063\u3066\u3057\u307e\u3046\u305f\u3081\u3001\u7701\u7565\u3057\u307e\u3059\u3002<!-- notionvc: ace42267-7fe8-49d1-b572-df5865636b2e --><\/p>\n\n<h2>5. \u7d50\u679c\u306e\u5206\u6790<\/h2>\n<div class=\"hcb_wrap\">\r\n<pre class=\"prism line-numbers lang-plain\" data-lang=\"Plain Text\"><code># \u30d0\u30fc\u30f3\u30a4\u30f3\u671f\u9593\u3092\u53d6\u308a\u9664\u304f\r\nburn_in = 500\r\npost_burn_in_samples = samples[burn_in:]<\/code><\/pre>\r\n<\/div>\r\n<p>\u30d0\u30fc\u30f3\u30a4\u30f3\u671f\u9593\u306e\u9664\u5916\u5f8c\u3001\u6b8b\u308a\u306e\u30b5\u30f3\u30d7\u30eb\u3092\u4f7f\u7528\u3057\u3066\u3001<span data-token-index=\"1\" class=\"notion-enable-hover\">\u03b2<\/span>0 \u3068 <span data-token-index=\"3\" class=\"notion-enable-hover\">\u03b2<\/span>1 \u306e\u4e8b\u5f8c\u5206\u5e03\u3092\u5206\u6790\u3057\u307e\u3059\u3002<!-- notionvc: a041f197-4494-462b-ba10-8775e5e0ef96 --><\/p>\r\n<div class=\"hcb_wrap\">\r\n<pre class=\"prism line-numbers lang-html\" data-lang=\"HTML\"><code>plt.figure(figsize=(12, 6))\r\n\r\n# beta_0\u306e\u4e8b\u5f8c\u5206\u5e03\r\nplt.subplot(311)\r\nplt.hist(post_burn_in_samples[:, 0], bins=30, density=True)\r\nplt.title('Posterior of $\\\\beta_0$ After Burn-in')\r\n\r\n# beta_1\u306e\u4e8b\u5f8c\u5206\u5e03\r\nplt.subplot(312)\r\nplt.hist(post_burn_in_samples[:, 1], bins=30, density=True)\r\nplt.title('Posterior of $\\\\beta_1$ After Burn-in')\r\n\r\n#\u4e8b\u5f8c\u5206\u5e03\u306e\u5e73\u5747\u3092\u6c42\u3081\u308b\r\nmean_beta_0 = np.mean(post_burn_in_samples[:, 0])\r\nmean_beta_1 = np.mean(post_burn_in_samples[:, 1])\r\nprint(mean_beta_0, mean_beta_1)<\/code><\/pre>\r\n<\/div>\r\n<p><img decoding=\"async\" src=\"https:\/\/since2020.jp\/media\/wp-content\/uploads\/2023\/12\/Untitled-2-300x62.png\" alt=\"\" width=\"300\" height=\"62\" class=\"alignnone size-medium wp-image-4191\" srcset=\"https:\/\/since2020.jp\/media\/wp-content\/uploads\/2023\/12\/Untitled-2-300x62.png 300w, https:\/\/since2020.jp\/media\/wp-content\/uploads\/2023\/12\/Untitled-2-768x159.png 768w, https:\/\/since2020.jp\/media\/wp-content\/uploads\/2023\/12\/Untitled-2-960x201.png 960w, https:\/\/since2020.jp\/media\/wp-content\/uploads\/2023\/12\/Untitled-2.png 968w\" sizes=\"(max-width: 300px) 100vw, 300px\" \/><br \/>\r\n<img decoding=\"async\" src=\"https:\/\/since2020.jp\/media\/wp-content\/uploads\/2023\/12\/Untitled-1-300x62.png\" alt=\"\" width=\"300\" height=\"62\" class=\"alignnone size-medium wp-image-4189\" srcset=\"https:\/\/since2020.jp\/media\/wp-content\/uploads\/2023\/12\/Untitled-1-300x62.png 300w, https:\/\/since2020.jp\/media\/wp-content\/uploads\/2023\/12\/Untitled-1-768x158.png 768w, https:\/\/since2020.jp\/media\/wp-content\/uploads\/2023\/12\/Untitled-1-960x201.png 960w, https:\/\/since2020.jp\/media\/wp-content\/uploads\/2023\/12\/Untitled-1.png 977w\" sizes=\"(max-width: 300px) 100vw, 300px\" \/><\/p>\r\n<p>&nbsp;<\/p>\r\n<p>MCMC\u3067\u30b5\u30f3\u30d7\u30ea\u30f3\u30b0\u3057\u305f\u4e8b\u5f8c\u5206\u5e03\u306e\u5e73\u5747\u306f\u3001<!-- notionvc: 5ee8311c-c8e4-4bd4-be42-272d23897758 --><\/p>\r\n<div class=\"hcb_wrap\">\r\n<pre class=\"prism line-numbers lang-plain\" data-lang=\"Plain Text\"><code>print(mean_beta_0, mean_beta_1)<\/code><\/pre>\r\n<\/div>\r\n<p>\u306b\u3066\u51fa\u529b\u3057\u3066\u3044\u307e\u3059\u3002\u305d\u306e\u7d50\u679c\u306f\u3001<!-- notionvc: 9da2e7d4-029e-4c20-b01e-80941dadf4d4 --><\/p>\r\n<p><img decoding=\"async\" src=\"https:\/\/since2020.jp\/media\/wp-content\/uploads\/2023\/12\/%E3%82%B9%E3%82%AF%E3%83%AA%E3%83%BC%E3%83%B3%E3%82%B7%E3%83%A7%E3%83%83%E3%83%88-2023-12-01-11.25.34-300x151.png\" alt=\"\" data-name=\"image\" width=\"111\" height=\"56\" class=\"aligncenter\" \/><\/p>\r\n<p>\u3068\u306a\u308a\u307e\u3057\u305f\u3002<\/p>\r\n<p>\u306f\u3058\u3081\u306b\u3001\u771f\u306e\u30d1\u30e9\u30e1\u30fc\u30bf\u306f\u3001<\/p>\r\n<p><img decoding=\"async\" src=\"https:\/\/since2020.jp\/media\/wp-content\/uploads\/2023\/11\/%E3%82%B9%E3%82%AF%E3%83%AA%E3%83%BC%E3%83%B3%E3%82%B7%E3%83%A7%E3%83%83%E3%83%88-2023-12-01-10.55.07-300x104.png\" alt=\"\" data-name=\"image\" width=\"190\" height=\"66\" class=\"aligncenter\" \/><\/p>\r\n<p>\u3068\u8a2d\u5b9a\u3057\u307e\u3057\u305f\u306e\u3067\u3001\u30b5\u30f3\u30d7\u30ea\u30f3\u30b0\u306b\u3088\u308a\u3001\u771f\u306e\u30d1\u30e9\u30e1\u30fc\u30bf\u306b\u6bd4\u8f03\u7684\u8fd1\u3044\u3082\u306e\u304c\u5f97\u3089\u308c\u305f\u306e\u3067\u306f\u306a\u3044\u3067\u3057\u3087\u3046\u304b\u3002<\/p>\r\n<p>\u4ee5\u4e0a\u3001Python\u3092\u7528\u3044\u3066MCMC\u3067\u30b5\u30f3\u30d7\u30ea\u30f3\u30b0\u3092\u884c\u3046\u65b9\u6cd5\u3092\u3054\u7d39\u4ecb\u3057\u307e\u3057\u305f\u3002\u4eca\u56de\u306f\u3001\u30a2\u30eb\u30b4\u30ea\u30ba\u30e0\u3092\u3088\u308a\u6df1\u304f\u7406\u89e3\u3059\u308b\u305f\u3081\u306b\u3001\u30e9\u30a4\u30d6\u30e9\u30ea\u306f\u7528\u3044\u305a\u306b\u95a2\u6570\u3092\u5b9a\u7fa9\u3057\u307e\u3057\u305f\u304c\u3001\u30e9\u30a4\u30d6\u30e9\u30ea<strong>PyMC3<\/strong>\u3092\u4f7f\u3046\u3068\u7c21\u5358\u306b\u5b9f\u88c5\u3067\u304d\u307e\u3059\u3002<\/p>\r\n<p><!-- notionvc: f85238b1-cb92-4b64-9455-a27cad579222 --><\/p>\r\n<p><!-- notionvc: 3bcb3d94-120c-4f48-83b0-8e64fa842950 --><\/p>","protected":false},"excerpt":{"rendered":"<p>\u3053\u306e\u8a18\u4e8b\u3067\u306f\u3001MCMC\uff08\u30de\u30eb\u30b3\u30d5\u9023\u9396\u30e2\u30f3\u30c6\u30ab\u30eb\u30ed\u6cd5\uff09\u3067\u30b5\u30f3\u30d7\u30ea\u30f3\u30b0\u3092\u884c\u3046\u624b\u9806\u3068\u3001Python\u3092\u7528\u3044\u305f\u5b9f\u969b\u306e\u9069\u7528\u4f8b\u3092\u3054\u7d39\u4ecb\u3057\u307e\u3059\uff01 \u30b5\u30f3\u30d7\u30ea\u30f3\u30b0\u624b\u9806\u306b\u3064\u3044\u3066 MCMC\u3092\u4f7f\u7528\u3059\u308b\u969b\u306e\u57fa\u672c\u7684\u306a\u624b\u9806\u306f\u4ee5\u4e0b\u306e\u3088\u3046\u306b\u306a\u308a\u307e\u3059\u3002 1. [&hellip;]<\/p>\n","protected":false},"author":16,"featured_media":3651,"comment_status":"open","ping_status":"open","sticky":false,"template":"","format":"standard","meta":{"content-type":"","swell_btn_cv_data":"","footnotes":"","_wp_rev_ctl_limit":""},"categories":[1246],"tags":[],"class_list":["post-4101","post","type-post","status-publish","format-standard","has-post-thumbnail","hentry","category-data-infrastructure"],"_links":{"self":[{"href":"https:\/\/since2020.jp\/media\/wp-json\/wp\/v2\/posts\/4101","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/since2020.jp\/media\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/since2020.jp\/media\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/since2020.jp\/media\/wp-json\/wp\/v2\/users\/16"}],"replies":[{"embeddable":true,"href":"https:\/\/since2020.jp\/media\/wp-json\/wp\/v2\/comments?post=4101"}],"version-history":[{"count":1,"href":"https:\/\/since2020.jp\/media\/wp-json\/wp\/v2\/posts\/4101\/revisions"}],"predecessor-version":[{"id":4193,"href":"https:\/\/since2020.jp\/media\/wp-json\/wp\/v2\/posts\/4101\/revisions\/4193"}],"wp:featuredmedia":[{"embeddable":true,"href":"https:\/\/since2020.jp\/media\/wp-json\/wp\/v2\/media\/3651"}],"wp:attachment":[{"href":"https:\/\/since2020.jp\/media\/wp-json\/wp\/v2\/media?parent=4101"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/since2020.jp\/media\/wp-json\/wp\/v2\/categories?post=4101"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/since2020.jp\/media\/wp-json\/wp\/v2\/tags?post=4101"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}