Skip to content

Instantly share code, notes, and snippets.

@flying-sheep
Last active December 4, 2018 15:32
Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save flying-sheep/5172a1f1f9757d4ece239e742ba41b08 to your computer and use it in GitHub Desktop.
Save flying-sheep/5172a1f1f9757d4ece239e742ba41b08 to your computer and use it in GitHub Desktop.
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "code",
"execution_count": 1,
"metadata": {},
"outputs": [],
"source": [
"import scanpy.api as sc\n",
"import matplotlib"
]
},
{
"cell_type": "code",
"execution_count": 39,
"metadata": {},
"outputs": [],
"source": [
"ad = sc.datasets.moignard15()"
]
},
{
"cell_type": "code",
"execution_count": 46,
"metadata": {},
"outputs": [],
"source": [
"sc.pp.pca(ad)\n",
"sc.pp.neighbors(ad, 5, method='gauss', knn=False)"
]
},
{
"cell_type": "code",
"execution_count": 47,
"metadata": {},
"outputs": [],
"source": [
"sc.tl.louvain(ad)\n",
"sc.tl.diffmap(ad)"
]
},
{
"cell_type": "code",
"execution_count": 63,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAaQAAAEFCAYAAACsDJN+AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAAIABJREFUeJzs3Xd0XNW1wOHfnqYZjXq3LMly7wU3MGCaaaa3BzwgtJBACiSUNAihJSThUZJAIIQQIHSS0DFgerXBBmzce5Ely+p1iqac98eMZcmSLduSPCN5f2tpWXNu2zMIbZ1zzz1bjDEopZRSsWaJdQBKKaUUaEJSSikVJzQhKaWUiguakJRSSsUFTUhKKaXigiYkpZRScUETklJKqbigCUkppVRc0ISkDngiYot1DEopTUhqPxKRfBH5r4hUisgGEbkm2j5HRO5ps9/zIvLP6PeXishnInK/iNSLyEoRmbUH1xosIh+LSKOIvCsifxWRp6LbikXEiMh3RWQz8H60/TQRWSYidSLyoYiMbnM+IyLD2rx+XER+G/3+KBHZIiI3ikiViGwUkQvb7HuSiCyPxlIqIjd0+8NUqh/SvwzVfiEiFuA14BXgf4EC4F0RWQVcDnwrIm8AA4BpwMQ2hx8M/AfIAs4CXhSRwcaYmt1c8hngM+BYYDowB3h1p32OBEYDYREZATwLnAF8CFwLvCYiY4wxLXvwFvOi8Q0EDgHmiMhCY8wq4FHgXGPMJyKSDgzeg/MpdcDRHpLaX6YB2caY240xLcaY9cAjwPnGmHLgKuAJ4M/AxcaYxjbHVgB/MsYEjDHPA6uAk3d1IREpil7vN9FrfUrHZARwqzGm2RjjBc4D3jDGvGOMCQB3Ay7g0L14jzcbY/zGmI+AN4Bzo+0BYIyIpBhjao0xX+/FOZU6YGhCUvvLICA/OhxWJyJ1wI1AbnT764AVWBVNIG2VmvarAG8C8ndzrXygxhjjadNW0sl+bdvyo+cFwBgTjm4fuJvrtFVrjGneRYxnAycBm0TkIxGZsYfnVOqAoglJ7S8lwAZjTFqbr2RjzEnR7b8DVgADROR/dzp2oIhIm9dFQNlurrUVyBCRxDZthZ3s1zbJlRFJmgBEr1cIlEabPEDb8+XtdK50EXF3FqMxZoEx5nQgB3gZeGE3sSt1wNKEpPaXL4EGEfmFiLhExCoi40RkmogcAVwGXBz9ul9E2vZMcoBrRMQuIv9D5L7PnF1dyBizCVgI3CoijmiP5NQu4nsBOFlEZomIHbge8AOfR7cvAi6Ixn0ikftPO7ster2ZwCnAv6OvLxSR1OhQYAMQ6iIWpQ5ImpDUfmGMCRFJCpOADUAV8A8ikxj+BfzYGFMaHa57FHisTa/oC2B49JjfAecYY6q7uOSFwAygGvgt8DyRBLOr+FYBFwH3R69zKnBqmwkNP4m21UXP/fJOpygHaon0ip4GrjLGrIxu+w6wUUQaiNwru6iL2JU6IIkW6FPxTEQuBa4wxhzezfM8D6w0xtzSI4G1P/dRwFPGmIKePrdSBxLtIal+KToUOFRELNEhttPp2KtRSsURfQ5J9Vki0rSLTbOBNOBFIBPYAvzAGPPN/opNKbX3dMhOKaVUXNAhO6WUUnGhV4bssrKyTHFxcW+cWiml+q2vvvqqyhiTHes4YqVXElJxcTELFy7sjVMrpVS/JSKbut6r/9IhO6WUUnFBE5JSSqm4oAlJKaVUXNCEpJRSKi5oQlJKKRUXNCEppZSKC5qQlFJKxQVNSEoppeKCJiSllFJxQROSUkqpuKAJSSmlVFzQhKSUUiouaEJSSikVFzQhKaWUiguakJRSSsUFTUhKKaXigiYkpZRScUETklJKqbigCUkppVRc0ISklFIqLmhCUkopFRc0ISmllIoLmpCUUkrFhX6dkJoCfv6zcTEV3sZYh6KUUqoLtlgH0Jsu/eQ5Pij7hmEpBXxx2g2xDkcppdRu9Ose0kT3RtYctYR/jv0IE/bGOhyllFK70a8T0tXDc7EKFDg9mHBDrMNRSim1G/16yC417UJ8ljBW+xAsttxYh6OUUmo3+nVCEksyrrQfxjoMpZRSe6BfD9kppZTqOw7YhGSMiXUISiml2ujXQ3a7csXzi3h0+WIcNuGdC07miKGZsQ5JKaUOeAdkD+m5NcugeAktBd9y48efxzocpZRSHKAJ6eIpRWBgSmozP5n0AY8ufJeKZk+sw1JKqQPaATlk9+DsmVxeVUxOzeWkusNUbp7HghvWUDjiUMZd/V8sNnusQ1RKqQPOAdlDApiaVUhLRWT1BldpNYM91dgWvcaPX7wzxpEppdSB6YBNSABpub9n9ROV1L+2prXtyjdupeLDf8QwKqWUOjAd0Akpa/hMpt34LaE2ba6ByVQ//0NWPXZlzOJSSqkD0QGdkLbznX0vISDlxKHk/mgaA248nPDXjzP/+V/p80pKKbWfaEIC/paej99iwz4wGQCxCI68ZFLm/IG5Px8T4+iUUurAoAkJGJCYwst546h9fjktWxvxb6qjZUMtAAUVK7XAn1JK7QeakIA7Jp/I0Mse5sRxF2BLd5EwKI2UE4YCIMDGq9IIBfyxDVIppfo5TUiAiHDZiOncfuxV1OACIOwNkjA8A8egFJII8+JPB8c4SqWU6t80IbVx8fCDGVz8DE89FyDU4CfrsklkfX8KzvHZjGnaypznb491iEop1W9pQmrjlY1LePHBHxJeE8aSlgBEek+uCbnY0pzkz7mFJn8wxlEqpVT/pAkpqsLbyNcPX8SYhc8zy/sZaz7YRvPCMgLVHhLH5pB73SG4i1P5zX/ejnWoSinVL2lCikq2OzmzfBkAQYS3c0fRuCQFizsVALFZSDliEJfPPQUTDMQyVKWU6pc0IUW5bHayTvo5DakDuH7c6Tw77hyChwzF6tyRfCxJDtyT86jaui6GkSqlVP+kCamNgv+5k0P+UsYjlz9F4Kssln5di7cmke2LNTgKUsg4ZwxLXj8vtoEqpVQ/pAmpE5nuBK4JzmH8srdoeOJ9tmckEwwDUFhTSlifS1JKqR51QNZD6kqy08YVF3yHqgdfIVjpofKhhZBoJ1ztxZJgJbC1iUDFOhIG6rJCSinVU7SHtAs5B5/FgO//C4BAaSP1axsJ1ngJbG3CAL7SZbENUCml+hlNSLuRfthFWFwpANRmFtFsiVSSFaD04YtjGJlSSvU/mpC6MPTO5WT84Bm+HH0c7vCOGXeVxqn3kZRSqgdpQuqCPWMgL6QU4F03H4HIcB1WskN1bLh9RqzDU0qpfkMnNeyBmblDOGf4EaTb7MxOz8e5+BUA/CWLMeEQYrHGOEKllOr7NCHtgYMyB7L60vuwiPDlzZNJjrZXZU8G0U6mUkr1BP1tuoesFgsPr5rHh8HI4qobnemcFLyBt1ZWxDgypZTqH7SHtBfuWfoxG4cfy38HTGSl52BsTSmMyE6KdVhKKdUvaA9pL9ww7kiGpufxk2mn8FrFvdzh/ymvl82PdVhKKdUvaELaC1eOmsHKs37BcdtWkl2zkVmVq7G/dAsmHI51aEop1edpQtoHqTMuoCF5AABHrf+MD954hpXbGmMclVJK9W2akPbB6rmP4WqsoAUbTbZUvvOejzF3fcjKbU2xDk0ppfosTUj7wLb0ZeyEKLNkc+3Yf5NofAwJbmZtlSYkpZTaVzrLbh8UXHgvW+bcx8Cpl/NSfgIVd/0ERBjqmg7kxTo8pZTqkzQh7YPkibMZPXE2AI2L5xB2TSeYfBZVX2+lcEiMg1NKqT5KE1I3JU88CcuQekydULMqSGGsA1JKqT5K7yH1gET3eiRYjq3xjViHopRSfZYmpB6QYFuJc9u12BpfiXUoSinVZ+mQXTctqdnKoPP+D+fwsdhHlOCt+zPO1KsRXXRVKaX2iv7W7IZbvnmbCa/cw+EfPo17aj5GVtLS/CrhwNpYh6aUUn2OJqRuWNtQBcCGxhqszkMQaw5WxwQs9kExjkwppfoeHbLrhvumn8aIlGxm5Q/H7hiMPe/ZWIeklFJ9lvaQuiHHlcwtBx3P4bmDYx2KUkrtNRFp2un1pSLyQPT7W0WkVEQWRb/+0NvxaA+ph32zpZ6P11dz6bRCUl32WIejlOon5PrXLgDuBIqAzcCN5p5Tn+nly95njLm7l6/RSntIPSgUNhz90Of89JVlXPfqsliHo5TqJ6LJ6BFgECDRfx+JtvcbmpB6kEUg0+0AIDfRQTgYinFESql+4k4gcae2xGh7d7jaDMktAm7fafu1bbaf0M1rdUmH7HqQiPDlT2ayYm0Fya8v4NtFKxl11bE4s1NiHZpSqm8r2sv2PeU1xkza/kJELgWmttmuQ3Z9WabbwWgTJNjsJ+Rtobm0JtYhKaX6vs172d4naULqYb7qRja9vBAsQsbEQaSP0+VWlVLddiPg2anNE23vNzQh9TB/VSOmJQhhQ+rIAVhs1liHpJTq46Kz6b4HbAJM9N/v7YdZdvuVGGN6/KRTp041Cxcu7PHz9gXGGLZ9ugoTCJF35GjEqjlfKbVnROQrY8zUrvfsn/S3ZQ+r8wSY8M4mhry9kXfWVsU6HKWU6jM0IfWgxkVv8N7bL1LZ3IIvGOb/PlyHMWFM2Bvr0JRSKu7ptO8e0vj1K5T8+Qwm5LpZdflRrGzIITnzFpoqvks4uAVX+q9wJB4T6zCVUipuaULqKdbIA7GusdmkpHnJTduELfEFgp7IrMyg/2tNSEoptRs6ZNdDkifOpugX7+FZVI1vXQ0tGxMQXNGtFhKSzo1pfEopFe+0h9SDksYcQ0LOFKof/RDnMBdZM76L1V6A1TEKq727D1QrpVT/pj2kHiaWyHNH4eZaxOIiIekMbI5RMY5KKaU6EhEjIve0eX2DiNwa/b5t+YmlInJab8ejPaQelv+9J6if/wzJk09vbXv6qy1c9+oyLppSwD2njY1hdEqpvkoeu6Fj+YnL7u7ug7F+4CwR+b0xprPnVO4zxtwtIqOBT0QkxxgT7uY1d0l7SD3MnjGQrJN+RkLeCGo//idrfz6cda/eS0VTCw/P24Sv8Rk8tf9HOFQf61CVUn1ENBl1LD8Rae+OIPB34Nrd7WSMWRHdN6ub19stTUi9xJggDYsfIFCznjNq/8uMQek8dlYIf8OjBDxv0dL8UqxDVEr1Hb1VfgLgr8CFIpK6qx1E5GAgDFT2wPV2SYfseom39v9IOyud5GMPxZLs4Bn7KtzWYKSDDLyzJsSpU2Ibo1Kqz+it8hMYYxpE5F/ANcDOT/FfKyIXAY3AeaY31pprQ3tIvcRX+jkAliQ7IlC+5Vne/Mdf2VqfxEebCrh6Tl6MI1RK9SG9XX7iT8B3AfdO7fcZYyYZY2YaYz7poWvtkiakXuDbsoyKv79D/VtrqX9jDbUbvYxIbeHkk5KoebmOX39+GdcfNSbWYSql+o5eLT9hjKkBXiCSlGJGE1IvsGcUEG7w0/TxZiq+qOP3a47D6rJjy3Qx6ugj+faGo7h65uBYh6mU6iOis+k6lp/o/iy7tu6hlyctdEXLT/SCxm9eo3nFR1TO/QtWE6DZkcjmi7/HCWMG4874MSJaI0kp1dGBXn5CJzX0sOq5f2bb0z8FdnQ/3S0eRldksnnc2YwRK8YYwoF1WGwFiMUZu2CVUiqOaELqQU1L5rYmIwCxOiDUAsCq5I84qvlDSpqLyUyaSqD5P1hsQ0jOfSRW4SqlVFzRhNSDxLr94xTsOUPIOvmXVG1cxkML1/HrMc1YBJLNRkItSQCEg+sJtazG6hgRu6CVUipO7HZSg4iMF5H5IlIiIn8XkfQ2277s/fD6FveYYxh048cU3/wZQ+9azesFB7HZYuWSxjcIBSL36sKSgiPx+NZjwmFdsUEppaDrWXYPAbcC44HVwKciMjS6zd6LcfVZ7pEzsQyaztUfvMvFnzxL8+ePYZcw1lAQAItpwNfwN5wpV+NKvwm7c1qMI1ZKqfjQVUJKMsa8ZYypM8bcDfwYeEtEDiEy9VB14uqXlvLgx6VYTJiUg3JIOXEoVlckf4sAxoO/6TkstoGxDVQppeJIVwlJ2q5vZIz5ADgbeJLI4n6qE8GwgcYsflFZwkGn5JM8M/JRhbwBtk+zN+FKmquux5iWWIaqlDqAiUioTXmJf4tIYrT9JhFZJiLfRrcfvD/i6WpSwx+B0cD87Q3GmG9FZBZwc28G1pc9cOY4znAsZ/g7rxJqnIYlyYGIYHHaEBHCngCWRDsmIAQbarCn6jJCSqndW36JdCg/MeYJ090HY73GmEkAIvI0cJWIzANOASYbY/wikgU4unmdPdJVD+lFYF0n7T7g6p4Pp39wJ9g4aMPT4Glh273zqX1lJQAiEtlBYO3TGyn/49uUPtTd1eOVUv1dNBl1KD8Rbe8pnwDDgAFAlTHGD2CMqTLGlPXgdXapq4T0F2BmJ+3HAff1fDj9R+qRkSWhwv4Q4Xp/u20Wl53CU/IRlx1HsR1jgrEIUSnVd/Rm+QlExAbMBpYAc4FCEVktIg+KyJE9cY09imN3SweJyHJjTKergIrIMmNMp+VPD/Slg7ZrqV5JTd2NOB11iEUIB8NYbDv+BggHQljsVhzus3Gl/TCGkSql4sGulg5afomEifSMdmbGPGH2eU1SEQkRSUIQ6SFdb4xpkcj6ZjOBo4ErgV8aYx7f1+vsqa7uIXX2AWynC7N2wZLkxeWrZ/vH6PEbktp84hZ7ZE27Fs9bJCSfj8WaEYMolVJ9wGY6n0jW3fITrfeQ2jLGhIAPgQ9FZAlwCfB4N6/Vpa6SSoWITN+5UUSm0cuVA/sDq2MCDveZbAkMp7kFkty7WFTVNONv/M/+DU4p1Zf0avmJtkRkpIgMb9M0icjq4r2uq4T0M+AFEblVRE6Nft1GpG7Gz3o/vL5NxIor7cdkL80j0RJubTfG4NtQS0uNl+0jpmX1lfRyMUalDhihliDhliD+mibK3l1CyZxvaGnwEvK14KtoiHV4ey06m65D+YkemGXXmSTgCRFZLiLfAmOILJDQ67osPyEiucAPgXHRpmXAA8aYil0do/eQ2tvw28MIBZeTedF4LC47xpgdM+6iXp2TRGnqD/nluSfEKEql+qagx4+/pgl/TTMJWUlUfLmOmgXrIr+2hU4f4benJ1J06hS2fbyS9PGF5MyIj/UkD/TyE3tcD0lEsgGMMV0O1WlCas+z+jO2PnEVWAVzzqlkZMzH2iYfmVCYuldWUfWtjykPVyGWPb89d+9H63hzRQV3nTKGgwpSuz5AqT7EhMJ4K+vZ+sEyrO4Eqheuh+CO31mWBBth/77NUnXmpuDb1oBYLUy+49yeCrlbDvSEtNtJDRL5M/4W4EdEhvckOivjfmPM7fshvn4hccRhDP1dZCLLdfOf4ZbM+e22i9VC6mkjCVSsoPqdu8g87ueIxcK2Rj/vrK7kpNE5ZCR2fC7NHwxx/avLAfjD+2t5/uIpvf9mlOohNYs3EWj0kjNjBP56D1s/WRHp2YS7Pna7fU1GtiQnWVOHsvX9ZaSPL9ync6ie19Usu58ChwHTjTEbAERkCPCQiFxrjNFnkfZSiiOLhXUupqZ527WLVci5aizB6ldo/HoEKVPP4ozHFjB/Uy3Hjchi7pUzOpwrwWblfw8ayJwV2zhvUv7+egtKdVvD2nI2PD8PgC1zFvXadaxJCSTmpWGxWrGnOEkelkvG+B2T1XIPG9lr11Z7r6uEdDFwnDGmanuDMWa9iFxE5OEpTUh76Wfjj+TeZWEGBf5Ltn19a/v2e0rWDBcPvPcYRz3xYy51TGY+V2LbzRDeMxdN7vWYleopoUCQZfe/RaCqqUfPK3YrJhjC4rCRMamY1JEDSBqUhc2V0KPXUb2rq4Rkb5uMtjPGVIqIlp/YB257AjdPOhFv3RZamtdjTHQF8CgRYVZiCckNW5nJG/z3h3cxa8LwXZ9QqTgRDoTY8MI8gs1+Bp97CI40d/vtwTCLbun5xxtyjxxNwQkTMWGDWHb36KSKd10lpN0tRa3LVHeDM/ViTLiGgPftDtscdljqHsDgaSdw8rh12B0DAZ2woOJbc0k1dcu2AFC7dAu5h0eGwzzltax+9ENCHv/uDt81h5XkQdmkjMwnUNfMwOMmYLFbCYdChFtC2FyR+6uajPq+rhLSRBHpbNK+AM5eiOeAIeLAlfZjvJULEdmK1b1j0sKQaZn8J30Mw3Kr8NX9iaBzAe5MnUOi4pu7IIOU4XkEmnykjYnU+vLXNLLiLx3/6NpjdgtTbv0fwoEQi+98mbA/QPXiTRTNPoiMSYOwuHbxsLnaY9FlghYCpcaYU0TkFOAOIhPZ7MCfjTEPR/e9CPg5YAWCwALgBmNMXU/EstuEZIzR/9q9SCyJJCXfzro7jyT3R8OxpewY7/7Z0HJCJvIXX7REiVJxzeKwMfyyowAw4TCerbVUL9rQvZMGwngrG0hId2NMZPpdqNHHhhfm4a9tYsDRkeU0y95dQtPmKgpPmYwrp3+OJnx143Mdyk9MufP8nngw9ifACiAleivm70Qmsm0RkQSgGEBETgSuBWYbY0qjiewSIBfokYSk69HFWEL+KMY8sI2gu33SERFsFvigZAau9OtiFJ1Sey/UEmTtk5+y4v63qfh0dbfPt/y+OdSvLCN1ZPuZpLXfRpZxCzT62Pr+MhrXbqPis+5fLx5Fk1GH8hPR9n0mIgXAycA/ok3JRDoq1QDGGL8xZlV0201EekOl0W0hY8w/22zvtq6G7NR+YMJeEqydP6B8VPE5iOyX2lhKdZunrJYVf51L65pYPbQalre8joITJmKxWvBVNxGobyYv2juyuRNIHT2Q5k1VpI/rt88U7a78RHd6SX8iMgSXDGCMqRGRV4FNIvIe8DrwrIl0T8cCX3fjWl3ShBQPJAGrfSjhwJoOs+7K7zoMuXYhKQWjOiw3pFSsNazbRu23m0gfX0Tjhgq2fbxyRzLqhiHfmYnNaaelrhlfRQO5h4/C6rQz+NyOz+OJRRj2nc7KtvUrRXvZ3qXovaIKY8xXInLU9nZjzBUiMh44FriBSP27S3c6djzwJJFEdqMx5vl9jaMtTUhxQMRCUvYDmHADDyy4hksKtwKRRVgzzh3Nkr8eQ3p1I8XXz2Gpazy+YIgjh2bFOGqliEzzbvRRtWB91ztH2dwJBJt3zLhLyE4mdfgAKj6PDLeljh5I+uiBPR5rH9cb5ScOA04TkZOITFJLEZGnjDEXGWOWAEtE5ElgA5GEtAyYDHwQ3T5JRB4AXN2IoR29hxQnRGxYrBnMyh/Tpk1IKExl3I9Hk/ODsaxf8CIz7v+Uox6cx5srtsUwWnWgMaEwZe8soey9pZhwZHKBv6aJcCC0y2OShuS0e21JdDDplrMZe93JuPLTIo1WC2N/MpuBJ0wkcWA6YrOSc8iwXnsffViPl58wxvzKGFNgjCkGzgfeB65q21uifemJ3wN3R+87bddjyQi0hxR3xhRcQ2NtBuWVc8hNbgYiicmRm0RiIISJrrISDGupCtW7Pilfz/rGai4aOoX6ZVvY+sEyALxba3Gkumlp8hD2BXYcYJXIPaOwiQyv/c8hlH+0nJZ6L56tteQcMhxrQuR5+pyDh7P59a/JnFSMWCyIBUb/6ARMOLxXiwsfKKbcef4zX934HPTOLLu2BPi5iDwMeIFmosN1xpg50UW234zOsKsDlgLdmNe/08V7owaPrvbdfR+v+piJSbe1aws1B6hJOIX13vM5ftQ+Dx0r1aVNTTUM+88fCJow904/jauyJrL8L2+2m6TgHpxN84Y2i/+LgDHYU12kjyukcsF6TCDEsEuOIHXEgP3/JvogXe1bxaUjRh7BkvnDyLEvBrcdZ7IDq9tONm9TWDCKbtzLVGqXgh4/K//2LoFmH6PzklhiayAzaKN2aQmu3DS85TseN2nXO4LWyQyBRl+76dctNT27bp3qvzQhxbHxhzxM5uO/xdiq+FFeGd8ZUk+Gw0qSY1SsQ1P9lLe8Dn9VIwCPbRqJcVhJCfvZum4pCDgykgg2egkHQrgLMjBhQ7DJ126SgiM1kZbaZtyFmaSPLyRr2tBYvR3Vx2hCinNvHXsVJ716C+fmbSDblkwoFKSufDWZRfFR4VL1L0nF2WROHRIphBc2iC9IyB/tCZlIb2fgiZOwOqxsnvMNBCMTHEb98FgS8zNBINjsp2ljJSkjBmB16K8Ytef07mGcm1aQxd+SgxSNTsFiFew2oPH3zC+d3+WxSu2tQIM30rOZWITYrNiSnAw6ezojrjia5KE5pIwcQPb0IYRDodZkBFC3rBSxCCKCPclJ+rhCTUZqr+lPTB8w4/DvY8zPWh+MtaU5GRK6iffensUR0y7BnqHPbKg9Z4zBX9lIQmYSYrVQ+vZiGjdUUHTqFDa8MB9f5Y71lINNIfyVjaSPKyR5SG5ru7sgE3FYMS2Rad82t661rLpPE1IfMGDQWBq2tl+lIcEKE7JfZuOdjzL87nUxikz1RSWvfU3l/DWkjsyn+NxDKP9oBQArHnwHW2L7ZaqsLgeJA9PbtYUDQdY98xmmJcSAo8eSMiKPpEHZ+y1+1X/pkF0fYIyfZuvhBAIG0+b5I0duMtZB+jeFgpCvhW2frcJTWrOL7QEqPl+Np6wW77bITDlvRT1NmyuxJTlbp2wHvS2kT9yxHtzw7x5NQnpSu3OZUJiQN3JfKRwMaTLq40TEKiLfiMjr0denRF8vFpHlInJltP1WESkVkUXRrz9E220icqeIrGmz7aZ9iUV/m8W5ULCUpoqrcJsAf/40i+8XrcY1NLN1e/iEYRgTQAv4HthK5iyieuF6LAl2Jt18ZoeHS7e8+Q1VCyLbR/3gWKq/3gjhMOue+KT9icKGkC8IRMqC25wdf66sTgcjvns0nrIasqYM6a23pHZSXzqrQ/mJ1IHv7bfyE1H3GWPu3un43wJ5wHhjjE9EkoHr9yUQ7SHFuXBwKxgPEODy48/Glt++1ku228fcxY/GJjgVN+zuSC0tm9vRfnVewIQNzaW1AFgTbDizUig4cSJB746izykjB+DMSQHAU1bH8CuOYexPZpOQ0b53tF3SoCxyZozAohMX9otoMupQfiLavs/2svxEZ8cnAt8DrjbG+KLHNBpjbt2XePRVJLOqAAAe/UlEQVSnKc7ZnVNxpv4IY3ykJJ3DquZXGED7dew+/fhr6l86nzOvuQt7ZuSB2XCLj2B9OY7s4hhErXqLCYUxYYPF3r52Zv5xE0gZMQBXbmqHVeF9VQ14yyIJKX1sAf6qRqoXbSRz8mBsiQkkD84hdVQ+5R+toPTtxQQbvVgd1l0mIxUT8VB+AuDaaNVYgF8AW4HNxpjGbsTQSntIfUBC0lk4ky9AxEpBy5Vs+KScYHTKbfOics7/6hHGbXyeLQ+eD0RmUW28YwZrbxhM9Vv3xTJ01Q0hX0u7lRGCHj9L732Dxb99kaZNVe32FYuQPDgHW+KOqsNbP1zG2n99jDEGZ24qYrVgT3ez4d/zKP9wOSWvf03B7EmkjooUvsuaNoTMKYMZcMxY3AWZqLjSq+Un2rYbY64AZgFfEik/8c82m+8zxkyKfnVYw05ELoveQyoRkb0uTqU9pD6m7v2HyRtiwWazEGpuIXFCLglD0ql+cjHhYOSPFBNswV+2HADf5kWxDFfto3AwxPL736altpnCUyaTc+gI/NVNtNRGFtxt2lxF0qBICRITClP27hIA8o8dj1gtBBq9lM2NtNlTEgk2+TChMHVLSkjITMZTWoszM7n1+HAojC0xgeKzD47Bu1V7IB7KT3RmLVAkIsnRobrHgMdEZClg3cUxu6QJqY9JP/IKmqt/DYC47IhFsKUkkPuj6YQ9NowJYbEnUHDNy3hWvE/mift0b1HFmAmGCDR4AfBXR9aCcxdmkn/ceAINXrKnRSYTeLfVs/m1r2haXwFAYn466eOLsLkTSB6WS3NJNWmj87EnJVC1YB1Z04aSeVAxeUeOwZWTQsgfYMVf59JS28ywi2eSMlwXQY1TNxK5h9R22K7b5SeAXwFES07cQLT8hDHmw+hubctPdHYOj4g8CjwgIldGJzVYgX0qc60JqY9xjzmGxPBcnv3qQU7Ob99jFmeAkr+cRtFP3iB54mySJ86OUZSqu6xOB8MuOYLmkmpyZgxvbR9w9FhMKIxYLYR8Lax4cC4mWpPI6nLgGhB5ZkgsFkZcfnTrcakj88k/dnzr68QBkXpEgdqm1rXrGjdWakKKU6kD33umvnQW9M4su7Z2WX5iN24C7gCWikhj9LgngLK9vriWn+ibaj3NmOrTsLa5C2iMoe7lVdA0EltqHgO/9zhi26c/VFSc2vLmIrZ9spK8I0eTMXEQy//yFgCODDdjf3JS62SH5i3VgOAuyOjynNs+W4WvsoGBx09odw9K7X8HevkJndTQR6UmeLC0n0yFiJB+5ihIWkXD/Gfxrv8yNsGpXlO3ojTy78oyXHlpFJw0iYyDihn1g+Nbk1HTpkpWPvQOKx+a22HyQ2dyDxvJoDOmaTJSMadDdn2UCXt3ftykVdLhhUhqPglFYzrfQfVZRadOofLLtWQfEhnGyz18FN6Keja+MI+kQVkMOGZcpKx4dOAjHAjGMFql9o4mpD7Kai/C7r6AQHPHIWR7tpu0Y6B+1fVYfbNJPfjcGESoekPK8DxShue1a6v4bBUNa8ppWFNO1vRhpAzLY9jFR0T2H5bX2WmUiks6ZNeHuVIvpax51/eIPMs+ofTB8/BtWUbYhHe5n+rb0sYWYnU5SB2V3zrsljoqv/X5IqX6Ck1IfZiIleUNJ+1yu2fpNiyJ6Zwy778kP3UT75et2Y/R9X1La7eypTnyYOqjX2xm4G3vcMc7q7s4av9LHTGASTefxbCLj0B2vrGoVB+iCamPO3vq1bjSbiDoN+w8Y9J9WBGeM27jvaYGPMEAc8vi75dpvHpp0xLGv3wPY176P0qb63l43ibKGnw8+NnGWIemVL+lCakfcLhns+nzwwhWNLdrd0/Mhfdv5pH6DZw/eBLXjDm8y3NV/Pdmtvz1PAJ1W3sr3D6htLkegKZAC/UBLzfOGsbUwlRuO2Fkr1436G1h9aMfsOqR9wg2+3v1WkrFG53U0E9MvugOfvUvH7/I+ap1cU0RYfAPJ+G85UVmORwMOPy83Z7Dt2UZVa/+FgDHgJHknHV7r8fdE/yhIOsaqxidmtthYdF9ddWoGVhEKHCnMSYtjzFpcMb43n9otHFtOY3rIovn1q/eSuZBxb1+TaXihfaQ+pFfH3ZmhzaLzUr6BePwfvEcv3juJr732b9Z31jd6fGO7ME4i6dgcSaTNO743g63x8ye+whjX7qb6758tcfO+eWSrRQ+VYY8vpqWXfRUGjdUULVwPSa0ZxNGqhaup+zdJYRbdj0VO3lYHslDckgqzu4wm06p/k57SP2Ie+ihbF7jJjXR067dnusmCLSs/JAP6ip4t3Q1G86NFHQMVJdQ+rcLsabmMvDKpxhy20KMMT3W09gfVtZXRv+t6LFzfv35WmaYEPhClKwuZ+hB7de1bKn3sPrRDyBsCHr85B0xerfn85bXsenFyIPKlgQ7eTNHdbqfzeVgxBXH9MybUKqP0R5SPzOg4OYObdY0Jzbg+yULeHXBExxbsbJ1W8OXL+BZ/QmNC/7TurJDX0pGAC/PupQbJ8zib4ee3WPnPOKEcSxNcLI+M5XB4ws6bBerpXVlBKtr98szecvr2PTSgsj+FsEVLYSnlGpPe0j9SNjfTMNnb2Id1b5oqMVmxTEklZb1kRv1v3I5CPs9mHCI5MlnUPf5k6zHyrmLPuCuxCxOKtz9X/vxZnp2EdOzI2VhVtRtIxAOMSGje8/gjEx3MfraE7CnuDrdbk9yMuaaEwnUe0kqzm5t95bXIXZra2kHgIr5a2kuiQyTjv7R8SQO7Hp9OaUORNpD6kcqX/kttR/d2+mSQlnfnYwtz40BfB8/xqqrc1nz03zC/iYKb13ICUOOYZmngZu+fpPD3niAFzcu2e/xd9fimjImvHwPk165j4/K1+3zeZo2VbL0ntdZes/r+Ko6FsL808frOfvxBZSEpF0yaly3jeX3v8XyP72Jd1t9a3vGhELsKS7Sxha0rsatlOpIe0j9iD2zkJbNDTStayRpaHK7bSJC8qwh1D69BDAYfxMG+NPch8it2czDZcu5Y8SxrGtw0BhsoaHFx1nF4zu9Tm9bWVfBzxa+xtF5w7hu3JF7fFx9i5dgdEWKGr9nl/s1rttGxfw1ZE8f2mm5BX9NM4QN4XCIQIMHZ9aOz7K6uYVrX1kGQEaig4fPGMumlxZgwoaUEXlgIgXvSt74GrFYcKS4KDp9KhN+efoevw+lDlTaQ+pHMmb9kMwbv+SaRR1n2wE4R2Qgjsh9j3JHMq9nj+QRr4+Dl73JYbWbOLdsMQdnDyIzIZHLhk/rViwNLT4eX7OADdEZfWsbqvj1V2+yqLq0y2P/uOR9Xi9ZwfULXqPO793jax6RN5SXjrmUZ468kMUrHbh+8Qa3z+34MPCmVxZSt2wLm17qvERKxsQiCk4+iEFnTSN5SG67bWkuO0cPyyTdJpydaqV26RZqFm+idslmECFtfKRqc+PabTSs3krVwvU0b6nZ4/eg1IFMe0j9TPbwSXyv8VmMoePQXciQMCwd34oq8loaOblyFZVDD8U6+hiaNn7NxBkX8Mfjv4fV0v7vlNKHL6Z5xfvkX/EYy1Om8uqybXz/kCK2BasoSkon15XMzq6a91+eXf8NAxNT2Hzur7ns0+f5dNsGXti4mNVn/7LD/tte+CUNC18i83/v5oyicbywcTFH5Q0l1eHcq/d/xqBxAIx64X18wTDPflPKb44f0W6fsD8Q+TcY6vQcYrGQe1j7B2DDwRBBTwuOFBfv/+BQVj/6AY0fLaV6aC7O7JRID2lIDu6BGXhKqrE67YR8AewpibjyUvfqPSh1oNKE1M+EAyuYMWAxQIekJAlWMi+aQMN7G2h8bwMCXDb/MVwHnUnIW8f0t/5I/cyLmfXBE2xsrOGs4vH8afzR1H/+JAD1nz/F6WWG8kY/L5UtZLllEbmuZNad/Uvc9va1dJzWyI9WqaeBiz9+jpEp2Xy6bQMjU7PZmQmHqZ5zFxjDC09fzyG/+oDm7/y+W5/DXaeM4YHPNnD14YM7bHMXZlK3vBR3YeYencuEw6x86B28W+soOn0q2QcPIxRNaiYcZuy17dcTHP/z07oVu1IHKk1I/YzVPgyrYyxltaXkJNZ1uo9Y23edmr55CYBQwza2/nQAN7iz+CxjMA+1eBmdmsvFp95E05K38Jcu5/bm+7nKfBd3Ygh8UOv34AsFOySkBw85i4/K17O+sZotW5byvNPKT6ccx4ixszrGY7HQeNRVlH/xAv8aMAFrVQmTMgfu8j3+eN6LPLnua/508GlcNnx6p/ucNi6P08Z1/mDp4PMOxVNWS2L+nk0wCAdCrZMUPKWR4behFx5O/aoy0kbvOk6l1N7RhNTPiCWRxIybCdT8ApGOCckYgwnvumy9hAIc1LCVgxq28m1KPh+Vr+OnZ90GVhtVL9/GYcDiq3/O4Imz+cfqIUzKzCfT6e5wHqfNzkvHXMLjaxdy2Tt3U736ExIyi3Dcu6nT60675AFum3g601t8XDh08m7f4+NrF9IcbOGZ9d/sMiHtjsVuJWlQ1h7vb02wM+R/D6NpY0XrA7CO1ESypw/b62srpXZNE1I/5Km5nYKkjr/4tz/w6hqXQ+O7G9ptM0AYsG4/hzWBdYnplG9ZwSfbNjBjylk0zHsae/YQCidMw2JP4CdjZ7YeH/LUs/lPp9NSthxH3kgGXvkkE7KLuXf6aWxd+Ta1qz/BnrvrX+AWsXDbQSfs0fu7d/ppPLf+G26eeOwe7d8T0scWkD624wOySqmeo7Ps+iPZ/UQAcduB1irXkTZ2JKOAxc41o2dz+eYFzBI4KGMgH1oclF73FoN+9jaWTiYaeFZ/infVR4QaK/Gu+ZS6Tx5r3ZY87RwcA8fhHn10h+PCYcP9n2zgn19sbm3buYzGzr4/8hDeO/EqHljxGVnP/IbXS5bvdn+lVN+gPaR+xlv/L0L+rzvdtn2NunBdZLHQ7XeSTJvvAezhAA8vfRkrcE7VWp6acjLzX76D8oRk6s/7HadFZ7K15R51JO4Js/Gum4/VmUzKlB1Tz2veupeW0qVUvrSCrFNvarc00XOLSrnm5aUADM5M5ItNtfz6rVVcfXgx953e8TrbNQb8/HvjtwC8sGExpxSO6frDUUrFNU1I/czSkvcZkd5xhh1Ehux862upefLb9u1tvt+enLZ3nR0tTUx74AyO8TUAcMerGbw37XTuO/g0/Bu/wZE3HKsrBYsziUHXz+HbmjKcVjvONrPp0mZehm/j16RMO7vDOnnF6YnYLILdKuSnOPnvkq2Ewob/LN7ampCa/UH8oTAZiTvWjEtxOPnd5Nm8t3UN1+/Fw7NKqfglXQ2P7IupU6eahQs7f+hQ9a5z/v4XHj35lV1uN+EwVY8uomVD5zPwtguxYwgPdiSqFoTjZ3yfx8NQ/MXDSM5IBv1uKYkOG++XreHYt/+OzWJh0enXMiZtxyy3Wr+HhVVbODJvCA5r+7+DNtZ4cFgt5Kc6+WBtFfd8uI7vHlzEmeMHsLXBx8R7PqLBF+S9q2Zw2OC9XwfOX9uMp6wGf00ziXlpWtZBxS0R+coYMzXWccSK9pD6kYB3Hn+b/S6BkAW7tfMaPWKxkHhQXpcJqW0yahELHqudtKAfOwZXyMa3KxZRDDRVbWTgrW/x6TVHUuVvxmAIhEPU+X3tznf83L+zsGoLlwybyuMzz2+3rTgjsfX7o4dlcfSwHTPg1ld7qGxqAWBxWcNeJ6RwIMTKv84l6InWNLIIE35+2i4XTVVKxY4mpH4k4P0Qu6V5t/uYUJimL7pevidEpEfks9h4duBBNFgdXF6ykJfyxrE50ckfRh7C+mrhy5RB1DVuY8bcu7hg6ET+NfN/SbI7ODS3uN35qnyRuNY0VHV6vXuXfsTqhkqm2ifx1vI6bpw1jMkFaRw2OIN7TxtDZXMLl08v3JOPYcd7aAmy8qF3diQjIvWGLA79sVcqHun/mf2II+lMAt73iUzgbm/7hIaWTfU0lnlwsutJDYbIPaRViZncMfw4wiL8delLpIb8rEnKxB4OUZ9k5R9JhzAjcTyJ3jU0hrz8ffV8Hjz0LN5dXcWwO99j9qgc7j8rskDr5MwCNjbVsrCqBH8oSEKbYbsVddu4fsFrAPyzupzA1kEs3FyHxSLcf+Y4rj1y6D59Hi21zfiiD7Smjh5I3lFjSEhzY3Xa9+l8SqnepdO++xGrfTiw6+JvnkXlLJ+zBZcJtUtAOy95J9GvUZ5qHl/8PM8uepb0YGQI7vZV7/D2/H+Q5W/GJhbE3YgnFOmBZNtTsYiFR78oYV21hwc+24gvEFkv7uBovaJhKVnYd1orr8CdyoiUbBwWK4dlFwOwuc7LhhoPX729hG9//zJbP1y215+HKzeV/GPHkTFpEIPOmEZSYSb25L1bG08ptf9oD6kfCQc2AA273F7/wUayKj14xYrL7FhYdOceUts2e/RpJb9YcZoQVgy5LU0M9tSwIMHNN1tXMbHGwxLbCI4vngHAjw4rZnVlE7NH5+CMVlX9+fijObVwDEXuNCzSPiEl250sO/MG/NEliAKhML97dw2vLC3nxHATgUYflfPXMuCosXv9mQw4ZtdTx5VS8UUTUj+ytfJpkjsZroPIlG+rzUIY2iUjgBJHMgNbGrGyIxHtnKB+N+wYEkyIYc1VVCUk8XVqARh44tv/MqahjGVDzue4k64F4IihmXxzfcep2KPTcju0bWezWLFZIsnLbrVw6wkjufWEkdQuLaH8oxVkH9z5sF3t0hJC/gCZkwf3udLrSqn2NCH1I7bAl7CL2yNVqwbxqq+JMc5VvDRgPKvdWZxSvowi60yGWmYSan6PkPGCxUWN9w3yWmoB2OxM4bYRxzM/Y1DkGmHh9Y3juWydnS9GZZHreRyAxooS0lz7fm/Gu62etf/6GHuKi+GXHYU1OvEgfVwh6eM6n8zQtLmK9c98BoDFaiFjUvE+Xz8cDFE691vEIgw8bgJi1dFspfY3TUj9SILN12m7CUPd+7N4znYKC8eugKRaEPg8YxAfrSoEi5NQ8smt+/+jaDg2z/sc0RBiQcYUvkjb8WNykDeJrFDkAdXk5TVckfQbDg98xVvFw6n89n1+MeGYfYr9y09Wk1TbTEttM97yOpKKul781OZ0IDYLJhTGltS9e0N1S0uo+HQVAElFWaSN0XXrlNrfNCH1Ew3bvtdhmA3AG4R3P5nKEmPIF8PszaOYYA3ylLOZfPca3PXPEki7EqzJYMIEJMjX7gCVqYdwpX8CBzVZaHaW8HT6NmhIZ3RtHj4MViM8GnSw2jaY1bk2yFrPbxa8yaGfeKkt8/BIWgb/+tFM0tusrgDQvKUGi8OKKyeVqgXrCDR6CU8czHlfbOXXFgtDstxMzG4/MSPQ5GPzywuwuZ0UnjoZiy0ytOfMSWHsT08iHAjhyu1eEbzEgRlYXQ7EIrgGpHXrXEqpfaMJqR8Ih+oxwfUdK8QCTywbyy9XTOS3QTu/sQdalxQ6J+BC6iZSmTmIVYGnsCUcRWrTPIrq5jM66VAqsiexfflVSxicWwdjrS7kmgQPAK+HbXy+/cfHmwJhC4das0ksbSQRyK2sY/YjXzD30qmkpEYeQq1fVcbaJz5GrBaK/+cQNr20AIBsmxXcTkJeP7m1jWx4fh7DL91xD2rDC/NoXLsNAGd2MrmHj2rdlpCR1COfoTM7hQm/Oh0R0eE6pWJEE1IfZ0yYutKzsFp3bgcMzF05EnzJiDMAJtSatCTan0ojjeEJV7DK4WFIiwN7ykWcXf4ZLw8YxyWFKxjvS+KaiiKuMBZesPpZGRYGWwyZYwvwXTidj9fXcPNbKxnqHsnDZ0xg08tf8dXyct72WzivtJw1f3wF96AsrE47yYNzIrGFwlgSbFgTHYR8AdLy01lywxBK/vkBpryOkLel9X1s+2xVazICSMjqOK3dGENLnQdHaiJiaZ+VTThMsMm/RyszbO95KaViQxNSHxbwL8NT9csOyQigeUMhW945lo+aUxickcilF0xgwz/eBxPC4MMVMmBJBBEywnZmeFIh5SJCQJplEtbwatYFHfw6OIoUfCBwmS0yO2+lO5HvHFxE1acryR6cx7xNtXy+EWYW5/C98w8lYUUZzz31CfboOonNmyKrM3i31UO0Wm3Z298y8qpjsTrsOKLJIvWSI6hbvqXd/ZvtFVoRGPmD40kq6Lh00KaXFlC9cD2ZkwdTfM7BhIMhNjw/j0CjD7EKTRsqGTBrHPmzdAq4UvFME1If5q25E/B0aF++fiwFTZeQctkgbl9Xy4mjsvn9M19yqYmswfAP90auaI6UazAY6ozQJJBmIFlgQDiF99ZP4k6LjfHFaVBTDhJdQRwYZUKsffITCBtch3hJsFrwBcNkuSP3i+ybKnBGk5ElwYYzKxlPaS2Buh2xesvr2PLGNwy/ZMfQnCM1kZwZI9q9l4EnTMTmTiB5aG6nyQh2JK3m7f+WVFO3bAsAEn0Oqrmkei8/XaXU/qYJqY/yeRZjwuUd2r9cVUzGuzOoYA0jxw7ickeQPz2/gCNqamkSeD9k5T/hceRkeBjus/JediI/LgmSHjbtHj5KMjbuDIGUVOLMSyVtXCHOrGRqFm8ma0oxm1/5ikCjj5yB6Sz92ShqvQGmFqYR9PipnL+29TxitzL88qNZfMeLra9NKAxhQ7glRN2KUpo2VZE3cxQ2d0KH9+NITaTw5N2XNC8++2Cqvor0kADcBRmkjBxAoNFH7ozhNG+pIeewEbs9h1Iq9rT8RB9jTJjmqusJtXzbYVupJ53m8t9jeesrLHYb2YcMY9vHKwmzY42orJMnM6jNL+eQP8DSe94g2BSZMm5PS8SRmkig0UtLzY6FWsVqYdJvzsJij/wNE/T4CTT6Wme3VS1YR6DJR8bkwSy75w1MMETykBxyjxiNp6wWITJkN/CECbQ0eKlbtoWMiUWsePAdCBuyDx5G0ekH7Kr7SgFafkJ7SH2N8XSajHxBCyOGPINtuAPvyBysCXaaN0dX1nbYWOF0MSLLTcGU4nbHWRPsjLnmRGqXllDy6lcE6jzkHjoCA5S+uTg6OwLsyU6kzc0qW2ICtsRIj6a5pLp1xpzN5WDUlbPwVTWSPr6QDc/No3ZpCVaXg0k3nxUJJ81NUlEWJhQmISMJf1Ujrjydaq3UgU4TUh/TWHVzu9fGQNgykpyiPyMSWSnBlRPptTjGFzF2QBo2VwLTOhkO286e5CR9bCFl7y4h5G3BYrdR8vrXYAwWu5VwIERLvQd/bRPOzOSOx6e4sLochPwBErJTSByYQeLAyP0eR1piu3/bEquFMT8+gUCzj4T0npm+rZTquzQh9SFryxaRZb5tv1K3JJGR/+Auj3F2Mk26rZA/QM2iTSQNymLcdScT8gVoWL01cp8HSCrOpmFNOVanA5vT0ek5HKmJjLv+ZMKBEI7U9oln4OxJpE8owpndeRwWh40EhyYjpZQmpD5j0eZlDLJcv9PDrxkkD3i0W+fd8uZiqr5ci9VpZ+JNZ2JLTCB9QhEN67ZhdTkoOm0K/pqmyBDdbnpZ24fvdiYiuAsyuxWjUurAoAmpD6hvXMwgy3Xs9MwnKfnPIdK9hzm3F6uzJNjZnu1siQkMvfDw1n22DwEqpVRv0oTUB3iqbyRxp4W0E9J/1+1kBDDw+AmkDMvFlZvWYZUDpZTanzQhxTlf/VMk2iNTso0Bi2MiSZm3YLH2TK9FLELKsLweOZdSSnWHJqQ4Fg7V4296csdrGUZazr0xjEgppXqPJqR4ZjxAZP04i20YqbkPxzYepZTqRZqQ4pjFNoDEjNsJB0txJJ0e63CUUqpXaUKKc3bXobEOQSml9gutRKaUUiouaEJSSikVFzQhKaWUiguakJRSSsUFTUhKKaXigiYkpZRScUETklJKqbigCUkppVRc0ISklFIqLmhCUkopFRc0ISmllIoLmpCUUkrFBU1ISiml4oImJKWUUnFBE5JSSqm4oAlJKaVUXNCEpJRSKi5oQlJKKRUXNCEppZSKC5qQlFJKxQVNSEoppeKCJiSllFJxQROSUkqpuKAJSSmlVFwQY0zPn1SkEtjU4ydWSqn+bZAxJjvWQcRKryQkpZRSam/pkJ1SSqm4oAlJKaVUXNCEpJRSKi5oQlL9hoiERGSRiCwTkcUicp2IWNpsny4iH4vIKhFZKSL/EJFEERklIvNExC8iN8TyPSh1ILPFOgClepDXGDMJQERygGeAVOAWEckF/g2cb4yZJyICnA0kAzXANcAZsQlbKQU6y071IyLSZIxJavN6CLAAyAJuAzDG/GY3x98KNBlj7u7lUJVSndAhO9VvGWPWE/kZ//927RAHYRgMoPD7BRYH4QwEhyDchYQzcA0EagLPTdCgCJdAEiRbEZsAs6BYM95nKlbxT72k7RiYAeduJ5LUxiCp76LrASR9xyCpt5ojuxK4AVdg3u1EktoYJPVSRIyAPVCk+qK0ANYRsXjbs4qISVczSvrkowb1RkSUwAUYAE/gAOxSSlXzfQlsqe+UKuAIbIAhcGrWCngA05TS/df/IP0zgyRJyoJHdpKkLBgkSVIWDJIkKQsGSZKUBYMkScqCQZIkZcEgSZKy8AJmL5YEX2aMVQAAAABJRU5ErkJggg==\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"del ad.uns['exp_groups_colors']\n",
"matplotlib.style.use('seaborn-colorblind')\n",
"sc.pl.diffmap(ad, color='exp_groups')"
]
},
{
"cell_type": "code",
"execution_count": 62,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAaQAAAEFCAYAAACsDJN+AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAAIABJREFUeJzt3Xd4lFXax/HvmZoe0oAkEBJ6rxKaosICAhZWXayr6Bbb6i5rl7WLq66u7oru2l0VX1wRXRCxBwsqAQTpPQlpkEJ6MpMp5/1jhpiQkABJmMnk/lwXVzJPvWeA/HKe5zznKK01QgghhK8ZfF2AEEIIARJIQggh/IQEkhBCCL8ggSSEEMIvSCAJIYTwCxJIQggh/IIEkhBCCL8ggSSEEMIvSCCJTk8pZfJ1DUIICSRxCimlEpRS7ymlCpVSGUqpW7zLP1JKPVVvu3eUUq96v5+nlFqjlHpWKVWmlNqplJp6HOdKUUp9rZSqUEp9rpR6Tin1lnddslJKK6V+o5Q6AHzpXX6+UmqbUqpUKbVaKTWo3vG0UqpvvdevK6Ue8X5/llIqRyl1j1KqSCmVqZS6ot62s5RS27215Cqlbmv1hylEAJLfDMUpoZQyACuA/wGXAT2Az5VSu4Brgc1KqZVAPDAWGFFv93HAUiAWuBBYppRK0VofbuaUbwNrgF8AqcBHwPKjtjkTGAS4lVL9gf8D5gCrgfnACqXUYK117XG8xe7e+hKB8cBHSqn1WutdwCvAXK31N0qpKCDlOI4nRKcjLSRxqowF4rTWD2mta7XW+4GXgEu11geB64H/AP8ArtJaV9TbtwB4Rmvt0Fq/A+wCZh/rREqpJO/57vOe61sahxHAA1rrKq11DXAJsFJr/ZnW2gE8CQQDE0/gPd6rtbZrrb8CVgJzvcsdwGClVITWukRr/eMJHFOITkMCSZwqvYAE7+WwUqVUKXAP0M27/kPACOzyBkh9ubrhKMBZQEIz50oADmutq+sty25iu/rLErzHBUBr7fauT2zmPPWVaK2rjlHjRcAsIEsp9ZVSasJxHlOITkUCSZwq2UCG1rpLvT/hWutZ3vULgR1AvFLqsqP2TVRKqXqvk4C8Zs6VD0QrpULqLevZxHb1Qy4PT2gC4D1fTyDXu6gaqH+87kcdK0opFdpUjVrrdVrrC4CuwAfAf5upXYhOSwJJnCrpQLlS6k6lVLBSyqiUGqqUGquUmgxcA1zl/fOsUqp+y6QrcItSyqyU+hWe+z4fHetEWussYD3wgFLK4m2RnNdCff8FZiulpiqlzMCtgB34zrt+E3C5t+5z8Nx/OtqD3vOdAZwLvOt9fYVSKtJ7KbAccLVQixCdkgSSOCW01i48oTASyACKgJfxdGJ4A/iD1jrXe7nuFeC1eq2itUA/7z4LgYu11sUtnPIKYAJQDDwCvIMnYI5V3y7gSuBZ73nOA86r16Hhj95lpd5jf3DUIQ4CJXhaRYuB67XWO73rfg1kKqXK8dwru7KF2oXolJRM0Cf8mVJqHvBbrfXprTzOO8BOrfX9bVJYw2OfBbylte7R1scWojORFpIISN5LgX2UUgbvJbYLaNyqEUL4EXkOSXRYSqnKY6yaCXQBlgExQA5wg9Z646mqTQhx4uSSnRBCCL8gl+yEEEL4hXa5ZBcbG6uTk5Pb49BCCBGwNmzYUKS1jvN1Hb7SLoGUnJzM+vXr2+PQQggRsJRSWS1vFbjkkp0QQgi/IIEkhBDCL0ggCSGE8AsSSEIIIfyCBJIQQgi/IIEkhBDCL0ggCSGE8AsSSEIIIfyCBJIQQgi/IIEkhBDCL0ggCSGE8AsSSEIIIfyCBJIQQgi/IIEkhBDCL0ggCSGE8AsSSEIIIfyCBJIQQgi/IIEkhBDCL0ggCSGE8AsSSEIIIfyCBJIQQgi/IIEkhBDCL0ggCSGE8AuBHUi1lbD9A6gq9HUlQgghWmDydQHt6oMbqNj1FdaufbBcl+braoQQQjQjoANpg6MPixL/TJwu5mGHG6s5sBuEQgjRkQX0T+iMvvPQykiBoSsVdqevyxFCCNGMgG4hzRzZA4cqoEdUELFhFl+XI4QQohkBHUihVhOXpSb4ugwhhBDHIaAv2QkhhOg4Om0gaa19XYIQQoh6AvqS3bH8dvEP5KxfTrBJMf+6W5jct5uvSxJCiE6vU7aQDm3+lI+7PMf7YYtIW7XE1+UIIYSgkwbS3NT+uLVin2Uw0V26885nX1FVUebrsoQQolPrlJfsfn3RpTChP89+YqPUGUa3nCLCl79FdEJvUidPx2A0+rpEIYTodDplCwmAhNFEGO0ABJuNVFu7kVNcRdnSm31cmBBCdE6dN5CAW6b15bzIPQyx/jz46hdqMhl79/iwKiGE6Jw6dSDFxvfkogsvAtx1y0rcoXz34xY2/fC17woTQohOqFPeQzpa35Qe7MvIYbu9K/sd0YDGlbmfEP0R/cbPRCnl6xKFECLgdeoW0hGjchdjdDsocwd5lygqtJXNudV89d7LPq1NCCE6CwkkgPDuJJd+y6igPMINdroYaog21gBQpLtAVYGPCxRCiMAngQRw9l8YlTqJ2XvupsZtotQdzK7aWM86pVi64mPcLpm+Qggh2pMEEoBSMOpKwmfdT5iuAMCs3BQ4Qyh2BoEphI/ffs7HRQohRGCTQKrHNOpyHrh0AoOq1mDFSXpND76v6UlubRjVQQn89PUHvi5RCCEClvSyq0fv+JDs/blUmSOx6COjNSjynOFEmWzsya+hj91BmNXs0zqFECIQSQvpiKpCNn7zMbtqYhgeY2GY2kNPUykhqpZDrnBWVyVTrMP5xwef+bpSIYQISBJIR1jCyIw6HQClXfSqWE9w+X6qaz0PzboxsK82mgGGStxud3NHEkIIcRIkkI4wBzOgf3/CDHbG57xAr6ofya82gjmobhO720S2I5L8wmIfFiqEEIFJAqmeISNGc86FV+C86n+Mzr+RBNNhDDWHwTu7bJkO4id7PP/30WrfFiqEEAFIAqkJ0aFBnNenC/16JjI5pgzqRg7yXKoLDgrB5XL5rD4hhAhE0suuCeFBZv44aywb16UTZnRwesgBal2KUJMLpzYQabRTVVVFRESEr0sVQoiAIYF0DH16JWFCs27dOroY7WgcKGUCgwKtKSsrk0ASQog2JJfsmtGrVy9MJk9mdzFWY3TbPCuUIj19rQ8rE0KIwCOB1IIZM2YwYcxw+pOBy/Bzj7tqu0PuIwkhRBuSQGpBcHAwiQVfcri42DPmndZot5Ngi5kvv/zS1+UJIUTAkEA6HkkTGFj+NSlVG4iNjUEZPJfxysrK0N4u4UIIIVpHOjUcj/gRBN+6mTHKwOqV73Ikx6sNwb6tSwghAoi0kI6XwQgbXiX0UDoAwfYCrv+qlI+35/m4MCGECAzSQjoR3y9izOEskku/5fHiUVgsZ9O/q3T9FkKItiAtpBMx4WYMMb2xTrqR6WPG89Hgn0jZ+6avqxJCiIAggXQiTrsW/rCe7PAxlOtQ8sJHsWl/kXRsEEKINiCBdBKSkpIwKc8zSPutg1mxdhs7D5b5uCohhOjYJJBOwrptu7G5DDjcGptb8dulPzH4kQ8llIQQohUkkE5C9eECTAZFkU3zfFYwQUZFYohib2GFr0sTQogOS3rZnYQxo0ayefsuuvdP4PzYCNau+RYFTEwM8XVpQgjRYUkgnYT4+Hji4+MByM/PJy83h507d6K1m7lz5/q4OiGE6Jjkkl0rxcfHk5+fT3l5OWlpab4uRwghOiwJpDYwfPhwQkND6du3r69LEUKIDksu2bWBfv36YbVasVgsvi5FCCE6LGkhtdahbYwc3I+UgcPYbhrIG9/n4pYHZYUQ4oRJILVG2qPw74lY3pxNjrMLOw/Z+HJnMQeKa3xdmRBCdDgSSK1Rst/ztTSLkT3CiA41M6BbKAldgprfTwghRCNyD6k1ZvwVovtC7zPpFRvK3+cO8nVFQgjRYUkLqTVC4+CsuyBpgq8rEUKIE6aUqjzq9Tyl1CLv9w8opXKVUpu8fx5r73qkhdTGNmYf5uu9Bcwb35vIYOl1J4RoG+oPiy8HHgWSgAPAPXrRFW+382mf1lo/2c7nqCMtpDbkcrs5+5+f86f3NvDnZT/6uhwhRIDwhtFLQC9Aeb++5F0eMCSQ2pBBKWJCrQDEhZhwOp0+rkgIESAeBY4eLDPEu7w1gutdktsEPHTU+vn11s9o5blaJIHUhpRSpN92Dh9ePZrST17mhhtuIC8vz9dlCSE6vqQTXH68arTWI4/8Ae47av3T9dZ/0spztUgCqY3FhFnpUnuY8vJyqqqq2L9/v69LEkJ0fAdOcHmHJIHUxg4ePMgrr7yCwWBg0qRJjBs3ztclCSE6vnuA6qOWVXuXBwwJpDZ28OBBbDYbbrebkSNHYjabfV2SEKKD8/am+x2QBWjv19+dgl52p5TS7TDu2mmnnabXr1/f5sftCLTWrFy5EofDwfnnn4/RaPR1SUKIDkIptUFrfZqv6/AVaSG1sdLqWq791sblX9n4YvchX5cjhBAdhgRSG8rPz2f5uu0UVtqxOd387fMduLXG7nD7ujQhhPB7MlJDG8nLy+O7777D5bJwxYgBFNfU8scpSSx4fzeHyu38fnIS43t38XWZQgjht6SF1EaUUgAcdIZhMlnoFh5GTomb/DI7bg3b8yp8XKEQQvg3CaQ2Eh8fz+TJk+lpqSDGWM3ASCdBJk+HBoOCmUPjfFyhEEL4N7lk14a6du1KUtcuBBfmEN0lmkljRtI9wkrvuBDiZY4kIYRolrSQ2tiRS3cOhwOr2cAvBsfSO+7oIaiEEML3lFJaKfVUvde3KaUe8H5ff/qJrUqp89u7HgmkNjZ27FiGDRvGxIkT65YtXpdBt7vf49ZlG3xYmRCiQ3sw8nIejMzkwUi392tbjPRtBy5USsUeY/3T3jHufgW8qpRq18yQQGpjwcHBDBgwgPDwcDIyMvj4449Zt3k7BRU2Xvh2Lx9uLuCVb7OpsMlI4EKI4+QJn0bTT7RBKDmBF4H5zW2ktd7h3fZYwdUm5B5SO3G6Nd9vy8JYU8UZMRbSU2KZNagnSzccBCAqxMyFo7v7uEohRAfR3PQTrR0+6Dlgs1LqiWNtoJQaB7iBwlaeq1kSSO3klW+z+b4wmiAVgV0buW+igWos7D3kaRmlZxZKIAkhjld7TT+B1rpcKfUGcAtQc9Tq+UqpK4EK4BLdHmPN1SOX7NpJ9qHDANi0EY3i+617Cdr5HtplJ6+8kjc37PFxhUKIDqS9p594BvgNEHrU8iPzIZ2htf6mjc51TBJI7aCsrIwBOpOBlkIGWwuIUeXkmnrxoZrKzNgCDlaVcNe0Ib4uUwjRcbTr9BNa68PAf/GEks9IILWDkJAQgg1O+lpLSDSV8VlWKU5MVGsLUYl92XzPbG4+a4CvyxRCdBT3lzU5/YR3eVt5inbutNASmX6iHeTl5VFQUMDuvXsxABUOWFGWyJzhPbhyfA8MBuXrEoUQfqizTz8hnRra2J49e/jpp58AT99MgHAzPJiqGJRYicGg0Fpz4LCN7pFWrCZppAohBEggtamDBw/WhRF4R23wtkBXr9vMU5tS6ROxnr494/lkezE9ooJ4ZE5/X5UrhBB+RX49b0MGw88fZ2hoKKNHjSI+KYU12SVsMI8FZWBfhZn9RZ57kzklNjKLjr5PKYQQnVOzgaSUGqaU+kEpla2UelEpFVVvXXr7l9exdO3alTPPPJOzzz6bmefMoHdFOmE1uYxL7IIBT0spzKKY1De6bp8Ku8tX5QohhF9pqYX0L+ABYBiwG/hWKdXHu87cjnV1WHFxcYRGdGH98mfhg+s4mJ2BUSmM3kFXK2s1S9blc+X4eK4/syfDEsN9XLEQQviHlgIpTGv9sda6VGv9JPAH4GOl1HigXZ/Y7chufncdj3yVg0ZxwBnFdlssDozetQqbw81HmwvpFmH1aZ1CCOFPWgokpZSKPPJCa50GXAS8iWdwP9EEp0vzv9qR/LXLQ6zRI8lwei7RGXBxJMcPVzt5fNV+HC63DysVQnRmSilXvekl3lVKhXiXL1BKbVNKbfauH3cq6mmpl93jwCDghyMLtNablVJTgXvbs7CObNHcsUzrFYw+tJ+sKid2bQQUbgyAwoQLJ0YsRo2jthZzsEzeJ4Ro3tKlSy/HM5hqEp4hg+65+OKLW/tgbI13egmUUouB65VS3wPnAqO11nbv1BSWVp7nuLTUQloG7GtiuQ24ue3LCQyhVhNJhkqsys1ZoZkMthzyrvHcR9JAP2MO40x7WZe+1md1CiE6Bm8YNZp+wru8rXwD9AXigSKttR1Aa12ktc5rw/McU0uB9E/gjCaWTwOebvtyAkdycjIAZlyEGhr2pHNhJMvVlVpt4EC1BadbbscJIZrV3PQTraaUMgEzgS3Ap0BPpdRupdTzSqkz2+Icx6OlQDpda73s6IVa68XA5PYpKTD06NGDM6fO4Ova/qyzJXiX/ny/qBYLX1f34ouDoby7Pt83RQohOor2mn4iWCm1CViP5zLgK1rrSmAM8Hs88x+9o5Sa18rzHJeW7iE1N+iaPFTbgsIaRXntzx9hsLuGGsPPo7tr70f49e7DzBwaR5cQ6UkvhGjSAZruSNba6Sfq7iHVp7V2AauB1UqpLcDVwOutPFeLWgqVAqVU6tELlVJjaeeZAwPBwO6h/GJQDMPDDmM5Kozqq3G4+WSbfJxCiGNq1+kn6lNKDVBK9au3aCSe0cXbXUstpNuB/yqlXgc2eJedBlwFXNqOdQUEg0Fx5fhEfgwuZcuPjnprNJGqBps2YsfzLNK+wkq01p7x74QQrWKz2VBKUVZWxjfffENNTQ2zZ8/GarVSUlJCYmKir0s8IRdffPHbS5cuhbbvZdeUMOBZpVQXwAnsxXP5rt21OP2EUqobcCMw1LtoG7BIa11wrH06+/QTR0tLS2PXwWrW2xJwYsTTz65h8NRW5jOwVxT3nd+oQSqEaEZFRQWFhYUUFBSQkJDAp59+yurVq3G73SjlGV0faPB9bGws8+bN48MPP2TcuHHMmDHDl2+hjkw/0QKt9SHgfqVUnPe1XFs6QcOGDcPh+JFeupAacxdWZFrQDVpCmtNiFN3tOWg99oRaSX//YgertufxxJxRjOoZ3fIOQnQgTqeTgwcPsmzZMiIiIkhLS8Ph+PlqQ3BwMDU1Ncfcv/4v3PW/LyoqYsmSJeTk5LBnzx6/CaTOrtlAUp6fjPcDN+G536SUUi7gWa31Q6egvoAQGxvL9OnTAahecSfL1RVHbaHYYu9KjNXB1h27GDpoAEopDpXX8NnOg8wakkB0aONhhuwOF7e+/yMAj322jXeubaqHvhD+6bvvvqOkpIQZM2ZQXFzMqlWr+OKLL3C5jn/A4ebCqDldunTh7LPPZtmyZYwfP/6kjiHaXkstpD8Bk4BUrXUGgFKqN/AvpdR8rbU8i3SCQsK7kJK7nQzr4AbLNYq08p78sLaSv0TmkZiYyJwXv+aHzCKmDezOp3+Y2uhYVrORy8b04qPteVwyWkZyEh3H1q1bWbRoEQCLFy9ut/NERkbSs2dPzGYzUVFRDBs2jHHjfh4FZ+bMme12bnHiWgqkq4BpWuuiIwu01vuVUlfieXhKAulETbyFe1nE34p6sqOi/kjfnst0NZjZsXIRm2KGMCM2nB8ywWQ4dmfIt685vX3rFaINORwO7rzzTg4ePNhmx1RKYTabcTgcBAUFMWnSJEaNGsWAAQMICTn6WVLhz1oKJHP9MDpCa12olJKHZk6GJRTDWXfSe30+O7YUemaUbXDPSPGTeTS9tGJomIN3501g2pAevqpWiONWW1vL888/T3l5OTfddBMxMTEN1judTq6++uo2P+95553HpZdeitvtbjBJpuh4Wgqk2pNcJ1pwwchulNU4+HZvaaN1tVgIrs6guusoYsIiMSr5Tyb83759+0hP98zbuXbtWmbNmgVAVlYWjz32GOXl5Sd1XKvVyoABAxg9ejSFhYX86le/wmKx4HK5sNlshIZ6nu+TMOr4WgqkEUqppv4VKUCGqG4Fi8nAFeMTyTxUxqEKB456fxW7SKF/8Vd8UTMcW34uW3MruGVqsu+KFeI49O7dm+HDh1NWVsaYMWMAKCgo4O677z7pY1osFl577TVqa2u54YYbqKmpYc2aNVxxxRVMmjSpLozEyVNKGfEMHZSrtT5XKXUu8DCejmxm4B9a6xe8214J3AEY8TyjtA64TWvd+Dfrk9BsIGmtjc2tF60TbDZyx/RefPHlalaVJFJbbxLeFZFXc2TupCCz/OYn/J/VauWuu+4CwO12k5WVxTfffNOqY9bW1pKXl0dcXFxdt+3S0lKee+45CgsLmTNnDgBLly5l9+7dXH311R3uodfjdfnllzeafuLtt99uiwdj/wjsACK8t2JexNORLUcpZQWSAZRS5wDzgZla61xvkF0NdAPaJJDkJ52PRURE8Ms559NoAl6lQBlwajvXTJJ7SKLjsNlsPPnkk9x9992sWrWqVcdSSnHbbbexceNGRo5sOOTa999/D0BZWRnLli1j69atrT6fv/KGUaPpJ7zLT5pSqgcwG3jZuygcT0OlGEBrbdda7/KuW4CnNZTrXefSWr9ab32rtfhgrGh/docbp7I0OSn8ghl9MBvl9wbRMWRmZrJgwYK61kxLI8G05Mj+Bw4c4NJLL8VoNFJQUEBRUVFd6yg8PJwxY8awe/fuBl26A0xz00+0ppX0DJ5LcOEAWuvDSqnlQJZS6gvgQ+D/tNZuYAjwYyvO1SIJJD9gNil6RgeTVVzTqNfd2jVfE/uLM4mN6iLj3Am/s23bNr777jsmTpzI9u3bWbFiRatDCOC2224jJCSEwsJCcnNzmTVrFiEhIdx0002NtjUYDNx6662tPqefa/PpJ7z3igq01huUUmcdWa61/q1SahjwC+A2PPPfzTtq32HAm3iC7B6t9TsnW0d9Ekh+wKAU957bl0qbk+WvPcmXIbO9azQba+Kxr/yMuFAzZ5xxOvurFDaHizP7dfNpzUIAPPfcc5SWlpKWlnbc+0RERDTocZeQkMCIESPqLreNGTOG0aNHAzBw4MC2Lbjjao/pJyYB5yulZuHppBahlHpLa32l1noLsEUp9SaQgSeQtgGjgTTv+pFKqUVAcCtqaECuBfkJk0F55kPqM6XeUkWZO5jvnQP4vCyRtTszmfDUJ5z1j89ZtS3XZ7WKzsflcvHuu+/y3nvv4XZ7Jpo8dOgQtbXHfvpj8OCGo5GEh4fz6quv8tRTT9Grl+dnq8lk4oknnuCSSy4hJSUFs9nMtGnT2u+NdFxtPv2E1vpurXUPrXUyntkbvgSur99aouHUE38FnvTedzqizcIIpIXkdy6echqGDfms2HKIYPORXneKaoLYXe65ogfItOei/WV9ByWZMPwS0tPTef/99z2Ls7KIiYmhtLSU6uqff0YajUa01rjdbkJCQrjxxhtZvnw5xcXFZGVlMW3aNIKCPE+LTJ8+nf/85z+cfvrpGAwGLBYLCxculIdbj+Htt99+29t/oT162dWngDuUUi8ANUAV3st1WuuPvINsr/L2sCsFtgKftNnJ2+J679Fk+onW+8/3maTtbPgImBknfePDGZIYwrnDArNrq/ATpQfg2VHgdsL0R8npeQF33nlng/tDgwcPZvv27XWvDQYDbrebmJgYxo4dy+rVq7Hb7dxxxx2MGDHCF++iw+ns00/IryJ+6uoJyQyJrCaEGizYAHBgYkd+DSFmeSZZtI+KigpuvfVWfn/7wxxweKYzKXNZSU9Pp2fPng22rd86Auou5ZWUlPDxxx9js9nQWlNQcMyp04RoQC7Z+bHbLxxPzl9H0MV5mE8j5vJlyHnYLVH0jm3Ty7ZC1MnOziY/Px+ABRnjCbaa6fXVIbZt+xKlFN26daO0tBS73U7v3r1xuVyUlZU16KQQExNDYWEhffv2Zdy4cZx99tm+ejuig5FA8nOVl3+E5Y2z2BwxhUrdBVxucorKSI6VUYxF2xs4cCBnn302aWlpuFwuKqs948WB55mgQ4cOcfnll2O1Wlm8eDG1tbUopXjkkUdITk5GKUVFRQU7d+5k+PDhdfeMhDgeEkh+bmCvRF4Z+HcyK3vWzXr+xve5DCCDuIETfVucCDiHDx+mb9++2O121q1bR0hICNdddx0VFRUsW7YMs9nMlClTWL16dV0PO6016enp9O7dG/B0605NTfXl2xAdlARSBzBx3Di++SKfI4lUi4U718Adhe/T67RzCA6WS3ji+GmtycvLo1u3bphMJpYsWcKOHTuYN28ezz33HHl5eXXblpWVkZeXR2pqKgsWLKhb3rdvX6xWK3a7HfCEkBCtJYHUAfROiAUaTmjmNlj4595EzqtcLbNeihPy+uuv89lnnzFq1ChuuOEGli9fDsC9995LWFhYg23DwsJISUlpsKy2tpZnnnkGu93OL3/5S4YPH86AAQNOWf0icEkvuw7A4XQzKVFh0g7qD3hXTRAHquR3CuHp8bZq1SoyMjKOuf6TTz4hMzOTnJwcAHJycti7dy+RkZF1XbarqqqYOPHnS8H33HMPcXFxDY7ldDqpqqoCPDPAShh1bEopo1Jqo1LqQ+/rc72vf1JKbVdKXedd/oBSKlcptcn75zHvcpNS6lGl1J566xY0d85jkZ9mfu5QuZ37l+/B6YLfuN/hHeMcSvn5t9iNZSE4XW5MMgBrp/bWW2+xevVqgoODeemllxo9XLp48WLS0tIIDg7mwQcf5JtvvsHlcvHEE0802M7lclFTUwN45iJqagrwkJAQFixYQEZGBmeddVa7vSfR0LzXNjeafuL1a4afsuknvJ7WWj951P6PAN2BYVprm1IqHDipwQXlp5ifK6yoxeZw43RrnKf9hkqsDdbbDOHcu3z7MfYWncWRezjh4eGNBuF1u911LaegoCASEhK47LLLGjxHNGrUqLp5hDIyMrj33nt5/PHH6dq1a5Pn69+/PzNmzMBqtTa5XrQtbxg1mn7Cu/ykneD0E03tHwL8DrhZa23z7lOhtX7gZOqRFpKfG5oYzhXjErA73UwcGscXW78h09GlwTbrd2fyj3e8uW76AAAdz0lEQVRz+N3syXW/0R49vbMIDC6XC5fLhcViabB87ty5jBw5ksTExEaBlJ+fT2ZmJgCpqank5+ezZs0aJk+eTFhYGIMGDWLUqFEsX76cJUuWUFpaisVioVs3GcDXj/jD9BMA872zxgLcCeQDB7TWFa2ooY4EUgcwbXBs3fc3TB/Iqyu/Z5fuCcpEvLGc2b1rUErxww8/MGXKFLTWpKWlUVpayvDhw+nfv78Pqxcnq7q6muLi4roREiorK1mwYAHl5eXcfffdDf5eDQZDo5GxP/jgA/bs2cNll11Gz549yc/PJy4ujueff56MjAw2btzIX//617rtp0yZQn5+PjExMfTp0+fUvElxvPxl+okGl+yUUsOPOuY1eC4BxgATtdbZJ1KTBFIHk5mRQZkpFpwmzLjId4XzWWUwY4PyCK51AZ5LNEeenC8rK/NlueIkORwO7r77bgoLC7nqqqs455xzOHjwIIWFhQDs3r27LpCOjMQNcPHFF2MymSgtLeW///0vANHR0ZSVleF0Ovnhhx+Ij48nIyOD7t271+3vcDgICwvjuuuu88G7FcfBH6afaMpeIEkpFe69VPca8JpSaitgPNGCJJA6mJSUFIz7MsAJDq1AKWoxs8bWi4hCzUy3xmg0MmHCBAoLC6V11EE5HA4OHz4MeKZ5AM+zP3PnzuXw4cNMmeKZpiQ7O5v//Oc/dYOcJicnM378eCIiIhg2bBh79+5lzJgxREZG8uWXXzJlyhTOOOMMzj//fBISErDZbCxYsIDCwkJuu+02hg8f3nRBwtfuwXMPqf5lu1ZPPwHcDeBtId2Gd/oJrfVq72b1p59o6hjVSqlXgEVKqeu8nRqMgOVY+zRHRvvugKrsLv617FO22o4e8Vszr08FZ00+3Sd1iba1ZcsW9u3bx/Tp0xv0dnO5XBiNRqqrq7nxxhvrRkwIDQ3loYceIj4+/rjPkZ+fXzfb6pw5c5g7d27bvglxQpob7bsde9nVD6TLgHeAPvw8/cQftdbrlVIPAJVH97Lz9sx7GLgYqPDutxL4m9b62BNmNVWHBFLHVFNZxg3/3Q+qfqtYM9RawNhEC0FBQYwdO1bmlgkwixcvZuXKlZx//vlMnDiRu+66C4Bu3brx+OOP13V22LdvH0qpuuF8mrNq1Sry8vKYO3cu4eHh7Vq/aF5nn35CLtl1UDYVAurosFFstXdD5R6il6WQPn36EBsb2+T+omP68ccfAdi4cSOXXnopV1xxBdnZ2VxxxRV1YbRr1y4eeughAO6///4WL9vKSB/CX0ggdVA1Dhd1o60eZa+9CyFBZoLDZHyxQDNv3jw+//xzpk+fDsDs2bPJzc1l0aJFDBgwgAsvvJDa2tq6ifSam2JcCH8jl+w6sHfX57FyS9Ex14/ubuCXQyMbTawmAstLL71EWloaAP/+97+JiIhg48aNgOeBV9FxdPZLdnKDoQO7cHQ8oc7Dx1xfdriItWvXerp+1z3XJgJNamoqoaGhjBo1qm5w1FGjRkkYiQ5HAqkDMxoUsV2OPXRLd0M5ZrMZ67uXwV8TIeOrU1hdACjYDuW5ALzy3V4SFyzj4VVbfFxUYyNGjOCll17i9ttvl04sokOTf70d3IO/msBvJvVA4ab+SOAAe+1RjOzVhaDcNeCohn1f+qbIjmjHCvjXBHhuHJTn8cKaveSV1fD8N7t9XZkQAUsCKQCc0T+aeQPthGBvsPyQjmD11ny2Df0LDL0Ixl3f4rG2bt3KDz/8UDfic6dVke/5WlsJ9nLumT6E05KieXBW+z44WllZycKFC3n44YfrRtsQorOQXnYB4swJ40jMuptHai7l5953ijTHQEKsJip7jGdsaLdmfwMpKytj586dgGfU6CFDhrR32W3DaYeSDIgdAKrpnocn7LRrPceKSIS4gcyJgzkj2r9zyNatW9m2bRsAmzdv5vTT5SFn0XlICymARExtagoSA+tr4snOPUTh0ttgxS2eH95NCA0NJSoqCpPJ1LFGel58ETw/Dj456VFUGvlow25+uWg9lz39KRUVlU1us3PnTlavXo3L5TquY65evZqlS5fWTfvdlKFDhzJ48GAGDhwow/iITkdaSAGka1wsZg7gOOqvtcxlBbeTguLDBJWtJXJfGvzJc3O+urqa9PR0rFYrqampTJ06Fa11oykM/FrRHs/X4ra7v/Pe52sIdlSAo4J1W3YwZeLYBuuLi4tZuHAhLpeLiooKzjvvvGaPd+DAAV588UXAMyfRueee2+R2YWFh/OUvf2mbNyFEByMtpABz87S+jZbZMIHBxK64c/ms78NkR4yuW5eTk0NRURG5ubl1g3l2qDACuPRtOP1WmP10mx3y9xedQ3VkD4KSBnHG2JGN1pvN5rqREVqacyo7O5uXX34Zq9WKwWCgR48ebVanEIFEWkgBxOl0UlmYC1ofdS/FQIEjmK5mT0eFgiG/I97pRGtNQkICWVlZBFXlELVkNpzzCPSb7ps3cLISx3j+ABTuArcDug1t1SF7RwXx2qN3ExUV1eT6iIgIHnvsMYqLixvMQ5SdnY3ZbK6b2gHg008/Ze/evQAsXLiQlJSUVtUmRKCSFlIA2bFjB1t372vyxn66rQdlDjOgycg+yIoVK1i5ciVOp5NpZ5/JGT/Nx3R4F3z5MLw6HXYsP/VvoLUOboF/T4R/nw6Z3570YXbt2sWf//xn/vznP5Ofn99o/TNpO7nopa8pI6hBGG3bto277rqLO+64g5ycnLrlEyZMIDo6mtTUVHr1ampKGyEESAspoISEhBBlsBGpqijTR19GUuyqjSXV7PkBe+RGfOnnfyeLCMp6zWd03puEFe8HRyXYK2DQ+af4HXgV7YbP7oXkM2DCH45/P3s5uJ2e72tKjrnZtm3b+Oyzz5gyZUqTHQcKCwvrpgovKSlpMJ1DcaWd+e9tACA61MKii0bx0ksv4Xa7GTFiBFprnE4nb7zxBkajkejoaK699loWLVp0/O9DiE5KWkgBpE+fPqSePpmP9hU0ub7IHYpTe1pP1trD9Cz5gag9/2UPyRSEDWF/9FnQ4zQIjoaRV7SuGHs5bFoMJZme14f3eVpfBze3vO+aZ2D3x/DpArCVHv85e02CSxbDha/wwP5eBM9fwkNNjKzw6quvkp6ezssvv9zkYSZOnMivf/1rfve73zF48OAG67qEmDm7fzdCDS5GB1WQnp7Od999xw8//IDBYGDcuHGAp/v2Tz/9RFpaGvv27Tv+9yBEJyYtpACT0DWWsckJTfaUU2gKnKHEmyqxW6LJNo+je0gt3WIiKT1cTLeB42HyNWBoOPNweno6hYWFnHbaaWTVGFm+JYffT+pHUs1OiOwJYV0bF/LhfNi6FMITYP42+N9NcOB72PY+3Pxjo823bNlCbm4OI4YNIX7AbM92yaeDNfLEPoCBnt5rSz5Ygc3h4v/WZ3LfzGENNjny0K/D4WjyEAaDodGUDA6Hg8rKSqKiovjyll+wcOFCvln6MSVDh5KYmIjL5WLw4MGkpKSwb98+goODqampITo6mqSkpBN7D0J0UhJIAWZ/YTUWg3d8u6M6N7gw8KMtgX6WYgZYi0Ep1lknk2ANxa7LSK/oxszqEkyL50DpARh4Ho5fLOTAgQOAp+vyBe8f4GC5jR673+T6ypchtCvcsgksR10iNAV5vlbkwfvXQUxfTyDF9mtUs9aaXbt2eer/7BXiz7kI7ml87+ZEPHHBKBZ9vYubzxzQaF3fvn1Zv349ffs27pHYFLfbzX333UdWVha/+c1vmDp1KjabDfBc+vzb3/7WYPt//vOfrapdiM5KAinAJMUE07drCDsPVmBUhiZnTHIftTQvLw8Au93O8k9W0yX4XLq7tjJo01uY4wYwcOAUDh48SFlZGdcPMPHQOugXVAWVeC6pOe2NA2n2U5D1LZRkUlWQRebpfydhyDVEJTdsrYCnm/nA7sHkZe6lb9GnkJsC3Zt5KPSj22DzOzDjrzDqyiY3OX94D84f3nT36j/84Q9kZmaSnJx87HPUU1tbS3Z2NgD79+9n6tSpzJ8/n40bNzJmzJjjOoYQomUyH1IAOlxVyx1Lt+N0H32L0PN37WkhHXvaiiPOyPo73eJ7wKVvs23HLnbs2AFA8vBUxiTHoH58A7oP81xaa8qhrbDpbb42TqKgtIaQkBBmzZrV9LbaDasf89x7mnofmEOOXdijCeCogt5nwa//1+L7aAvr1q1j586dnHfeeXTp0uWUnFN0Pp19PiRpIQWg59MONBFGcGSMu3xHeBOBpD2hoDz3j0xuG+G2PNizDbK+IzFxONnZ2Z55d/okooxGGH9D3d4Oh4M1a9ZQXl5OeHi4Z46ebkNhxqOEb9xIQem+5h8gVQY4+ziH/pmxELa+B5PvOL7t28DYsWMZO3ZsyxsKIU6a9LILQBZT83+tdu3ttNCgdazqwkjhYlzmInbFnkN51/EQP5wuhes4Z2AYZ5xxBkajsdExi4qKKCoqora2luLiYjIzM+vWJSYmEh4eTlxcXKP93G7Ns6t38er3P/dEa7HVPuYauGoFpL8IT6R4euQJITo8aSEFmPc3HmR7ftODgXou2SmCDd5ndY50eDiq84PGyJqU+aCMZOHi/E1L2LfuU0Ich0mcdTsMnN3oyHFxcXTv3p3i4mLMZjOJiYl16/bs2UNFRQU7duxg0KBBDXr/LdmQyS1LPZd3U2JCWZtZzF8+/ImbzxzA0xc1c3+mtgK2f+D5ftv70P+cFj8bIYR/k0AKMO/+mIvFYD7GAKmKKEMVqSFH9WCrv11dOHlaWU5t4MOccGrjLwdg9Oq36JXxDcZzHqWktIywsDDMZjMmk8kzVcKhrZ4edvXusyQnJ1NSUkKPHj0a1ZQcE4bJoDAbDSREhvDepo243JqlGw/UBVKV3Ynd6SI6tN7suNYImHKvZxbcE3l4Vgjht6RTQ4AZ8OBKJiQ1N2+PmwnBucSYWpiAz+1q+DySN6iU28msPXey5rRnKa2oJSgklCm/mEaIxeQJhzcuAIMJrv8W4n4eVoeaEsjb6OkAYbQ0OFVmcSUWo4GELiGk7T7IU1/s4DcT+/LLET3JL6thxF9XUm5z8MXNU5nUp4lnnlpQWFhIZmYmBQUFJCUlMWxY455+QvgD6dQgAsam7HImpyRhd7gwGI51H8lAjjOi5UCqF0YGtwOjy47DHIZWRopVFJv255IcF0dVRRkJ9yxjza3TGVJdDGjP4Ka2sobHe+uXnkAacTnM+VeDVckxYXXfn92/O2f3/3lg0v1FFRRWeuYP+im39IQDqba2lr/85S9UVFQAYDQa+ec//3nMQVOFEL4jgRRA1maU4nDpZsIIwE0v07HHeauj3YDCqGvpc/hLTK4q9sTMJLn0WxIdGUwreBGLeyAhlRn8wj2RPq/fAUPnwJwXPM8k9RzX8HjVxZ6vxccYRuf7RVC8l/dj5rF4Wzn3zBjK6J7RTOrTlb9fOJrCSjvXTuhzPB9DHZvNxn333VcXRuCZKsJqtTazlxDCVySQAsi0QTGs3V+Ku8mrsJ4ODTFGG9GGCtzacsxODZ7ed4oIWzYj898Crfmh1y04TKGE23LQysgIvR93USbfdLucV5zvE1RbAj++Duf+nU93HuLGB/7HzMEJPDvX21U6fqRn9Ie8jZ4HaU31QqFwl2fcOmCHvZD3ymexPuswBoPi2V+dxvwpg07q8ygqKqobdXvMmDFccMEFxMXFERLSzDNOQgifkW7fAaRXTAihlsZdso9IMJUz1rADt8HaMICO7vygFChFeXASXyffwdd9FlBr9FxW+zHxGj7q9zg2UwQYTJxh2YfZ4enVVxrcA5SBV77fy76iShZ9vRubwzu995H5imJ6g9Hc8HwRCZ6hhYwW6DEegAMlVWQUV/La0uXceOONfPDBByf8efTo0YOLL76YSZMm8dvf/pa+ffsSGXmCY+MJIU4ZaSEFkJwSG5V21zHX9zMXYzKGYXDZcBuDfl7RaEK/essMnn8iBu3ArSygDNjMUZRb4omr3oUzdyM/Wsayu1xT2Pc33AncNLk/uwsqmDk4gSCzNyAn/Qn6z4RIT2g1YA2HG9eC0849llBud7lZ+PFW/rclh/jCnyguLeXzzz9nzpw5J/yZXHjhhSe8jxDCNySQAsiSdTkcu8+kwq0MoFTDMAKCawuoscR5guJIEB0VUCPz38JtsFBuTSTIWUZ01W5cysDXKXdSEpREFeFcN+tsACb37cbGu5oYIiiu8UCndQwmsHj+OZqNBh6YPZwHZg8nPT2W5cuXM3Xq1CZ3S09Pp7q6mjPPPLPjTb0uhGhAAimA7Mg/ds85w6GthIftIDSqO8ml39LFlkNG5Hh+zKpgbWYVycnJmM1mnA4Ho/sn4Az2jKoQWlvA6Lw36FblGcfO4Vbcsn8qZa5ZjBw/lhhjTyKBkvJKuoRYjnn+luTk5PC3v/2NqKgo7rrrLoKCPKGZmppKampqk/vs3r2bZ555BgCz2cykSZNO+vwOh4N33nkHo9HI3LlzmxyNQgjRviSQOgHtdmPbsoo/hY3m9dg36G/ehVFpulZtZ3HBHFwuV4NJ5CZVrSQlJYUdrhRMhRuINe2ou9u4syaaMpenQ8Ln63azK6EbqdEuptesgm8PwOnzT6rGJR+vprCwkMLCQg4cOED//v1b3Cc0NNQTok5nq+8Npaen89FHHwHQr18/Tjut0z4KIoTPSCAFiHs/2N30pHxuG733vsnGoFCCXDZ+f+iXJDiLSAy38+vwbQwYMoINGzbUTVZnNsIAayFRBdks3G/CoeOpjUthdnQGH9kHs7ysN6BxodgW3Je8Ks1F7jR+F7aS2s9W8ejHB/kpu5RDvc9i5e0XEhXSsIv1/v37sVqtJCYmkpaWRklJCaedMZUHNtYy1hRLn7jwBsMOAZSXl/PKK68QHh7OvHnzMJk8/2wTExN54okncDgc9OjR9FQTxyslJYXQ0FCMRiO9evVq1bGEECdHAikAVNQ4yS6xNXkPJa+ogE01iQS7ihhfu8Xb+Rso1iwp7k5U/nf8amRX9mblEpuYQlD8YCguhaJ1dbMmObSB28ov4B+2Kfyq4lMADlgSOGTtDhrWOpKp1mbyzcns3pNNMFB+YBczn0tjxW8nEhcVAcCmTZt44oknMBqN3HTTTbz00ksAGC1BhEVE4q5UVOdnsGjRIu6888669/Dcc8+xZYtnKvKEhIQGU1h069atTT7DhIQE/vUvzwO7RwJPCHFqSbfvDs6tNTcv2d5oudYard1sKa5hk7MnZrPnh+zPkeX5rsRuYPmGXGoqy8lzRLJ+627WmycQaarl4V7f8tuum1la1J+sfDeDq/dw2BCOSxlJTU3F9sxlfHHzVOg5jtt6/5ek274mdcIkqkK7kmWOx7HtK/540/U88MADPPHEE3WT3LlcLoKCgggLC8NoNNInOYkt98wmNcHTtbyqqqquylWrVtWFEUB8fHyT77WwsBC3293483G7KSk5jgeB8QSRhJEQviP/+zqwvQVVPPnJ/ibX6cK91G5ZwcHIqaTEhrPoit/z6MMPYnLZsJiM2N2Guh/glW4rW6rjYKun44LtkJ05SYos1ZVllum4ycaMmyH2DABcscnMm9CXj1Z+SI9h4/g+s4jvMmBUnzz+dPNNpG7cSvBTT6C0pwv67t27AU/HBZPJhFKKJUuW8OCDDxIUFFQ3jM9dt/2Z9evXN7h/k5HhOafBYODhhx8mJSWl0Xt96aWXWL16NZMnT+b666/H4XDw3HPPUVpaitFoZMeOHVx88cXSBVwIPyeB1IH9+6sD2JyNO3rXVhxkVg87A3/5EMMzyjlnUAJ/fuEd4lwOXBg5q2sxq/JivFtrapWFWmXG4q7FgpMCu5Xf7Z3Ogah+TO7RhcM52SilcGuNAiLsh3nqqadwu92cVVKO1WTC5nARG+a5X5S1cwsGt2eKi5CQELp3787+/fspKiqqq/HAgQO89dZb3H777XXLYmJimDFjRoP3cumllxIeHs6QIUOaDCP4ObSOfN23bx/p6ekAWCyenn979uw5wU9XCHGqSSB1UDvyyyiqdDRaXliYS+SmN/kIGJc6lmG1Gdz/4jdYsjbiwEi2pTuvm0ZxTvf1FFVBZEp/dmzLxOKqbnAcuzbR7XAGWdVBJCUlkZqaSnx8PN999x2TJ0/mtddeo7S0lP69erB11hhKqms5rVcMFRUVfPrpp3XHMZvN3HXXXfz+978HwGq14nA4cLvd2Gw2NmzYwO7du5k9ezYRERGN3k90dDS//vWvm/0srrvuOr766ivOOOMMAHr37s3IkSMpLS3lnHPOYe/evcycOfNEP2IhxCkm0090MG6tefzj/ew6WNVonVVXMTXRwf/eeAGr1cq0adNYsWLFzx0ZgMuvuoZzz5lWt4/NZmP+/PmUlXlG546NjSUmJobS0lIOHTpUt53RaOSVV16pa3FUVFRQVlZW17stLS2NsrIyJk+ezPz583E4HAwZMoRzzz2XzMxMlFJkZ2dzySWXUFJSwrp165g4cSL33nsvLpeLadOmcc0117TPhyZEByHTT4gOpabW3WQYGbDz/FVjMZpMTBrah6CgoLrLVNbgEGojE+kXF8aUyQ0fHg0KCuKxxx4jPT2d1157jaKiImbMmIHBYGDx4sV195mioqIa3PAPDw8nPDwcgL1799b1mAsNDeWBBx4gPz+f8ePH8+yzz7J27VrCwsJ48cUXAU/o9evXD6fTSdeuXcnPzycpKantPywhRIcigdTB/PPzjEbLekZZuf+8oRiNnk6TR57jiYmJoVevXoSGhjZ5OeyIyMhIUlNTeffdd6mqqsJqtfLGG2/gdruxWq3Y7XaKi4spKCige/fujfaPjo4mNDQUm81GQkICKSkpdfd7YmNj62o5mslk4tFHH6W8vJy4uLgT/zCEEAFFAqkD+WpXAbsKGt7rCTErHp5z7DHimuomXZ/NZmPNmjX069ePp556iurqajZv3ozT6emUMGDAADZv3kxoaOgxp22Ijo7m6aefpra2lujo6AbrLr/8csaPH09CQkKT+1qtVgkjIQQggdRhfLr9IIt/ONTg4ddIq2LhRSc3V9ARixcv5osvviAkJIQXXniB8PBwJkyYwLZt2wgJCeHaa6/l0KFDhIWFNdvKCgsLa3K5Uoo+fU5sYj0hROckgdQBbDpQyttrDzUaieHpS4diMLRuhOsjrZ7g4OC644eFhfGnP/2pbpujh/IRQoj2IIHUATz9RRaKhsHzpylJrQ4jgEsuuYShQ4eSlJTUwtTnQgjRviSQ/NyKnw7WhZHWmsHdQ7lpSjJhQW3zV2cwGBg2bFibHEsIIVpDAsmPVdqcfLCpsO51jygrd87q68OKhBCi/Ugg+bEahxuX2/PgclJ0EA9d0PIcQUII0VFJIPmxuHALf5yazKFyO1MHNX6ORwghAokEkp8blXTsrtZCCBFIpFuVEEIIvyCBJIQQwi9IIAkhhPALEkhCCCH8ggSSEEIIvyCBJIQQwi9IIAkhhPALEkhCCCH8ggSSEEIIvyCBJIQQwi9IIAkhhPALEkhCCCH8ggSSEEIIvyCBJIQQwi9IIAkhhPALEkhCCCH8ggSSEEIIvyCBJIQQwi9IIAkhhPALEkhCCCH8ggSSEEIIvyCBJIQQwi9IIAkhhPALEkhCCCH8gtJat/1BlSoEstr8wEIIEdh6aa3jfF2Er7RLIAkhhBAnSi7ZCSGE8AsSSEIIIfyCBJIQQgi/IIEkAoZSyqWU2qSU2qaU+kkp9WellKHe+lSl1NdKqV1KqZ1KqZeVUiFKqYFKqe+VUnal1G2+fA9CdGYmXxcgRBuq0VqPBFBKdQXeBiKB+5VS3YB3gUu11t8rpRRwERAOHAZuAeb4pmwhBEgvOxFAlFKVWuuweq97A+uAWOBBAK31fc3s/wBQqbV+sp1LFUI0QS7ZiYCltd6P5994V2AosMG3FQkhmiOBJAKd8nUBQojjI4EkApb3kp0LKAC2AWN8W5EQojkSSCIgKaXigH8Di7TnRuki4Gql1Lh621yplOruqxqFEA1JpwYRMJRSLmALYAacwJvA37XWbu/6CcATeO4puYGvgflABLDe+9UNVAKDtdblp/o9CNGZSSAJIYTwC3LJTgghhF+QQBJCCOEXJJCEEEL4BQkkIYQQfkECSQghhF+QQBJCCOEXJJCEEEL4hf8Hy5eVv6Ex4tAAAAAASUVORK5CYII=\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"del ad.uns['exp_groups_colors']\n",
"matplotlib.style.use('tableau-colorblind10')\n",
"sc.pl.diffmap(ad, color='exp_groups')"
]
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 3",
"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.7.1"
}
},
"nbformat": 4,
"nbformat_minor": 2
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment