{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# 9. Mixture Models and EM" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "import matplotlib.pyplot as plt\n", "import numpy as np\n", "from sklearn.datasets import fetch_openml\n", "%matplotlib inline\n", "\n", "from prml.clustering import KMeans\n", "from prml.rv import (\n", " MultivariateGaussianMixture,\n", " BernoulliMixture\n", ")\n", "\n", "np.random.seed(2222)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## 9.1 K-means Clustering" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "# training data\n", "x1 = np.random.normal(size=(100, 2))\n", "x1 += np.array([-5, -5])\n", "x2 = np.random.normal(size=(100, 2))\n", "x2 += np.array([5, -5])\n", "x3 = np.random.normal(size=(100, 2))\n", "x3 += np.array([0, 5])\n", "x_train = np.vstack((x1, x2, x3))\n", "\n", "x0, x1 = np.meshgrid(np.linspace(-10, 10, 100), np.linspace(-10, 10, 100))\n", "x = np.array([x0, x1]).reshape(2, -1).T" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAARMAAAD8CAYAAABUzEBbAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/MnkTPAAAACXBIWXMAAAsTAAALEwEAmpwYAAA+/UlEQVR4nO2dd5hU1fnHP+/U7bssvTdBAVFQFBVrLAFLrFFMzM+SaBI1xpLEGmti1xhji7HHFkERFLBgASsiWOi9LR2W7Tv9/f1xZ5eZnTuwZWYb5/M8++zce8+999zdme+c8j3vK6qKwWAwNBVHS1fAYDC0D4yYGAyGlGDExGAwpAQjJgaDISUYMTEYDCnBiInBYEgJKRETEXlORLaKyIKYfYUi8qGILI/+7pDk3LEislREVojIDamoj8FgaH5S1TJ5ARhbZ98NwEeqOgj4KLodh4g4gceBccBQ4HwRGZqiOhkMhmYkJWKiqrOA4jq7TwdejL5+ETjD5tRDgRWqukpVA8Dr0fMMBkMbw5XGa3dV1U0AqrpJRLrYlOkJrI/ZLgJG211MRC4DLgPIzs46eN/99klxdQ0GQw3z5v64XVU7N+ScdIpJfRCbfbb+flV9Gnga4OBRB+rs2e+ns14Gw16N29V9bUPPSedszhYR6Q4Q/b3VpkwR0DtmuxewMY11MhgMaSKdYjIFuDD6+kJgsk2ZOcAgEekvIh5gfPQ8g8HQxkjV1PBrwFfAviJSJCK/Bu4FThSR5cCJ0W1EpIeITANQ1RBwJfA+sBh4Q1UXpqJOBoOheUnJmImqnp/k0PE2ZTcCJ8dsTwOmpaIeBoOh5TAOWIPBkBKMmBgMhpRgxMRgMKQEIyYGgyElGDExGAwpwYiJwWBICUZMDAZDSjBiYjAYUoIRE4PBkBKMmBgMhpRgxMRgMKQEIyYGgyElGDExGAwpwYiJwWBICUZMDAZDSjBiYjAYUoIRE4PBkBKMmBgMhpSQVjERkX1F5PuYnzIRubpOmWNFpDSmzK3prJPBYEgPac2bo6pLgRFQmwp0AzDJpuhnqnpqOutiMBjSS3N2c44HVqpqg5P7GAyG1k9zisl44LUkxw4XkR9EZLqIDGvGOhkMhhTRLGISTbD1M2CCzeF5QF9VPRD4F/B2kmtcJiLfisi327ftSFtdDfFoeAfq/xoNrmr8NTSElZfe0J5prlzD44B5qrql7gFVLYt5PU1EnhCRTqq6vU65uFzD6a7w3o6qohWPQvVkEA9oCHUNQAruRxwF9btGpAQtux8CXwKKuoYgeTcgrn7prLqhhWiubs75JOniiEg3EZHo60OjdTJNjxZGfdOh+h0gAFoB+CC0BC25qX7nawTd+YeokISAMIQWojt/j0bK9nS6oQ2SdjERkSys9KBvxez7nYj8Lrp5DrBARH4AHgXGq6ppebQ01f8DfHV2RiD0I5HSv7HHf1Hge4hsxhKSGhQ0gFZPT2lVDa2DtHdzVLUK6Fhn31Mxrx8DHkt3PQwNJFKR/Jj/E/AdDpkJ2V9jzi8CW8HxQ3hNU2tnaIUYB6zBHs8RgDPJQT9a/fbuz3fuk+RABrj2a3y9DK0WIyYGWyT7QpDs3ZTw7/4C7iHgGgx4Y3Y6wZGDZJyYghoaWhtGTAxxqIYBEGcnKHyZeDGowQvOAWjlc2j1e6jWHVsBEUE6PAyZZ4Dkg2SB9ydIh2cQR1Zan8HQMkhbHOs8eNSBOnv2+y1djXaDqqJVr0PVy6Cl4OwBWVfiyDwa9X+Flt4CRIAgkAEoIIAPJBMkAyl4CnH1bMnHMKQQt6v7XFUd1ZBzTMvEgFa9BFXPWkICEN4I5Xeg/m8Q7+FI4X8h6xeQcQq4DsQSlmhrRKshUoqW3d24e0cqiVRNIFJyPZGyR9CQWW3RVmku05qhlaIagqpXIKGr4kcrn0G8hyKuHkjOpQBEtp2K1UKJJWJ5SCJVDerCaKQM3flrCBdjjcE4Ud87kP83xHt4E57K0BKYlsneTqQcouMkCYSL0nprrXoFwjvYNZgbBvxo+d2oRtJ6b0PqMWKyt+PIBXHbH3P2RMPF8fsyTgI8dS8CruENH1j1zwRs1uxoNYTXN+xahhbHiMlej4JzX5v9AqFl6I5ziBT/Hxpcae3N/jW4+lizMzis344OSN5NaKQS9c9GAwvr17JINvWsGr2+oS1hxkz2crT0ZgjNtzuC1e0IQ2gVWnIldJyIOLKhw7MQ+AZCK8DRHTKOQavfgYrHQVxABCQbzbkW8Y7BiouViGSeg5Y/RLxt3wmuQYizc8qf1ZBeTMtkL0ZDayEwF9uuRkLhIPhnACDitGZ5sn+FZJ6ABn6Ain8BftDK6AzPdii7Cd02lohvlv01M8aC+6CYHQKOTkj+XYm3j1Si4W1mLKUVY1omezOhlVZLQvfgZgXAh4Y3I3X2RqomQ8UDuzmvGspuIlI1GnQLSAGSdS6ScTTqmwXBuTFlFSIlEFwO0ZaJRqrQsnsh8BngAEc2mv0nHJlHN+hRDenHiMnejLN38pmcukgm4hpqOWSDS9HIdvB/Cf5363d+aHbtSy1bgoYvAN80Em35frTycSTjCKts6a0QnEftdHTEb3lgnP9CPEPrd29Ds2DEZC9G3INQ12AILSG+q+PAemvU7POAswfq6AA7zoRIFYnhCRqCDypfID48QQzRKWkNb4XgdyR2wwJo1auI52/1uptqBALfoIG5Vjcq40TEWdjYyhuSYMRkL0cKHkDLHwH/R0DIWtGb80cIzAbfdCBiTQd7joKSq2iaiMQSwrLm21zPEY1YEd6apBumVoiDeqAaRHdeC+El1lgOXrTqP5D3AOId2YT6G+pixGQvRxzZSP7NqN4ARJAaz4lnGORcgkbKrbU5Va+RtCVhh6NbNDjS7o75sBYSxopFBmRfbL109QO1u6cL3CPqVQ2tfgdCi9klWn4rRlP5X8EzOelMk6HhmNkcAxCdobExr2nJX6JjFvUVEheSeyeOThOR3OsgbshWkNy/RI/dBjhBOoDkYL0Vo2UDS9FIGeLIgazxWC2YGhzW+E3WL+tXHd972LZ+NGBNbRtShmmZGGpRVaLheAGIVE9L4kFJhgvJvRXJ/AkAknmmdd3yh63t3L8gmadFj50ICFr+N3AOgfBSrLERH/inoqE5UPgykv0bcPZCq1+zZnpcw8E9GAJfoZ4jEWchGiqyVj2HloFrMJJ1fswK5mRvcSV58CdDYzBispejGkArngTfO6B+1DUMyb0WXAOhoiHRNHcJiWo5Ig+jen2toIALyTwN1WpE7kf1WiTzBCxBuY/4QdagtWbH9wmSeRJkjkMyxxGpessyxgW+RhHgn2jWL6D69ei4SsQKeu1/HwoeQ9z7IpmnouXLSGidSJ71jIaU0RwBpdeIyPxoHuFvbY6LiDwqIitE5EcROcjuOob0oKW3QvWU6KphhdACtOQKNDjf6gokxUnsN7vk3hgjJOOBpxG5GFUfknlmjJBcGD32C6ts5vFI7tU21/dZdaipZ2iDJST4oz8+63fVC9GB1RozWwS0Gi2/x9rMGAueMVhdJY9l05ccJP/uuFaYoek0V8vkuLp5cGIYBwyK/owGnoz+NqQZDW2AwBwSvB4aBN+H7PqA1sUJHf5rzQBVvQoEwHskACLTgO+j5WZFBeV5QKNC8kX02DxEPgDOBu9RUO5IvJ9W7nrt/yRJfZIE9wqtiHbbHEjBHWhwGQS/t8ZoMo5CJMP+PEOjaQ0DsKcDL6nF10CBiHRv6UrtFYTXRdfS1CUEoXVRq7vNimJnX8SZF13Z6wPCaNndVmwUzgOujyk8C5GjETmaXUICcDNwtpXtr+webIUiZrGfde0GWunDG3ddyj3Yct5mnmiEJE00h5go8IGIzBWRy2yO9wRi15sXRffFYdKDpgFn3yRTr25w74vk3wauQYmHw2utBFv+j6ltGQRmoqW3RT/0VxEvKBuAjTHbNwOXW0JS+lcI2K3d8SLuXVHsxXs0DWtIm8HV5qY5xGSMqh6E1Z25QqyvqFjsOq4JbVdVfVpVR6nqqE6dO9qcYqgvqkHU/7U1++E+mISg0eK2vsUdueA5jMTWSRjCm0j4NwVmoqV3RDeuAk6wuftY4HKrHqV3Rtfc2BFBXbty2It7AGSdF62rA0ssvOAahe3b2NHVimVraDaaIwnXxujvrSIyCTgUiP0qKgJ6x2z3Iv5rzJBCNLgYLfkTlm9ErZaJa7g1NavV4DoAyb0acXZBg8sh8BWJYRqxzsULVMXsE8R7qHVUqxFZZ3PeWmtQVjIQ70g08An24x5BKPkt2uH52mleR85lqPdY1DfDmib2HAbeMVByDYSWYo39eEFcSP5dZoC1mUmrmIhINuBQ1fLo65OAO+sUmwJcKSKvYw28lqrqpnTWa29FNYiWXAe7csVbhBYgHR6v7VaohogU/yb6AU2WvcBP/NtHan0ku2Ztltmct7h2UDbeh2JzH61CK59H8m/ZtS+8GXxTrPKBT6CyAPLvRSLFlifG0Qm8x1txVwzNSrpbJl2BSdFvCBfwqqq+V5NnOJomdBpwMrAC62vu4jTXae8l8C32TtYgWv3OLjEp+Ut08d+e2HUtybmmjpDEDraejzV4+r/o9qx4QdEQWvFP+1v4P0eDixH3EMucVnYHcbNPWg0lf4RObyPeQ+pRZ0O6SKuYqOoq4ECb/bG5hhW4Ip31METR6iQHIrXTsBqphOA3DbywEzLGASDyJomzNpdHX/cD7ou+noXI28B4yDg5GlzJbramAi35A+o+BJy9sKK/1SVgTXEniWiv4S0QXArOblYUN9P9SQvGAbs34TnIPn6JZCLeY6zXoTVYY+INSc4WRisegdzrEbkAWA08ReysDYDIVdHy9wFXAuNRjVirlnc37as+Syxcxdi3rKJBleru1bDlrvXNiAbNDlszWAUPI478BjyfoT60Bp+JoZkQRwHkXMauGRGsjHyuYZZxDMDZEfsJtj3gm4qW3xcNq/hXYCKx07/x08YTgRstISm7B/zT63MDiOwEMhMPaQTcCQ1gtHoS+D4GAtFwkj4IrURL6xcHxdAwTMtkL8ORdR7qPsBamh+pQDKOBe/RSNS8Js5uqGs4hH6wOXsPLRbfVOto7vWIHB7jI7Gmf7X0Nsi/I3qsIUISJVIGrr4QWk3tuIlkQMZpiMtmGrh6IokrhkMQ/BaNVFirkg0pw4jJXoi4hyDuIcmPF9yLltwOoZqlVC7Ivtzypfin7v7ivqlWaMes8WjlMxD4fNexwExrLVD2JVYCLv+HDax5NTj3scZnfB9a3bPMM8CbJB5s0jEiR3QtkhGTVGLExJCAOHKRwofQSInVGnD2QMRleVT8M0iM21oH/3uo/z37Y4FZqK3jtT6Ewf8hkncdknX2not7johGi6szzuIs3BXNzZAyzJiJIQENLCBSehdaepe1OC46OCruIVHHbENxg3MwUN/EWrt7W0Z20+KIR7J/DY48djl8XUAGknOTmdFJA6ZlYogjUvUmVDyBFV9E0eAPUP02mnsrVL8VDWLUQByFSId/II58Itt+Gr8a2P4Eks7uODqC5NbrtuLsBIWvoNWTIfA9uHojmWcjrl4Nqb2hnrRJMQlGAhRVrafAnUOOu0NLV6fdoJGKqJDEdmN8EFoLOy/GGnxNFr7RleSYG7J+gTjyrZkeya2HmOwmRKTrAIhsBWfXPVzDQhy5SPYFkH1BvcoboCK4s1HntcluTnXAwaKNmXyxdiVFVSbBdcoILkgSksCPtT4n2YfcgdWVsDlX3IhYng6t+JeV6a8pBD5Bd/wC9X2x57KGBlFUtZ6iqvWUhm2m3+tBmxQTV9hBx9IMNhfls3DjFoqq1lMR3NloRTVEkRwaZlarOS8XCh4g6bJ/7xg0Um5FdGtIhHtbgoAfLb8TVbsFiIbGUFS1ntUlJRSVFbK+pJ6J2erQJrs5Xo+Lffp0orA4k1lFqwhpkJzsahzsoH9BBb2yeu/5IoZE3EOj3ZBqGiQq6kMkG829EcrviWndOJD8exBHJhpc2oBUpPUksABM7ptGE9uqX11SgtsxgEF59es+2tEmxaSGwsJsjmYAK7fuJAiUhDqwrWAtDMIISiMQcUDBw2jJtaDlgERDFAyF0CKSTgmLA0JLcWSejHqPsAY7xQ2eEbvSZzi7WuEgbfFEhaaqzn4XluO13OYcTdIlM9SHoqr1rCxx4o0uK3A7CpskJNDGxQQsQSks3LXc/MNF8M3ytazuXAJAQYaX4YWDW6h2bQ9x9YWOE6zxk0gFeIaDZFhBjwKfYb/QTqxgRIA4siCaJziuhKMA9Z4QzRwYK0peKHgCETdadhOEt4OIFfQ57w6IlKDld0ZNZrEX9FotKUO9qAjupCRYUbtdIyRNFZBY2ryY1OXEofswZ1kG/ujf7Qd+hCEYQWkAIg7wHBC/r+BvRAKLoeRK4sXACY4O4Nlzd0Py/oJW5IPvbUscnH2RnGsRz75WgcLXrLiyGgZXX0QcqCoExkbNZ4A4AUHy7jPZ+OpJRXAnq8u3oY5OZDoL2FFZjddBSoUEQKwIAG2LEQeO1I/e/6ReZecsK2IzP3LgkF3eAiMsjUcDi6zEWeHNgIJrOJL/V8TZuf7XUAXCteuB6nVOaC0E5wJ5kHGkCQq9B2InI1aXb6PYl0NX767Qyp0zdx88qlte7lxVHdWQe7a7lkldDhnciznLYOl8q5lcJlsp2Wc+R/UY3sI1a5uIZyjS8VU0vAPEY8WJbeg1RGjoW09cfa1FfoY9EtsSAdgZFZI9CUhTafdiApag1LBiXQE/rFgGzKd/QQGAMb81AnGatS2tCfuWiPU/6urdc0skFewVYhLLPn06wTpYsGwTW6PxcXJyVzKm70AjKIY2SUu1ROqS7oDSvYGXgG5Yiy2eVtV/1ilzLDAZKzwXwFuqWjfodEqp8agAFFdUs6DIxxesZFgPa0DKtFQMrZ26LZHKQE8K3FZIheZqidQl3S2TEHCdqs4TkVxgroh8qKqL6pT7TFVPTXNd4qiZTi4szLZaKkUg0SheGhUW41UxtEZqzGbi6MwOXxVVgZ4pn5lpDOkOKL0J2BR9XS4ii7Gy9dUVkxalputDqbW9oDwf2AI9qFV701IxtAZqbO9uxwAsz4+3VQgJNOOYiYj0A0YCs20OHy4iP2Al3/qTqi60Of8y4DKAXj1Tv4R8nz6dal8bm76hNZFq23u6aBYxEZEc4E3gatW6GaCYB/RV1QoRORl4G0hIcKuqTwNPg+UzSWd9a2z6xWXVUAarKyLGpm9oEdJhe08XaRcTsRZnvAm8oqpv1T0eKy6qOk1EnhCRTqraxLXqTSPWpr8PnYxN39AsNIftPV2kezZHgGeBxar6cJIy3YAtqqoicihWWIQd6axXYzA2fUO6qRGS4kB2Wm3v6SLdLZMxwK+A+SLyfXTfTUAfqM3sdw7wexEJAdXAeG2lHv9Y89ucZfDDYktQajDCYmgosVO8JcEKVpY46ertSDhsDf63xBRvY0n3bM7n7CGjk6o+BjyWznqkA2PTNzSVBLNZlbNFzGapYq9zwKaS+JZKBstXrMPY9A27ozXY3tOFEZMUccjgXrAMVi3baWz6Bltai+09XRgxSSGHDO7FwGJLOIxN3wCt0/aeLoyYpBhj0zfU0Fpt7+nCiEkaMTb9vZfWbHtPF0ZM0kwymz7sxCXFpqXSjmgrtvd0YcSkGYmz6QOrKzx8U2ls+u2BtmR7TxdGTJoZY9NvH7Rl23u6MGLSwhibftujrdve04URk1aAsem3fpLZ3gs92RTm5bVgzVoPRkxaGcam3/pob7b3dGHEpBVibPotT3u2vacLIyatnLo2/RDl9OhawsjOxqOSLtq77T1dGDFpA9S16Vt5f4roX2CN2pqWStPZm2zv6cKISRshwaa/bBPVvYxNPxXY2d7bWiyR1oARkzaIsemnDjvbuxGSxmHEpI1ibPqNZ2+3vacLIybtAGPTrz/G9p4+miM6/Vjgn4ATeEZV761zXKLHTwaqgItUdV6669XeMDZ9e4ztvflId3R6J/A4cCJQBMwRkSl10oOOw8qTMwgYDTwZ/W1oAsamv0tIYlsixvaePtLdMjkUWKGqqwBE5HXgdOLTg54OvBSNSP+1iBSISPdoalFDE9gbbfr2tnfjEWkO0i0mPYH1MdtFJLY67Mr0JJqjuIZ0pwdt7+wNNn1je29Z0i0mdmku6ubEqU+ZZk0P2l5pjzZ9Y3tvPaRbTIqA2OmEXljJyRtaxpBiamz6W5b7qCrIpCSypc3Z9I3tvXWRbjGZAwwSkf7ABmA88Is6ZaYAV0bHU0YDpWa8pHk4ZHAviosrASiu8LYJm76xvbde0p3RLyQiVwLvY00NP6eqC0Xkd9HjTwHTsKaFV2BNDV+czjoZ4tmTTb815f0xtvfWTdp9Jqo6DUswYvc9FfNagSvSXQ/DnrGz6X8RFZQaWkpYjO299WMcsIY47Gz6M9XqWjS3Td/Y3tsWRkwMSWlJm76xvbc9jJgYdkusTb+wOJO5m9Nj0ze297aPERNDvSkszOZgurNye2pt+sb23j4wYmJoELEtFWi8Td/Y3tsfRkwMTaIxNn1je2+fGDExNJlkNv2CDC8AHTIya41lYGzv7RUjJoaUUtemXxGpIiMnSO/Ou5ZTGdt7+8SIiSHlxNr0oYCVW3dS5c6kIDcLMC2R9ooRE0NaiB2kLa6oJs+bRaERkHaNo6UrYDAY2gdGTAwGQ0owYmIwGFKCERODwZASjJgYDIaUYMTEYDCkBCMmBoMhJRgxMRgMKSFtpjUReQA4DQgAK4GLVbXEptwaoBwrFl9IVUelq04GQ2MJhEJUBAN0yMjEymhrqEs6HbAfAjdGg0rfB9wIXJ+k7HGquj2NdTEYGkUwHOaRr7/k3eVLAPA6XfTOz2dndTW98/K5eOTBjOjWvYVr2TpIWzdHVT9Q1VB082usfDgGQ5vi/i8/Y+qKpQTCEQLhCOWBAIu2bWNTRQXfbNzANe9PZebq1S1dzVZBc42ZXAJMT3JMgQ9EZG40BagtInKZiHwrIt/u2GEaMYb0U+kP8MHK5fhD4aRlfKEwD8/+HCvJwt5Nk7o5IjID6GZz6GZVnRwtczMQAl5JcpkxqrpRRLoAH4rIElWdVbeQSQ9qaG6KfdU46jE+UlztozIQIMfrbYZatV6aJCaqesLujovIhcCpwPGaRLpVdWP091YRmQQcCiSISVtAVdmwfBP+6iB9h/bC5Xa2dJUMTaBrdjaRerQ4nA4hw2UW4KdzNmcs1oDrMapalaRMNuBQ1fLo65OAO9NVp3SyafVW7r/oMYo37UQcDpxOB7//x0WMOunAlq6aoZF4XC5OHzyECYsXJi3jdTk5ffBQXE7zxZHOMZPHgFysrsv3IvIUgIj0EJGaDH9dgc9F5AfgG2Cqqr6XxjqlhXA4wp0/f4jNq7birwrgq/BRWVrFo1c8w6bVWxPKl24rY/bUeSz8cinhcMT2mvM/X8J9Fz7GTafczZuPvEtlma0eG9LMNYePoXNWlu0xt9PBTwcO4g+Hjm7mWrVO0tYyUdV9kuzfiJVbGFVdBbT5r+6FXyzBV+FLGIQLB0N89MosLrjlnNp9E//xLm8/Nh2324WiZGZncPPr19Br0K7pxWnPzOD1+yYTqA4AsH7JBj59/Uvu/eAWsvPs39j1IRJRRDA+iQYgIjz3s7O45ZMZLNq+FacIOW4vlx18CIf26MmLP3zH2JdfJEyEI3r35ZrRR9Ape1cQqPlbNvPOsiUEQmGOHzCQI/v0bbd/f9PRSwFl2ytsR/PDoQg7NpbUbv84axHvPPE+IX+IkN+aNfdX+rnvV//i0a/+johQXeHj9XvfJuAL1p4X9Ico2VbKBy98yplXnRx3j02rt/LfOyew8IslZGZncNJFx3H6FWNxunY1OtcsXM9zN73K8u9W4/a4OPa8MVzw17PxZHhS/JdoP6gqVYEgWR43nbKzeerU0ymurqI6EKJHXi4Al0x5i5XFxQQjVuty5prVzN+ymTfOGY/X5eLJb2fzv4XzCYYjKDBr3RpG9+zF3cef1C4FxYhJCtj30IGEbKYPvVkeRv5k/9rtD178FH+0tVGDKpTvrGDVD2sZOKIfq+evw+V2xYkJWIIyb8aPcWJSvKWEm0+5m+pyq1Xkrwrw9r+msXHlZq589BIAtm/Ywe1nP4ivwkpFEfAF+eT1z9m6bjs3/PcPKfsbtCW+XL+W/8ydw4bycvoXdOD3h4yOM569uWgBT8/7lqpgkEyXi0tGHsx5w4ZTmJkFmVaZeZs2srakpFZIAMKqVAQC/HP2V3y6djUlPl/cfatDIWZvKGLOhiIO7dU8+ZqbE7M2JwV07tWR439xJN6sXd/07gw3Xft24bCf7VodUBXN2VsXEaG60nrj5XXMIRgM2ZbL75IXt/3ecx8TqA7EtYoCviBfT53L9o3FVpnnPyXkjxewoD/Ewq+WsmnVlgY8ZftgxqqV3PTxhyzZsYPyQIAft27h6venMnfjBgCmLFnCY3O+pszvJxSxTGr/nvsNkxYvirvOl+vWUh1K/D9Vh0K8s2xJgpDEHp+5Zk3Kn6s1YMQkRVx053h+9/BFDDlsEP0P6MO5fzqNOyf/BY9nV+PvsFNH4clM7FoE/UF679cTgF6De+By2c8M9B7cM257+dxVhIKJLSKPx03R0o0ArFmwjlAwcZDX5XKxaVXi4HB7RlX51zdfJpjQ/KEwj835GoBnv/8WX53jvlCYZ7+fW7td4Q8wcfEC23sIVgtld5QHLXEv9/v5fO1avtu0kUjEfiC+LWG6OSlCRDj81IM5/NSDk5Y5dvwRfPL6Z2xYsZlAdTD2ZK4ecws3vnIV3Qd2I+gP2p4/840v+WbaPKorfBx80oF06duZZXNWJswIBYMhuvbrAsDAEf1Y8s1yQoH4D0goGKTHIDu/YfslGA6zvcq+dbhqp9WS215lP2tWXF2NqiIifLR6JYL9mEd93JQz167iH19l8vbShbgdTlQhy+3mkXGnMLBDYb2epTViWibNiMfr5o7JN3DCL4/G4dz1pw8FQlRX+Lj/osep2FmBI0nLpHhLCRtWbKZ4cwmfvPoZ3300H6cnvqzL42LIofvQvb8lJj+96FjcHjex433uDDcHHjuMbn07p/4hWxnVgSDBsCWkbqeTbJfbtpwAs4vW0Ss3z/a4y+Hgg5XLAdjpq8YftrfYd8/JSSIzuwiEI0xcNJ9AOEJlMEhVKMj26iqunj61TbdQjJikmYA/iK/KX7vt8bjYtmEHERt/STgYZvvG4rixlzhivvZCoQhVJZV4Mz24M9wg4HI7OfLMQ7n2md/Xlivs1oG7plzPsDH74XI5ycrNZOzFx3HVk0mXQbULFmzZwgVvTeDEl5/n+Jee4/ZPP6IqEOSXB4zA60x82/vDEa774D2y3C48NsdDkQh3zPyEI5/7DxMWzrdtgThFuGjEwXiTfBnEYicZVcEAP2zZXI+na52Ybk6aKN1WxpPXvciPMxehGiGvUz6HnTKS484/ki1rki9U1LDym3su4PE/PkvQF0JVcTiESCTx7RsKRSiPZs7zZHrYf8x+/PbB/0uYduw1uAe3vH5Nah+wFbOxrJyr3nu3doA0rMrHa1axtaKCx0/5GcuLt/PR6lUJ54VVWV1SwjlD9mfy0sVUBuO7mwqENUJxksHVsCqTFi/ghjHH8PL879lQVoYvFKpX1wesrnJFILDngq0UIyZpIBJRbvnZfWxbv0s0SreW8v7zn/LBCzMRh31DOOgL8Nq9kwj6goQC4dpZGq2HJSFQHWDBF0tY+f0a9hnZPyXP0VaZsGg+oTrdhWA4wqId21hTUkKuJ7m/xh8O89naNWS4XAliUh+W7NjB7TM/JtttdafyvF58oSD+mJao22G1fIJ16xgJt+nYKKab0wjCoQjffvADEx56h0/f+DKuGwOw4PPFbN+ww/ZcVbXt4gBEFFb/uI6iZZviymi4ft9t4WCIJd+sqOdTtF9Wl+xM+KACuMXBhrIyOmZm43Ikf+tvqCgnkGRMpL5UBoNUh0KU+f04xEGmy4VDhK452fz16OMY3KlT3OJAt9PB5YccRm4bXnlsWiYNpKq8mtvPfICt67fjq/TjzfLy8t8mctfkG2oHPZfNXYXadEuSISK7FZn6Eg5FkrZ69ib279yV7zZvJFDn7xmIhBnQoZB9Cgt5ecH3hJL8uSOq+Gw8JI1BsbwlApy931CuPeJIRIRj+/Xn33O/4Y2FC2pXJk9ZupgxvfrQKz8/JfdubkzLpIFMfOgdNq7ajK/Sao34q/xU7qziiaueqy3ToWvD3gypDKwz6dFpVO3liwLPGjoMr9MdF4vE63JydJ/+9MjLpVtuLvf+5KSkMzuQ2AVpKgq8tXQxS6OBvXZUV/Hm4oUEIxHCqgTDEVbv3Mnl06YQbqMzOkZMGsgXk+ckeDZUldXz19Z+iI88czQttfTCV+ln2nMft8zNWwmFmZk8f/pZHNOnH1kuN52ysrjwgIO4/difAOALBumRm8db5/2CrN0ISl0O69WLTklWEI/q3pNuOTm7PT+iyjtLrFiyU5YuJlyn9apAVSDInA1F9a5Ta8J0cxpIfUSiqrya7PxsKkoq91i2Y49CDj7xAD5782uqK+JnCfI753HYKQfxxeQ5VOzc87XA8qy89Y932V60g0vv+xVOm2nOvYGeeXncfcJJcftUlafnzuG1BT/gQAhGIgkDtcm47vAxnDN0f4qrqrh82jusLS2pPXbpyFFcctDBlPv9XPXeuyzZnny2rmZQd2tFhW3rJ4Kyvbpttiz3zndaEzjyrNG4PPEaLA5h4Ij+ZOVl8e6/P+APo2+ol5B069eZu6ZczyV/P5+bXruazJzM2mOF3Qq4Y9JfuPhv53PrhOvI67j7b71YImHly8nfMv2ZGfV/sDZGud/P3I0bWFuys97nTFi0gFcX/IAvFKYqFCIYidRr2rZGSAAKs7J44uTT6JtfAOwSEoBcr5dHx57Kfp06Jb3WuH0GAXBIz95k2kRnC6tyQJe26Uw2YpIEVeWrd+dy+9kP8peT7mLSo9OorvBxzjWn0nvf7mRke3E4HWTkeMnrmMPl/7yE1fPX8cYDU2zXwtSlW7/O3DrxTxR2KyASVgaN7M9Nr/2RzJxMCrsVcOvEP9GtX2ciYaXPfj255X/XNkhQAtUBpj/b/ro7qsqzc7/l1Nde4oYZ73Ph5Df5zTuT2Fltb5OP5aUfvt9tcGg7aoTEr8pvsOKJ1gjKnw8/kksOOpiwKldiRUzfnaAIsE/HQtaXlvLl+nWEIpE4t2yGy8mJ/QfQp6CgQXVsLRgxScJ/75jAU9e+wJLZy1m3qIi3/jmVW067B4fTwW1v/pmf/OIo+u/fh1EnjeCe6TfTtW8nZk38klBgz7MAsUKy/qsiXjr+RUrXlTJoZH9ufu3qWiHZ/OMWXjzuBbYv3dEoQakqszdXtWU+XbOalxd8TyAcoSIYxB8Ks3T7dm7++MM9nrvTt2fBieWiESNrheQSEaYDF7NLUM4aOswSEhEmAb8jXlDqRmjzuJwUV1dx0dsT+Wj1ytqWkQC9cvP40+FHcvPRxzWojq0JIyY27Ni0kw9emom/apcbMegPsX1DMR++NJM/H38HH70yi5U/rOGbafO47tjbWbu4CL8vaOtUjcXlcnLTa9fUCsmkX02ieHkxb5wzgdJ1pewzsl+tkLx53kR2rtzJxHMn1ArKdTFW+T0RDoX44KWZ7SoNw2vRbkosoUiEhdu2sL1y911LZwNHxWvW6RSLUOPe8WEJyieALyokU2rqAdSsLc71einI2NVtdTsd/GzQfjzx7TdUhUJxK4trpo9PHrRvmw6alDYxEZHbRWRDNP7r9yJycpJyY0VkqYisEJEb0lWfhrBs7ircnsT+rL8qwHvPf8LOTTtrhSbgC1Jd4ePxPzzH6JMPTr6uJgaH03rDREKR2g96+YYy3jhnAlsXbWPdF+t587yJ+Mus6edIOFLrW4ldBJjTIZtz/3I65/75dDyZ7rjFgzV1e/muibzx4BTaC6U+v+1+lzgoC9gfA2vVr3s3RjU7Hvn6SxZs3UJ3YCK7ssj5gAuA/WOEBOBS4Jbo63s/n8ny4l3GxVyPlynLlvB1kf1MTVnA1+CWU2sj3S2Tf6jqiOjPtLoHRcQJPA6MA4YC54vI0DTXaY/kd8xFbYbmnE4HJVtLCdm4ndYt2UBuYTYHnXggLk/yhV6hUJi7z3+Ekq1l9D2qD2e8cAbODEu4yjeU8fKJ/2XiuRNqhSSjQwY/n3AunYd0YsPyTTx0yRO119KIcuYfxnHWH0/m3vduoVPPxOXrgeoAU//9YYJLt4ZgMMSMl2dx+9kPcvcv/8nsqfNadUvm8N59bEXB6XDQJ8/e31Pm93Hh2xPxNdDVWhEMcs17U1mwdQu9sQQlNqJM7Ef/UuD26Ot7P5/J5KVL4q5V4vMlXWkMlt1/wsL5Dapfa6OluzmHAitUdZWqBoDXgdNbuE7sN3oQuQU5CW5Sp8dJZm5mkrPghVvfYNvabQk+lLpsXLmFO3/+YJygOGwEKKMgXkjuOvdhSraV1R6vLK3iu4+tID09BiafAXC6nGxbn2jvD4cj/H38I7x0xwSWzF7OjzMX8cQ1z/PsTa/ttv4tyf8dOJJ8b0btyl7BMqT9+YijkqabmLhwIWV+n20OnHyvl0N79LRdSQy7BGVbZQW9gT/blOkI3Ba99nPz5iYICVCv/DuvLpjP5+vW7rFcayXdYnKliPwoIs+JSAeb4z2B9THbRcSLfy3NmR7U4RBueeNaeg/ugSfTTWZOBtn5WVz12KWccMHRSc9bNmcFK75fU697bFy5hceuehaAvkf1oe9RfRLK7HfmEDoP6UQ4rDz46yfjhKSGZ65/mVf+/iYbVmymx8CutvcKBUN07F6QsP+7j+azesG62ij4YHXlZk740jZFR2ugMDOTl8/6Of93wEj279KFEwYM5IlxP+PEgbbJEACYs7EowVoPkOlyceexx/PI2FP49chR5LjtDWzn738gnbNzCKvyrs3xHcCs6FjHGfsNqZ02biiBcJjXFvzYqHNbA00SExGZISILbH5OB54EBgIjgE3AQ3aXsNmXLPPf06o6SlVHdeyYfB4/VXTt04n7Z9zKfe//lVsnXMe/v3+QUT89kDP+MC5hbKIxFHYr4Df3XgDA5h+3sP6rxL708qnLKF1XitMpXP7Pi+N8KDUUby7hnac+4PqT7mTgyH64MxLHegaM6EeWTYqM7z9ZiL8ysfsjYqXvaK3kZ2Tw64NG8Z/TzuTO405gaJcuuy3fPTfXNs2nonTKzkZE+NWBI3nrvF8mLACs8ZHUzNokc+5cQvy0cWMFZUNZaaPOaw006VOhqieo6v42P5NVdYuqhlU1AvwHq0tTlyIgNkx3L2BjU+qUaroP6Er/4X1qU316vG5OuODouFQSYBnX6rvILtZHUjNrE6pKXO5eubWydpYn1oeSgEIoEObtR6fTrW/iB2v1j2tZ+OXShP35nXNtU5g6nU5yO9R/Crq1c96w4bjrfAG4HA76FXRgQEyYxIVbt8Stizl36P5xQhI72DoeeIb4QdnYaeMnTj6NgoyMBte1Q2bybnRrJ52zObGBGc4E7CLwzgEGiUh/EfFg/Y9a/dTD+TeeSe8hPcnI9uLyuPBmeayVv/VYKexwOrj5tavjhCR2sPWsV85i3OMnxw3KxgrKNf9OHiEtHIqwcUVipK5AdZB3nvwgYf8xPz/CtpXldDoYecLwPT5LW2HfTp258pDDcDsdOERwirBfp048fFL8BGMgFI5rFu/fxeo2bhLhy5j9l2I1s8eROMszPfq6MCuL7ntYq2PHif0HNvic1kI61+bcLyIjsLota4DfgpUeFHhGVU9W1ZCIXAm8DziB51Q1eWLXVkJmTgb3TLuZRV8tY93iDcx4eWa9I707nA5yot/6FRvLCUZbJLGzNgBZHbN4+6K3CftC+Ep9VG2vIr9PfkK6i7okk7PijYm28y69O3L1U7/lsT88i6qiqmTmZnL9C1fi8dZ/AVxrZ9HWrTwxZzaRiBJRxSXCyuJirnrvXdaUlJDj8TB+2AEc269f3Hn/+PoLBnTowMDCjkwEzsH6Vrw9eryotJTe+fm1x8YAd6uCCE/M+ZrFu1mjk4w1pSW1gavbGtKapwGTMeLAkfrR+5+0dDUAKzzjFaNvrJfztYb+w/tw8+vXkJOfxfJpy/n45o8569WzrVmbFZvJzsuioEseaz9bx/Qrp/GzZ39Gj1E92Fa0gzt//nBcBLf64HQ6+OklP+H/bvu57fFgMMTK71bj9nroP7wPjhTHRFmxbjt5PXIpzM/ec+E0cOHbE1m2wz5YVQ0ZLifH9x/I1OXL4vZ3yMjgX+NOZWBhR7YCNZ3Iez+fyUerVvKPsaewf5eubAM6RUXgiTlf898ff0i4h7Dn6PUZLhd3HXsCR/btW9/HSwvd8nLnquqoPZfchRGTRrBl3XZeu2cSCz5bjCfTQ+m2UsLJIu0kIVZQgtUh3Jmu2unfrLxMbp3wJwq65NUea6yQgDU79dDMO2uDNzU3LSkmwXCYY198tl5Ts8mIFRSA+z6fydvR6d8ct7tWUICkQlJT1h8O7zFWypF9+vLAiWMbXd9U0BgxaWmfSZtj59ZSbhr3d2ZPnUtFSSXFm3Ym5K0BK1L8yOOH061/ZzKyvQmDs6vnr+Pv4/9BRWlVnJCUbCuL86E0VUgAnG4n30yb16hz2zpOEVxNbGnt9Pn4w/R3mbJ0CX/9eEatkMAuH8qkRQv5+2efJhWSmrL1EbVgE0NGthQmnkkDmf7sR/irfPGDrdGXLreTUDCMN8tDQec8rnj0EjKzM/jtyD/VRmaLZfX8ddxx5gMcdvooPnjhU0pjfCQbV27htjPv55hzj+DjVz9j+4biRtc56A8xe9o8Tr+iZb/tWgKHw8HYgYN5b+UyW69Jfdnp83HP5zNtj1UEg9z/1ef1us6esv15XU7G7TO4wfVrDRgxaQCqyldT5tqGGMjI8TLqpBEoMOzwwYw541C8mR7WLFhPKEnuYID1yzZS9NA7tjNBW9Zu440HJqek7rkdWma8ojVw9WFHsLWygu+2bMItjnp1NVqKEV27c8KAtjmjY8SkAXzyvy/YsdG+hRAORTjn2lPp1i9+XMLldRHZQ3T5GiFxeV1EQpEmB5auizfLw08v/klKr9mWyIyOa6wvLaWorJTuObn8ctKEJo2jpBqnCCcPGsyNRx7TJmdywIhJg3jrH1OTftD3HTUgQUgAeu7TjcJuBWxZs5U9vXdD/sZFRBeH0KFzPhk5Xjav2ZZQx+N/cRQHHd9+fCONpXd+Pr3z8wmEQjhEbMWkPjMuDSXZvWrwOB1cf8RRjDMhCPYeSraUJD129b9/Z7tfRPjz81eQ1ymPzJwMvFkeXF5XylJSdB/Yld/c80tueOUqthUlph11uZz1C1y7F+FxuRjds1eCdd7jdHDB8APIS3Humpq72Fn6AZzioGNWdpsWEjAtkwbRY5/urFuyIWF/Vm4mEx6YzPCjhjLyhOEJQZx77tONx7+5lx9nLaJ0Wxn7HjKQ2895kLJt5U2qj8MpnHPNqYw541C++3gBbo+bYJ3WTSgUZs2CdU26T3th3qaNvL7gR7ZVVjCiWw/Wl5ayvaqSmkiwQzt15dcHjaJnXgGPzP4iIQhTYwnVZGZM0jqpDoX4bvMmRvfqbXu8rWDEpAFccOs5PHjJEwR88eto/FV+3n/hU2ZO+MoKrzjhOjx1giu53M64rsYFN5/DMze+TKC64Skoa/BkeOjQtQCwQhAEA4nXcrkd9BueuCJ5b2PS4kU8+s2XtQKxsmQn+R4vtx3zE0r8fgYVFjKks9VNPX2/IWS4XTwz91u2VFUQCtcv8PSeSHYNr9OZNIVGW8J0cxrAAUcP5YaXrmLQQQPIzMmoXddS4zPxVfpZs2g9H/3XfgoxlqPPOYyrn/otAw/sZ12rgd0eESG3Qw77jbainXft24kDjx2GOyPeBh8OK4MPbpuzA6nCHwry6DdfxbU0guEIpQE/i7Zv5Wf77lcrJDX8dOAgJpx7PrMuupSzhwwjw5U84FVTcYpw4oDkIRTaCkZMGsiAEX259c3ruHXCdXgyE9evBKqDfPbW7Hpd66Djh/P3qTfSY2DX5LFjHda0s8PpwOF04HI7cXtd9BzUnVv+d22cCF315GUJLleNKE9c/TxrF7fNxE6pYGXxTtv4r8FwhC/Xr7c5I55rDx/D70eNpntuDlk26SnsSDY+0iEjg+65OWS4nGS6XHTOyuKRsaeQ34gVxq0N082pJ2sWrOepP73I+sUbEIcw9Ih9k075uhu4SK6qPHnsT5fbhdvj4u53b6RL3y6sW1xEZk4G3QckBkLyV/rZtGpLwv6gP8jb/5rOH5+4tEH1ai/kZ2QQUvtZuI5Ze17yLyKcO2w45w6zuqm3fDyDz9at3q0JzukQPOIgGI7UGtW8Lic3jDmGo/r2ZW1pCWFVBhR0aPMDrzWYlkk9KN5Swu3nPMiaBesJhyOEgmEWfrGEiM0AnTfLwwm/TB6NzY5Dxo60jSsC1nRxVWk1/3tgCi63kwEH9LUVEoBt67fjcid+P2hEWbcXt0x65uUxuLAjrjof2gyXk/P3P6DB17vtmOM4b9hwctxWPuPeuXm4nQ6yXC6yXG68Lid3HnM8L55+DicO2Ic+efkc3qs3/xp7Kkf26cPz383jd+9O5sJJE7n0nUks2to6o9o1FNMyqQcz/jsrYVVwOBTBneEiKyOTiGqtsBw6biRjzrSLA5Wc0y7/KV9N+ZbS7WUJg7tgjcnMm7HncH5d+na2dds6HELfYXv3IOy9J4zlLzPeY0XxDtziIKQRfjdqNIf2bPgMitvp5PJDDuPyQw6r3Vfu9/NV0TpEhMN79iHHa2UpuO3YeLPg/Z/PYvrKZbXjNwu3beOK6VN4/vSz6VdgF9m07WDEpB4ULd1oG2LA5XJx8d/G48nwULajnP1GD6L3vj0afP3cgmzun3ErH7/2OS/fOcHW3FafUJE5+VkcO/5IZr7xZVxcV5fXzZl/GNfgerUnCjMzeea0MykqLWWnr5p9OnQk05O6mC25Xi8nDRy02zKlPh/vrlhKsE73KBiO8NIP33HrMW3bpWzEpB4MOngA33+yIKHVEA6F6T+8D70GN1xA6pKZk8Epl57A8nmr+Xb6vLh0Gi63kyNOq99q8IvuPI8OXfKZ/uwMqkqr6Te8DxfdNb5RItce6ZWfT698+5QYe+K7TRuZsHABO33VHN23H2fsOzSpIBWVlrJ0x3a65+QypHNnRIQN5WV4HM4EMQmrsrQRgZRaG0ZM6sFx5x/JO0++TygYrnWYujPc7H/kkJQISSyX/P181i5az84tJYQCIVxuF116d+KCJIGN6uJ0Ojjrjydz1h9tc54ZGskbC+fz5Leza7sni3dsY/LSJTz/s7PiBCUciXDHzE+YuXYVbnESQemVl8+j406hR24uwUjiOJtDhEEdOzbbs6QLIyb1ICc/i7un38wrf3uT7z9egCfTzQm/PJoz0vCBzSvM4aFP7uDHmQvZuHILvfftwbAx+6U8+pmh/lT6Azw+ZzaBmDgj/lCYzZXlTF62mPExg7j/Wzi/dqYngPXFs7pkJ3fO/ISHf3oyJw0cxIerVsQlUPc4HfzqwJHN90BpIm1iIiL/A/aNbhYAJao6wqbcGqAcCAOhhkZ3ai469ShstqlVh0MYcdz+jDhu/2a5n2H3LNy+FY/DEScmYAnKrLVr4sTkzcULbXMhf7txA5X+ADeMOZqOmZlMXLSQqmCQQYUdue7wIxkYEyW/rZI2MVHV82pei8hDwO4Sghynqm2/02hol+R7vbZBjYTE1BTVyWLXCPgjYbIdHn43ajS/GzW6zQaOTkbafSZi/bXOBVpvzkmDYTcM7tiJzlnZCa5Wr8vJuUPjW49H9ulj67btlp1Dhzou1/YkJNA8prWjgC2qujzJcQU+EJG5IpI8KYzB0EKICI+MPYVeuXlkulzkuC1j2hWHHMaB3brHlf3twYfQISMTb3Qtj9vhIMPl4pajjmt34lGXJnVzRGQGYJcx+2ZVrYk3eD67b5WMUdWNItIF+FBElqjqLJt7XQZcBtCrZ6+6hw2GtNI9N5fXzzmP5Tu2UxrwM6xTV7JspoU7ZmXz2tnn8c6yxXy3ZRN98wo4e8gwuuXmtkCtm5e0proQERewAThYVffo5xaR24EKVX1wd+VaOtWFoWG0dN4cQ8NpjakuTgCWJBMSEckWkdya18BJ2KcRNRgMrZx0i8l46nRxRKSHiEyLbnYFPheRH4BvgKmq+l6a62QwGNJAWk1rqnqRzb6NwMnR16uAA9NZB4PB0DyYEAQGgyElGDExGAwpwYiJwWBICUZMDAZDSjBiYjAYUoIRE4PBkBKMmBgMhpRgxMRgMKQEIyYGgyElGDExGAwpwYiJwWBICUZMDAZDSjBiYjAYUoIRE4PBkBKMmBgMhpRgxMRgMKQEIyYGgyElGDExGAwpwYiJwWBICU0SExH5uYgsFJGIiIyqc+xGEVkhIktF5KdJzi8UkQ9FZHn0d4em1MdgMLQcTW2ZLADOAuKSZonIUKzI9MOAscATIuK0Of8G4CNVHQR8FN02GAxtkCaJiaouVtWlNodOB15XVb+qrgZWAIcmKfdi9PWLwBlNqY/BYGg50pXqoifwdcx2UXRfXbqq6iYAVd0UTRFqS2x6UMDfqXuH9pisqxOwvaUrkSba67O11+fat6En7FFM6plPOOE0m31NykOqqk8DT0fr9G1DUxe2Bdrrc0H7fbb2/FwNPWePYqKqJzSiLkVA75jtXsBGm3JbRKR7tFXSHdjaiHsZDIZWQLqmhqcA40XEKyL9gUFY6T/tyl0YfX0hkKylYzAYWjlNnRo+U0SKgMOBqSLyPoCqLgTeABYB7wFXqGo4es4zMdPI9wInishy4MTodn14uin1bsW01+eC9vts5rmiiGqThjIMBoMBMA5Yg8GQIoyYGAyGlNBmxKSp1v22gojcLiIbROT76M/JLV2npiAiY6P/lxUi0q4cziKyRkTmR/9PDZ5KbS2IyHMislVEFsTsa/BSlzYjJjTdut+W+Ieqjoj+TGvpyjSW6P/hcWAcMBQ4P/r/ak8cF/0/tWWvyQtYn51YGrzUpc2ISQqs+4bm51BghaquUtUA8DrW/8vQilDVWUBxnd0NXurSZsRkN/QE1sdsJ7PutyWuFJEfo83PtrySuj3+b2JR4AMRmRtd7tGeiFvqAiRd6lJDutbmNIrWYt1PN7t7TuBJ4C6sZ7gLeAi4pPlql1La3P+mgYxR1Y3RNWUfisiS6Lf8XkmrEpM0W/dbDfV9ThH5D/BumquTTtrc/6YhqOrG6O+tIjIJq1vXXsSkwUtd2kM3p77W/TZB9B9Xw5lYA89tlTnAIBHpLyIerIHyKS1cp5QgItkiklvzGjiJtv2/qkuDl7q0qpbJ7hCRM4F/AZ2xrPvfq+pPVXWhiNRY90PEWPfbKPeLyAis7sAa4LctWpsmoKohEbkSeB9wAs9Fl1q0B7oCk0QErM/Rq6r6XstWqXGIyGvAsUCn6PKY27CWtrwhIr8G1gE/3+N1jJ3eYDCkgvbQzTEYDK0AIyYGgyElGDExGAwpwYiJwWBICUZMDAZDSjBiYjAYUoIRE4PBkBL+H4/rAblne3BzAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "kmeans = KMeans(n_clusters=3)\n", "kmeans.fit(x_train)\n", "cluster = kmeans.predict(x_train)\n", "plt.scatter(x_train[:, 0], x_train[:, 1], c=cluster)\n", "plt.scatter(kmeans.centers[:, 0], kmeans.centers[:, 1], s=200, marker='X', lw=2, c=['purple', 'cyan', 'yellow'], edgecolor=\"white\")\n", "plt.contourf(x0, x1, kmeans.predict(x).reshape(100, 100), alpha=0.1)\n", "plt.xlim(-10, 10)\n", "plt.ylim(-10, 10)\n", "plt.gca().set_aspect('equal', adjustable='box')\n", "plt.show()" ] }, { "cell_type": "markdown", "metadata": { "collapsed": true }, "source": [ "## 9.2 Mixture of Gaussians" ] }, { "cell_type": "code", "execution_count": 4, "metadata": { "scrolled": true }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAARMAAAD8CAYAAABUzEBbAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/MnkTPAAAACXBIWXMAAAsTAAALEwEAmpwYAAAx/UlEQVR4nO2dd3hUVfrHv++kzkwm9NCkiCIIKgoYRBdFV1kEC7oqWFjbytpW8SfuWh4VRF1ddtfeUOxtLYviiggq6KqIJBQBEQiIVLHQSU++vz/OTKbdSWYm9yaT8H7mOc/Mvafcc2+4X855z3vOEZJQFEWpL67GroCiKM0DFRNFUWxBxURRFFtQMVEUxRZUTBRFsQUVE0VRbMEWMRGRZ0XkJxFZHnKutYjMEZE1/u9WMfIOF5FVIlIkIjfbUR9FURoeu1omzwMYHnHuZgAfk+wJ4GP/cRgikgbgMQCnAugD4HwR6WNTnRRFaUBsEROSnwHYHnH6TAAv+H+/AGCURdZ8AEUk15EsB/C6P5+iKE2MdAfLbk9yKwCQ3CoieRZpOgPYGHK8CcAgq8JEZByAcQDg9XoH9O7d2+bqKooSoLCw8BeS7RLJ46SYxINYnLP07yc5FcBUABg4cCALCgqcrJei7NeIyA+J5nFyNGebiHQEAP/3TxZpNgHoEnJ8AIAtDtZJURSHcFJMZgC42P/7YgDvWqRZCKCniBwoIpkAxvjzKYrSxLBraPg1APMB9BKRTSJyOYD7AJwiImsAnOI/hoh0EpGZAECyEsC1AD4EsBLAGyRX2FEnRVEaFltsJiTPjxH1W4u0WwCMCDmeCWCmHfVQFKXxUA9YRVFsQcVEURRbUDFRFMUWVEwURbEFFRNFUWxBxURRFFtQMVEUxRZUTBRFsQUVE0VRbEHFRFEUW1AxURTFFlRMFEWxBRUTRVFsQcVEURRbUDFRFMUWVEwURbEFFRNFUWxBxURRFFtwVExEpJeILAkJu0VkfESaoSKyKyTNHU7WSVEUZ3B03xySqwAcCdRsBboZwHSLpP8jeZqTdVEUxVkaspvzWwBrSSa8uY+iKKlPQ4rJGACvxYgbLCJLReQDEenbgHVSFMUmGkRM/BtsnQHgTYvoRQC6kewH4BEA78QoY5yIFIhIwc8//+xYXZUIfvwRmDULWL48+TIqKoDSUvvqpKQkDdUyORXAIpLbIiNI7ia51/97JoAMEWlrkW4qyYEkB7Zrl9B+ykoykMD11wPduwNjxgCDBgH5+cAvv8Rfxi+/AGefDXi9QE4OcNxxwMqVjlVZaVwaSkzOR4wujoh0EBHx/8731+nXBqqXEovnnwemTQPKyoBdu4DiYmDRIuCss+LLX10NnHAC8N//mpZJVRUwfz5w7LHA9u2OVl1pHBwXExHxwGwP+p+Qc1eKyJX+w3MALBeRpQAeBjCGJJ2ul1IHDz4I7NsXfq6qCvj8c+Dii03LpTY+/RTYsMEISQDSiNMLL9heXaXxcXRoGABIFgNoE3HuyZDfjwJ41Ol6KAmyY0fsuDffBEaMAEaPjp2mqMi0TiIpKQFW6HbSzRH1gFWsOe00ID3G/zUlJcATT9Se/4gjANN7DcfrBY4+uv71U1IOFRPFmttvB3JzY8eXlNSePz8fOOooIDs7eC49HWjRArjwQnvqqKQUKiZKOFVV5rtjR+DbbwG3OzqN2w0cfjgwaRLw0kvGOBuJCPDhh8DVVwNt2gA+H3DeeUBBgRnZUZod0hRtnQMHDmRBQUFjV6P5QAL//Cdw333Ar78CPXqY41GjgJkzgXPPBSorgfJy002pqgLS0oyBNicH8HiAL78EDjqose9EsQkRKSQ5MJE8jhtglSbAPfcYIQmM3qxbZ7oi06cbQ+vy5cCzzwJbtgCbNwOffBJ0Qtu717RMLr0U+OyzxK+9Zw/w3HPAnDnGp+Waa4DevW27NaXh0JbJ/k5FhemG7NkTHZefDyxYEH4uLw+w8kDOyDD+I4l0YbZvBwYMAH76yQhSejqQmRkcLVIajWRaJmoz2d/ZsSPcFySUNWviLyeZ/5Tuv9+0dgI2l8rKYCvHalhZSWlUTPZ3WrcGsrKs4w4+GNgWMQPiggui07tcxlU+UcPq9OnGDhPJvn3A6tWJlaU0Oiom+ztVVaarEYnLBSxeDHTrZkZuli0z5++6y9g0cnKMEdbnM12f558Hdu82IzhffRVfyyLW0HNVlSlXaVKomOzv/P73ZiQmEtJ0O8rKjAH2+OONWOTmAoWFwL//DUyeDEydCqxfb0Z9OnQww7+nnAJ06QK8805wqNmK664zo0OhpKUZ/5TOne28S6UhINnkwoABA6jYwMqVpNtNGumoPXg85FNPWZfz0UdkZqZ1vpwccvp063zV1eTIkaSISetykV26kJs2RafdvZvcuJGsqrLt9pXYAChggu+ltkz2Z5YtM6Mw8VBcbFogkUydCpx8srXtAzBDx2edBQwfDvTta2YSv/OOiZs+HZg7N2i8ra42yxYsWRKef/RooG1b4JBDTItlutXKn0pjo34m+zOHHGK6MvGQk2PWNKmqMt2cLVvM8gLTpsWX/8MPg78LC4G//tXYWSK9Z0tKgAkTgJEjzfF55xm/loBYlZQAF11kRCg/P75rKw2Cisn+TL9+QP/+wMKFxjYSIC3NtFgCjmlZWcYrNi8POOAAYzuxcqGPl337jL0l1pB0UZH53rjRiEZo3QAjKFOmGH+UeKiuBmbPNqLUoYNxyGvfPvn6K5ZoN2d/Z+ZM043IyjIjOPn5wBdfALfcYjxSu3YFbrgBeOYZY1j98cf6CUmAigrjhm9Fhw7me9Mm62FrMn4fmPJy4KSTzJSAKVOA224zbv/z5iVVbaUWEjWypEJQA6wDVFaSZWXR53fsIE86iczIiM9QGwhdu8YXF2kA9njIp5821965k8zOjs6fkUFed1189/XYY6bMyDLy8sw9K5ZADbBK0qSlGVf2SE47zXQ1YnVJIklPB15/HfjhB+Dxx8PXNBEBnnrKxL3yirlm+/ZAy5amVSRiQmGhcbVv0QK48cbw4WOXyxzfdFN89XnxReuWVEkJsHRpfGUocaFiogQJjKoEeOEF0+WJPB+L9HQjEoEV2K66CnjssaBIPPUUMG6cibvgArN8webNZpQnI8NcZ98+M6kwP9/YSiZPNqJ02GGm+zNqlBGYDz4IeucWFZlrHXOM+V67NlinWKNV1dWxF39SkiPRpkwqBO3m2EhJiekyeL3Gz2PwYHLxYtMFaN06/m5Nejr573+bMnfuJG+4gSwuNsePP05OnWp+79tn4nbtMsevvWZ8USLL83rJl18Or+vjj5tuj8djQnY2OXGiyZ+WxhpfFa+XLCw0eZ5/3hxHlt+li/FzUSxBEt0cx198AOsBLAOwxKqCAARmIekiAN8A6F9XmSomNnLGGdF2iZwc8rPPrG0NoeKRnh48fvFFU97OneTRR5tzp5wSFBTSCMlJJ5m4Y44JCsqzz1pf4+qrg3mLiqztJwGHt8jQr5/JV1VFjh5t7iUri/T5yJYtg2KjWJLKYtK2lvgRAD7wi8oxABbUVaaKiU2sXRvbwHnllebls3pR09LIVavIu+4Ktgp27zZlRgpDQFBChSRSgHbsMC2KyOtcdFGwrvfdl7gROLTlsWgR+cAD5CuvmLootZKMmKRCp/FMAC/6b+ArEWkpIh1Jbm3sijV7Vq0yQ6+Ru+1VVADffWeGVD/+ONq7tXdvoFUrM7O3uNjYHy65xMzXufRSYOtWMwQLmEWPevc2r/fGjcEy7r8fGDvWXCvWkgOhk/3Ky+O33QT4/nvjHwOY+T5HHZVYfiUxElWfRAOA72G2AC0EMM4i/r8AfhNy/DGAgRbpxgEoAFDQtWtX25V4v2TdOuuWSWYmOWECuX07OWhQdHx6OtmnT3RL4eyzyYoKU/Y998RuMdx/v0lTXk6eeaZ1GrebfOaZYF2XLYt/HlGgjkVFDf5ImwtI0W5OJ/93HoClAI6PiH/fQkwG1FamdnPqSVkZ+cEH5FtvmYl2kYLi85EbNpi0kyZZdy+83nCbSSCcc07wOqedFh0/alQw/rzzYotBZia5fHl4vW+91QiKy2W6Vm43OWxY0PgaGnr0UANrPUhGTBwfGia5xf/9E4DpACInVGwC0CXk+AAAW5yu137LwoVm5fnRo4HLLgM++shs2dmihRkqHTrUDAd36WL8MALbe0ZSXR2+jQVghn+HDTO/i4vNWrKRrFsX3CZj6FDrvXUA060ZPDh8mPeee0zdJkwwXaQXXzSTBgcPDnrTut3mXt58M3bZijMkqj6JBABeAL6Q318CGB6RZiTCDbBf11WutkySpKzMerjX4yEXLgymKy8nBwywNopGGmpDR1VCh38jja1WRlnSDPfGGpEByLFjw+9h+nSyRQsyN9cYf7t3J7/5hpw92wwTT50aHCVSkgap1s0B0AOma7MUwAoAt/nPXwngSv9vAfAYgLUwQ8hR9pLIoGKSJDNnmpcw8oV1ucg//SmYbtiw2kXEKjz6qMlrJSSXX05eemlsQXnwwdjltmhBfv21SbdmTWzX+PLyBn2UzZ2UExOngopJkrzxhrGHWL20Y8aYNLt3195SsAppaeTevSb/U0+FxwWMrWS0UXbatOA1a2sFeTzGUDthgrX9xucj338/9n3/8INp0SxapHaUOFExUWrn55+tR2+8XvLNN02a+fMTFxPAtDwCq6BNmBAuJOXl0aM8N99sjquqyD/8oe7yPR7j6GYVl5NDvvBC9P1WVpKXXGLuOTfX3OeAAeQvvzj7nJsBKiZK3fzrX+bFDLQEvF7y5JODL/v69dajI4kKyty55jsw/Bs6bByIi1dIAqFHD2vX+OxsM8wdyUMPRXeLMjLIESMcfMDNAxUTJT6+/pocN44891zT9Qm85AGOP976Za7LIBspKJF+JKGCkqiQBOwnAwaEC4TXG3s5gp49rcvJzDRu/0pMVEwUe9ixgxw+PDj/JjvbGFgvuyy+LtDYseSSJeTpp0fHjRpl4i64IDEhAUxdLr+cfOQRMyHxlFPI//wnth2kQwfrcrKzyS1bGvSRNjWSERPdHlSJzc8/m3VFDjrI+KAsXGgWhA74iTQG2dnAzp2xNw4L5Y9/NMsoRK5ze+CBxn9F/VBiotuDKvYwf75xCvvDH8xm5IGX8eijzXydRMnMNGvNxruxVlpa7Ljq6uAG63UxebLZR9ntNscZGca57dlnVUgcIBUm+impxKOPmpXjS0pMp+Czz4AnnjCLHj32mFkFLVE6dDALOrdpY7xTd++uPb3LFXvzro4dzSTDeOjYEVi5Enj6abNaXK9ewLXXmm1PFdtRMVGC7NpllkMMnUVcXGxeyCOPNK2CWFtjpKdbx2VmGvf3Nm1M/lat6haT2paIHDLEzD7u2rXO2wFgrveXv5igOIp2c5Qg8+dbrwNbWmrmysQSEpfLdB+slkjMyjIbaAFmlfut9VxZ4o03gEMPBd57r37lKLajYqIEadHCdG0SpVUr4P33Y6+pesYZwI4dZve/WDv/xUt5uWktXXRR9H46SqOiYqIEGTTICEOixsniYiNE06aZFkpurgktWxqR8XrNbGGrVk99mD/f3vKUeqFiogRxucw2ngccYEZecnPNUOzQocERkVj5CgqA8883K8a/+irw1lvm95AhJk23brFbJdnZ1iM9GRmxja2k/eKk1AsVEyWc3r3NBuUffGBEYetWYNYsYMSI2N0Yl8uIBWD2JB450uz+F/qyt20LjBkTLUpuN/DllyYcfLBp2Xi9ZoPyuXPNMG7ovjmh+QYNsuWWFXvQ0RwlGpcLOO648HNvvWWc1o4/Pny0Jy3N7EE8dGjd5U6dakTliSdM16h3bzPcHFibdfVqs+1nRYUxsrpcpgUydqxxPgOMoLlcwIwZtfujKA2OesAqibFgAXDxxab1QppV2l55BejUKf4ySDMyFGuDLCu++85sPN66tTHoxtqnWLGFZDxgtWWiJMagQebF/vFHM+wbrwNZKCKJCQlgWjG9eyd+LaXBUDFRkqNDh8augZJiqAFWURRbcFRMRKSLiMwVkZUiskJErrdIM1REdonIEn+4w8k6KYriDE53cyoB3EhykYj4ABSKyByS30ak+x/J0xyui6IoDuJoy4TkVpKL/L/3AFgJoLOT11QUpXFoMJuJiHQHcBSABRbRg0VkqYh8ICJ9Y+QfJyIFIlLw888/O1lVRVGSoEHERERyALwNYDzJyPnniwB0I9kPwCMA3rEqg+RUkgNJDmzXrp2j9VUUJXEcFxMRyYARkldI/icynuRuknv9v2cCyBCRtk7XS1EUe3F6NEcATAOwkuS/YqTp4E8HEcn31+lXJ+ulKIr9OD2acxyAsQCWicgS/7lbAXQFAJJPAjgHwFUiUgmgBMAYNkUff0XZz3FUTEh+DrOXcG1pHgXwqJP1UBTFedQDVlEUW1AxURTFFlRMFEWxBRUTRVFsQcVEURRbUDFRFMUWVEwURbEFFRNFUWxBxURRFFtQMVEUxRZUTBRFsQUVE0VRbEHFRFEUW1AxURTFFlRMFEWxBRUTRVFsQcVEURRbUDFRFMUWGmJ1+uEiskpEikTkZot4EZGH/fHfiEh/p+ukKIr9OL06fRqAxwCcCqAPgPNFpE9EslMB9PSHcQCecLJOiqI4g9Mtk3wARSTXkSwH8DqAMyPSnAngRRq+AtBSRDo6XC9FUWzGaTHpDGBjyPEmRO81HE8a3R5UUVIcp8XEapuLyD1x4kmj24MqSorjtJhsAtAl5PgAAFuSSKMoSorjtJgsBNBTRA4UkUwAYwDMiEgzA8Af/KM6xwDYRXKrw/VSFMVmnN7Rr1JErgXwIYA0AM+SXCEiV/rjnwQwE8AIAEUAigFc6mSdFEVxBqf3GgbJmTCCEXruyZDfBHCN0/VQFMVZ1ANWURRbUDFRFMUWVEwURbEFFRNFUWxBxURRFFtQMVEUxRZUTBRFsQUVE0VRbEHFRFEUW1AxURTFFlRMFEWxBRUTRVFsQcVEURRbUDFRFMUWVEwURbEFFRNFUWxBxURRFFtQMVEUxRYcW7ZRRKYAOB1AOYC1AC4ludMi3XoAewBUAagkOdCpOilKspSWArt2Ae3aAS79L9gSJx/LHACHkTwCwGoAt9SS9kSSR6qQKKlGeTlw9dVAq1ZA9+5A27bAoEHAgQcCw4YB//tfY9cwdXBMTEjOJlnpP/wKZj8cRWlSXHUV8PzzpmVSWgrs2AF8/TWwfj0wZw4wfDgwfXpj1zI1aKgG22UAPogRRwCzRaRQRMbFKkC3B1Uamt27gVdfBUpKYqcpLgauuw5g1B6U+x/1spmIyEcAOlhE3UbyXX+a2wBUAnglRjHHkdwiInkA5ojIdyQ/i0xEciqAqQAwcOBA/dMpjrNtW3z2kW3bjD2lZUvHq5TS1EtMSJ5cW7yIXAzgNAC/9e+PY1XGFv/3TyIyHUA+gCgxaQoQxEqsRDGK0Q/9kIGMxq6SUg+6dgWqq+tOl54OeL3O1yfVcaybIyLDAfwVwBkki2Ok8YqIL/AbwDAAy52qk5OswRocikORj3ychJPQHu0xI2onVKUpkZUFXHFF7WncbmDcOCBD/99w1GbyKAAfTNdliYg8CQAi0klEAjv8tQfwuYgsBfA1gPdJznKwTo5QhSqciBOxGquxD/uwB3uwAztwPs7HGqyJSr8N2/A23sZczEUVqizL/Bgf4zSchnzkYzImYyd2OnwXihUPPQR06mQdl5UFXHghMGVKw9YpZSHZ5MKAAQOYSszmbProIyI+GczgBE4ISzuRE5nNbOYylz762Imd+C2/DUvzAB+gh56acrKZze7szh3cUa96VrGK1ayuVxn7I5s3k7/5DZmVRXo8ZKdO5LRp5IYN5JVXkj6fOX/eeSZtKF98QV5+OXnBBeS775JVVY1zD4kCoIAJvpeNLgzJhFQTk5f5MnOYEyUmIDiGY2rSzebsMJEAQaGwG7vVvOS7uZtuuqPKyWY27+bdUddezdU8nafTSy87sAPv5t2sYEVYmsVczGN5LF100U03r+E1LGaxsw+liVNVRe7aFf7y//gjWVREVlebMGCAERgzlkOmpZGdO5P79pl8N99Mut2kiIn3esmzzzZ5Ux0Vk0ZiPdczm9lRAuClly/whZp0ozjKUnBymMOv+TVJch7nsQVbWKY7hseEXXczN7MlW1IoNWncdPNCXliT5gf+ENVqymY2T+WpDfNwUpD33yf79ydbtSKPPZb89NPw+MceI9u0ITMyTJoHHogWgLlzyZycoJAEgtdLjhtHtmsXHReI//DDhrrT5FExaUSu5bX00hv2Uh/BI1jK0po0QznUUiRymcuP+TFJcgVXWLZMQHAUR4Vd86/8KzOZGZUum9ncwA0kyQmcYJnGTTdXcVXDPaAU4fXXTZck9AV3u8lPPjHxzzwTHe/xGIEJZcIEa7EItFBixQHkVVc1/H0nSjJiorMMbOJhPIzn8ByGYigGYAAmYRK+xJfIQlZNmnNxLjzwROUtRSkOx+EAgD7og/QYI/aH4bCw4wVYgHKUR6XLQhZWYAUAYAmWWKbJQAZWY3X8N9gMIIEJE4yjWSglJcBNN5nfkyZFxxcXA5MnB4937QIee8z6Gi4XUGVtUw/LDxhv2vfeAz79tO48TYJE1ScVQiq2TOKhhCXsz/5RdpOAQfZLfslf+atlSwIEO7MzD+WhPIAH8Cpexct4GdOZbtnqWM3VJMmbeTOzmGXZeiliUSM/kYalpIR0uaxbC9nZJk2sVoVIsKszdWp06yU0XW2tksC1rr/efOfmGgNux47ksmWN9miigHZzUp9SlvIG3sA0pkW94G3Yhqu5OkpsQo21gd8ZzGAe86LSZjGLwzis5nqbuIm5zI2yq0R2mZore/aQZWXmd3U12bKl9QvudhtbRq9e1vEZGeTLL5ty7r47tmh07x6foFiJWseOZGVl4z2rUJIRE+3mOEwpSrEP+2qOs5CF9Vhv6V9SjnJswAZ4Ye1OSQSdiCtQgR3YAQ88cMMNgSATmbgAF+A/+E9Nus7ojPmYj5NwEjKQgZZoiWtxLV7H6zbeZerx1VfAEUeY2b4+H3DRRcCePaY743ZHpy8pAUaOBHJyrOMrKoCxY41z2sMPW8/FSUsDbrvNOn8kVp61e/cCn39ed96UJVH1SYXQFFomP/JHjuAIpjOdLrrYkR35Z/6Zi7mYR/AIy5ZHLnM5m7P5Nt+mh56a1oSLLsv0oR8PPTydp6sfCcl166JHWrKyyOOPN0O2550Xu8XgdpM33mi6H3W1LqxC//6mBXPEEaYO8bRSAiE31/iipALQbk5qUMUqHsgDLV96oVh2cQJdl4EcyL7sG2YLiZXeSlAWcEFj336jc8MNplsS+bJ6POSKFWbotraXumdP0+VIRkwCdpPcXCMmbdoY20hofGZmuH9KqC1l+/bGfnqGZMREuzlJUIlKzMAMTMREPI/nw7oxgHGF/wE/WOYlGNOFvhrVKEABVmAFKlFZcz5W+kgqUIH/QVfrWbHCdEsiycgA1q4FOnasfS7N2rVAWVny1yfN8gV79wLbt5uJgDk5phvUtSvw4ovAkUeGTw7Mzgbuu890y5oqKiYJshu70R/9cSEuxCRMwrW4Ft3RPWwOznzMRzXimG7qRyAA4heNWFSgAi79k2LwYDNvJpKyMuCww4BLL61dTKqrjX3FDkgjKsXFwJVXmkWVRo8GPvvMrOCWmWnEhgSeeQYoKrLnuo1Cok2ZVAiN2c0Zz/FRQ61CCfNOfZpPx9UtceLTlm25kzsb7fmkAtu2ka1bh4+YuN3k6NHBNB98kLxdJNngcpEFBeb6P/xgvGEju0edOpEVFdb31ZBAuznO8zpeRxnC28AEUYhC7ILxRroAF9S0NhqaPdiDh/Fwo1w7VcjLAxYuBM46y4zkdOpkRlleftnEFxcDPXqYVoLPF3+5v/ud6SJZcdJJQLduteevrgamTTO/n346uitGmhbRRx/FX6dUQsUkQeIRiV3YhVaIr/PbBV1wNa6GD9H/qtujPa7BNWiDNnHXrwxlmIRJuByX17vb1JTp0QN46y1ju9i82YhJWhpw++1mhfkBA4AOHUwXJB4eeQSYNQtYtAjo3Ts8btIk4OOPgcWLTbm1sXu3+d640SxWHUl1NbB1a3x1SjkSbcqkQmjMbs4ETojq5rjo4nE8jiT5D/4jpgdr5OdgHsxN1ZtIkl/xK+YytyauMztzDdeQJJdxGduxXULdHQ89nMIpjfacnGb7djOfZuXK+PM89FBsz9XawiOPhJezdSvZu7eJmzQpul4DBsQuKzDJ7+WXo7s5ge7YqhSYMgUdGraPalbzDb7BIRzCfuzHu3k3d3M393AP+7M/c5jDNKbRRx/bsz3Xci0LWRjTe7U2IamsMm6PAUEJFZJAXDKC0oVdHH9ODU11NTlxohlGbdHCvHyDBpE//VR33g4dkheSkhLyrLPI2bPN8datwcl/FRVmvZK33zbHsQRFhNyyhVy92qTPzAy363g85MUX2/3EkkPFxEZu4A1hs4Czmc1DeSiLWcy93Mv/4/8xn/m8iBdxE40oXMfr4nIwCxWSed/P46GPHsp129eRJBdwQY2QFGwuYK9HenH5tuUkExeUFmzh+HNqaN56K/p/9IwM8oQT6s4ba15OrHDbbSZfSQn5u9+xpuUQEBTSCEnACS49PVxQOnUKL8/tJhcvNnNxQucAiZAHH0w+91zqLJ6kYmITG7kx5vok/+Q/2Z3da4TGTTdzmculXMoreEWdL3gGM/h99fckjZB47vEQE8GuD3StERTSCEnL+1oSE8G8KXk1gvI5P49bTNx083E+3qy8YgcPtn7xs7OjVzmLxMqRrbbw3HMm38aNZLdu4aIwcyZZXBztTTshZGG9fv2C57OyyOuuI4cPt75Wx46pIyRkcmLi5ILSE0Vks3/91yUiMiJGuuEiskpEikTkZqfqkwjzMd9yZfl92IdH8Ag2Y3ONo1oJSrAbuzEWY3EOzok5ryaUNKQBACqqK4yiA9iwawOGvjAUS39cik++/wQnv3QydpbuBABUVVehisaYGro8QVu0xd3+jweemnIDlKAEEzABd+COhJ9BqvLrr9bn09PNlP5Y/PSTcQxLhOuvN3N8DjgAmDs3OFpTUgKMGGF293vjjWD68eOD68GOGwcsXRqMa9XKjODMirHC8fbtQJPfDipR9Yk3AJgIYEIdadJg9iHuASATwFIAfeoq2+mWyTzOs1zTNZ3pli2WwKeQhRzDMZZT/kM/vdiLW6u3kiTnrJ1D991uYiIsQ5v723Dpj0tJkt/yW3Zgh5pyWrFVTavjO37HHuwRs4Wyl3st77WMZXyKT/F4Hs9hHMa3+FZKt2TGj7duYbRsSZaXW+f59VfT5ahr0SKrkJtLzp9vylm3juza1Trd+PHB611xRXR8PF2sQLcqFUAqdXPiFJPBAD4MOb4FwC11le20mFSxit3YLcr+4aGHecyLKRLH8TgO4qBahSSWoGRNzooSktb3t44pJIHPf/nfmnrHEhMffVzO5VH3WclKDuGQMKOxl15eySsdfb71Yds20yUIzHcRMYbL116Lneeuu6LnxwRC27bkKaeYrkttgrLJmLj4/PPR8e3aBbsoEycmLlih3af33nP+GcZDKorJegDfAHgWQCuLNOcAeCbkeCyAR2OUNw5AAYCCrl27OvIAQ1nHdTych9NDD330sRVb8R2+w9t5e1xiEc/nt/xtzfVGvjIySkyunXktSfPSH8JDLMvozM68iTdxJVdyOIdbpslmtqVX7Lt813Ih7NDFlVKRX34xL+3gweSYMeSCOuY2Dhli/fLm5Jih2upq8m9/M6NDVukCw78VFeTIkdZpZs0yaUKHjZMJJ57o6KOLmwYXEwAfwWyaFRnOhNkTJw3GMe4eAM9a5D/XQkweqeu6DelnsoqrWMhCltO0oUtYEvcs3to+ocO/BZsL6L3HGyUmHf/RscYoG+mHEvoRCrOYxTt5p+X6sUM4xPLeruJVluV56OFTfKphHnADMHasdTfD4yGXhzTYduyI7kKFCkldSxeEDhsnKygN8P9kXCQjJvUywJI8meRhFuFdkttIVpGsBvA0zLafkWwC0CXk+AAAW+pTJ7s5BIegP/rXGGSzkY0/4U9RBlqX/xMPndEZ8zAPB+NgFG4pxMkvnYx9Ffui0m3duxVDXxiK73d8j0EYhNmYjVzkRqUjiDKU4V7ci4NwUFR8IQoxF3OjzuchD5nIjDqfhrSEvG5TnfHjo42vGRlAnz5A377Bc199BVQGJ2vjuuuAO+4w5y68MNzYetllwNtvhxtlzzwTmDPHeNbOnWsMtInSvn3ieVIFJ0dzQmcxnAXrbT8XAugpIgeKSCaAMUDq76l5H+7D4TgcOchBFrLghRcuuOKaKZyGNMzBnDAhCYzatHG3wawLZ+HVs1+FO90s1xUY5QkIypt4M2bZFajAd/gu6nwxijEF0dvOXYJLLBevTkc6RmJknffSVOjfH/j7342gpKWZMHAgMHNmeLrSUtM+CHDMMeZ70yZg3rzg+fHjzRybs8+OHuX5j3+Ruw4dgO7dE6/r6NGJ50kZEm3KxBsAvARgGYzNZAaAjv7znQDMDEk3AsBqmFGd2+IpOxUWR6pmNT/hJ3yID7EP+8Td9cliFrdVbyNJTl85nRl3ZUSN2pDhozy+e32cv9EMKSzl0lrLt1pgGgQP5+GW9/Ee32MLtmAuc5nDHHZmZy7iIucfYAOyYIFxdEtPN12JjAxz3K+fOdeuHXnvveR334V3Odq1I7/5xpSxYgWZlxc+arPG9FK5bp3xQ7nkkqAh9i9/Sa6bc/nlqbFJF1LJAOtkSAUxCfAjf6xzKDjy05/9ub3aLKn19rdvs+M/OtYIyUquDBvlyZuSxy82fEGSXF+9nt3ZPaFrBQTm//h/Me+hjGX8H//HhVzIKqaQ55RNHHVU3S+xx0Neemn0+VBB2bIlWOYVV4QPG2/dWreQxDM87PWSM2Y07POxQsWkgVjHdTyP57E1W7MLuzCDGQm/4KGCsq98H8ng8G/osHEgLlkhAc1ExFQenXGSsrLE3ehrExQyfNnHUEEha2+RtGhh5uPUdb3TT2/wxxSFikkDsIVb2IqtwnxQQreRCHwymcmRHMme7Mkc5ljO2QkVlEg/klBBqY+QgGZo+D7e12jPrDGprIztY5KooDzzjFlgKTIuN5d8/HHrlk1kiMdxbtiwuu/LaVRMGoBYW3IGBAQ0jl8H8SBu53ZWsIJt2Cbmi96XfXkX72J7to+KO4gHcTInsyu7Ji0kgc/RPLrRnllj88c/2iMoDRE8nuD+PI1JMmJivQ+lYglBvIE3LLfb9MGHM3AGAOBEnIjzcT488GAxFlumD7ACKzAREy1HgtZiLW7H7bbUvTVa21JOU+TBB81iRJ99ZtZcLSmxXpgoFRgypOmO6KiYJMCzeBYbsdEyrhKVmIiJOBgHh53PRnadK54FhCQb2ahAhe0rpHnhxZ/xZ1vLbEp4vWaC3Zo1ZsHm7t2Bww9Prf1909OBiy8Gpk41+xU3RVRMEuAu3BW2BUUox+LYKCEBgN7ojc7ojCIUgWCt5ZeiNKl6ueBCB3SADz4UoShKjK7AFc3KbyRZevY0obQ09gbjIqbDYScul/UOfgHcbuCJJ4A//MFcv6nSRDWwcdiK2ItzxnImEwhmYAbykAcffPDCi2xk27YlRS/0wpN4Eh/iQ/yAH6KEJAMZuv1FBNnZZnHo9PTo8zfdBLSx2fk3Lc2IRKwWR1qacXJrykICqJgkxKE41PJ8S7TE7bgd7+Jdyy5Kb/TGRmzEa3gND+NhLMEStEO7etcnDWm4E3fiClyBjdiILERvFlOBCizG4npfqzkwb55xeR840LRQevQwm2N5vSYccwwwcSLwt78BHo99162oCJpYrdi719hzmjyJWmxTITTWaM5szracSBfwOs1hDgdzMEtZWmdZL/LFuNeLjfXJYQ7ncR5Jci3XWq61kslM3sgbnX40Kc8TT4QvJp2VZZYymD7dDPkuXBie/pVXzFKKWVmJ7RecTHC7oxetbmygQ8POM5dzOZiDmctcS9d1Dz18iA/FVdb7fJ/5zGcuc+NaOzb0IxR2Y7cwj9VRHBUldi66+CbfdOpxNAmKi61Xgs/KIm+9te7811yT3Kr28YacHLOsQiqhYtIA7OVelrGMhSy0XI0tGZ+Oo3l0TNFw0UUffUzzfzKZyWxmsy/7ci3XhpVTylIewSMsBW4pl8a4evNnwYLYu/cdeWTd+auqzDYZBx5oFoOORyBied3m5ZHduxtx83rJzp3JL75w/hkkSjJiojaTOFmMxeiP/miJlvDBh1twS8whXCvbRW0EdgK0IhOZyEIWlmM5SlCCL/EllmIplmM5eqBHWNq92IvVWB1VRilKcS/uTahOzYm2ba03MgeM4bMuXC6zHMG6dWYTrdGj615PNiPD2F3SQpbl9XjM0O+6dWbHwa++AjZsAI49Nv57SWkSVZ9UCA3dMtnMzVGtkAxmWHrCeunlS3wpofL/wr/UunFXOtN5Ds+ps5wCFsRcQOlQHprs7TcLBg8OzhoO9TYN3bYiXsrKzBycFi2Me3zAtuLzmeB2my05Vq0iL7qIPOQQ8tRTTQukstIsI9mmjanPoEF1rxTXGEC7Oc5wO2+3nBnsppst2II++uj2f8ZybMIzb3/lr+zO7pbG3cAnm9l1lrOd2y2NsC66eAEvSPb2mwXbtpkX1+02XR6Ph3zwQfvK376dfPVVsxbtzp2x0/3pT9H2F6+X/PZb++piByomDnE2z7Z8wX308SW+xLf4Fh/n45aLNsfLbu7mA3zActJgoMUTD1fz6qhRIg899apbc2LNGtNC2LOn4a/9yy+mBRNpR0lLM0tLphLJiInaTOJgMAbDDXfU+UpUoj/64/f4Pa7CVeiLvha548MHH8ZjPM7FuVFLQmYiE6MR34SNh/EwbsWtaIu2SEc68pGPT/BJverWnDj4YGOjyMlJPO+nnwLnnguccALwr3/Vvul5UZFZ5vHrr4P+JevWAVkW5rSqKrPpeZMnUfVJhdDQLZPt3M485oWtpuamm6fxNNuv9TN/Zi/2oo8+ZjGLOczh4TycO7jD9msp8RO56bnbbRaNjmzhVFSQ559vZin7fKYLc+SR5M8/m2A1e9nlIi+8sHHuKxbQbo5zbOAGjuEY5jKX7dmed/AOlrHMkWtVsYof8AM+yAc5h3Oa5epnTYldu6z31XG7yQceCE87ZUq0TSQjgxwxwsRfdll0WR4PuWxZg99WrSQjJmLy2Y+I/BtAL/9hSwA7SR5pkW49gD0AqgBUkhxYV9kDBw5kQUGBbXVVlNr46CPg9783w8KRnHBC+GLTPXoA338fnS4z02z/6fEAd94JPPoosGcP0K+f+X3ccY5VPylEpDCedzEUx2YNk6zp5IvIP4FanCmAE0n+4lRdFKU+tG5tPetXJHprin3RO5bUpC0tBXJzgXvuMYFs+pP7QnHcACsiAuA8AK85fS1FcYKjjgI6d46e9et2A3+OWCbm9NOjZyMDQNeuQLuIuZ3NSUiAhpk1PATANpJrYsQTwGwRKRSRcQ1QH0VJCBGzuFLPnmZ2cW6uEZK//x34zW/C006ebETD7R/8y8w0eZ57rvmJRyT16uaIyEcArBySbyP5rv/3+ai9VXIcyS0ikgdgjoh8RzJqQrZfaMYBQNeuXetTbUVJmO7dgZUrgSVLgO3bgfx8wOeLTtexo0k3bZpZVqBXL+Dqq4MbdTVnHDPAAoCIpAPYDGAAyU1xpJ8IYC/Jf9SWTg2wiuIsyRhgne7mnAzgu1hCIiJeEfEFfgMYButtRBVFSXGcFpMxiOjiiEgnEQns8toewOcishTA1wDeJznL4TopiuIAji4oTfISi3NbYPYXBsl1APo5WQdFURoGnZujKIotqJgoimILKiaKotiCiomiKLagYqIoii2omCiKYgsqJoqi2IKKiaIotqBioiiKLaiYKIpiCyomiqLYgoqJoii2oGKiKIotqJgoimILKiaKotiCiomiKLagYqIoii2omCiKYgsqJoqi2EK9xEREzhWRFSJSLSIDI+JuEZEiEVklIr+Lkb+1iMwRkTX+71b1qY+iKI1HfVsmywGcDSBs0ywR6QOzMn1fAMMBPC4iaRb5bwbwMcmeAD72HyuK0gSpl5iQXElylUXUmQBeJ1lG8nsARQDyY6R7wf/7BQCj6lMfRVEaD6e2uugM4KuQ403+c5G0J7kVAEhu9W8Rakno9qAAykSkOW7W1RbAL41dCYdorvfWXO+rV6IZ6hSTOPcTjspmca5e+5CSnApgqr9OBYluXdgUaK73BTTfe2vO95VonjrFhOTJSdRlE4AuIccHANhikW6biHT0t0o6AvgpiWspipICODU0PAPAGBHJEpEDAfSE2f7TKt3F/t8XA4jV0lEUJcWp79DwWSKyCcBgAO+LyIcAQHIFgDcAfAtgFoBrSFb58zwTMox8H4BTRGQNgFP8x/EwtT71TmGa630Bzffe9L78CFkvU4aiKAoA9YBVFMUmVEwURbGFJiMm9XXdbyqIyEQR2SwiS/xhRGPXqT6IyHD/36VIRJqVh7OIrBeRZf6/U8JDqamCiDwrIj+F+m4lM9WlyYgJ6u+635R4gOSR/jCzsSuTLP6/w2MATgXQB8D5/r9Xc+JE/9+pKfuaPA/z7oSS8FSXJiMmNrjuKw1PPoAikutIlgN4HebvpaQQJD8DsD3idMJTXZqMmNRCZwAbQ45jue43Ja4VkW/8zc+mPJO6Of5tQiGA2SJS6J/u0ZwIm+oCIOZUlwBOzc1JilRx3Xea2u4TwBMAJsPcw2QA/wRwWcPVzlaa3N8mQY4jucU/p2yOiHzn/19+vySlxMRh1/2UId77FJGnAfzX4eo4SZP72yQCyS3+759EZDpMt665iEnCU12aQzcnXtf9JoH/DxfgLBjDc1NlIYCeInKgiGTCGMpnNHKdbEFEvCLiC/wGMAxN+28VScJTXVKqZVIbInIWgEcAtINx3V9C8nckV4hIwHW/EiGu+02Uv4vIkTDdgfUA/tSotakHJCtF5FoAHwJIA/Csf6pFc6A9gOkiApj36FWSsxq3SskhIq8BGAqgrX96zJ0wU1veEJHLAWwAcG6d5ag7vaIodtAcujmKoqQAKiaKotiCiomiKLagYqIoii2omCiKYgsqJoqi2IKKiaIotvD/Cl8EPY2eEpQAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "gmm = MultivariateGaussianMixture(n_components=3)\n", "gmm.fit(x_train)\n", "p = gmm.classify_proba(x_train)\n", "\n", "plt.scatter(x_train[:, 0], x_train[:, 1], c=p)\n", "plt.scatter(gmm.mu[:, 0], gmm.mu[:, 1], s=200, marker='X', lw=2, c=['red', 'green', 'blue'], edgecolor=\"white\")\n", "plt.xlim(-10, 10)\n", "plt.ylim(-10, 10)\n", "plt.gca().set_aspect(\"equal\")\n", "plt.show()" ] }, { "cell_type": "markdown", "metadata": { "collapsed": true }, "source": [ "### 9.3.3 Mixtures of Bernoulli distributions" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ "/var/folders/9s/lky4p_js2czgsr4_5962ffbw0000gn/T/ipykernel_10929/1003235212.py:6: DeprecationWarning: `np.float` is a deprecated alias for the builtin `float`. To silence this warning, use `float` by itself. Doing this will not modify any behavior and is safe. If you specifically wanted the numpy scalar type, use `np.float64` here.\n", "Deprecated in NumPy 1.20; for more details and guidance: https://numpy.org/devdocs/release/1.20.0-notes.html#deprecations\n", " x_train = (x_train > 127).astype(np.float)\n" ] } ], "source": [ "x, y = fetch_openml(\"mnist_784\", return_X_y=True, as_frame=False)\n", "x_train = []\n", "for i in [0, 1, 2, 3, 4]:\n", " x_train.append(x[np.random.choice(np.where(y == str(i))[0], 200)])\n", "x_train = np.concatenate(x_train, axis=0)\n", "x_train = (x_train > 127).astype(np.float)" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAABGoAAADOCAYAAABmdeTtAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/MnkTPAAAACXBIWXMAAAsTAAALEwEAmpwYAAAekUlEQVR4nO3de29U1/0u8E0INxuMMeYO5u4UEgKlVaOoL6Bvt6+gqqKqqhRVbSCQAOUOxhjbGAzmYlLOHz8dHemc+T4cNthZwOfz5zzaw56ZvfZasxj5WfP69esOAAAAgF/fZ7/2CQAAAADwP2zUAAAAADTCRg0AAABAI2zUAAAAADTCRg0AAABAIz5P4Zo1a1RC8Ul7/fr1ml/7HAYxNvnUGZvQJmMT2mRsQpuqsekXNQAAAACNsFEDAAAA0AgbNQAAAACNsFEDAAAA0AgbNQAAAACNsFEDAAAA0IhYzw3A+7NmTd2M+fq1dkoAAMAvagAAAACaYaMGAAAAoBE2agAAAAAaYaMGAAAAoBE2agAAAAAaYaMGAAAAoBHqud9SqtdNVO8CiepuAPh0+E5BJV0ba9euLbPPPqt/g/Hf//73nc7p/5auw19++aXM+l73fc/lQ+YXNQAAAACNsFEDAAAA0AgbNQAAAACNsFEDAAAA0AgbNQAAAACNsFEDAAAA0IiPup471X+tX7++zDZu3Fhmqf4rZZ9/Xr/Vz58/H/h4qjZLWapmW15eLjP1wHyK0njpc0yqPzSO+BSl+bavapz1rR9937WlXdfv3tJ1K3Mu8LZWszr6Q6nrVafN20rXTJoj0vfGDRs2lFmq7q7m4i1btrz1MV2Xz7/6btt1Xbe4uFhmz549K7MXL16UWRpjL1++LLPW+UUNAAAAQCNs1AAAAAA0wkYNAAAAQCNs1AAAAAA0wkYNAAAAQCNs1AAAAAA04oOo505VY6mOa+vWrWW2bdu2MhsaGiqzHTt2lNnu3bvLbHR0tMyq1/fkyZPymCtXrpTZ3Nxcmd2+fbvMUiWaakFat27dujJL95BNmzaV2ebNmwc+nuoWR0ZGyixVGabs8ePHZZZqB1OV4cLCQpm9evWqzNwLPnx9a6OTNG+mMZaOS/N0Nb+nmtFUF5qksZLm6VRBmubpNKZVd7NaUj3wL7/88l6fL82p6ZpP97K+3x3Seaa5MZ1nlfV5H2lP3wruVLO9cePGMkvzZvq+uXfv3oGPHz16tDzmyJEjZZau+TTvp/F369atMrt8+XKZXbhwocwePnxYZtX83spa1y9qAAAAABphowYAAACgETZqAAAAABphowYAAACgETZqAAAAABrRVOtT9Zexh4eHy2PSX8U+cOBAmR0+fLjMDh48WGanTp0qs/RXs7dv315mVUtF+ivVMzMzZfb999+X2V//+tcy+/nnn8tsamqqzJaXl8sM3qf0l/WT9evXl1m6h1TNMvv27SuPmZiYKLN0/ukvzM/OzpbZ9evXex3Hx61vs1MaD33bm1IDY5qnf/Ob35RZ1USR/q10H0iNUKm96caNG2V28eLFMktz+LVr18rs6dOnZZbaolh9feerSmqI6XseaUykRpd0n6jWhGmOS+eYxnQ6Ln13SOM9nefS0lKZpQanamz2bW1spZHmU5KutdQw1ndOTQ2iad48ceJEmZ05c2bg46m1eGxsrMzS9+UkNaulRsQ///nPZZbO829/+1uZVXNqK/OpX9QAAAAANMJGDQAAAEAjbNQAAAAANMJGDQAAAEAjbNQAAAAANMJGDQAAAEAjVr2eO1WGVtWDqb4sVZQdP368zL7++usyO3bsWJml6t3x8fEyS6+het379+8vj3ny5EmZ7dy5s8x27NhRZqm6+7vvviuz6enpMkvV3eoFeVvpmknVmOm4VBNY1f3t2rWrPCbVHKbzSFWcz549K7NUr5rqBdO5GJsfvjTXpms+zVV9quy7ruu++OKLMjt9+nSZ/fa3vy2zLVu2DHw8jc3NmzeX2fz8fJml8VfVhHddrhdP1d3p80l14FWt6atXr8pjWH2pyjeNv3SvT+O9Gitdl6+1vveQyuPHj3v9W6Ojo2WWxnTK0mdw//79MkvvZbpPVBXAaS3fp+6blZOu0ZSlmvtUwb1v374yS99v//CHP5RZNSelMbZt27YyS2uCdevWlVn1vfdNx3377bdllu4v6bvv7OzswMdbGWN+UQMAAADQCBs1AAAAAI2wUQMAAADQCBs1AAAAAI2wUQMAAADQCBs1AAAAAI1Y9XruVPtaVXKl+q9UV5jqxlK1YFVx2XW5Li/VVKd6tkqqGU0V3Nu3by+zkydPltnNmzfL7MqVK2W2uLhYZgsLC2UGqyXV06aazuq4VKFaVXp3Xb5fVRWBXZeru1MlYaoLVdn7cUtzTpr/0jU6PDxcZqmWfnJyssz++Mc/9vr3qixd8+l1p/rOJNX1piyta9I95C9/+UuZ/fvf/x74+KNHj8pj+qxNeDcrMf6GhobKLFXhpurd9O+l86zmpFR3myp5U812WientXBay69Zs6bMZmZmej1nVc+9vLxcHmOObku6V6Yxliq401oyjb+JiYlex1VjM9XE37t3r8yeP39eZmltfezYsTIbHx8vs3S/Su9Jes5072mBX9QAAAAANMJGDQAAAEAjbNQAAAAANMJGDQAAAEAjbNQAAAAANMJGDQAAAEAjVr2eO6lq6lL9V6qtvXXrVpml2uhUO7h169YyS9WDqYKvqmdLtYMnTpwos6NHj5ZZkioQ0+tONaOphjK9J/C2UnVi35rLqqYzVf2lKsZUAzg9PV1mi4uLZZZqFVNdKB+3zz6r/x8mXYdp/ktzy9dff11mf/rTn8qs73ip6m4fPnxYHpPGSpqr0jhKFaSp5jdVcKf1ydmzZ8useu0XL14sj0nrFlZfug5TBXe6nlKlbaqwTmOzGn9d13WvX78e+PjOnTvLY549e1ZmqeZ+z549ZTY6OlpmqbI+3XfSZ5Dek6WlpYGPp7VJWtOw+tJ1kcZt+ozTPJ3mlrReTOdSSd/j0nyUasnTvSWNozTe03P2rdmuxlmav6t73ErwixoAAACARtioAQAAAGiEjRoAAACARtioAQAAAGiEjRoAAACARtioAQAAAGjEqtdzpyqyqsIs1UfOzs6WWaq0TRWkqTYs1Z6l80z/XlWr2LcG8PHjx2XWt3bwfdeewWpKVXp96mlTFeqhQ4fKbGZmpsxSdXAa030qQfl4VHNqmqtStmPHjjL76quvyuybb74ps3379pXZvXv3yixV6P70009vfUySjkv1xnv37i2ziYmJMktzY3q/0hrqzp07Ax+/efNmeYx64JVT1bum2tc0NsfHx8ssXYepFnvz5s1llmrpqyr4rqvn23StpSy9J2mMpXVrem2pcrhvnX01ltLrbqUemDdLVdrp+1/6/FNl/dzcXJmla7sa7y9evCiPSevPNDb7fm9M4+/+/ftlltbQ6X71/Pnz/78T+5X4RQ0AAABAI2zUAAAAADTCRg0AAABAI2zUAAAAADTCRg0AAABAI2zUAAAAADRi1eu5U6VcVUWXKrdSVV6qSxseHi6zVHuWasPSc6bjqsrFVAm6devWMkt1i6nKN73P6bhUbaZCkBakmtmUVbWmhw8fLo9J9capJndxcbHMrl+/XmZp3PLxq+qaU41zmltOnTpVZkeOHCmz48ePl1map6empsrshx9+KLOq1jvVdz548KDMUk1qmr8nJyfLbHZ2tsx+//vfl9nQ0FCZpfM8dOjQwMe3b99eHtO3zpz+0vWUsi1btpTZ2NhYmaVa73QvuHHjRpmlOt9qDKbxkGp+R0dHy2zTpk1lltbJad5MdcrpXjYzM/PWx1kjfzjSWnF5ebnM0rW9EtXQqda9qu5OFdxLS0tllr5vbty4scySVC/e577zpudsnV/UAAAAADTCRg0AAABAI2zUAAAAADTCRg0AAABAI2zUAAAAADTCRg0AAABAI5qq565qJ1PVWKrKS5VoqZovVYqlKsOqyrfrum5iYqLMqlrTqmqz63JdYXpP5ufny2x6errMFhYWyixV1kEL+l6jVfVgqiJOYzNV4d66davM0v0q1fXy6UrzX5pb0rX9u9/9rsxSpXSq2U5ZGhPnz58f+PjDhw/LY1LNaKozT+9lGtOperXvZ5Dm92p9ks6flVOtd9etW1cekyq40xjbvXt3mY2MjJTZs2fPyiyt+9KcVK0lU7VuWlunmu302vpWd6d1cqoxTnOxefrDlz7DdK9P4z3dz9O1nb5vpvtEVVOdarbTd/Dx8fEyS2Msjc3Z2dkym5ubK7N0v0rjNs39LWj77AAAAAA+ITZqAAAAABphowYAAACgETZqAAAAABphowYAAACgETZqAAAAABrRVGdjVWWYKr371u6muspUz51qzyYnJ3tlVc1a32rBVDX29OnTMks1jWvXri2zvtJ5qvzmfUr1gik7fPjwwMcPHjxYHvPgwYMyS7WDqYp4cXGxzGCQVJu5a9euMktz1fbt28ss1fXevn27zC5cuFBmFy9eLLOqhjvNcX3nlfXr1/d6zvS67969W2YnTpzodS5VltY05uF3k+aPap2ZPsNUzz08PFxmaUyn9Vsat2mdnOay6jmrauCu63+/SjXFGzZsKLMrV66U2fLycpk9evSozNJ4qa6TVN2c6qDT9yJWRhrrSRpHfeut070gVd3fv39/4ONpjkj3j7GxsTJL32HTvJOqtKempnodl6T7Ugv8ogYAAACgETZqAAAAABphowYAAACgETZqAAAAABphowYAAACgETZqAAAAABqx6vXcq1kpl+q/Uj1iqhs7efJkmR05cqTMdu/e/dZZqmZLWaoIHB0dLbNUwZbq5VJl+dLSUpm1XonGxyNd26nm8PTp0wMfTxWq9+7dK7NU15vquaFSzXM7duwojzl69GiZ7d+/v8zS3JJqa8+fP19mly9fLrPp6ekye/78eZlV+lZRp38r1eum6uNUAby4uNjr36vqiFNNcXpPeLO0pq2yNB+lz/fAgQNltmfPnjJL691Ud5vWfanOt1rTpuuwmmu7ruuOHTtWZocPHy6zubm5Mktr4fR+pdf94sWLMqs+81T5nO5Jfa473k0am6kSO43NM2fOlNnk5GSZpXk6XYfVnPTgwYPymHTNp+9/6f16+vRpmaX5Np1nGi99sxaYoQEAAAAaYaMGAAAAoBE2agAAAAAaYaMGAAAAoBE2agAAAAAaseqtT6sp/eXz9Fesx8fHyyy1vTx+/LjMFhYW3vrfS38VO/017fS601+3Tq1V6a9wp2aL+fn5MkutT9VfwvfX7OkjtTik8V41SqTrMI2HCxculFmfFhuomllSC0W616fjkp9++qnMrl69Wmap7ezVq1dvfR6pqSbp2wqR2ps2b95cZps2bSqz9BrSvadqi0rrFq1PK6dap6XPd2RkpMxSe+jevXvLLLWEpmsjrdEmJibKrBq3qVE13ZNSs1N6T1J7Wppv0+eT1t7puGpNm84jNUJZC6+M9BmmhtDt27eXWRorJ06cKLPUFpXmlpmZmTKrrt++c1VqmHr58mWZpXkz3XfS9+L0HfzZs2dl1vpYMkMDAAAANMJGDQAAAEAjbNQAAAAANMJGDQAAAEAjbNQAAAAANMJGDQAAAEAjPvh67lQtuWHDhjJL9dypnm1paanXcdPT02VWVYoNDQ2Vx+zatavM0utO71d6T/bv319mU1NTZZYq0ebm5sqsqixMlW59K9Zar2b71KRKyr7HpTGR6hH37ds38PFUSfjjjz+WWRor0Ec176Ta3TS3pHksVc//5z//KbP79++XWarg7lOZ3ff5kjRvpurSVHV+6NChMtu5c2eZnT9/vswePnw48PFUAcy7SddG36r4ShrTfeveUwVwqiNeu3ZtmVXXW6q2PnDgQJmNj4+XWbonpfXno0ePeh2X1slpnVGdZ6oiThnvZv369QMfT+Mo3Zf37t1bZql6PlXWV+fYdbmmOtViV9KaIFWPj4yMlFm6t6T7x+LiYpmlCu67d++WWVqzV/fHtF5IY/N9f6f0ixoAAACARtioAQAAAGiEjRoAAACARtioAQAAAGiEjRoAAACARtioAQAAAGjEB1/PnaT6rFS/t7CwUGapditViqV60qoOLtWezczMlFk6bvPmzWWWauJSvWeq7k7vc/p8qvcyVaylc+xby8rq61vPnaoxU519qk6sqn5T/WGqKU7H9X3dfNqqCt1UjZnqblNtZqoEraqh33RckqqPK33v9em4VJN69OjRMjt79myZperu9Lrn5ubKrKounZ+fL4/pU+XK/5E+q+Xl5YGPp2strWPS3JLWdilL4z3Veqe1cFVhndaYqR44/Vu3bt0qs1TXm9bQaZ1ZfaZdl6vCq+es1hhd9/5rfj81aU1V3dPT95mJiYkyO3bsWJmdOnWqzNJcnKrgL126VGZpLby0tDTw8b713OvWrSuzdG9J96R79+6VWVpnpHtx+r5ZZen50thM/1YfflEDAAAA0AgbNQAAAACNsFEDAAAA0AgbNQAAAACNsFEDAAAA0AgbNQAAAACN+CDquVPFWqr/2rJlS5mlus1Uu5WqLFMFaarrqs4z1QCmis6tW7eW2eHDh8ss1cTt2LGjzFJlZPoM0r9XfQbp806Vl6kmnNWXxnSqxEt1oanu7/jx42WWqrurSs379++Xx6TxkO4tfasF+bRVYyLNH6n2M11r6T6axma6b6e5OM23faq70zmm7Isvviizr776qsy+/PLLMtu9e3eZpVrhVEd8/fr1gY8/ePCgPIZ3k+qVq2s0rd/u3LlTZmncVp991+V63bRuSvN0mm937tw58PF0zU9OTpZZqgCenp4usytXrpRZGhOLi4tl9ujRozJL99x0f6mkayt9bin7lKS5ZWRkZODjqWY7ZSdOnHjrf6vr8mec1pKpwjo9Z1XdPTY2Vh6T5u89e/aUWao6T2MsVXCntXd63Un1nvR9PvXcAAAAAB8pGzUAAAAAjbBRAwAAANAIGzUAAAAAjbBRAwAAANAIGzUAAAAAjWiqnruqAkw1nENDQ2WWqqGHh4d7PWeqe0vVfOk5+zxfVSfWdV03OjpaZps2bSqzVP2YPoNU6/b48eMySxWV1bXw/Pnz8phU5cqHI1VLpuswVcinWtA0lubn5wc+fvXq1fKYJ0+elFmqOYQ+quraVGmbarZT/Wwam2m+PXDgQJmleS7d06sszbXbtm0rs3SOhw4dKrNvv/2213FPnz4ts1TBPTU1VWZVtXPfz5R3U9W7pvc8VdP++OOPZdan/rnr8pyU1otpLO3bt2/g40ePHi2PSe9JOseff/65zK5du1ZmN27cKLNUkZ7uV6nO9/Xr1wMfT1W+xuabpWsjzQVVrfTJkyfLY9L1m+716TNOddPpvp1ed5r7t27dOvDxNH+n75Rp3kzfs9O1ndbkfb8fpOfsc59+3xXciV/UAAAAADTCRg0AAABAI2zUAAAAADTCRg0AAABAI2zUAAAAADTCRg0AAABAIz6Ieu5UO5jqstNxqdps8+bNZZYqgJNUM1rV/aVqxF27dpVZqiJOVdpJqtlOFWYLCwtllipiqyxVoqVqttWsUuPNqqrKrsvjtqoW7Lo8JrZv315maZzdvn174OOp2jNJr626/71JOi69z3wcXrx4MfDxdM9+8OBBmaVKzZTt37+/zGZnZ8ssXb+Li4tlVtWJptrgVDN6+PDhMkvVq6medOPGjWVW3Vu6rusuXbpUZtevXy+z6nNNtcGsvrRmSmMzXU9p/ZOkKt+0Fk5zS7VOXlpaKo958uRJmaV72c2bN8ssVR+n9WcaL+mz65Op4F45aW03MjIy8PE0Rxw5cqTM9u7dW2Z912FpvFfz/puyag2a5s302tJ36bm5uTKbmpoqswsXLpRZqlxPFdzpvay+g6f71WqOW7+oAQAAAGiEjRoAAACARtioAQAAAGiEjRoAAACARtioAQAAAGiEjRoAAACARjRVz13VhqUK7qpiretyle/BgwfLLFW67dy5s8xSPWKq565qCVMlWqpCTRXc6b1M5z8zM1NmqYL7+fPnZZbqGKs67b41dyqM25KutVSxl+qBUy19uu6T+fn5gY+n2r5U7ZnquVPNYV99K7+TPuPF+Fs51dxy586d8phr166VWaogTePvm2++KbM0F589e7bMUj13qtusTExMlFk6x1Rdmmo6UwX3uXPnyuz8+fNlduPGjTKr6lBVAP86+tzbqrVP1+X1VLrHpnGb5p103aSa3KruNq0VU11vOo9Us52yJM3T6fNJn0G15jE23016/9K1Xa0J0zGpyr6qpO+6XC+f5pb0/TbdW/pch2k+Ta87jdt79+6V2ffff19mly5dKrP0XTTVkqd1eWUl1uR9+EUNAAAAQCNs1AAAAAA0wkYNAAAAQCNs1AAAAAA0wkYNAAAAQCNs1AAAAAA0YkXquVP1bqq7qirF+tavjY6Ollmq6UzV16lKLdUKJ1Ud+Pj4eHlM39rBhw8fltn169d7ZamCNP17qeawqp5L10KqX1MB3JY0blPtZ6orTBWCqXpwenq6zB49ejTw8dnZ2fKYdP+rxnrX5fOv6k7fJI2XNCbe93gx/lZONRek2sx//etfZbZr164yO336dJmlOfXQoUNlluaIdP1WWarZ7lu7myqx5+fny+zvf/97maXPINVzV/ekrutXQcrKed/3vXSNpnknzR/puLSmTWOzes50Hql6PNXuPn36tNdzJuk9SWMsrWvSZ0d/fddN1bzz8uXL8pjFxcVeWZr/9uzZU2Zzc3O9/r30up88eTLw8TQPp5rtc+fOlVmaN9P3xjSnpvGe7knpO3P1nH3ucV33/se6X9QAAAAANMJGDQAAAEAjbNQAAAAANMJGDQAAAEAjbNQAAAAANMJGDQAAAEAjVqSeO1VaVRXcSarUTM+XataqirKuy1W+Y2NjZbZjx44yW79+fZlVdcTptaVqtrt375bZzZs3y+zixYtlliqMU3V3Os9Ul1Z9Pn0qvfl1pOs31VimrG8Fad+a+KqGO90/UiVoujeme0R6bSlL/x4fh6ouNl3zV69eLbPvvvuuzNLYTLW1qbp7cnKyzPrU5Kbxl+pCL126VGaXL18us3/+859ldufOnV7nkj67tK7hw7cS65i+c2qqPk7HVXW3w8PD5TGPHz8us4WFhTJL4z3N7X3nzfT5pHsgKyN9jqnW/f79+wMf/8c//lEe03ccpXrrJF1r6dpO46Wq/E7z39TUVJn98MMPvc4jzWPpfU7jPVlaWiqz6hpK7/9qrq39ogYAAACgETZqAAAAABphowYAAACgETZqAAAAABphowYAAACgETZqAAAAABqxIvXcqZa3T/VgqvGan58vs1Sz/eDBgzK7cuVKmaW66WTLli1l9ujRo4GPz8zMlMdcu3atzFIlaHpt6bhUF5oq8FLWpx5R3fCHI90Hkr71l6kK8PPP61tdqvu7devWW/9b1Xjuunz+qT4wjSO19J+2an5Mc1VVTdp1+dpO88eZM2fK7MsvvyyzQ4cOlVm631dzf1oTTE9Pl1mq4L5w4UKZPXv2rMzSvJnuOyq4GSTNqakeuO9cnKS5rKrnTtd1ytK9LN0j0nFp/KX5NknrDGN6ZaT1T3UdJqmmOn0nO3fuXJmNjY2V2datW8ts3bp1vc4lrSWrqu30/S/Vi6dxlO4Radz2rcVOtd6pzrzPeawmv6gBAAAAaISNGgAAAIBG2KgBAAAAaISNGgAAAIBG2KgBAAAAaMSa9FeN16xZ0+tPHvf9C/PVcekc019Z37RpU5nt3r27zEZGRsostTcNDQ2V2fDwcJlVf2E+NUbMzc2V2ZMnT8osNXr0bdtJ7U19/tJ217XT7vT69ev3X5fwHvQdm61I94j0l+5Tk1saf+n6Xb9+fZlV12Hf9or0F/KTdP6fKmPzf6S2lz7HrESW5uI0h6dsw4YNZVZJLTCp6WUlxvTHzNhcGSsxVtLaNM3FKRsdHR34eFrXpTVBGmNp3ZoaGFMjUHrOtKZN51kdt9rNMp/S2OzTkpbGUZrj0joyrVv7Nrmlay1d29VaMl3XK9Hgyv+rGpt+UQMAAADQCBs1AAAAAI2wUQMAAADQCBs1AAAAAI2wUQMAAADQCBs1AAAAAI1YkXruVvStFEu1g+k5+9ZbV8e94bMps776PmcrVdor4VOqMvzQ9R2bfa77dMzHPB5aYmxCm4zNLM0fK1GFm55zaGiozNauXdvrOassrYNTLfKLFy/KbHl5+a3Po+vyPJ2qij/0ymFjE9qknhsAAACgcTZqAAAAABphowYAAACgETZqAAAAABphowYAAACgETZqAAAAABrxUddzw7tSZQhtMjahTcZmf33ruT8EH3q19cfA2IQ2qecGAAAAaJyNGgAAAIBG2KgBAAAAaISNGgAAAIBG2KgBAAAAaISNGgAAAIBGfP5rnwAAAHzqVFgD8L/5RQ0AAABAI2zUAAAAADTCRg0AAABAI2zUAAAAADTCRg0AAABAI2zUAAAAADRijSpAAAAAgDb4RQ0AAABAI2zUAAAAADTCRg0AAABAI2zUAAAAADTCRg0AAABAI2zUAAAAADTifwHH19WzmkBLvQAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "bmm = BernoulliMixture(n_components=5)\n", "bmm.fit(x_train)\n", "\n", "plt.figure(figsize=(20, 5))\n", "for i, mean in enumerate(bmm.mu):\n", " plt.subplot(1, 5, i + 1)\n", " plt.imshow(mean.reshape(28, 28), cmap=\"gray\")\n", " plt.axis('off')\n", "plt.show()" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "anaconda-cloud": {}, "kernelspec": { "display_name": "Python 3 (ipykernel)", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.9.7" } }, "nbformat": 4, "nbformat_minor": 1 }