Skip to content

Instantly share code, notes, and snippets.

@yuu-ito

yuu-ito/コーヒーの注文待ちをシミュレーションしてみた.ipynb

Last active December 3, 2016 04:10
Show Gist options
  • Save yuu-ito/a72dc0a2b45ed51c952bddc4984d390d to your computer and use it in GitHub Desktop.
Save yuu-ito/a72dc0a2b45ed51c952bddc4984d390d to your computer and use it in GitHub Desktop.
Download ZIP
コーヒーの注文待ちをシミュレーションしてみた http://techlog.voyagegroup.com/entry/advent2016
Display the source blob
Display the rendered blob
Raw
コーヒーの注文待ちをシミュレーションしてみた.ipynb
{
"cells": [
{
"cell_type": "markdown",
"metadata": {},
"source": [
"# コーヒーの注文待ちをシミュレーションしてみた\n",
"\n",
"これは [VOYAGE GROUP Advent Canlendar 2016](http://techlog.voyagegroup.com/entry/advent2016) の3日目のエントリです。\n",
"\n",
"こんにちは、最近RよりもPythonを書いてお仕事している yuu-ito です。\n",
"\n",
"弊社ではバリスタさんが美味しいコーヒーを淹れてくれる「ガーデン」というサービスがあるのですが、\n",
"\n",
"オススメです。もし弊社に商談などで寄る際はぜひ注文してみてください。\n",
"\n",
"https://www.instagram.com/officecafe_garden/\n",
"\n",
"今回は最近業務で使い始めたPythonでコーヒーの注文待ちのシミュレーションしてみます。\n",
"\n",
"Python シミュレーションフレームワークである[SimPy](http://simpy.readthedocs.io/en/latest/)を使って実装します。\n",
"\n",
"コードはこの記事を読んでわかる通り[Jupyter](http://jupyter.org/)で以下に書いていきます。\n",
"\n",
"(Jupyter notebookって本当に便利ですよね！)"
]
},
{
"cell_type": "code",
"execution_count": 1,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"import matplotlib.pyplot as plt\n",
"import seaborn as sns\n",
"import pandas as pd\n",
"import numpy as np\n",
"import simpy\n",
"\n",
"%matplotlib inline"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"STATE_FILE = \"./state.csv\"\n",
"def log_state(m):\n",
" # 実行結果をファイルに追記\n",
" with open(STATE_FILE,\"a\") as f:\n",
" f.write(m+\"\\n\")\n",
"\n",
"def get_state_results():\n",
" return pd.read_csv(STATE_FILE)\n",
"\n",
"def reset_state():\n",
" # ファイルのクリア\n",
" with open(STATE_FILE,\"w\") as f:\n",
" f.write(\"name,time,state\\n\")\n",
"\n",
"class Person(object):\n",
" def __init__(self, env, name, freq, coffee_garden):\n",
" self.env = env\n",
" self.name = name\n",
" self.freq = freq\n",
" self.coffee_garden = coffee_garden\n",
"\n",
" def needs_coffee(self):\n",
" # 一定の確率でコーヒーが飲みたくなる\n",
" if np.random.randint(SIM_SIZE) <= self.freq:\n",
" log_state(\"%s,%d,コーヒー飲みたい！\" % (self.name, self.env.now))\n",
" return True\n",
" else:\n",
" log_state(\"%s,%d,仕事しよ...\" % (self.name, self.env.now))\n",
" return False\n",
"\n",
" def run(self):\n",
" # シミュレータが実行するメソッド。毎ループごとにコーヒーを飲みにいくか判断する\n",
" while True:\n",
" if self.needs_coffee():\n",
" with self.coffee_garden.request() as req:\n",
" yield req\n",
" yield self.env.timeout(1) # コーヒーを淹れる時間\n",
" log_state(\"%s,%d,おいしかった！！\" % (self.name, self.env.now))\n",
" yield self.env.timeout(1)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"ざっとコード貼ってしまいましたが、\n",
"\n",
"\n",
"大事なとこは **Personクラス** でもつ **env** 変数と **run()** メソッドです。\n",
"\n",
"このあとメインとなるシミュレーションの実行処理を書くのですが、その処理（process）が\n",
"\n",
"上の **Personクラス** のオブジェクトから **run()** メソッドを実行するように設定します。\n",
"\n",
"- http://simpy.readthedocs.io/en/latest/simpy_intro/process_interaction.html\n",
"\n",
"以下がシミュレータの準備のコードになります。"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"SIM_SIZE = 10 # ループを回す回数\n",
"FREQ = 1 # Personがコーヒーを飲みに行きたくなる頻度\n",
"NUM_OF_PEOPLE = 5 # Personの人数\n",
"\n",
"# simpy のEnvironmentオブジェクトを生成\n",
"env = simpy.Environment()\n",
"\n",
"# コーヒースタンドオブジェクトの生成\n",
"coffee_garden = simpy.Resource(env, 2)\n",
"\n",
"# Personオブジェクトを生成。runメソッドをプロセスとして設定\n",
"for name in [\"person_%03d\" % i for i in np.arange(0,NUM_OF_PEOPLE,1)]:\n",
" env.process(Person(env, name, FREQ, coffee_garden).run())"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"SimPy 特有の処理の書き方ですが、\n",
"\n",
"複数のプロセス（今回だと**Personオブジェクト**）が共有するモノに対して処理を行う場合、\n",
"\n",
"またそのモノに容量の制限があるオブジェクトを用意したいときは[simpy.Resource](http://simpy.readthedocs.io/en/latest/simpy_intro/shared_resources.html)オブジェクトを利用します。\n",
"\n",
"今回だとバリスタさんとコーヒースタンドの場所がそれにあたります。"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"それでは実際に実行してみましょう。"
]
},
{
"cell_type": "code",
"execution_count": 4,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"reset_state()\n",
"env.run(until=SIM_SIZE)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"csvファイルが吐かれているはずです。読み込んでデータを可視化してみましょう。"
]
},
{
"cell_type": "code",
"execution_count": 5,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"text/html": [
"<div>\n",
"<table border=\"1\" class=\"dataframe\">\n",
" <thead>\n",
" <tr style=\"text-align: right;\">\n",
" <th></th>\n",
" <th>name</th>\n",
" <th>time</th>\n",
" <th>state</th>\n",
" </tr>\n",
" </thead>\n",
" <tbody>\n",
" <tr>\n",
" <th>0</th>\n",
" <td>person_000</td>\n",
" <td>0</td>\n",
" <td>仕事しよ...</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1</th>\n",
" <td>person_001</td>\n",
" <td>0</td>\n",
" <td>仕事しよ...</td>\n",
" </tr>\n",
" <tr>\n",
" <th>2</th>\n",
" <td>person_002</td>\n",
" <td>0</td>\n",
" <td>コーヒー飲みたい！</td>\n",
" </tr>\n",
" <tr>\n",
" <th>3</th>\n",
" <td>person_003</td>\n",
" <td>0</td>\n",
" <td>コーヒー飲みたい！</td>\n",
" </tr>\n",
" <tr>\n",
" <th>4</th>\n",
" <td>person_004</td>\n",
" <td>0</td>\n",
" <td>コーヒー飲みたい！</td>\n",
" </tr>\n",
" </tbody>\n",
"</table>\n",
"</div>"
],
"text/plain": [
" name time state\n",
"0 person_000 0 仕事しよ...\n",
"1 person_001 0 仕事しよ...\n",
"2 person_002 0 コーヒー飲みたい！\n",
"3 person_003 0 コーヒー飲みたい！\n",
"4 person_004 0 コーヒー飲みたい！"
]
},
"execution_count": 5,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"dat = get_state_results()\n",
"dat.head()"
]
},
{
"cell_type": "code",
"execution_count": 6,
"metadata": {
"collapsed": false,
"scrolled": true
},
"outputs": [
{
"data": {
"text/html": [
"<div>\n",
"<table border=\"1\" class=\"dataframe\">\n",
" <thead>\n",
" <tr style=\"text-align: right;\">\n",
" <th>name</th>\n",
" <th>person_000</th>\n",
" <th>person_001</th>\n",
" <th>person_002</th>\n",
" <th>person_003</th>\n",
" <th>person_004</th>\n",
" </tr>\n",
" <tr>\n",
" <th>time</th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" </tr>\n",
" </thead>\n",
" <tbody>\n",
" <tr>\n",
" <th>0</th>\n",
" <td>仕事しよ...</td>\n",
" <td>仕事しよ...</td>\n",
" <td>コーヒー飲みたい！</td>\n",
" <td>コーヒー飲みたい！</td>\n",
" <td>コーヒー飲みたい！</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1</th>\n",
" <td>仕事しよ...</td>\n",
" <td>仕事しよ...</td>\n",
" <td>おいしかった！！</td>\n",
" <td>おいしかった！！</td>\n",
" <td>None</td>\n",
" </tr>\n",
" <tr>\n",
" <th>2</th>\n",
" <td>コーヒー飲みたい！</td>\n",
" <td>コーヒー飲みたい！</td>\n",
" <td>コーヒー飲みたい！</td>\n",
" <td>仕事しよ...</td>\n",
" <td>おいしかった！！</td>\n",
" </tr>\n",
" <tr>\n",
" <th>3</th>\n",
" <td>おいしかった！！</td>\n",
" <td>おいしかった！！</td>\n",
" <td>None</td>\n",
" <td>仕事しよ...</td>\n",
" <td>仕事しよ...</td>\n",
" </tr>\n",
" <tr>\n",
" <th>4</th>\n",
" <td>コーヒー飲みたい！</td>\n",
" <td>コーヒー飲みたい！</td>\n",
" <td>おいしかった！！</td>\n",
" <td>仕事しよ...</td>\n",
" <td>仕事しよ...</td>\n",
" </tr>\n",
" <tr>\n",
" <th>5</th>\n",
" <td>おいしかった！！</td>\n",
" <td>おいしかった！！</td>\n",
" <td>仕事しよ...</td>\n",
" <td>コーヒー飲みたい！</td>\n",
" <td>仕事しよ...</td>\n",
" </tr>\n",
" <tr>\n",
" <th>6</th>\n",
" <td>仕事しよ...</td>\n",
" <td>仕事しよ...</td>\n",
" <td>仕事しよ...</td>\n",
" <td>おいしかった！！</td>\n",
" <td>コーヒー飲みたい！</td>\n",
" </tr>\n",
" <tr>\n",
" <th>7</th>\n",
" <td>コーヒー飲みたい！</td>\n",
" <td>仕事しよ...</td>\n",
" <td>仕事しよ...</td>\n",
" <td>仕事しよ...</td>\n",
" <td>おいしかった！！</td>\n",
" </tr>\n",
" <tr>\n",
" <th>8</th>\n",
" <td>おいしかった！！</td>\n",
" <td>仕事しよ...</td>\n",
" <td>仕事しよ...</td>\n",
" <td>仕事しよ...</td>\n",
" <td>仕事しよ...</td>\n",
" </tr>\n",
" <tr>\n",
" <th>9</th>\n",
" <td>仕事しよ...</td>\n",
" <td>仕事しよ...</td>\n",
" <td>仕事しよ...</td>\n",
" <td>仕事しよ...</td>\n",
" <td>仕事しよ...</td>\n",
" </tr>\n",
" </tbody>\n",
"</table>\n",
"</div>"
],
"text/plain": [
"name person_000 person_001 person_002 person_003 person_004\n",
"time \n",
"0 仕事しよ... 仕事しよ... コーヒー飲みたい！ コーヒー飲みたい！ コーヒー飲みたい！\n",
"1 仕事しよ... 仕事しよ... おいしかった！！ おいしかった！！ None\n",
"2 コーヒー飲みたい！ コーヒー飲みたい！ コーヒー飲みたい！ 仕事しよ... おいしかった！！\n",
"3 おいしかった！！ おいしかった！！ None 仕事しよ... 仕事しよ...\n",
"4 コーヒー飲みたい！ コーヒー飲みたい！ おいしかった！！ 仕事しよ... 仕事しよ...\n",
"5 おいしかった！！ おいしかった！！ 仕事しよ... コーヒー飲みたい！ 仕事しよ...\n",
"6 仕事しよ... 仕事しよ... 仕事しよ... おいしかった！！ コーヒー飲みたい！\n",
"7 コーヒー飲みたい！ 仕事しよ... 仕事しよ... 仕事しよ... おいしかった！！\n",
"8 おいしかった！！ 仕事しよ... 仕事しよ... 仕事しよ... 仕事しよ...\n",
"9 仕事しよ... 仕事しよ... 仕事しよ... 仕事しよ... 仕事しよ..."
]
},
"execution_count": 6,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"# データフレームを整形して見やすく。\n",
"dat.pivot(\"time\",\"name\",\"state\")"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"意図通りに動いているぽいので、シミュレーションのシナリオを整理して設定してみます。"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## シミュレーション・シナリオ\n",
"\n",
"- 社内の就業時間(8時間)中に注文できる。\n",
"- 毎日行く人もいると思うけど社員全員ではないはず。1日 0〜2回のランダムで\n",
"- 社員数は300人。だいたい弊社こんなものだったはず。\n",
"- バリスタさんは3人常駐してくれている。\n",
"- 注文から淹れるのに3分ほど時間がかかるとする。\n",
"\n",
"1ループを1分としてシミュレーションを行うとすると、8時間＝60分 x 8 = 480ループになります。\n",
"\n",
"バリスタさんも休憩したりで、全員がいない場合もあると思いますがとりあえず3人固定で。\n",
"\n",
"待ち人数については毎単位時間ごとに並んでる人をモニタリングする処理を追加します。"
]
},
{
"cell_type": "code",
"execution_count": 7,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"# スタンドの前で待つ人のモニタリング\n",
"def monitor(resource, env):\n",
" while True:\n",
" # バリスタさんにコーヒーを淹れてもらっている人数\n",
" log_state(\"Brewing,%d,%d\" % (resource._env.now, resource.count))\n",
" # キャパ越えで待ってる人\n",
" log_state(\"Waiting,%d,%d\" % (resource._env.now, len(resource.queue)))\n",
" yield env.timeout(1)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"それぞれのPersonの動きには注目するのは不要になったので、結果を返す処理は除いておきます。"
]
},
{
"cell_type": "code",
"execution_count": 13,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"class Person(object):\n",
" def __init__(self, env, name, freq, coffee_garden):\n",
" self.env = env\n",
" self.name = name\n",
" self.freq = freq\n",
" self.coffee_garden = coffee_garden\n",
"\n",
" def needs_coffee(self):\n",
" return np.random.randint(SIM_SIZE) < self.freq\n",
"\n",
" def run(self):\n",
" while True:\n",
" if self.needs_coffee():\n",
" with self.coffee_garden.request() as req:\n",
" yield req\n",
" yield self.env.timeout(COFFEE_TIME) # コーヒーを淹れる時間\n",
" yield self.env.timeout(1)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"パラメタを再設定して改めてシミュレータを実行してみます。"
]
},
{
"cell_type": "code",
"execution_count": 18,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"SIM_SIZE = 60*8 # ループを回す回数\n",
"MAX_FREQ = 2 # Personがコーヒーを飲みに行きたくなる最大頻度\n",
"NUM_OF_PEOPLE = 300 # 人数\n",
"NUM_OF_BRISTA = 3 # バリスタさんの人数\n",
"COFFEE_TIME = 3 # コーヒーの淹れてくれる時間\n",
"\n",
"# 人ごとにコーヒーの飲みたくなる頻度を設定\n",
"def get_random_freq():\n",
" return np.random.randint(MAX_FREQ)\n",
"\n",
"# simpy のEnvironmentオブジェクトを生成\n",
"env = simpy.Environment()\n",
"\n",
"coffee_garden = simpy.Resource(env, NUM_OF_BRISTA)\n",
"\n",
"# Personオブジェクトを生成しrunメソッドをプロセスとして設定\n",
"for name in [\"person_%03d\" % i for i in np.arange(0,NUM_OF_PEOPLE,1)]:\n",
" env.process(Person(env, name, get_random_freq(), coffee_garden).run())\n",
"\n",
"# コーヒースタンドのモニタリング処理\n",
"env.process(monitor(coffee_garden, env))\n",
" \n",
"reset_state() \n",
"env.run(until=SIM_SIZE)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"やっとここからが本題です。(コードばかりで全く文章書いていないですが...!?\n",
"\n",
"1日のうちに待っている人数がどのくらいの頻度で起きていたのでしょう。\n",
"\n",
"ヒストグラムを見てみます。"
]
},
{
"cell_type": "code",
"execution_count": 20,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"text/plain": [
"<matplotlib.axes._subplots.AxesSubplot at 0x7f373d02de50>"
]
},
"execution_count": 20,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAhQAAAFoCAYAAAAPcmLCAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAAPYQAAD2EBqD+naQAAGfVJREFUeJzt3X+QZWV95/F3920Gx4HrTA3LDEbZRMGvgorUGJlJBLHG\nwC5Y/ohZFZNaiLsb+SFF4SqoqwViZZOQMEUCS0k0iZBas+UGjUaFYcGAhSK4CDoO8kU3mEFxhl8T\nekIGh77T+8e5Fy/t9Mzc83T36dv9flV1dd1znnPv9z7zzOlPP/c5p0cmJyeRJEkqMdp0AZIkafgZ\nKCRJUjEDhSRJKmagkCRJxQwUkiSpmIFCkiQVM1BIkqRiBgpJklTMQCFJkooZKCRJUrGxQRpHxIeA\ntwIvBXYC3wAuzMz7+9rcApzQd9gkcHVmnt3X5oXAJ4ATgR3AtcAHM3N3rXchSZIaNegMxfHAFcBx\nwBuAA4AbI2JpX5tJ4M+BVcBq4DDggt7OiBgFvkIVZtYCpwNnAJfUegeSJKlxA81QZOYp/Y8j4gzg\nYWANcFvfrn/NzEemeZqTqWY4Xp+ZjwKbIuKjwB9GxMWZOTFITZIkqXmlayiWU81IPD5l+29HxCMR\nsSki/vuUGYy1wKZumOjZCDwPOLqwHkmS1ICBZij6RcQIcDlwW2be27frfwL/BDwEvBK4FHgJ8Fvd\n/auBbVOeblvfvu/UrUmSJDWjdqAArgKOAn69f2Nmfqrv4eaI2ArcHBG/kpkP7OM5J/f3xScnJydH\nRkb2u1hJkvSMGf8BWitQRMSVwCnA8Zn50300v6P7/QjgAWAr8KtT2qzqfp86czGtkZERxsd30ul4\nYcj+arVGabeX2m8DsM/qsd8GZ5/VY78NrtdnM23gQNENE28GXpeZW/bjkGOpZh56weN24MMRcUjf\nOoqTgCeAe/dw/LQ6nd1MTDiABmW/Dc4+q8d+G5x9Vo/91rxB70NxFXAa8CbgyYjozSw8kZlPRcSL\ngHdRXRb6GHAMsAG4NTO/1217I1Vw+OuIuJDqstKPA1dm5tOlb0iSJM29Qa/yOBNoA7dQLbrsfb29\nu38X1f0pNgLfB/4Y+N9UAQSA7s2r3gh0qG6MdS3waeCiem9BkiQ1bdD7UOw1gGTmj6nufrmv53mQ\nKlRIkqQFwL/lIUmSihkoJElSMQOFJEkqZqCQJEnFDBSSJKmYgUKSJBUzUEiSpGIGCkmSVMxAIUmS\nihkoJElSMQOFJEkqZqCQJEnFDBSSJKmYgUKSJBUzUEiSpGIGCkmSVMxAIUmSihkoJElSMQOFJEkq\nZqCQJEnFDBSSJKmYgUKSJBUzUEiSpGIGCkmSVMxAIUmSihkoJElSMQOFJEkqZqCQJEnFDBSSJKmY\ngUKSJBUzUEiSpGIGCkmSVMxAIUmSihkoJElSMQOFJEkqZqCQJEnFDBSSJKmYgUKSJBUba7qAuq77\n/BfZseMpOrt3N13KQP7NypWsW7u26TIkSZpRQxso/uQvb+a5q45uuoyBPWfHbQYKSdKCM7SB4jkH\nreC5z1vVdBkDO+DpZU2XIEnSjHMNhSRJKmagkCRJxQwUkiSpmIFCkiQVM1BIkqRiBgpJklTMQCFJ\nkooZKCRJUjEDhSRJKmagkCRJxQa69XZEfAh4K/BSYCfwDeDCzLy/r82BwAbgHcCBwEbg7Mx8uK/N\nC4FPACcCO4BrgQ9m5nD9pS9JkgQMPkNxPHAFcBzwBuAA4MaIWNrX5nLgVOBtwAnA84HrejsjYhT4\nClWYWQucDpwBXFLrHUiSpMYNNEORmaf0P46IM4CHgTXAbRHRBt4NvDMzb+22+V3g+xHxmsy8EziZ\naobj9Zn5KLApIj4K/GFEXJyZE6VvSpIkza3SNRTLgUng8e7jNVQh5eZeg8xMYAuwrrtpLbCpGyZ6\nNgLPA4bv75FLkqT6f748IkaoPt64LTPv7W5eDezKzPEpzbd19/XabNvD/t6+79StaRiMjMDYWDNr\nYVut0Wd9177ZZ/XYb4Ozz+qx3wY3W31VO1AAVwFHAa/dj7YjVDMZ+7I/bYZaa6zFihXLGq2h3V66\n70Z6FvusHvttcPZZPfZb82oFioi4EjgFOD4zH+rbtRVYEhHtKbMUh/LzWYitwK9OecpV3e9TZy4W\nnM5Eh+3bn2zktVutUdrtpYyP76TT8YKa/WGf1WO/Dc4+q8d+G1yvz2bawIGiGybeDLwuM7dM2X0X\nMAGsBz7fbf8S4HCqS0wBbgc+HBGH9K2jOAl4AriXBW5yEiYmmh30nc7uxmsYNvZZPfbb4Oyzeuy3\n5g16H4qrgNOANwFPRkRvZuGJzHwqM8cj4i+ADRGxneoeE38GfD0zv9VteyNVcPjriLgQOAz4OHBl\nZj5d/pYkSdJcG3RlxplAG7gFeKjv6+19bc4HvgT8bV+7t/V2dm9e9UagQzVrcS3waeCiwcuXJEnz\nwaD3odhnAMnMnwHndr+ma/MgVaiQJEkLgNfZSJKkYgYKSZJUzEAhSZKKGSgkSVIxA4UkSSpmoJAk\nScUMFJIkqZiBQpIkFTNQSJKkYgYKSZJUzEAhSZKKGSgkSVIxA4UkSSpmoJAkScUMFJIkqZiBQpIk\nFTNQSJKkYgYKSZJUzEAhSZKKGSgkSVIxA4UkSSpmoJAkScUMFJIkqZiBQpIkFTNQSJKkYgYKSZJU\nzEAhSZKKGSgkSVIxA4UkSSo21nQB0mzbtWsXmzdvqnVsqzVKu72U8fGddDq7Z7iyfTv66FewZMmS\nOX9dSRqUgUIL3ubNm7hgw+c4eOXhTZcykB2PbeHS98Gxx65puhRJ2icDhRaFg1cezvLVRzZdhiQt\nWK6hkCRJxQwUkiSpmIFCkiQVM1BIkqRiBgpJklTMQCFJkooZKCRJUjEDhSRJKmagkCRJxQwUkiSp\nmIFCkiQVM1BIkqRiBgpJklTMQCFJkooZKCRJUjEDhSRJKmagkCRJxQwUkiSpmIFCkiQVGxv0gIg4\nHvgAsAY4DHhLZn6xb/9fAadPOeyGzDylr80K4ErgjcBu4DrgvMx8cuB3IEmSGldnhmIZcA9wDjA5\nTZvrgVXA6u7XaVP2fwZ4GbAeOBU4Abi6Ri2SJGkeGHiGIjNvAG4AiIiRaZr9LDMf2dOOiHgpcDKw\nJjPv7m47F/hyRLw/M7cOWpMkSWrWwIFiP50YEduA7cBXgY9k5uPdfeuA7b0w0XUT1WzHccAXZqkm\nSZI0S2YjUFxPtSbiAeDFwB8AX4mIdZk5SfURyMP9B2RmJyIe7+5b0EZGYGysmbWwrdbos74vFsP8\nflut0cbGS4nFOtZK2Gf12G+Dm62+mvFAkZmf7Xu4OSI2Af8POBH4h70cOsL0azIWjNZYixUrljVa\nQ7u9tNHXn2vD/H7b7aWNj5cSw9z3TbHP6rHfmjdbH3k8IzMfiIhHgSOoAsVW4ND+NhHRAlYA22a7\nnqZ1Jjps397MxSyt1ijt9lLGx3fS6exupIYmjI/vbLqE2sbHdzY2Xkos1rFWwj6rx34bXK/PZtqs\nB4qIeAGwEvhpd9PtwPKIOLZvHcV6qhmKO2a7nqZNTsLERLODvtPZ3XgNc2mYTzLD/m817PU3wT6r\nx35rXp37UCyjmm3oXeHxoog4Bni8+3UR1RqKrd12fwTcD2wEyMz7ImIj8MmIOAtYAlwB/I1XeEiS\nNJzqrMx4NXA3cBfVmofLgG8DHwM6wCuprtRI4JPAt4ATMvPpvud4F3Af1dUdXwK+Bryn3luQJElN\nq3MfilvZexD5d/vxHP8M/M6gry1JkuYnr7ORJEnFDBSSJKmYgUKSJBUzUEiSpGIGCkmSVMxAIUmS\nihkoJElSMQOFJEkqZqCQJEnFDBSSJKmYgUKSJBUzUEiSpGIGCkmSVMxAIUmSihkoJElSMQOFJEkq\nZqCQJEnFDBSSJKmYgUKSJBUzUEiSpGIGCkmSVMxAIUmSihkoJElSMQOFJEkqZqCQJEnFDBSSJKmY\ngUKSJBUzUEiSpGIGCkmSVMxAIUmSihkoJElSMQOFJEkqZqCQJEnFDBSSJKmYgUKSJBUzUEiSpGIG\nCkmSVMxAIUmSihkoJElSMQOFJEkqZqCQJEnFDBSSJKmYgUKSJBUzUEiSpGIGCkmSVMxAIUmSihko\nJElSMQOFJEkqZqCQJEnFDBSSJKmYgUKSJBUbG/SAiDge+ACwBjgMeEtmfnFKm0uA/wwsB74OnJWZ\nP+zbvwK4EngjsBu4DjgvM5+s+T4kSVKD6sxQLAPuAc4BJqfujIgLgfcC7wFeAzwJbIyIJX3NPgO8\nDFgPnAqcAFxdoxZJkjQPDDxDkZk3ADcARMTIHpqcB3w8M/++2+Y/AtuAtwCfjYiXAScDazLz7m6b\nc4EvR8T7M3NrrXciSZIaM6NrKCLiV4DVwM29bZk5DtwBrOtuWgts74WJrpuoZjuOm8l6JEnS3Bh4\nhmIfVlMFg21Ttm/r7uu1ebh/Z2Z2IuLxvjYL1sgIjI01sxa21Rp91vfFYpjfb6s12th4KbFYx1oJ\n+6we+21ws9VXMx0opjPCHtZb1Ggz9FpjLVasWNZoDe320kZff64N8/ttt5c2Pl5KDHPfN8U+q8d+\na95MB4qtVMFgFc+epTgUuLuvzaH9B0VEC1jBL85sLDidiQ7btzdzMUurNUq7vZTx8Z10OrsbqaEJ\n4+M7my6htvHxnY2NlxKLdayVsM/qsd8G1+uzmTajgSIzH4iIrVRXb3wXICLaVGsj/ke32e3A8og4\ntm8dxXqqIHLHTNYzH01OwsREs4O+09ndeA1zaZhPMsP+bzXs9TfBPqvHfmtenftQLAOOoAoAAC+K\niGOAxzPzQeBy4CMR8UPgR8DHgR8DXwDIzPsiYiPwyYg4C1gCXAH8jVd4SJI0nOqszHg11ccXd1Gt\nebgM+DbwMYDMvJQqIFxNNeOwFPj3mbmr7zneBdxHdXXHl4CvUd23QpIkDaE696G4lX0Ekcy8GLh4\nL/v/GfidQV9bkiTNT15nI0mSihkoJElSMQOFJEkqZqCQJEnFDBSSJKmYgUKSJBUzUEiSpGIGCkmS\nVMxAIUmSihkoJElSMQOFJEkqZqCQJEnFDBSSJKmYgUKSJBUzUEiSpGIGCkmSVMxAIUmSihkoJElS\nMQOFJEkqZqCQJEnFDBSSJKmYgUKSJBUzUEiSpGIGCkmSVMxAIUmSihkoJElSsbGmC5Ck+WLXrl1s\n3rypkddutUZpt5cyPr6TTmd3rec4+uhXsGTJkhmuTNo/BgpJ6tq8eRMXbPgcB688vOlSBrbjsS1c\n+j449tg1TZeiRcpAIUl9Dl55OMtXH9l0GdLQcQ2FJEkqZqCQJEnFDBSSJKmYgUKSJBUzUEiSpGIG\nCkmSVMxAIUmSihkoJElSMQOFJEkqZqCQJEnFDBSSJKmYgUKSJBUzUEiSpGIGCkmSVMxAIUmSihko\nJElSMQOFJEkqZqCQJEnFDBSSJKmYgUKSJBUzUEiSpGIGCkmSVGxspp8wIi4CLpqy+b7MPKq7/0Bg\nA/AO4EBgI3B2Zj4807VIkqS5MVszFN8DVgGru1+v7dt3OXAq8DbgBOD5wHWzVIckSZoDMz5D0TWR\nmY9M3RgRbeDdwDsz89butt8Fvh8Rr8nMO2epHkmSNItmK1AcGRE/AZ4Cbgc+lJkPAmu6r3lzr2Fm\nZkRsAdYBBgpJkobQbASKbwJnAAkcBlwMfC0iXk718ceuzByfcsy27r4Fb2QExsaaWQvbao0+6/ti\nMczvt9UabWy8lBjWsTZs9U41rOOlxLCOtSbNVl/NeKDIzI19D78XEXcC/wS8nWrGYk9GgMmZrmU+\nao21WLFiWaM1tNtLG339uTbM77fdXtr4eCkxbH0/bPVONezjpcSw/9stBLP1kcczMvOJiLgfOAK4\nCVgSEe0psxSHUs1SLHidiQ7btz/ZyGu3WqO020sZH99Jp7O7kRqaMD6+s+kSahsf39nYeCkxrGNt\nmMcKDO94KTGsY61JvT6babMeKCLiIODFwDXAXcAEsB74fHf/S4DDqdZaLHiTkzAx0eyg73R2N17D\nXBrmk8yw/1sNW/3DPFZg+Pp7Ji3m9z5fzMZ9KP4Y+Huqjzl+CfgYVYj4X5k5HhF/AWyIiO3ADuDP\ngK97hYckScNrNmYoXgB8BlgJPALcBqzNzMe6+88HOsDfUt3Y6gbgnFmoQ5IkzZHZWJR52j72/ww4\nt/slSZIWAK+zkSRJxQwUkiSpmIFCkiQVM1BIkqRiBgpJklTMQCFJkooZKCRJUjEDhSRJKmagkCRJ\nxQwUkiSpmIFCkiQVM1BIkqRiBgpJklTMQCFJkooZKCRJUjEDhSRJKmagkCRJxQwUkiSpmIFCkiQV\nM1BIkqRiBgpJklTMQCFJkooZKCRJUjEDhSRJKmagkCRJxQwUkiSpmIFCkiQVM1BIkqRiBgpJklRs\nrOkCJEkaRrt27WLz5k1NlzGwVmuU9etPmPHnNVBIklTD5s2buGDD5zh45eFNlzKQHY9t4R4DhSRJ\n88fBKw9n+eojmy5jXnANhSRJKmagkCRJxQwUkiSpmIFCkiQVM1BIkqRiBgpJklTMQCFJkooZKCRJ\nUjEDhSRJKmagkCRJxQwUkiSpmIFCkiQVM1BIkqRiBgpJklTMQCFJkooZKCRJUjEDhSRJKmagkCRJ\nxQwUkiSp2FiTLx4R5wDvB1YD3wHOzcxvNVmTJEkaXGMzFBHxDuAy4CLgWKpAsTEiDmmqJkmSVE+T\nH3mcD1ydmddm5n3AmcC/Au9usCZJklRDI4EiIg4A1gA397Zl5iRwE7CuiZokSVJ9Ta2hOARoAdum\nbN8GxNyXM3dGRuC73727kdceHR3hoIOew7/8y1Ps3j3ZSA1N+MEPkh2PbWm6jIHteGwLP/jBwbRa\nw7d2eljH2rCOFRju8VKiybE2rONltmoemZyc+//sEXEY8BNgXWbe0bf9UuC1mflrc16UJEmqrako\n+yjQAVZN2X4ovzhrIUmS5rlGAkVmPg3cBazvbYuIke7jbzRRkyRJqq/J+1BsAK6JiLuAO6mu+ngu\n8OkGa5IkSTU0soaiJyLOBi6g+ujjHqobW/3fxgqSJEm1NBooJEnSwrC4ri+SJEmzwkAhSZKKGSgk\nSVIxA4UkSSpmoJAkScUMFJIkqViTN7baq4g4B3g/sBr4DtU9Kr61l/b/AbgE+GXgfuCDmXn9HJQ6\nrwzSbxFxOvBXwCQw0t38VGY+dy5qnQ8i4njgA1R//fYw4C2Z+cV9HHMicBlwNLAF+P3MvGaWS503\nBu2ziHgd8A9TNk8Ch2Xmw7NW6DwSER8C3gq8FNhJdUfgCzPz/n0ct6jPa3X6bbGf1yLiTOAsqjED\nsBm4JDNv2MsxMzLO5uUMRUS8g+qEfRFwLNUPxo0Rccg07dcBnwE+CbwK+Dvg7yLiqLmpeH4YtN+6\nnqAKH72vfzvbdc4zy6huqnYO1QloryLil4EvATcDxwB/CnwqIn5jFmucbwbqs65J4Eh+Ps4WTZjo\nOh64AjgOeANwAHBjRCyd7gDPa0CNfutazOe1B4ELqQL/GuCrwBci4mV7ajyT42y+zlCcD1ydmdfC\nM4nrVODdwKV7aH8ecH1mbug+vigiTgLeC5w9B/XOF4P2G8BkZj4yR/XNO93UfgM88/dk9uUs4B8z\n84LeU0TEa6n6/v/MTpXzS40+63kkM8dnp6r5LTNP6X8cEWcAD1Od8G+b5rBFf16r2W+wiM9rmfnl\nKZs+EhFnAWuB7+/hkBkbZ/NuhiIiDqAaLDf3tmXmJHATsG6aw9Z19/fbuJf2C07NfgM4KCJ+FBFb\nImKx/fZTx1oW+ViraQS4JyIeiogbI+LXmi6oYcupZm0e30ubRX9e24P96TfwvAZARIxGxDup/k7W\n7dM0m7FxNu8CBXAI0OIX/4z5Nqqpqz1ZPWD7hahOvyXV7MWbgN+mGg/fiIhfmq0iF4Dpxlo7Ig5s\noJ5h8FPgPcDbgN+kmpK9JSJe1WhVDenO6lwO3JaZ9+6lqee1PgP026I/r0XEyyNiB/Az4CrgrZl5\n3zTNZ2yczdePPPZkhP3/vLZO+4Vq2n7IzG8C3+w9jojbqabEfo9qHYb2T2/a3/G2B90FdP2L6L4Z\nES+m+pjo9GaqatRVwFHAr9c4djGf1/ar3zyvAXAf1Rqv5VRB/tqIOGEvoWKqWuNsPs5QPAp0qP4C\nab9D+cUU1bN1wPYLUZ1+e5bMnADuBo6Y2dIWlOnG2nhm7mqgnmF1J4twnEXElcApwImZ+dN9NPe8\n1jVgvz3LYjyvZeZEZv5jZn47M/8b1QL986ZpPmPjbN4Fisx8GrgLWN/b1p3qWk91ydCe3N7fvus3\nmP4zowWnZr89S0SMAi+nmqLWnu1prJ3EIhprM+RVLLJx1v2h+Gbg9Zm5ZT8OWfTnNajVb1OP97xW\n/ayf7iPZGRtn8/Ujjw3ANRFxF9VvMudTLSr5NEBEXAv8ODM/3G3/p8CtEfE+4MvAaVQLFP/LHNfd\ntIH6LSI+SjU1+EOqqbELqC6v+tScV96QiFhG9ZtL72OLF0XEMcDjmflgRPwB8PzM7E3NfwJ4b0T8\nEfCXVP8Rf4vqt6dFYdA+i4jzgAeorod/DtX/y9dTnbQWhYi4iuq89CbgyYjo/Ub4RGY+1W1zDfAT\nz2s/V6ffFvt5LSJ+H7ieaq3SwVTrSF5H9YvPrP78nHczFACZ+Vngv1LdaONu4JXAyX2XAb2AvgUj\nmXk7VSf8HtX18b8JvHkfC3cWnEH7DVgB/DlwL9VAOghYN8DnbAvBq6n66i6qzwwvA74NfKy7fzXw\nwl7jzPwR1aW4b6Aaa+cD/ykzp66SXsgG6jNgSbfNd4FbgFcA6zPzlrkpd144E2hTvf+H+r7e3tfm\nhXhem2rgfsPz2irgWqp1FDdRhYOTMvOr3f2z9vNzZHJysa7vkSRJM2VezlBIkqThYqCQJEnFDBSS\nJKmYgUKSJBUzUEiSpGIGCkmSVMxAIUmSihkoJElSMQOFJEkqZqCQJEnFDBSSJKnY/wfUS2Zlgx8x\nyAAAAABJRU5ErkJggg==\n",
"text/plain": [
"<matplotlib.figure.Figure at 0x7f373d10cb10>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"dat = get_state_results()\n",
"res = dat.pivot(\"time\",\"name\",\"state\")\n",
"res.Brewing.hist()"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Brewing(コーヒーを淹れてもらっている状態)の人数が左から 0,1,2,3 となり\n",
"\n",
"一番右が、バリスタの方3人が同時にコーヒーを淹れている状態です。\n",
"\n",
"50に達していないので、1日(60x8=480)のうちの1割程度です。"
]
},
{
"cell_type": "code",
"execution_count": 21,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"text/plain": [
"<matplotlib.axes._subplots.AxesSubplot at 0x7f373af7f890>"
]
},
"execution_count": 21,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAhQAAAFoCAYAAAAPcmLCAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAAPYQAAD2EBqD+naQAAGYdJREFUeJzt3X2UZHV95/F3P8A4DDbMCWEGVogxylclLmFHyMxJQNyJ\nsAq7PrARTHYjEkx4kGOIDzysLgrrRpPjLBp04SRucPBo1gOLCirDghGPMCvsiDjh4avZPICSGYKM\nNE7QoWt6/7i3tShmmL71677V1bxf5/TpqXt/Vf2t71RXffpXv3trZHp6GkmSpBKjgy5AkiQNPwOF\nJEkqZqCQJEnFDBSSJKmYgUKSJBUzUEiSpGIGCkmSVMxAIUmSihkoJElSMQOFJEkqNt5kcERcDFzc\ns/n+zHxpvX8JsA44BVgCbADOzsyHu27jEOAK4DjgcWA9cEFm7uzzPkiSpAFrFChqfw2sBUbqy1Nd\n+y4DXg2cDEwCHwOuBY4BiIhR4EvAQ8Bq4GDgamAH8J4+apEkSQvASJMPB6tnKF6bmf9qF/smgH8C\nTs3M6+ptAdwHrM7MOyLi1cAXgIMy85F6zO8DHwR+PjOnem9XkiQtfP2soXhRRHw/Iv5fRHyqfgsD\nYBXVjMctMwMzM4EHgDX1ptXA5pkwUdsA7Acc3kctkiRpAWgaKP4PcBpwAnAm8IvA1yJiGbAS2JGZ\nkz3X2Vrvo/6+dRf76RojSZKGTKM1FJm5oeviX0fEHcA/AG8Efrybq40As3lfZfbvvQDT09PTIyMj\nex4oSZJ6zfkLaD+LMn8qMx+LiO8ALwRuBvaOiImeWYoD+dksxBbgqJ6bWVF/7525eEYjIyNMTj5B\np+PBIW0YGxtlYmKpPW+RPW+fPW+fPW/fTM/nWlGgiIh9gV8CPglsojriYy0wsyjzMOBQ4Pb6KhuB\niyLigK51FMcDjwH3Nv35nc5OpqZ8ALbJnrfPnrfPnrfPng+/pueh+BPgeqq3Of4F8H6qEPGXmTkZ\nEZ8A1kXENqpzTHwUuC0z76xv4iaq4HB1RJwPHARcClyemU/OxR2SJEnta7oo83nAp4H7gb+kOkx0\ndWb+oN5/HnADcA3wVarzTZw8c+X65FUnAR2qWYv1wFU8/WRZkiRpiDQ6D8UCM71t23anyFoyPj7K\n8uXLsOftsefts+fts+ftq3s+54sy/SwPSZJUzEAhSZKKGSgkSVIxA4UkSSpmoJAkScUMFJIkqZiB\nQpIkFTNQSJKkYgYKSZJUzEAhSZKKGSgkSVIxA4UkSSrW6OPLF5K/WP9pHn/8CTo7h+vDzQ5asYK1\nrzxu0GVIkjSnhjZQ/Pl132TZypcNuozGlmy8xUAhSVp0hjZQ7P2cfVmyz/6DLqOxvf55yaBLkCRp\nzrmGQpIkFTNQSJKkYgYKSZJUzEAhSZKKGSgkSVIxA4UkSSpmoJAkScUMFJIkqZiBQpIkFTNQSJKk\nYgYKSZJUzEAhSZKKGSgkSVIxA4UkSSpmoJAkScUMFJIkqZiBQpIkFTNQSJKkYgYKSZJUzEAhSZKK\nGSgkSVIxA4UkSSpmoJAkScUMFJIkqZiBQpIkFTNQSJKkYgYKSZJUzEAhSZKKGSgkSVIxA4UkSSpm\noJAkScUMFJIkqZiBQpIkFTNQSJKkYgYKSZJUzEAhSZKKjZdcOSIuBD4AXJaZf1hvWwKsA04BlgAb\ngLMz8+Gu6x0CXAEcBzwOrAcuyMydJfVIkqTB6HuGIiKOAt4K3N2z6zLgROBk4FjgYODaruuNAl+i\nCjOrgTcDpwGX9FuLJEkarL4CRUTsC3wKOAP4Ydf2CeB04LzMvDUz7wLeAvxaRBxdDzsBeDHw25m5\nOTM3AO8FzomIohkTSZI0GP3OUHwMuD4zv9Kz/eVUMw+3zGzIzAQeANbUm1YDmzPzka7rbQD2Aw7v\nsx5JkjRAjWcEIuJU4FeowkOvFcCOzJzs2b4VWFn/e2V9uXf/zL7et1AWlZERGB8fvrWwY2OjT/mu\n+WfP22fP22fP2zdfvW4UKCLieVRrJF6VmU82uOoIMD2LcbMZM9TGxsdYvnzZoMvo28TE0kGX8Kxj\nz9tnz9tnz4df0xmKVcDPA5siYqTeNgYcGxFvA/4NsCQiJnpmKQ7kZ7MQW4Cjem53Rf29d+Zi0elM\nddi2bfugy2hsbGyUiYmlTE4+QafjwThtsOfts+fts+ftm+n5XGsaKG4GXtaz7SrgPuCDwPeBJ4G1\nwHUAEXEYcChwez1+I3BRRBzQtY7ieOAx4N6G9Qyd6WmYmhreX5pOZ+dQ1z+M7Hn77Hn77PnwaxQo\nMnM7PS/6EbEd+EFm3ldf/gSwLiK2UZ1j4qPAbZl5Z32Vm+rbuDoizgcOAi4FLm/4NookSVog5mJl\nRu+6h/OAG4BrgK8CD1GdkwKA+uRVJwEdqlmL9VSzHBfPQS2SJGkAis/7kJn/uufyT4Bz66/dXedB\nqlAhSZIWAY/TkSRJxQwUkiSpmIFCkiQVM1BIkqRiBgpJklTMQCFJkooZKCRJUjEDhSRJKmagkCRJ\nxQwUkiSpmIFCkiQVM1BIkqRiBgpJklTMQCFJkooZKCRJUjEDhSRJKmagkCRJxQwUkiSpmIFCkiQV\nM1BIkqRiBgpJklTMQCFJkooZKCRJUjEDhSRJKmagkCRJxQwUkiSpmIFCkiQVM1BIkqRiBgpJklTM\nQCFJkooZKCRJUjEDhSRJKmagkCRJxQwUkiSpmIFCkiQVM1BIkqRiBgpJklTMQCFJkooZKCRJUjED\nhSRJKmagkCRJxQwUkiSpmIFCkiQVM1BIkqRiBgpJklTMQCFJkooZKCRJUjEDhSRJKmagkCRJxQwU\nkiSpmIFCkiQVG28yOCLOBM4Cnl9vuge4JDNvrPcvAdYBpwBLgA3A2Zn5cNdtHAJcARwHPA6sBy7I\nzJ0ld0SSJA1O0xmKB4HzgVX111eAz0fES+r9lwEnAicDxwIHA9fOXDkiRoEvUQWZ1cCbgdOAS/q+\nB5IkaeAazVBk5hd7Nr0nIs4CVkfE94HTgVMz81aAiHgLcF9EHJ2ZdwAnAC8GXpmZjwCbI+K9wAcj\n4n2ZOVV6hyRJUvv6XkMREaMRcSqwD7CRasZiHLhlZkxmJvAAsKbetBrYXIeJGRuA/YDD+61FkiQN\nVqMZCoCI+GWqAPEcqjUQr8/M+yPiSGBHZk72XGUrsLL+98r6cu/+mX13N61n2IyMwPj48K2FHRsb\nfcp3zT973j573j573r756nXjQAHcDxwB7E+1VmJ9RBz7DONHgOlZ3O5sxgy9sfExli9fNugy+jYx\nsXTQJTzr2PP22fP22fPh1zhQ1Osc/ra++M2IOBp4O/BZYO+ImOiZpTiQn81CbAGO6rnJFfX33pmL\nRakz1WHbtu2DLqOxsbFRJiaWMjn5BJ2OB+S0wZ63z563z563b6bnc62fGYpeo1SHiG4CpoC1wHUA\nEXEYcChwez12I3BRRBzQtY7ieOAx4N45qGXBm56Gqanh/aXpdHYOdf3DyJ63z563z54Pv6bnofgA\n8GWqw0efC/w28Arg+MycjIhPAOsiYhvV+oqPArdl5p31TdxEFRyujojzgYOAS4HLM/PJubhDkiSp\nfU1XZqygOhHV/cDNVEd2HJ+ZX6n3nwfcAFwDfBV4iGqdBQD1yatOAjpUsxbrgauAi/u9A5IkafCa\nnofijD3s/wlwbv21uzEPUoUKSZK0SHicjiRJKmagkCRJxQwUkiSpmIFCkiQVM1BIkqRiBgpJklTM\nQCFJkooZKCRJUjEDhSRJKmagkCRJxQwUkiSpmIFCkiQVM1BIkqRiBgpJklTMQCFJkooZKCRJUjED\nhSRJKmagkCRJxQwUkiSpmIFCkiQVM1BIkqRiBgpJklTMQCFJkooZKCRJUjEDhSRJKmagkCRJxQwU\nkiSpmIFCkiQVM1BIkqRiBgpJklTMQCFJkooZKCRJUjEDhSRJKmagkCRJxQwUkiSpmIFCkiQVM1BI\nkqRiBgpJklTMQCFJkooZKCRJUjEDhSRJKmagkCRJxQwUkiSpmIFCkiQVM1BIkqRiBgpJklTMQCFJ\nkooZKCRJUjEDhSRJKmagkCRJxQwUkiSp2HiTwRFxIfB64MXAE8DtwPmZ+Z2uMUuAdcApwBJgA3B2\nZj7cNeYQ4ArgOOBxYD1wQWbuLLkzkiRpMJrOUBwD/Cnwq8BvAHsBN0XE0q4xlwEnAicDxwIHA9fO\n7IyIUeBLVGFmNfBm4DTgkr7ugSRJGrhGMxSZ+ZruyxFxGvAwsAr4ekRMAKcDp2bmrfWYtwD3RcTR\nmXkHcALVDMcrM/MRYHNEvBf4YES8LzOnSu+UJElqV+kaiv2BaeDR+vIqqpByy8yAzEzgAWBNvWk1\nsLkOEzM2APsBhxfWI0mSBqDRDEW3iBihenvj65l5b715JbAjMyd7hm+t982M2bqL/TP77u63pmEw\nMgLj48O3FnZsbPQp3zX/7Hn77Hn77Hn75qvXfQcK4OPAS4Ffn8XYEaqZjD2ZzZihNjY+xvLlywZd\nRt8mJpbueZDmlD1vnz1vnz0ffn0Fioi4HHgNcExmPtS1awuwd0RM9MxSHMjPZiG2AEf13OSK+nvv\nzMWi05nqsG3b9kGX0djY2CgTE0uZnHyCTseDcdpgz9tnz9tnz9s30/O51jhQ1GHitcArMvOBnt2b\ngClgLXBdPf4w4FCqQ0wBNgIXRcQBXesojgceA+5lkZuehqmp4f2l6XR2DnX9w8iet8+et8+eD7+m\n56H4OPAm4N8B2yNiZmbhscz8cWZORsQngHURsY3qHBMfBW7LzDvrsTdRBYerI+J84CDgUuDyzHyy\n/C5JkqS2NV2ZcSYwAXwVeKjr641dY84DbgCu6Rp38szO+uRVJwEdqlmL9cBVwMXNy5ckSQtB0/NQ\n7DGAZOZPgHPrr92NeZAqVEiSpEXA43QkSVIxA4UkSSpmoJAkScUMFJIkqZiBQpIkFTNQSJKkYgYK\nSZJUzEAhSZKKGSgkSVIxA4UkSSpmoJAkScUMFJIkqZiBQpIkFTNQSJKkYgYKSZJUzEAhSZKKGSgk\nSVIxA4UkSSpmoJAkScUMFJIkqZiBQpIkFTNQSJKkYgYKSZJUzEAhSZKKGSgkSVIxA4UkSSpmoJAk\nScUMFJIkqZiBQpIkFTNQSJKkYgYKSZJUzEAhSZKKGSgkSVIxA4UkSSpmoJAkScUMFJIkqZiBQpIk\nFTNQSJKkYgYKSZJUzEAhSZKKGSgkSVIxA4UkSSpmoJAkScUMFJIkqZiBQpIkFTNQSJKkYgYKSZJU\nzEAhSZKKGSgkSVIxA4UkSSpmoJAkScXGm14hIo4B3gWsAg4CXpeZX+gZcwlwBrA/cBtwVmb+Tdf+\n5cDlwEnATuBa4O2Zub3P+yFJkgaonxmKZcC3gHOA6d6dEXE+8Dbg94Gjge3AhojYu2vYp4GXAGuB\nE4FjgSv7qEWSJC0AjWcoMvNG4EaAiBjZxZC3A5dm5vX1mN8BtgKvAz4bES8BTgBWZeZd9ZhzgS9G\nxDszc0tf90SSJA3MnK6hiIhfBFYCt8xsy8xJ4BvAmnrTamDbTJio3Uw12/Grc1mPJElqR+MZij1Y\nSRUMtvZs31rvmxnzcPfOzOxExKNdYxatkREYHx++tbBjY6NP+a75Z8/bZ8/bZ8/bN1+9nutAsTsj\n7GK9RR9jht7Y+BjLly8bdBl9m5hYOugSnnXsefvsefvs+fCb60CxhSoYrOCpsxQHAnd1jTmw+0oR\nMQYs5+kzG4tOZ6rDtm3DdzDL2NgoExNLmZx8gk5n56DLeVaw5+2z5+2z5+2b6flcm9NAkZl/FxFb\nqI7e+DZARExQrY34WD1sI7B/RBzZtY5iLVUQ+cZc1rMQTU/D1NTw/tJ0OjuHuv5hZM/bZ8/bZ8+H\nXz/noVgGvJAqAAC8ICKOAB7NzAeBy4D3RMTfAH8PXAp8D/g8QGbeHxEbgD+LiLOAvYE/BT7jER6S\nJA2nflZmvJzq7YtNVGsePgx8E3g/QGb+MVVAuJJqxmEp8OrM3NF1G78F3E91dMcNwNeozlshSZKG\nUD/nobiVPQSRzHwf8L5n2P9D4D80/dmSJGlh8jgdSZJUzEAhSZKKGSgkSVIxA4UkSSpmoJAkScUM\nFJIkqZiBQpIkFTNQSJKkYgYKSZJUzEAhSZKKGSgkSVIxA4UkSSpmoJAkScUMFJIkqZiBQpIkFTNQ\nSJKkYgYKSZJUzEAhSZKKGSgkSVIxA4UkSSpmoJAkScUMFJIkqZiBQpIkFTNQSJKkYgYKSZJUzEAh\nSZKKGSgkSVIxA4UkSSpmoJAkScUMFJIkqZiBQpIkFTNQSJKkYgYKSZJUzEAhSZKKGSgkSVIxA4Uk\nSSpmoJAkScUMFJIkqZiBQpIkFTNQSJKkYgYKSZJUzEAhSZKKGSgkSVIxA4UkSSpmoJAkScUMFJIk\nqZiBQpIkFTNQSJKkYgYKSZJUbHzQBUiSNIx27NjBPfdsHnQZjY2NjbJ27bFzfrsGCkmS+nDPPZt5\n97r/xXN/7tBBl9LI4z94gG8ttkAREecA7wRWAncD52bmnYOsSZKk2Xruzx3K/itfNOgyFoSBraGI\niFOADwMXA0dSBYoNEXHAoGqSJEn9GeSizPOAKzNzfWbeD5wJ/DNw+gBrkiRJfRhIoIiIvYBVwC0z\n2zJzGrgZWDOImiRJUv8GtYbiAGAM2NqzfSsQ7ZfTnpER+Pa37xp0GY2Njo6w777P4Uc/+jE7d04P\nupxnBXvePnvevmHu+Xe/mzz+gwcGXUZj81XzQjvKYwSY1SPqrz7zvpF5rmWevHbQBUiS5sDatcdy\nzjmDrmLhGNQaikeADrCiZ/uBPH3WQpIkLXADCRSZ+SSwCVg7sy0iRurLtw+iJkmS1L9BvuWxDvhk\nRGwC7qA66mMf4KoB1iRJkvowMj09uEUwEXE28G6qtz6+RXViq/87sIIkSVJfBhooJEnS4uCnjUqS\npGIGCkmSVMxAIUmSihkoJElSMQOFJEkqZqCQJEnFFtpnefxURJwDvBNYCdxNdY6KO59h/G8ClwDP\nB74DXJCZX26h1EWjSc8j4gzgd4BfrjdtAi56pv8jPV3Tx3nX9U4FPg18LjPfML9VLi59PLfsB/xX\n4PXAcuAfgD/IzBtbKHdR6KPnfwCcCRxK9VEN1wAXZuZPWih36EXEMcC7qD7V+yDgdZn5hT1c5zjg\nw8DhwAPABzLzk01+7oKcoYiIU6ju2MXAkVQPwA0RccBuxq+henL9M+BXgM8Bn4uIl7ZT8fBr2nPg\nFVQ9Pw5YDTwI3BQRB81/tYtDHz2fud4vAH8CfG3ei1xk+nhu2Qu4meqF7Q1Un4b8VuD7rRS8CPTR\n898C/qge/2LgdOAU4AOtFLw4LKM6WeQ5zOIDNyPi+cANwC3AEcBHgD+PiFc1+aELdYbiPODKzFwP\nEBFnAidSPbD+eBfj3w58OTPX1ZcvjojjgbcBZ7dQ72LQqOeZ+R+7L9czFidTfR7Lp+a92sWh6eOc\niBil6u9/Bo4F9mun1EWjac9/F9gfWJ2ZnXrb8H1e9WA17fka4OuZ+T/ryw9ExGeAo9sodjGoZ89u\nhJ9+TtaenAX8bWa+e+YmIuLXqf7v/vdsf+6Cm6Go/yJYRZWUAMjMaaq/Etbs5mpr6v3dNjzDeHXp\ns+e9lgF7AY/OeYGLUEHPLwYezsy/mN8KF58+e/5vgY3AxyNiS0RsjogL62CnPeiz57cDqyLiqPo2\nXgC8Bvji/Fb7rLaaOXgNXYi/FAcAYzz9Y8y3Ur3/tisrG47XU/XT814fopoG7n1Qatca9zwifg14\nC3DG/Ja2aPXzOH8B8JtUz5WvBi4F3gFcNE81LjaNe56Zn6EKzl+PiB3Ad4G/yswPzWehz3K7ew2d\niIgls72RhRgodmeEWbwXVDBeTzerHkbEBcAbqRb+7Jj3qha3XfY8IvYFrgbempnbWq9qcXumx/ko\n1RPr72XmXZn5War38s9qq7hFarc9rxcHXkS1KPNIqrUrJ0XEe1qrTlD9H0GD19GFuIbiEaBD9Qmk\n3Q7k6QlqxpaG4/VU/fQcgIh4J9Unxq7NzHvmp7xFqWnPfwn4BeD6rvdERwHqv+IiM/9unmpdLPp5\nnP8jsKOepp9xH7AyIsYzc2ruy1xU+un5JcD6rrf17qkD9ZXAf5mXKrW719DJJn8kLrgZisx8kuoQ\nxLUz2+on0LVU763tysbu8bVX1du1B332nIh4F/CfgBMy8675rnMx6aPn9wEvozqK6Yj66wvAV+p/\nPzjPJQ+9Ph/ntwEv7NkWwD8aJvasz57vA+zs2bYTGJnlAkM1t6vX0ONp+Bq6EGcoANYBn4yITcAd\nVCtN9wGuAoiI9cD3MnPmfcyPALdGxB9SLdx5E9VCoLe2XPcwa9TziHg31V8Sb6JahT2Tbn+Umdtb\nrn1Yzbrn9V8J93ZfOSJ+CExn5n2tVj3cmj63/HfgbRHxEeBy4DDgQuCyluseZk17fj1wXkR8C/gG\n8CKq55rP98wUaTciYhlVEJ4JYC+IiCOARzPzwYj4I+DgzHxzvf8Kqsf5h4D/QRUu/j3VYthZW3Az\nFAD1+5TvoHoQ3QX8S6q/gv+pHvI8uhb0ZOZGqhe236M69vYNwGsz8ylPwNq9pj2neg95L6oTzjzU\n9fWOtmoedn30XIX6eG75HtVfakdRnT/hMuC/US1C1iz08Ti/lOq8FZcC91CdX+jLVGsqNDsvp+r1\nJqo1EB8Gvgm8v96/EjhkZnBm/j3Voby/QfUaeh7wu5nZaJH9yPS0gU+SJJVZkDMUkiRpuBgoJElS\nMQOFJEkqZqCQJEnFDBSSJKmYgUKSJBUzUEiSpGIGCkmSVMxAIUmSihkoJElSMQOFJEkq9v8BRcAa\nt4CLqH0AAAAASUVORK5CYII=\n",
"text/plain": [
"<matplotlib.figure.Figure at 0x7f373d10cd10>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"res.Waiting.hist()"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"次はWaiting(前に先約がいて待ちの状態)のヒストグラムです。\n",
"\n",
"待ち人数が1人である割合はごくごくわずかですね。基本的に待たずに注文ができているようです。\n",
"\n",
"**バリスタさん3人がいれば、弊社の社員数であればほぼ待たずにコーヒーを注文できるようです。いいね！**"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"ただ、このままだと実行する度に結果にばらつき起きるので\n",
"\n",
"複数回実行した結果から、**各待ち人数の発生するばらつき具合**も確認しておきます。"
]
},
{
"cell_type": "code",
"execution_count": 36,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"from collections import Counter\n",
"\n",
"def get_hist_data(simulate_data):\n",
" # スタンド前での待ち人数を集計した結果を返す\n",
" d = simulate_data.pivot(\"time\",\"name\",\"state\")\n",
" return dict(Counter(d.Brewing)), dict(Counter(d.Waiting))\n",
"\n",
"def simulate():\n",
" i = 0\n",
" # 人ごとにコーヒーの飲みたくなる頻度を設定\n",
" def get_random_freq():\n",
" return np.random.randint(MAX_FREQ)\n",
"\n",
" SIM_SIZE = 60*8 # ループを回す回数\n",
" MAX_FREQ = 2 # Personがコーヒーを飲みに行きたくなる最大頻度\n",
" NUM_OF_PEOPLE = 300 # 人数\n",
" NUM_OF_BRISTA = 3 # バリスタさんの人数\n",
" COFFEE_TIME = 3 # コーヒーの淹れてくれる時間\n",
"\n",
" # simpy のEnvironmentオブジェクトを生成\n",
" env = simpy.Environment()\n",
" \n",
" coffee_garden = simpy.Resource(env, NUM_OF_BRISTA)\n",
" \n",
" # Personオブジェクトを生成しrunメソッドをプロセスとして設定\n",
" for name in [\"person_%03d\" % i for i in np.arange(0,NUM_OF_PEOPLE,1)]:\n",
" env.process(Person(env, name, get_random_freq(), coffee_garden).run())\n",
" \n",
" # コーヒースタンドのモニタリング処理\n",
" env.process(monitor(coffee_garden, env))\n",
" \n",
" reset_state() \n",
" env.run(until=SIM_SIZE)\n",
" return get_hist_data(get_state_results())"
]
},
{
"cell_type": "code",
"execution_count": 44,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"brew_stats = []\n",
"wait_stats = []\n",
"for i in np.arange(0,200,1):\n",
" b, w = simulate()\n",
" brew_stats.append(b)\n",
" wait_stats.append(w)"
]
},
{
"cell_type": "code",
"execution_count": 45,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"text/plain": [
"<matplotlib.axes._subplots.AxesSubplot at 0x7f371d073510>"
]
},
"execution_count": 45,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAgkAAAFoCAYAAADdImiaAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAAPYQAAD2EBqD+naQAAIABJREFUeJzt3X2QXNV9p/Gnu2f0wkgtRoPQSyAIjH0wJPFSOFjULviF\nLGywE3B5129JxcS1CRjHcVxlG5NNCgfXrrOkTLkwUcBsNphUki0n3hi/gVyG2Ak2Ml6MHSSs42Aj\nUBAahpmJZjQjwah79o/bg4bhAN09c6d77jyfKlXX9D3d9zRHor/zu+ecW5qenkaSJGmucqc7IEmS\nupMhQZIkJRkSJElSkiFBkiQlGRIkSVKSIUGSJCUZEiRJUpIhQZIkJRkSJElSkiFBkiQl9bTSOIRw\nJfA+YGvjqd3AdTHGuxrHVwI3AO8AVgI7gKtijE/Neo+TgZuBNwDjwO3Ax2KM9fl8EEmStLBarSTs\nA64Gzmn8uQe4I4Tw6sbxTwNvBt4GXABsAb4w8+IQQhn4Glk42Qa8B7gcuK7tTyBJknJRmu8NnkII\nw8CHycLAEPDOGOPfN44F4EfAthjj/SGEXwa+BGyOMT7daHMF8MfAhhjj0Xl1RpIkLZi25ySEEMoh\nhHcCxwH3kVUWeoC7Z9rEGCPwOHBe46ltwEMzAaFhB7AOOKvdvkiSpIXXckgIIfxcCGEceAbYDrw1\nxrgH2AQ8G2Mcm/OSwcYxGo+DiePMaiNJkrpASxMXG/YArwGOJ5t7cHsI4YKXaF8Cmrmm0dJ1j+np\n6elSqdTKSyRJUqapL9CWQ0Jj3sBPGz9+P4RwLvBB4PPAihBCdU414USOVQsOAL845y03Nh7nVhhe\n0sjIBOVysUNCpVKmWl3N2NhhajUXfyx1jmfxOKbFspzGs7+/r6l27VQS5iqTLXd8ADgKXAjMTFx8\nFfCzwHcabe8Dfj+EcMKseQkXAQeBh1s5ab0+Tb0+v0mXS0WtVufo0WL/hV1OHM/icUyLxfE8ptV9\nEv47cCfZUsi1wK8BrwcuijGOhRD+HLghhDBKtgfCjcC3Y4zfa7zF18nCwF+GEK4GNgOfAG6KMU4t\nxAeSJEkLo9WJixvJNj/aA3yDbEXDRTHGexrHPwR8Bfg74JvAfrJ5CwA0Nkx6C1Ajqy7cDtwGXNvu\nB5AkSfmY9z4JnTI0NL40O96Cnp4y/f19jI5OWPoqAMezeBzTYllO47lhw9qmJvV57wZJkpRkSJAk\nSUmGBEmSlGRIkCRJSYYESZKUZEiQJElJhgRJkpRkSJAkSUmGBEmSlGRIkCRJSYYESZKUZEiQJElJ\nhgRJkpRkSJAkSUmGBEmSlGRIkCRJSYYESZKUZEiQJElJhgRJkpRkSJAkSUmGBEmSlGRIkCRJSYYE\nSZKUZEiQJElJhgRJkpRkSJAkSUmGBEmSlGRIkCRJSYYESZKUZEiQJElJhgRJkpRkSJAkSUmGBEmS\nlGRIkCRJSYYESZKUZEiQJElJhgRJkpRkSJAkSUmGBEmSlGRIkCRJSYYESZKU1NNK4xDCNcBbgTOA\nw8B3gKtjjD+e1eabwAWzXjYN3BJjvGpWm5OBm4E3AOPA7cDHYoz1tj6FJElacK1WEs4HPgO8Dvgl\noBf4eghh9aw208BngY3AJmAz8NGZgyGEMvA1soCyDXgPcDlwXVufQJIk5aKlSkKM8ZLZP4cQLgee\nAs4B7p11aDLGOPQib3MxWSXijTHGp4GHQgh/CPxxCOHjMcajrfRJkiTlY75zEo4nqxyMzHn+10II\nQyGEh0II/2NOpWEb8FAjIMzYAawDzppnfyRJ0gJpqZIwWwihBHwauDfG+PCsQ38FPAbsB34BuB54\nFfCfG8c3AYNz3m5w1rEfNnP+crlEuVxqr/NLRKVSft6jljbHs3gc02JxPF+o7ZAAbAfOBP797Cdj\njP9r1o+7QwgHgLtDCKfGGB99mfecbvbk69f3USoVOyTMqFZXv3wjLRmOZ/E4psXieB7TVkgIIdwE\nXAKcH2N88mWaf7fxeDrwKHAA+MU5bTY2HudWGF7UyMjEsqgkVKurGRs7TK3mwo+lzvEsHse0WJbT\nePb39zXVruWQ0AgIlwKvjzE+3sRLziarEMyEifuA3w8hnDBrXsJFwEHg4cTrk+r1aer1pgsPS1qt\nVufo0WL/hV1OHM/icUyLxfE8ptV9ErYD7wJ+FZgIIcxUAA7GGI+EEE4D3k22xHEYeA1wA/CtGOOu\nRtuvk4WBvwwhXE22RPITwE0xxqn5fiBJkrQwWp2dcSVQBb5JNjFx5s/bG8efJds/YQfwI+BPgL8l\nCxUANDZMegtQI9uM6XbgNuDa9j6CJEnKQ6v7JLxkqIgx/ivZLoov9z77yIKCJEnqUq7zkCRJSYYE\nSZKUZEiQJElJhgRJkpRkSJAkSUmGBEmSlGRIkCRJSYYESZKUZEiQJElJhgRJkpRkSJAkSUmGBEmS\nlGRIkCRJSYYESZKUZEiQJElJhgRJkpRkSJAkSUmGBEmSlGRIkCRJST2d7sBysnfvo4yNHWy6faVS\nplpdzdjYYWq1etOvq1bXsXXrqe10UZKk5xgSFsnw8DDbtp1Nvd78l327KpUKu3Y9wsDAQO7nkiQV\nlyFhkQwMDLBz54MtVRImjxzlnx8d4RdOXc9xq5ofqmp1nQFBkjRvhoRF1OolgH8dOsRn77mfS950\nBidtWJNTryRJSnPioiRJSjIkSJKkJEOCJElKMiRIkqQkQ4IkSUoyJEiSpCRDQhfr7Slz8sa19PY4\nTJKkxec+CV3sZzasYftH38To6ARHj+a/U6MkSbP5K6okSUoyJEiSpCRDgiRJSjIkSJKkJEOCJElK\nMiRIkqQkQ4IkSUoyJHSxJ4YOcdX19/DE0KFOd0WStAwZErrY1NE6+wbHmXIjJUlSB7S042II4Rrg\nrcAZwGHgO8DVMcYfz2qzErgBeAewEtgBXBVjfGpWm5OBm4E3AOPA7cDHYox+G0qS1CVarSScD3wG\neB3wS0Av8PUQwupZbT4NvBl4G3ABsAX4wszBEEIZ+BpZQNkGvAe4HLiurU8gSZJy0VIlIcZ4yeyf\nQwiXA08B5wD3hhCqwHuBd8YYv9Vo85vAj0II58YY7wcuJqtEvDHG+DTwUAjhD4E/DiF8PMZ4dL4f\nSloMe/c+ytjYwabbVyplqtXVjI0dplZrrWhWra5j69ZTW+2iJM3LfG/wdDwwDYw0fj6n8Z53zzSI\nMcYQwuPAecD9ZNWDhxoBYcYO4M+As4AfzrNPUu6Gh4fZtu1s6vXFuUJWqVTYtesRBgYGFuV8kgTz\nCAkhhBLZpYV7Y4wPN57eBDwbYxyb03ywcWymzWDi+MwxQ4K63sDAADt3PthSJWFwdJKbv7ibKy87\ni439x7V0vmp1nQFB0qKbTyVhO3Am8B+aaFsiqzi8nGbaAFAulyiXS802X3AHRiY58ky+V0YOjB5+\n3mOeVq3sYdP61r64lrvTT39FS+2fHJnk537ay787++fZ7H/rQqhUys971NLmeL5QWyEhhHATcAlw\nfoxx/6xDB4AVIYTqnGrCiRyrFhwAfnHOW25sPM6tMLyo9ev7KJU6ExL2Dx3io9u/s2jn2/5/H1qU\n89zysQvZsmHNopxrOerv72P7R9/U6W4oB9Xq6pdvpCXD8Tym5ZDQCAiXAq+PMT4+5/ADwFHgQuDv\nG+1fBfws2XJJgPuA3w8hnDBrXsJFwEHgYZo0MjLRsUrC4NA4AFdeehZbTujL7TzlSpm+vpVMTDxD\nvcWJbq3Y//QEN9+xm8GhcVb3dK46U3Tzmbio7uSYFstyGs/+/ua+u1rdJ2E78C7gV4GJEMJMBeBg\njPFIjHEshPDnwA0hhFGyPRBuBL4dY/xeo+3XycLAX4YQrgY2A58AbooxTjXbl3p9mnq96asTC6pW\ny867sf84TsrxN++enjL9/X2Mjk5wNMcNlWY+T602net5lKnV6v53LhjHtFgcz2NavfByJVAFvgns\nn/Xn7bPafAj4CvB3s9q9beZgY8OktwA1surC7cBtwLWtd1+SJOWl1X0SXjZUxBifAT7Q+PNibfaR\nBQVJktSlnMIpSZKSDAmSJCnJkCBJkpIMCdIieWLoEFddfw9PDB3qdFckqSmGBGmRTB2ts29wnCmX\nVklaIgwJkiQpyZAgSZKSDAmSJClpPneBXNbWTY1T37+PI8/kdze/SqVM78hqDue8j3h9eJJ1U+O5\nvb8kaWkyJLRheuIQVzz2Raa2TzP3DldL1RWUmJ44F1jb6a5IkrqEIaENpb413HLKZfzum1/J5oF8\nKwmLcUeyJ4cnufGr/8Lv9S3f20QPjkxy5NlavucYnQSyu27O3FQrT6tWVNi4Pr+/n5KKz5DQpoO9\naylvOZlVm/L7zbunp8ya/j6mcr4LZHnlOAd7D+T2/t1ucGSSaz67c9HOd/MduxftXJ/87W0GBUlt\nMyRo2ZupIPzWr5zJloHm7rHejkqlxNrqasbHDudeSdg/PMGtX3449+qIpGIzJEgNWwb6OCXnylB/\nfx+jo73eq17SkuASSEmSlGRIkCRJSYYESZKUZEiQJElJTlyUKNYOmuAumpIWhiFBy14Rd9AEd9GU\nNH+GBC17RdtBE9xFU9LCMCRIFGsHTXAXTUkLw4mLkiQpyZAgSZKSDAmSJCnJkCBJkpKcuDgPjw3m\nuw69UikxPDGV+10D9w9P5PbeS0lRxhMcU0kLw5DQhlo9+x/8bXfu6XBPFtaqFZVOd6EjijqesHzH\nVNLCKE1P5/sbTV6GhsY72vGf7h+jUi7leo7B0UluvmM3V156Fhv781u/D9mXycb1+Z6jmxVtPMEx\nXQzHbv+d/7JW5W85jeeGDWub+h+elYQ2nbalmvs5KpVsDLec0MdJG9wUJ0+OpyS9kBMXJUlSkiFB\nkiQlGRIkSVKSIUGSJCUZEiRJUpIhoYv19pQ5eeNaenscpiJwPCUtNe6T0MWW05rd5cDxLB7HtFiW\n03g2u0+Cv9JIkqQkQ4IkSUoyJEiSpCS3ZZbatHfvo4yNHWy6faVSplpdzdjYYWq11q53Vqvr2Lr1\n1Fa7KEnz0nJICCGcD3wEOAfYDFwWY/zSrON/AbxnzsvuijFeMqtNP3AT8BagDnwB+GCM0fvbakkY\nHh5m27azqdcXZ3JTpVJh165HGBgYWJTzSRK0V0noA34A/G+yL/eUO4HLgZnZk8/MOf7XwEbgQmAF\ncBtwC/DrbfRHWnQDAwPs3PngolYSDAiSFlvLISHGeBdwF0AI4cWWUDwTYxxKHQghnAFcDJwTY3yw\n8dwHgK+GED4cYzzQap+kTmi1/L+clldJKoa8Ji6+IYQwGELYE0LYHkJYP+vYecDoTEBo+AYwDbwu\np/4sSU8MHeKq6+/hiaFDne6KJGkZymPi4p1klyEeBV4BfBL4WgjhvBjjNLAJeGr2C2KMtRDCSONY\nU8rlEuVyU3tBLFm1adg3OE5tOvstVEtbpVJ+3qOWPse0WBzPF1rwkBBj/PysH3eHEB4CfgK8AfiH\nl3hpiaya0JT16/solYodEoYnpgDo61tJf39fh3ujhVKtru50F7TAHNNicTyPyX0JZIzx0RDC08Dp\nZCHhAHDi7DYhhArQDww2+74jIxOFryRMTDzz3OPoqAs/lrr5TFxUd3JMi2U5jWezv3jmHhJCCCcB\nA8CTjafuA44PIZw9a17ChWSVhO82+771+jT1erFv31Bv/CWt1+pOdCuQmuNZOI5psTiex7SzT0If\nWVVg5tf400IIrwFGGn+uJZuTcKDR7n8CPwZ2AMQY94QQdgC3hhDeR7YE8jPA37iyQZKk7tHO7IzX\nAg8CD5DNIfgU8H3gj4Aa8AvAHUAEbgW+B1wQY5ya9R7vBvaQrWr4CvCPwBXtfQRJkpSHdvZJ+BYv\nHS7+UxPv8W+4cZIkSV3NdR5d7Pg1K3nXRYHj16zsdFckScuQIaGLHb92Je+++AyOX2tIkCQtPkOC\nJElKMiRIkqQkQ4IkSUoyJEiSpCRDgiRJSjIkSJKkJENCF3t2qsZjB8Z4dqrW6a5IkpYhQ0IX2//0\nBL/zJ//A/qe9A6QkafEZEiRJUpIhQZIkJRkSJElSkiFBkiQlGRIkSVKSIUGSJCUZEiRJUlJPpzug\nF7flhD5u+sgbWWWUkyR1gF8/XWxFb4VTNlVZ0VvpdFckScuQIUGSJCUZEiRJUpJzEhbR3r2PMjZ2\nsOn2lUqZanU1Y2OHqdXqTb+uWl3H1q2nttNFSZKeY0hYJMPDw2zbdjb1evNf9u2qVCrs2vUIAwMD\nuZ9LklRchoRFMjAwwM6dDy5aJcGAIEmaL0PCImr1EkBPT5n+/j5GRyc4ejT/CoQkSbM5cVGSJCUZ\nEiRJUpIhQZIkJRkSJElSkiFBkiQlGRIkSVKSIUGSJCUZEiRJUpIhQZIkJbnjoiQ1eBM26fkMCZKE\nN2GTUgwJkkR7N2EbHjvC337zJ/yXN7yCgeqqpl/nTdi0VBgSJKmhnZuwvf2t/9GbsKmwnLgoSZKS\nDAmSJCnJkCBJkpJanpMQQjgf+AhwDrAZuCzG+KU5ba4D/itwPPBt4H0xxkdmHe8HbgLeAtSBLwAf\njDFOtPk5JEnSAmunktAH/AB4PzA992AI4Wrgd4ArgHOBCWBHCGHFrGZ/DbwauBB4M3ABcEsbfZEk\nSTlpuZIQY7wLuAsghFBKNPkg8IkY45cbbX4DGAQuAz4fQng1cDFwTozxwUabDwBfDSF8OMZ4oK1P\nIkmSFtSCzkkIIZwKbALunnkuxjgGfBc4r/HUNmB0JiA0fIOsKvG6heyPJElq30Lvk7CJ7Mt+cM7z\ng41jM22emn0wxlgLIYzMavOyyuUS5XKqkFEclUr5eY9a2hzP4nlyZJL/dut3+Z23/Tyb1x/X6e5o\nnvw3+kKLtZlSicT8hTbaPGf9+j5KpWKHhBnV6upOd0ELyPEsjuGJKfYNjrNiZS/9/X2d7o4WiP9G\nj1nokHCA7Mt+I8+vJpwIPDirzYmzXxRCqAD9vLAC8aJGRiaWRSWhnZvHqDs5nsUzMfHMc4+joy7O\nWuqW07/RZkPtgoaEGOOjIYQDZKsW/hkghFAlm2vwp41m9wHHhxDOnjUv4UKycPHdZs9Vr09Trzdd\neFjSarW6W74WiONZHPXGF0ndMS0U/40e084+CX3A6WRf6gCnhRBeA4zEGPcBnwb+IITwCLAX+ATw\nr8AdADHGPSGEHcCtIYT3ASuAzwB/48oGSZK6RzuzM15LdungAbI5BJ8Cvg/8EUCM8XqyL/1byCoD\nq4FfjjE+O+s93g3sIVvV8BXgH8n2VZAkSV2inX0SvsXLhIsY48eBj7/E8X8Dfr3Vc0uSpMXjraIl\nFdLgyCRHnq3le47RSQD2Pz1BrZbvHKlVKypsdJmlFpkhQVLhDI5Mcs1ndy7a+W6+Y/einOeTv73N\noKBFZUiQVDgzFYTf+pUz2TKQ3/4FlUqJtdXVjI8dzrWSsH94glu//HDulRFpLkOCpMLaMtDHKZvW\n5vb+PT1l+vv7GB3tdcmcCsm9JyVJUpIhQZIkJRkSJElSkiFBkiQlGRIkSVKSqxskFdK6qXHq+/dx\n5Jn89hWoVMr0jqzmcM53DawPT7Juajy395dejCFBUuFMTxziise+yNT2aR7vdGcWyBWUmJ44F8hv\nSac0lyFBUuGU+tZwyymX8btvfiWbB/KtJFSrqxnLuZLw5PAkN371X/i9vjW5nUNKMSRIKqSDvWsp\nbzmZVTlvprSmv4+p0YlcN1MqrxznYO+B3N5fejFOXJQkSUmGBEmSlGRIkCRJSYYESZKUZEiQJElJ\nhgRJkpTkEkhJhfXYYL67FFYqJYYnphgfO0ytNp3befYPT+T23tJLMSRIKpxaPfvCvu3OPR3uycJa\ntaLS6S5omTEkSCqc07ZU+YPfeC2VcinX8wyOTnLzHbu58tKz2Nif386OkAWEjevzPYc0lyFBUiGd\ntqWa+zkqlSyEbDmhj5M2uGWyiseJi5IkKcmQIEmSkgwJkiQpyZAgSZKSDAmSJCnJkCBJbertKXPy\nxrX09vi/UhWTSyAlqU0/s2EN2z/6JkZHJzh6tN7p7kgLzvgrSZKSDAmSJCnJkCBJkpIMCZIkKcmQ\nIEmSkgwJkiQpyZAgSZKSDAmS1KYnhg5x1fX38MTQoU53RcqFIUGS2jR1tM6+wXGm3EhJBWVIkCRJ\nSQu+LXMI4Vrg2jlP74kxntk4vhK4AXgHsBLYAVwVY3xqofsiSa3Yu/dRxsYONt1+cHSSg4M/4UcP\nr2Kk/7imX1etrmPr1lPb6aJa0Op4ViplqtXVjI0dplZrvjpU5PHM694Nu4ALgVLj56Ozjn0a+GXg\nbcAY8KfAF4Dzc+qLJL2s4eFhtm07m3q99UsH//RXrbWvVCrs2vUIAwMDLZ9LzZnPeLaqyOOZV0g4\nGmMcmvtkCKEKvBd4Z4zxW43nfhP4UQjh3Bjj/Tn1R5Je0sDAADt3Prhov3kW8QulmzieCyOvkPDK\nEMITwBHgPuCaGOM+4JzGOe+eaRhjjCGEx4HzAEOCpI5ptWTc01Omv7/Pu0B2Kcdz/vKYuLgTuBy4\nGLgSOBX4xxBCH7AJeDbGODbnNYONY5IkqUsseCUhxrhj1o+7Qgj3A48BbyerLKSUgOlWzlMulyiX\nSy/fcAmrVMrPe9TS5ngWj2NaLI7nC+V1ueE5McaDIYQfA6cD3wBWhBCqc6oJJ5JVE5q2fn0fpVKx\nQ8KManV1p7ugBeR4Fo9jWiyO5zG5h4QQwhrgFcDngAfIVjpcCPx94/irgJ8lm7vQtJGRiWVRSWhn\nEo26k+NZPI5psSyn8ezv72uqXR77JPwJ8GWySww/A/wRWTD4PzHGsRDCnwM3hBBGgXHgRuDbra5s\nqNenqddbukKxZNVqdSfRFIjjWTyOabE4nsfkUUk4CfhrYAAYAu4FtsUYhxvHPwTUgL8j20zpLuD9\nOfRDkiTNQ2l6emn+Nj40NL40O94Cl+MUi+NZPI5psQyOTrL9i7u56rKz2NjCDppL0YYNa5u6Xu8U\nTkmS8IZdKYYESZKUZEiQJElJhgRJkpRkSJAkSUmGBEmSlJT7jouSJM3X4MgkR56t5XuO0UkA9j89\nQa2W7yr7VSsqbFzf/cssDQmSpK42ODLJNZ/duWjnu/mO3Ytynk/+9rauDwqGBElSV5upIPzWr5zJ\nloHm7jnQjkqlxNrqasbHDudaSdg/PMGtX34498rIQjAkSJKWhC0DfZyyaW1u739sB81ed9BscOKi\nJElKMiRIkqQkLzdIkrreuqlx6vv3ceSZ/Cb6VSplekdWc3jsMLVafpcb6sOTrJsaz+39F5IhQZLU\n1aYnDnHFY19kavs0j3e6MwvkCkpMT5wL5DfHYiEYEiRJXa3Ut4ZbTrmM333zK9k8kG8loVpdzVjO\nlYQnhye58av/wu/1rcntHAvFkCBJ6noHe9dS3nIyq3Je3bCmv4+p0YlcVzeUV45zsPdAbu+/kJy4\nKEmSkgwJkiQpycsNkqQl4bHBfFcEVColhiemFmXHxaXCkCBJ6mq1evaFfdudezrck4W1akWl0114\nWYYESVJXO21LlT/4jddSKZdyPc/g6CQ337GbKy89i439+d54ybtASpK0QE7bUs39HJVKFkK2nNDH\nSRu6f3niYnDioiRJSjIkSJKkJEOCJElKMiRIkqQkJy5Kkgpp795HGRs72HT7wdFJDg7+hB89vIqR\nFlY3VKvr2Lr11Ha62PUMCZKkwhkeHmbbtrOp11u/B8M//VVr7SuVCrt2PcLAwEDL5+p2hgRJUuEM\nDAywc+eDLVUS2r0LZLW6rpABAQwJkqSCavUSQE9Pmf7+PkZzvgvkUuLERUmSlGRIkCRJSYYESZKU\nZEiQJElJhgRJkpRkSJAkSUmGBEmSlGRIkCRJSYYESZKUZEiQJElJhgRJkpTU0Xs3hBDeD3wY2AT8\nEPhAjPF7neyTJEnKdKySEEJ4B/Ap4FrgbLKQsCOEcEKn+iRJko7p5OWGDwG3xBhvjzHuAa4EJoH3\ndrBPkiSpoSMhIYTQC5wD3D3zXIxxGvgGcF4n+iRJkp6vU3MSTgAqwOCc5weB0MwblMslyuXSQver\nq1Qq5ec9amlzPIvHMS0Wx/OFOjpxMaEETDfTcGBgTbETwizV6upOd0ELyPEsHse0WBzPYzoVl54G\nasDGOc+fyAurC5IkqQM6EhJijFPAA8CFM8+FEEqNn7/TiT5JkqTn6+TlhhuAz4UQHgDuJ1vtcBxw\nWwf7JEmSGkrT001NAchFCOEq4KNklx1+QLaZ0v/rWIckSdJzOhoSJElS93KdhyRJSjIkSJKkJEOC\nJElKMiRIkqQkQ4IkSUoyJEiSpKRuu3eDZgkhvB/4MLAJ+CHZPhLf62yv1I4QwvnAR8jufroZuCzG\n+KXO9krtCCFcA7wVOAM4TLZL7NUxxh93tGNqWwjhSuB9wNbGU7uB62KMd3WsU13CSkKXCiG8A/gU\ncC1wNllI2BFCOKGjHVO7+sg2DHs/Td7ETF3rfOAzwOuAXwJ6ga+HELwr0NK1D7iaLMSfA9wD3BFC\neHVHe9UF3EypS4UQdgLfjTF+sPFziewv8o0xxus72jnNSwihjpWEwmgE96eAC2KM93a6P1oYIYRh\n4MMxxr/odF86yUpCFwoh9JKl2btnnosxTgPfAM7rVL8kJR1PVh0a6XRHNH8hhHII4Z1k9xK6r9P9\n6TTnJHSnE4AKL7xt9iAQFr87klIaFb5PA/fGGB/udH/UvhDCz5GFglXAOPDWGOOezvaq86wkLC0l\nvJ4tdZPtwJnAOzvdEc3bHuA1ZHNN/gy4PYRwRme71HlWErrT00CN7O6Ys53IC6sLkjoghHATcAlw\nfozxyU73R/MTYzwK/LTx4/dDCOcCHyRb9bBsWUnoQjHGKeAB4MKZ5xplzQvJlltJ6qBGQLgUeGOM\n8fFO90e5KAMrO92JTrOS0L1uAD4XQngAuB/4ENlEmts62Sm1J4TQB5xOdskI4LQQwmuAkRjjvs71\nTK0KIWwH3gX8KjARQpip+B2MMR7pXM/UrhDCfwfuJFtBthb4NeD1wEWd7Fc3MCR0qRjj5xtLq64j\nu+zwA+CZ9IExAAAAiElEQVTiGONQZ3umNr0W+AeyOSXTZHtgAHwOeG+nOqW2XEk2ht+c8/xvArcv\nem+0EDaSjd1m4CDwz8BFMcZ7OtqrLuA+CZIkKck5CZIkKcmQIEmSkgwJkiQpyZAgSZKSDAmSJCnJ\nkCBJkpIMCZIkKcmQIEmSkgwJkiQpyZAgSZKSDAmSJCnp/wO4wM1FCKlVeQAAAABJRU5ErkJggg==\n",
"text/plain": [
"<matplotlib.figure.Figure at 0x7f371f80fe10>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"pd.DataFrame(brew_stats).fillna(0).boxplot()"
]
},
{
"cell_type": "code",
"execution_count": 46,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"text/plain": [
"<matplotlib.axes._subplots.AxesSubplot at 0x7f3737889a50>"
]
},
"execution_count": 46,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAgkAAAFoCAYAAADdImiaAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAAPYQAAD2EBqD+naQAAHIVJREFUeJzt3X2QZFWZ5/FvZjY0TUFWNw10w6CDbzwCjgyDYLEK6jLC\nKs6AS6zgOqFooIDIOsw6C7gYKIQzrBOyBCLRzKyx2LLqGhqK+NYsGOqKINgwLA30cXDkRZECupqu\ntvq1Kmv/uFlUUZx6yayXW2V+PxEV2XnvubdPPpFR+atzzz1ZGR4eRpIkabxq2R2QJEkLkyFBkiRl\nGRIkSVKWIUGSJGUZEiRJUpYhQZIkZRkSJElSliFBkiRlGRIkSVKWIUGSJGUtaaVxRFwOXD5u88aU\n0hHN/UuBq4EzgaXAOuDDKaWnx5zjJcAa4M3AVmAtcElKqdHma5AkSXOgpZDQtAE4Cag0nw+O2XcN\n8DbgDKAf+DzwDeAEgIioAt8DngR6gIOBLwG7gMva6IskSZojlVa+4Kk5knBaSunPMvvqwDPAWSml\nbza3BfAw0JNSujsi3gZ8GzgopfRss825wFXAASmlwfHnlSRJ5WhnTsKrIuK3EfGriLipefkA4BiK\nkYnbRxqmlBLwOHB8c1MP8MBIQGhaB3QDR7bRF0mSNEdaDQl3AWcDpwDnAS8DfhIRXcBqYFdKqX/c\nMb3NfTQfezP7GdNGkiQtAC3NSUgprRvzdENE3A08BrwL2DHBYRVgOtc0pn/dAxgeHh6uVCpTN5Qk\nSeNN6wO0nYmLz0spbYmIXwKvBG4D9oyI+rjRhAMZHS14Cjh23GlWNR/HjzBMqq9vgGq1vJBQq1Wp\n15fR37+doaHOvjHDWhSswyhrUbAOBeswaqHUYsWKrmm1m1FIiIh9gFcAXwTWU9zpcBIwMnHxMOCl\nwM+ah9wJfDwi9h8zL+FkYAvwUCv/d6MxTKPR0uDDnBgaajA42Nlv+hHWomAdRlmLgnUoWIdRi6UW\nra6T8A/ALRSXGP4I+BRFMPhqSqk/Ir4AXB0RmynWQLgWuCOldE/zFLdShIEvRcTFwEHAlcB1KaXd\ns/GCJEnS7Gh14uIhwJeBjcBXKW557EkpbWruvwj4DvB14EcU6yGcMXJwc8GkdwBDFKMLa4EbefEC\nTZIkqWStTlx89xT7dwIXNn8mavMERVCQJEkLmN/dIEmSsgwJkiQpy5AgSZKyDAmSJCnLkCBJkrIM\nCZIkKcuQIEmSsma0LPMfukcf/TX9/Vuy+6a7/na93s2hh75srrooSdKcMSRMYNOmTfT0HE2jMbO1\ntWu1Ghs2PMLKlStnqWeSJM0PQ8IEVq5cyV133TfhSELv5m2s+daDnHf6kaxasfeE56nXuw0IkqRF\nyZAwickuE/zmmd/TfecODj/itRxywD7z2CtJkuZHR4eE3r5t7Ng11N6xm7cB8OSzAwwNtf+V1Xvt\nWWPVfhOPREiSVJaODQm9fdu49B/vmvF51tz84IzP8fcf6jEoSJIWnI4NCSMjCB/8iyM4eGVXy8fX\nahX2rS9ja//2tkcSntw0wD/d8lDboxmSJM2ljg0JAN27t7J6Zx+rdu5o+dharUp95zL23jn5LZCT\naezcRvfurW0dK0nSXOvYkDA88HvOfexb7L5+mMdL7Me5VBgeOA7Yt8ReSJL0Yh0bEipd+3DDH5/O\nf+g5mIPaudxQrbB311K2DexkqNHe5YZnntvOV+74DX/d5d0RkqSFp2NDwlBjmC177Mv/WL8VKHHI\nf4992WvPWnn/vyRJE+jYkPDyg+tc9t7XUatW2jq+d/M21tz8IOedNvliSlPxFkhJ0kLVsSEBiqDQ\nrlqtCBcH79/lYkqSpD9IfgukJEnKMiRIkqQsQ4IkScrq6DkJU3n00V9P+i2QW3p/xcMP7UXfFN8C\nOdkXRUmStFAZEiawadMmenqOptGYfDXF//u/Jj9PrVZjw4ZH/LpoSdKiY0iYwMqVK7nrrvsmHEmo\n1arU68vo7598WeZ6vduAIElalAwJk5jsMsGSJVVWrOhi8+YBBgfb++4GSZIWMicuSpKkLEOCJEnK\nMiRIkqQsQ4IkScoyJEiSpCxDgiRJyjIkSJKkLEOCJEnKMiRIkqQsQ4IkScoyJEiSpCxDgiRJyjIk\nSJKkLEOCJEnKMiRIkqQsQ4IkScoyJEiSpCxDgiRJyjIkSJKkLEOCJEnKMiRIkqQsQ4IkScoyJEiS\npCxDgiRJyloyk4Mj4lLg08A1KaW/aW5bClwNnAksBdYBH04pPT3muJcAa4A3A1uBtcAlKaXGTPoj\nSZJmT9sjCRFxLPBB4P5xu64BTgXOAE4EDga+Mea4KvA9ioDSA7wPOBu4ot2+SJKk2ddWSIiIfYCb\ngHOA58ZsrwMfAC5KKf04pXQf8H7gDRFxXLPZKcCrgfeklB5IKa0DPgFcEBEzGtmQJEmzp92RhM8D\nt6SUfjhu++soRghuH9mQUkrA48DxzU09wAMppWfHHLcO6AaObLM/kiRplrX8l3tEnAX8KUUgGG8V\nsCul1D9uey+wuvnv1c3n4/eP7Bt/+SKrWq1QrVam1ee5UKtVX/DYyaxFwTqMshYF61CwDqMWWy1a\nCgkRcQjFnIO3ppR2t3BoBRieRrvptAFgv/26qFTKCwkj6vVlZXdhwbAWBeswyloUrEPBOoxaLLVo\ndSThGOAAYH1EjHxC14ATI+IjwL8DlkZEfdxowoGMjhY8BRw77ryrmo/jRxgm1Nc3UPpIQr2+jP7+\n7QwNdfZNGdaiYB1GWYuCdShYh1ELpRYrVnRNq12rIeE24E/GbbsReBi4CvgtsBs4CfgmQEQcBrwU\n+Fmz/Z3AxyNi/zHzEk4GtgAPTbcjjcYwjca0Bx7mzNBQg8HBzn7Tj7AWBeswyloUrEPBOoxaLLVo\nKSSklAYY90EeEQPAppTSw83nXwCujojNFGsgXAvckVK6p3nIrc1zfCkiLgYOAq4ErmvxEoYkSZpD\nszFzYvyf8xcB3wG+DvwIeJJizQQAmgsmvQMYohhdWEsxGnH5LPRFkiTNkhmvS5BS+rfjnu8ELmz+\nTHTMExRBQZIkLVCL4x4MSZI07wwJkiQpy5AgSZKyDAmSJCnLkCBJkrIMCZIkKcuQIEmSsgwJkiQp\ny5AgSZKyDAmSJCnLkCBJkrIMCZIkKcuQIEmSsgwJkiQpy5AgSZKyDAmSJCnLkCBJkrIMCZIkKcuQ\nIEmSsgwJkiQpy5AgSZKyDAmSJCnLkCBJkrIMCZIkKcuQIEmSsgwJkiQpy5AgSZKyDAmSJCnLkCBJ\nkrIMCZIkKcuQIEmSsgwJkiQpy5AgSZKyDAmSJCnLkCBJkrIMCZIkKcuQIEmSsgwJkiQpy5AgSZKy\nDAmSJCnLkCBJkrIMCZIkKcuQIEmSsgwJkiQpy5AgSZKyDAmSJCnLkCBJkrIMCZIkKcuQIEmSsgwJ\nkiQpy5AgSZKylrTSOCLOA84HDm1uehC4IqX0g+b+pcDVwJnAUmAd8OGU0tNjzvESYA3wZmArsBa4\nJKXUmMkLkSRJs6vVkYQngIuBY5o/PwRujojDm/uvAU4FzgBOBA4GvjFycERUge9RhJMe4H3A2cAV\nbb8CSZI0J1oaSUgpfXfcpssi4nygJyJ+C3wAOCul9GOAiHg/8HBEHJdSuhs4BXg18JaU0rPAAxHx\nCeCqiPhkSmlwpi9IkiTNjrbnJERENSLOAvYG7qQYWVgC3D7SJqWUgMeB45ubeoAHmgFhxDqgGziy\n3b5IkqTZ19JIAkBEvIYiFOxFMafgnSmljRFxNLArpdQ/7pBeYHXz36ubz8fvH9l3/3T7Ua1WqFYr\nrXZ/1tRq1Rc8djJrUbAOo6xFwToUrMOoxVaLlkMCsBE4ClhOMfdgbUScOEn7CjA8jfNOp83z9tuv\ni0qlvJAwol5fVnYXFgxrUbAOo6xFwToUrMOoxVKLlkNCc97Avzaf3hsRxwEfBb4G7BkR9XGjCQcy\nOlrwFHDsuFOuaj6OH2GYVF/fQOkjCfX6Mvr7tzM01Nk3ZliLgnUYZS0K1qFgHUYtlFqsWNE1rXbt\njCSMV6W43XE9MAicBHwTICIOA14K/KzZ9k7g4xGx/5h5CScDW4CHWvlPG41hGo2WBh/mxNBQg8HB\nzn7Tj7AWBeswyloUrEPBOoxaLLVodZ2ETwPfp7gVcl/gPcCbgJNTSv0R8QXg6ojYTDFf4VrgjpTS\nPc1T3EoRBr4UERcDBwFXAtellHbPxguSJEmzo9WZE6soFj/aCNxGcUfDySmlHzb3XwR8B/g68CPg\nSYp5CwA0F0x6BzBEMbqwFrgRuLzdFyBJkuZGq+sknDPF/p3Ahc2fido8QREUJEnSArY47sGQJEnz\nzpAgSZKyDAmSJCnLkCBJkrIMCZIkKcuQIEmSsgwJkiQpy5AgSZKyDAmSJCnLkCBJkrIMCZIkKcuQ\nIEmSsgwJkiQpy5AgSZKyDAmSJCnLkCBJkrIMCZIkKcuQIEmSsgwJkiQpy5AgSZKyDAmSJCnLkCBJ\nkrIMCZIkKcuQIEmSsgwJkiQpy5AgSZKyDAmSJCnLkCBJkrIMCZIkKcuQIEmSsgwJkiQpy5AgSZKy\nDAmSJCnLkCBJkrIMCZIkKcuQIEmSsgwJkiQpy5AgSZKyDAmSJCnLkCBJkrIMCZIkKcuQIEmSsgwJ\nkiQpy5AgSZKyDAmSJCnLkCBJkrIMCZIkKcuQIEmSsgwJkiQpy5AgSZKyDAmSJClrSSuNI+JS4J3A\nq4HtwM+Ai1NKvxzTZilwNXAmsBRYB3w4pfT0mDYvAdYAbwa2AmuBS1JKjZm8GEmSNHtaHUk4Afgc\n8Hrgz4E9gFsjYtmYNtcApwJnACcCBwPfGNkZEVXgexQBpQd4H3A2cEVbr0CSJM2JlkYSUkpvH/s8\nIs4GngaOAX4aEXXgA8BZKaUfN9u8H3g4Io5LKd0NnEIxEvGWlNKzwAMR8Qngqoj4ZEppcKYvSpIk\nzdxM5yQsB4aBvubzYyiCx+0jDVJKCXgcOL65qQd4oBkQRqwDuoEjZ9gfSZI0S1oaSRgrIioUlxZ+\nmlJ6qLl5NbArpdQ/rnlvc99Im97M/pF990/n/69WK1SrlZb7PVtqteoLHjuZtShYh1HWomAdCtZh\n1GKrRdshAbgeOAJ44zTaVihGHKYynTYA7LdfF5VKeSFhRL2+bOpGHcJaFKzDKGtRsA4F6zBqsdSi\nrZAQEdcBbwdOSCk9OWbXU8CeEVEfN5pwIKOjBU8Bx4475arm4/gRhgn19Q2UPpJQry+jv387Q0Od\nfVOGtShYh1HWomAdCtZh1EKpxYoVXdNq13JIaAaE04A3pZQeH7d7PTAInAR8s9n+MOClFLdLAtwJ\nfDwi9h8zL+FkYAvwENPUaAzTaEx74GHODA01GBzs7Df9CGtRsA6jrEXBOhSsw6jFUotW10m4Hng3\n8JfAQESMjABsSSntSCn1R8QXgKsjYjPFGgjXAneklO5ptr2VIgx8KSIuBg4CrgSuSyntnvlLkiRJ\ns6HVmRPnAXXgR8CTY37eNabNRcB3gK+PaXfGyM7mgknvAIYoRhfWAjcCl7fefUmSNFdaXSdhylCR\nUtoJXNj8majNExRBQZIkLVCL4x4MSZI07wwJkiQpy5AgSZKyDAmSJCnLkCBJkrIMCZIkKcuQIEmS\nsgwJkiQpy5AgSZKyDAmSJCnLkCBJkrIMCZIkKcuQIEmSsgwJkiQpy5AgSZKyDAmSJCnLkCBJkrIM\nCZIkKcuQIEmSsgwJkiQpy5AgSZKyDAmSJCnLkCBJkrIMCZIkKcuQIEmSsgwJkiQpy5AgSZKyDAmS\nJCnLkCBJkrIMCZIkKcuQIEmSsgwJkiQpy5AgSZKyDAmSJCnLkCBJkrIMCZIkKcuQIEmSsgwJkiQp\ny5AgSZKyDAmSJCnLkCBJkrIMCZIkKcuQIEmSsgwJkiQpy5AgSZKyDAmSJCnLkCBJkrIMCZIkKcuQ\nIEmSsgwJkiQpy5AgSZKylrR6QEScAPwtcAxwEHB6Sunb49pcAZwDLAfuAM5PKT0yZv8K4DrgHUAD\n+Abw0ZTSQJuvQ5IkzbJ2RhK6gH8GLgCGx++MiIuBjwDnAscBA8C6iNhzTLMvA4cDJwGnAicCN7TR\nF0mSNEdaHklIKf0A+AFARFQyTT4KXJlSuqXZ5r1AL3A68LWIOBw4BTgmpXRfs82FwHcj4mMppafa\neiWSJGlWzeqchIh4GbAauH1kW0qpH/g5cHxzUw+weSQgNN1GMSrx+tnsjyRJal/LIwlTWE3xYd87\nbntvc99Im6fH7kwpDUVE35g2U6pWK1SruYGM+VGrVV/w2MmsRcE6jLIWBetQsA6jFlstZjskTKRC\nZv5CG22et99+XVQq5YWEEfX6srK7sGBYi4J1GGUtCtahYB1GLZZazHZIeIriw34VLxxNOBC4b0yb\nA8ceFBE1YAUvHoGYUF/fQOkjCfX6Mvr7tzM01CitHwuBtShYh1HWomAdCtZh1EKpxYoVXdNqN6sh\nIaX064h4iuKuhf8HEBF1irkGn282uxNYHhFHj5mXcBJFuPj5dP+vRmOYRmPaAw9zZmioweBgZ7/p\nR1iLgnUYZS0K1qFgHUYtllq0s05CF/BKig91gJdHxFFAX0rpCeAa4LKIeAR4FLgS+A1wM0BKaWNE\nrAP+KSLOB/YEPgd8xTsbJElaONqZOfE6iksH6ynmEHwWuBf4FEBK6TMUH/o3UIwMLAPellLaNeYc\n/xHYSHFXw3eAn1CsqyBJkhaIdtZJ+DFThIuU0ieBT06y/zngr1r9vyVJ0vxZHPdgSJKkeWdIkCRJ\nWYYESZKUZUiQJElZhgRJkpRlSJAkSVmGBEmSlGVIkCRJWYYESZKUZUiQJElZhgRJkpRlSJAkSVmG\nBEmSlGVIkCRJWYYESZKUZUiQJElZhgRJkpRlSJAkSVmGBEmSlGVIkCRJWYYESZKUZUiQJElZhgRJ\nkpRlSJAkSVmGBEmSlGVIkCRJWYYESZKUZUiQJElZhgRJkpRlSJAkSVmGBEmSlGVIkCRJWYYESZKU\nZUiQJElZhgRJkpRlSJAkSVmGBEmSlGVIkCRJWYYESZKUZUiQJElZS8rugBa+Rx/9Nf39WybcX6tV\nqdeX0d+/naGhRrZNvd7NoYe+bK66KEmaA4YETWrTpk309BxNo5H/8J+uWq3Ghg2PsHLlylnqmSRp\nrhkSNKmVK1dy1133TTqS0Lt5G2u+9SDnnX4kq1bsnW1Tr3cbECRpkTEkaEpTXSb4zTO/p/vOHRx+\nxGs55IB95qlXkqS55sRFSZKUZUiQJElZhgRJkpTlnATR27eNHbuG2j9+8zYAnnx2gKGh4bbOsdee\nNVbtl5/0KEkqhyGhw/X2bePSf7xrVs615uYHZ3T833+ox6AgSQuIIaHD7dg1RPfurbz7DYdwwPJl\nbZ2jVq2wd9dStg3sZKjR+kjCM89t5yt3/GZGoxmSpNlnSOhwwwO/59zHvkX1sWF2t3mO3cCOGfRh\nOXAuFYYHjgP2ncGZ5t5kq0+68qSkPzSlhoSIuAD4GLAauB+4MKV0T5l96jSVrn244Y9P5z+d+ioO\nWtneUP90Phwn87tN27j2u//CX3ct7DUWZmP1SVeelLSYlBYSIuJM4LPAh4C7gYuAdRFxWErp2bL6\n1Ym27LEvT9SWU13a3l/xtVqFbUuXsXXp9rYmLj61dC+27LGwRxBg6tUnXXlS0h+aMkcSLgJuSCmt\nBYiI84BTgQ8AnymxXx1lZA7Bjd/fWHJPijscyjbVnR6Vvfane6/9s/t27LGN7lU76D7wFXRPEBIA\nHntq66R98E4PSQtFKSEhIvYAjgH+bmRbSmk4Im4Dji+jT53q5QfXuey9r6NWrbR9jt7N21hz84Oc\nd9rEf0FPZSF8MPb2beOqz/8f9hra1fY5VgHf/N8/mVE/dtT25JIL3lp6PSSprJGE/YEa0Dtuey8Q\n0zlBtVqhOoMPtpmq1aoveFzMDnvp8kn3P/ror9myZeIveNrat40tvb9i67N709XIf7B1dy/8CXu7\n+vuLSZy0t9bDbGlQYVf/61lyYHlzNB6895dsfua5Cfc/kjaw5blNE+6vVissWVJjcHCIxgR3vHQv\nX8kr4zUTnmPFAcs58s8Om36n58BUdYDJazEbdYDFUYuZ1gF8T4xYSHWoDA/P/y/EiDgI+C1wfErp\n52O2fwZ4Y0rp38x7pyRJ0guU9Wfws8AQxejsWAfy4tEFSZJUglJCQkppN7AeOGlkW0RUms9/Vkaf\nJEnSC5V5d8PVwBcjYj2jt0DuDdxYYp8kSVJTKXMSRkTEh4H/QnHZ4Z8pFlP6RWkdkiRJzys1JEiS\npIVr8d+/J0mS5oQhQZIkZRkSJElSliFBkiRlGRIkSVKWIUGSJGWVuZjSohYRFwAfA1YD91Os8XBP\nub2aXxFxAvC3FN/oeRBwekrp2+X2av5FxKXAO4FXA9spVg29OKX0y1I7Ns+aX/d+PnBoc9ODwBUp\npR+U1qkFoPn++DRwTUrpb8ruz3yKiMuBy8dt3phSOqKM/pQpIg4G/hvwNoqFA/8FeH9K6d5SOzYF\nRxLaEBFnAp+lePMfTRES1kXE/qV2bP51USyCdQGU/NWJ5ToB+BzweuDPgT2AWyNiWam9mn9PABdT\nhMZjgB8CN0fE4aX2qkQRcSzwQYrfEZ1qA8WCeaubP28stzvzLyKWA3cAO4FTgMOB/wxsLrNf0+FI\nQnsuAm5IKa2F5/+COhX4APCZMjs2n5p/If4Anv/ujY6UUnr72OcRcTbwNMUH5U/L6FMZUkrfHbfp\nsog4H+gBHi6hS6WKiH2Am4BzgE+U3J0yDaaUnim7EyW7BHg8pXTOmG2PldWZVjiS0KKI2IPil//t\nI9tSSsPAbcDxZfVLC8pyipGVvrI7UpaIqEbEWRTDqneW3Z+SfB64JaX0w7I7UrJXRcRvI+JXEXFT\nRLyk7A6V4C+AX0TE1yKiNyLujYhzpjxqATAktG5/oMaLv9K6l2IoTR2sOaJyDfDTlNJDZfdnvkXE\nayJiK8Ww6vXAO1NKG0vu1rxrBqQ/BS4tuy8luws4m2KI/TzgZcBPIqKrzE6V4OUU83UScDKwBrg2\nIv6q1F5Ng5cbZk+Fzr4ur8L1wBHAG8ruSEk2AkdRjKacAayNiBM7KShExCEUQfGtKaXdZfenTCml\ndWOeboiIuymG2d8F/M9yelWKKnB3SmnkstP9EXEkRXC4qbxuTc2Q0LpngSGKiThjHciLRxfUQSLi\nOuDtwAkppd+V3Z8ypJQGgX9tPr03Io4DPkrxy7BTHAMcAKwfM1enBpwYER8BljYvUXaclNKWiPgl\n8Mqy+zLPfseL5+U8DPz7EvrSEi83tKj5l8F64KSRbc1fBCdR3PqmDtQMCKcBb0kpPV52fxaQKrC0\n7E7Ms9uAP6G43HBU8+cXFH8xHtWpAQGen8z5CooPzU5yBxDjtgWLYPKiIwntuRr4YkSsB+6muNth\nb+DGMjs135rXFV9JcakF4OURcRTQl1J6oryeza+IuB54N/CXwEBEjIwybUkp7SivZ/MrIj4NfJ/i\nVsh9gfcAb6K4BtsxUkoDwAvmo0TEALAppdRRd3lExD8At1B8GP4R8ClgEPhKmf0qwX8H7miumfE1\nitulz6G4PXZBcyShDSmlr1Hc43oFcB/wWuCUDrzN53UUr389xXyMzwL3Uvwi6CTnAXXgR8CTY37e\nVWKfyrAKWEsxL+E2imH3k53dD3TufKVDgC9TvCe+CjwD9KSUNpXaq3mWUvoFxYJr7wYeAP4r8NGU\n0ldL7dg0VIaHO/W9K0mSJuNIgiRJyjIkSJKkLEOCJEnKMiRIkqQsQ4IkScoyJEiSpCxDgiRJyjIk\nSJKkLEOCJEnKMiRIkqQsQ4IkScr6/y6SNmeCmvzbAAAAAElFTkSuQmCC\n",
"text/plain": [
"<matplotlib.figure.Figure at 0x7f371d7b7290>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"pd.DataFrame(wait_stats).fillna(0).boxplot()"
]
},
{
"cell_type": "markdown",
"metadata": {
"collapsed": true
},
"source": [
"Waitingのほうは横軸が5の位置まであるので、何度か6人待ちの状態が発生していたようです。\n",
"\n",
"それでも、待ち0の箱が最大値の480に集中しているので\n",
"\n",
"待つ事になるのはほとんど起きなさそうです。"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## まとめ\n",
"\n",
"- PythonのシミュレーションフレームワークSimPyを使ってコーヒー注文待ちのシミュレーションをしました。\n",
"\n",
"- 3人のバリスタさんがいれば弊社の規模感の客数ならうまく捌けていそう。\n",
"\n",
"**「お昼休み直後、おやつの時間はよりコーヒーを飲みたくなる」** とか \n",
"\n",
"**「たまに仲間を引き連れて注文しにくる」** みたいなシナリオを考えてみても良さそうです。\n",
"\n",
"それでは良い年末を。 明日のエントリもお楽しみに！"
]
}
],
"metadata": {
"anaconda-cloud": {},
"kernelspec": {
"display_name": "Python 2",
"language": "python",
"name": "python2"
},
"language_info": {
"codemirror_mode": {
"name": "ipython",
"version": 2
},
"file_extension": ".py",
"mimetype": "text/x-python",
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython2",
"version": "2.7.12"
}
},
"nbformat": 4,
"nbformat_minor": 0
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment