chap3ex3.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Exercise 3\n",
    "\n",
    "저자는 “첫째 아이가 좀더 늦게 태어날 것 같은가?” 라는 질문으로 시작했다. 이 질문을 다루는데, 아이 집단간에 평균에 차이를 계산했다. 하지만, 동일 여성에 대해 첫째 아이와 다른 아이 사이에 차이가 있을 수 있다는 가능성을 간과했다.\n",
    "이 질문을 다루는데, 적어도 아이가 두명인 응답자를 선택하고 쌍별 차이(pairwise difference)를 계산한다. 이 질문 구성이 다른 결과를 도출해 내는가?\n",
    "\n",
    "힌트: nsfg.MakePregMap을 사용한다."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "%matplotlib inline\n",
    "import numpy as np\n",
    "import pandas as pd"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "defaultdict(list,\n",
       "            {1: [0, 1],\n",
       "             2: [2, 3, 4],\n",
       "             8195: [8952],\n",
       "             8196: [8953],\n",
       "             8197: [8956],\n",
       "             6: [5, 6, 7],\n",
       "             7: [8],\n",
       "             8200: [8958],\n",
       "             2732: [3085, 3086, 3089, 3090, 3091],\n",
       "             8202: [8959],\n",
       "             12: [10],\n",
       "             8194: [8950],\n",
       "             14: [11, 12],\n",
       "             2733: [3092, 3093],\n",
       "             8208: [8960],\n",
       "             8209: [8961],\n",
       "             18: [17],\n",
       "             21: [19, 20],\n",
       "             23: [21],\n",
       "             24: [23, 24, 25],\n",
       "             4100: [4597],\n",
       "             8219: [8965, 8966],\n",
       "             28: [26],\n",
       "             31: [27, 28, 29],\n",
       "             8224: [8974, 8975],\n",
       "             36: [31],\n",
       "             8229: [8980],\n",
       "             38: [32, 33, 34],\n",
       "             39: [35],\n",
       "             8233: [8985, 8986],\n",
       "             44: [39],\n",
       "             46: [43],\n",
       "             8239: [8993],\n",
       "             49: [44, 47],\n",
       "             51: [48, 49],\n",
       "             8248: [8994, 8995, 8996],\n",
       "             57: [52, 53, 54],\n",
       "             8250: [8997, 8998],\n",
       "             60: [55, 56],\n",
       "             8253: [9002],\n",
       "             8254: [9004, 9005, 9006, 9007],\n",
       "             63: [58, 59],\n",
       "             8256: [9009],\n",
       "             6547: [7226, 7227, 7228],\n",
       "             69: [60],\n",
       "             70: [62, 64],\n",
       "             71: [65],\n",
       "             72: [67],\n",
       "             73: [68],\n",
       "             6839: [7538],\n",
       "             77: [70, 71],\n",
       "             8270: [9017, 9018, 9019],\n",
       "             80: [72, 74, 75, 76],\n",
       "             81: [77],\n",
       "             8274: [9023, 9024],\n",
       "             86: [78, 79],\n",
       "             8279: [9026, 9027, 9028, 9029],\n",
       "             90: [82],\n",
       "             91: [83, 87, 88, 89],\n",
       "             92: [90, 91, 92],\n",
       "             95: [94],\n",
       "             8288: [9031, 9032],\n",
       "             8290: [9033],\n",
       "             8292: [9034],\n",
       "             101: [95, 96],\n",
       "             106: [97, 98, 99],\n",
       "             8300: [9042],\n",
       "             8301: [9043],\n",
       "             8303: [9045, 9046, 9047, 9048],\n",
       "             5732: [6301, 6302],\n",
       "             114: [100],\n",
       "             115: [101],\n",
       "             8308: [9050, 9051, 9052],\n",
       "             118: [104, 105, 106],\n",
       "             119: [108],\n",
       "             8313: [9058],\n",
       "             123: [110],\n",
       "             8319: [9062],\n",
       "             8322: [9067],\n",
       "             132: [113],\n",
       "             8326: [9068, 9070, 9071],\n",
       "             135: [114],\n",
       "             8328: [9076, 9077],\n",
       "             1388: [1554, 1555, 1556],\n",
       "             138: [115],\n",
       "             139: [116, 117],\n",
       "             2754: [3116, 3119],\n",
       "             142: [118, 119, 120, 121],\n",
       "             8336: [9080],\n",
       "             145: [124, 126, 127, 129],\n",
       "             8339: [9083, 9084, 9085],\n",
       "             149: [130, 131, 132],\n",
       "             150: [133, 134, 135, 136, 137],\n",
       "             152: [140, 142, 143],\n",
       "             153: [147, 148],\n",
       "             8347: [9090, 9091],\n",
       "             156: [149],\n",
       "             160: [154, 155, 156],\n",
       "             8355: [9105, 9106, 9108],\n",
       "             8357: [9109],\n",
       "             1393: [1559, 1560, 1561, 1562],\n",
       "             8360: [9112],\n",
       "             8361: [9116],\n",
       "             8362: [9118],\n",
       "             172: [158, 159],\n",
       "             173: [160],\n",
       "             8366: [9120, 9122],\n",
       "             181: [163],\n",
       "             8374: [9124, 9125],\n",
       "             183: [165, 166],\n",
       "             184: [167],\n",
       "             186: [168],\n",
       "             8223: [8970, 8971, 8973],\n",
       "             190: [170, 171, 172, 174],\n",
       "             1397: [1567, 1568, 1569, 1570],\n",
       "             193: [176],\n",
       "             8389: [9136],\n",
       "             8391: [9137],\n",
       "             8392: [9140],\n",
       "             8393: [9143, 9144, 9145],\n",
       "             8394: [9146],\n",
       "             8397: [9151],\n",
       "             8398: [9152],\n",
       "             8399: [9153, 9154],\n",
       "             209: [178],\n",
       "             210: [180, 181, 182],\n",
       "             4131: [4625],\n",
       "             213: [184, 185],\n",
       "             215: [186, 187],\n",
       "             8409: [9158, 9159],\n",
       "             218: [188, 190],\n",
       "             219: [191],\n",
       "             222: [193],\n",
       "             8416: [9162],\n",
       "             8418: [9163],\n",
       "             227: [204, 205],\n",
       "             8230: [8981, 8982, 8983, 8984],\n",
       "             230: [206, 207],\n",
       "             234: [208, 209],\n",
       "             235: [210],\n",
       "             238: [211],\n",
       "             8432: [9173, 9174],\n",
       "             244: [212, 213, 214],\n",
       "             245: [215, 216],\n",
       "             246: [220],\n",
       "             247: [222],\n",
       "             248: [223, 224, 225],\n",
       "             249: [228],\n",
       "             8442: [9180, 9181],\n",
       "             251: [229, 230],\n",
       "             252: [233, 234],\n",
       "             254: [239, 240, 241],\n",
       "             255: [246],\n",
       "             8448: [9186],\n",
       "             8449: [9187],\n",
       "             259: [251, 252],\n",
       "             8452: [9189, 9190],\n",
       "             8454: [9193, 9197],\n",
       "             263: [253, 256],\n",
       "             269: [258, 259, 260, 261],\n",
       "             8464: [9203, 9204],\n",
       "             275: [262],\n",
       "             277: [263, 264, 265],\n",
       "             278: [267, 268, 269, 270, 271, 272, 273],\n",
       "             281: [274, 275, 276],\n",
       "             8475: [9209],\n",
       "             284: [277, 278, 279, 280, 281],\n",
       "             285: [282],\n",
       "             8478: [9215, 9216],\n",
       "             287: [285, 286, 287],\n",
       "             288: [288, 289, 290],\n",
       "             8481: [9218],\n",
       "             8483: [9222, 9223],\n",
       "             8193: [8947, 8948],\n",
       "             1414: [1597, 1598, 1599, 1600],\n",
       "             8486: [9224, 9225],\n",
       "             8488: [9230],\n",
       "             8490: [9231],\n",
       "             299: [299],\n",
       "             300: [305],\n",
       "             302: [306],\n",
       "             8495: [9235],\n",
       "             309: [311, 313],\n",
       "             314: [314, 315, 316, 317, 318],\n",
       "             315: [321, 322],\n",
       "             8508: [9242, 9244],\n",
       "             319: [323, 324, 325, 326],\n",
       "             320: [327, 328, 329],\n",
       "             8513: [9245, 9246, 9247],\n",
       "             8514: [9248],\n",
       "             324: [330, 334],\n",
       "             325: [335, 336],\n",
       "             8518: [9252],\n",
       "             327: [338],\n",
       "             328: [339, 340],\n",
       "             8524: [9257, 9258, 9260, 9261],\n",
       "             334: [341, 343, 346],\n",
       "             336: [348, 349],\n",
       "             8529: [9266, 9267],\n",
       "             8531: [9271, 9272],\n",
       "             340: [350, 351, 352, 353],\n",
       "             8534: [9279],\n",
       "             343: [354, 356],\n",
       "             8538: [9280],\n",
       "             348: [359, 360, 361],\n",
       "             349: [362, 363],\n",
       "             350: [365, 367],\n",
       "             8545: [9284, 9285],\n",
       "             354: [368, 369],\n",
       "             356: [370, 371],\n",
       "             357: [374],\n",
       "             361: [375],\n",
       "             362: [378],\n",
       "             8555: [9294, 9295],\n",
       "             364: [380],\n",
       "             365: [381, 382, 383, 384],\n",
       "             8559: [9300, 9301],\n",
       "             368: [385, 386, 388],\n",
       "             8561: [9306, 9307],\n",
       "             372: [389, 392, 393],\n",
       "             374: [395],\n",
       "             379: [398, 399],\n",
       "             382: [402],\n",
       "             383: [404, 405],\n",
       "             8576: [9322, 9323],\n",
       "             385: [408, 409],\n",
       "             390: [411, 413],\n",
       "             1432: [1612, 1613],\n",
       "             403: [416],\n",
       "             8596: [9339],\n",
       "             8598: [9341, 9343],\n",
       "             407: [419],\n",
       "             408: [421],\n",
       "             8601: [9351, 9353],\n",
       "             8602: [9354, 9355],\n",
       "             413: [423, 428, 429],\n",
       "             8606: [9361, 9362],\n",
       "             415: [433, 435, 436, 437],\n",
       "             416: [438],\n",
       "             418: [440, 441],\n",
       "             419: [443, 444],\n",
       "             420: [447, 448],\n",
       "             8616: [9366],\n",
       "             425: [451],\n",
       "             427: [452],\n",
       "             8620: [9370, 9371, 9372, 9373],\n",
       "             429: [455],\n",
       "             8622: [9376],\n",
       "             8624: [9377],\n",
       "             8625: [9378, 9379, 9380],\n",
       "             434: [459],\n",
       "             8628: [9381, 9382, 9384],\n",
       "             437: [461],\n",
       "             438: [463, 464, 465, 466],\n",
       "             8265: [9013],\n",
       "             442: [468],\n",
       "             444: [470, 474],\n",
       "             445: [475],\n",
       "             8638: [9390, 9392],\n",
       "             8639: [9394, 9395],\n",
       "             8641: [9396],\n",
       "             451: [477, 478],\n",
       "             8644: [9398],\n",
       "             455: [481],\n",
       "             8650: [9402, 9403, 9404],\n",
       "             459: [488, 489],\n",
       "             8652: [9406],\n",
       "             8653: [9407, 9408, 9409],\n",
       "             8654: [9410, 9411],\n",
       "             8269: [9016],\n",
       "             8656: [9412, 9413, 9414],\n",
       "             465: [492, 493],\n",
       "             8658: [9417, 9418, 9419, 9420],\n",
       "             1443: [1633, 1635, 1636, 1637, 1638],\n",
       "             468: [496],\n",
       "             8661: [9421, 9422, 9423],\n",
       "             8662: [9424],\n",
       "             472: [497],\n",
       "             473: [498, 499, 500],\n",
       "             478: [501, 502],\n",
       "             479: [503, 504],\n",
       "             480: [505],\n",
       "             481: [506, 507, 508, 509, 510, 511, 512],\n",
       "             8674: [9429, 9430, 9431],\n",
       "             8675: [9432, 9433],\n",
       "             8677: [9434],\n",
       "             8273: [9020],\n",
       "             8681: [9441, 9442, 9443, 9444],\n",
       "             8682: [9446],\n",
       "             491: [514, 515],\n",
       "             494: [517, 518],\n",
       "             497: [520],\n",
       "             8691: [9449, 9450],\n",
       "             8693: [9451],\n",
       "             8696: [9452],\n",
       "             505: [521, 523],\n",
       "             506: [525],\n",
       "             8699: [9454],\n",
       "             509: [527],\n",
       "             510: [530],\n",
       "             512: [534],\n",
       "             514: [535, 536],\n",
       "             515: [540, 542],\n",
       "             8709: [9468, 9469],\n",
       "             519: [543],\n",
       "             520: [546, 547, 548],\n",
       "             8713: [9472, 9473],\n",
       "             522: [549],\n",
       "             523: [550, 551, 552],\n",
       "             2818: [3178],\n",
       "             8718: [9481],\n",
       "             8720: [9484, 9485],\n",
       "             530: [554, 556],\n",
       "             8724: [9491, 9492],\n",
       "             535: [559],\n",
       "             8728: [9497],\n",
       "             8731: [9499, 9500, 9501],\n",
       "             540: [562, 563, 564, 565],\n",
       "             541: [566],\n",
       "             542: [567, 570],\n",
       "             8735: [9504, 9505],\n",
       "             548: [573],\n",
       "             8743: [9510, 9511],\n",
       "             552: [577, 578],\n",
       "             553: [579, 580],\n",
       "             555: [582],\n",
       "             556: [583],\n",
       "             560: [584, 585],\n",
       "             562: [590],\n",
       "             8755: [9514, 9515],\n",
       "             564: [591, 592],\n",
       "             8757: [9517],\n",
       "             567: [594, 595],\n",
       "             8760: [9521, 9522, 9523, 9524],\n",
       "             1460: [1648, 1649],\n",
       "             571: [597, 598],\n",
       "             572: [599, 601],\n",
       "             573: [602, 603],\n",
       "             574: [606, 609, 610],\n",
       "             8767: [9531],\n",
       "             578: [614],\n",
       "             580: [616],\n",
       "             581: [617, 618, 620],\n",
       "             7096: [7748],\n",
       "             8775: [9537],\n",
       "             586: [624, 625, 626, 627],\n",
       "             587: [629],\n",
       "             589: [630, 631, 633],\n",
       "             590: [634, 635],\n",
       "             591: [636],\n",
       "             8784: [9546],\n",
       "             595: [639, 640, 641],\n",
       "             596: [642],\n",
       "             599: [644],\n",
       "             601: [647, 648, 649, 650, 651],\n",
       "             602: [656],\n",
       "             6927: [7612],\n",
       "             8797: [9552, 9553, 9555],\n",
       "             606: [657, 658, 659, 660, 661],\n",
       "             607: [662, 663],\n",
       "             8800: [9557],\n",
       "             609: [664, 665],\n",
       "             8802: [9561],\n",
       "             8808: [9562, 9563],\n",
       "             8811: [9564, 9566],\n",
       "             620: [668],\n",
       "             622: [669],\n",
       "             623: [671],\n",
       "             8817: [9568, 9569],\n",
       "             627: [672],\n",
       "             628: [673, 674],\n",
       "             6022: [6654, 6655, 6656],\n",
       "             632: [676, 677],\n",
       "             635: [680, 683, 684, 685, 687, 688],\n",
       "             637: [689, 690, 691],\n",
       "             8831: [9576, 9577, 9578],\n",
       "             642: [693, 694, 695],\n",
       "             8835: [9585, 9586],\n",
       "             8838: [9587],\n",
       "             8839: [9588, 9589, 9590, 9591],\n",
       "             652: [697],\n",
       "             653: [699, 700],\n",
       "             8847: [9604, 9605],\n",
       "             656: [703],\n",
       "             657: [705, 706, 707],\n",
       "             8850: [9610],\n",
       "             659: [708, 709, 710],\n",
       "             663: [713],\n",
       "             664: [714, 715],\n",
       "             665: [716, 717],\n",
       "             8859: [9616, 9617, 9618],\n",
       "             668: [719, 720],\n",
       "             8861: [9621],\n",
       "             671: [721],\n",
       "             8865: [9622],\n",
       "             8866: [9624],\n",
       "             676: [725],\n",
       "             8869: [9633],\n",
       "             678: [726],\n",
       "             8871: [9635, 9636, 9637],\n",
       "             681: [730],\n",
       "             8874: [9641],\n",
       "             684: [731, 732],\n",
       "             686: [739],\n",
       "             8880: [9642, 9643],\n",
       "             689: [741, 742],\n",
       "             693: [743, 744, 745, 746],\n",
       "             695: [747, 749],\n",
       "             10752: [11673, 11674],\n",
       "             8890: [9653, 9654],\n",
       "             8893: [9655, 9656],\n",
       "             702: [756, 757, 760, 761, 762, 763],\n",
       "             707: [766, 767, 768, 769, 770],\n",
       "             708: [772, 774],\n",
       "             710: [776, 777, 778],\n",
       "             711: [779],\n",
       "             712: [780, 781, 782],\n",
       "             714: [784],\n",
       "             8912: [9673, 9674],\n",
       "             8913: [9675, 9676],\n",
       "             723: [785],\n",
       "             724: [792],\n",
       "             8918: [9677, 9678],\n",
       "             8920: [9679, 9680, 9681],\n",
       "             8922: [9683],\n",
       "             8923: [9684],\n",
       "             732: [793, 794],\n",
       "             8925: [9687],\n",
       "             734: [795, 796],\n",
       "             737: [798, 799, 800],\n",
       "             739: [802],\n",
       "             8932: [9691],\n",
       "             8934: [9692],\n",
       "             743: [804, 805, 807],\n",
       "             745: [808, 809, 810, 811],\n",
       "             8938: [9701, 9703, 9704],\n",
       "             750: [813, 814],\n",
       "             752: [816, 817],\n",
       "             8945: [9708, 9710, 9711, 9712],\n",
       "             755: [819],\n",
       "             756: [822, 823, 824],\n",
       "             758: [826, 828],\n",
       "             759: [829, 830],\n",
       "             760: [832, 833],\n",
       "             8954: [9714],\n",
       "             8959: [9715],\n",
       "             8352: [9102, 9104],\n",
       "             8962: [9716],\n",
       "             772: [842],\n",
       "             8965: [9718, 9719],\n",
       "             775: [845],\n",
       "             776: [846, 847],\n",
       "             8972: [9721],\n",
       "             781: [850],\n",
       "             782: [851, 853],\n",
       "             8975: [9722],\n",
       "             784: [854, 855],\n",
       "             785: [857],\n",
       "             786: [859],\n",
       "             787: [861],\n",
       "             788: [863],\n",
       "             789: [864],\n",
       "             790: [866, 867, 868],\n",
       "             792: [869],\n",
       "             8986: [9734, 9735],\n",
       "             796: [870, 871, 872, 873, 874],\n",
       "             4229: [4732],\n",
       "             8994: [9738],\n",
       "             804: [877, 878],\n",
       "             805: [879, 880],\n",
       "             8998: [9741],\n",
       "             8999: [9742],\n",
       "             808: [885, 887],\n",
       "             9001: [9747, 9748],\n",
       "             810: [888],\n",
       "             811: [889, 890, 891],\n",
       "             814: [893, 894, 895],\n",
       "             816: [897],\n",
       "             9009: [9751],\n",
       "             820: [900, 901, 902, 903],\n",
       "             9013: [9753, 9754],\n",
       "             822: [905],\n",
       "             6875: [7575, 7576],\n",
       "             827: [909, 910, 911],\n",
       "             829: [912, 913],\n",
       "             830: [914, 915, 916],\n",
       "             832: [917, 919],\n",
       "             9026: [9765, 9766],\n",
       "             835: [924],\n",
       "             836: [925, 926],\n",
       "             9033: [9767, 9768, 9769, 9770, 9771],\n",
       "             843: [928, 929],\n",
       "             844: [931, 932, 933, 934, 935],\n",
       "             845: [936, 937, 938, 939, 940, 941],\n",
       "             847: [944, 945],\n",
       "             848: [947, 948, 949],\n",
       "             849: [950, 951, 952, 953],\n",
       "             9042: [9780, 9781],\n",
       "             9043: [9782, 9784],\n",
       "             852: [956, 957, 958, 959],\n",
       "             854: [960, 961, 962],\n",
       "             859: [963, 964],\n",
       "             861: [965],\n",
       "             862: [969, 970],\n",
       "             863: [971, 972, 973],\n",
       "             865: [978, 979, 980, 981],\n",
       "             867: [983],\n",
       "             868: [984, 985],\n",
       "             9061: [9800, 9802],\n",
       "             9062: [9807, 9808, 9809],\n",
       "             8337: [9081],\n",
       "             872: [987, 989, 990],\n",
       "             874: [991, 992],\n",
       "             9067: [9812],\n",
       "             876: [994],\n",
       "             877: [1000],\n",
       "             9070: [9816],\n",
       "             880: [1002],\n",
       "             881: [1003, 1004],\n",
       "             882: [1006, 1007, 1009],\n",
       "             9075: [9827, 9828],\n",
       "             884: [1011],\n",
       "             886: [1013],\n",
       "             9081: [9830, 9831],\n",
       "             890: [1016],\n",
       "             9083: [9832],\n",
       "             9085: [9833, 9834],\n",
       "             895: [1017, 1018],\n",
       "             9088: [9839, 9840],\n",
       "             897: [1021, 1022],\n",
       "             9090: [9841, 9842],\n",
       "             899: [1023, 1024],\n",
       "             900: [1025],\n",
       "             901: [1027, 1028, 1029],\n",
       "             904: [1031, 1032],\n",
       "             1516: [1700, 1701],\n",
       "             910: [1034],\n",
       "             9103: [9856, 9857],\n",
       "             9106: [9858],\n",
       "             917: [1035],\n",
       "             918: [1036, 1037],\n",
       "             8345: [9087, 9089],\n",
       "             921: [1038],\n",
       "             9114: [9868],\n",
       "             9116: [9871],\n",
       "             9118: [9873, 9874],\n",
       "             9120: [9875, 9876, 9877, 9878],\n",
       "             929: [1042, 1044],\n",
       "             930: [1045, 1046],\n",
       "             932: [1048, 1050],\n",
       "             9128: [9885, 9886, 9887],\n",
       "             8348: [9092, 9093, 9094],\n",
       "             9131: [9889, 9890],\n",
       "             943: [1056],\n",
       "             4123: [4623],\n",
       "             947: [1060],\n",
       "             948: [1062, 1063],\n",
       "             953: [1065, 1066],\n",
       "             954: [1067, 1068, 1070],\n",
       "             4207: [4706],\n",
       "             9149: [9894],\n",
       "             959: [1075],\n",
       "             9152: [9897],\n",
       "             961: [1077, 1078, 1079, 1080],\n",
       "             9154: [9899, 9900, 9903],\n",
       "             10214: [11077, 11078],\n",
       "             9157: [9904, 9905, 9907, 9910],\n",
       "             966: [1081, 1082, 1083],\n",
       "             967: [1084],\n",
       "             9160: [9917, 9918],\n",
       "             969: [1089, 1090],\n",
       "             971: [1091, 1092],\n",
       "             972: [1094, 1096],\n",
       "             9166: [9923, 9924],\n",
       "             9169: [9925, 9927],\n",
       "             9171: [9929, 9930],\n",
       "             980: [1098, 1102],\n",
       "             2894: [3239],\n",
       "             9174: [9931, 9932, 9933, 9934],\n",
       "             984: [1103, 1104],\n",
       "             985: [1105],\n",
       "             989: [1112, 1113],\n",
       "             990: [1115],\n",
       "             992: [1118],\n",
       "             9185: [9947],\n",
       "             9186: [9949, 9950],\n",
       "             996: [1119, 1120],\n",
       "             997: [1121, 1122],\n",
       "             999: [1125],\n",
       "             1000: [1126, 1127, 1129],\n",
       "             1004: [1130],\n",
       "             1005: [1131],\n",
       "             1007: [1133, 1134],\n",
       "             9200: [9953],\n",
       "             9202: [9955, 9956],\n",
       "             9203: [9957],\n",
       "             9204: [9961],\n",
       "             5762: [6337, 6338],\n",
       "             1015: [1139],\n",
       "             1016: [1144],\n",
       "             1017: [1145, 1147, 1148],\n",
       "             9210: [9966, 9967, 9968],\n",
       "             9213: [9969, 9970, 9972, 9973],\n",
       "             9215: [9975],\n",
       "             1025: [1149],\n",
       "             1026: [1151],\n",
       "             1028: [1152],\n",
       "             1029: [1156],\n",
       "             1031: [1160],\n",
       "             9225: [9979],\n",
       "             9227: [9980, 9981, 9982],\n",
       "             9228: [9983],\n",
       "             9229: [9987, 9988, 9989],\n",
       "             9230: [9990],\n",
       "             1040: [1163],\n",
       "             9233: [9996, 9997],\n",
       "             1042: [1166, 1167, 1168],\n",
       "             9235: [9999, 10000, 10001],\n",
       "             1044: [1172, 1173],\n",
       "             9244: [10006, 10007],\n",
       "             9247: [10009],\n",
       "             1058: [1177, 1178, 1179, 1180],\n",
       "             1060: [1182],\n",
       "             1061: [1184, 1185],\n",
       "             1062: [1186, 1187, 1188],\n",
       "             1063: [1192],\n",
       "             1064: [1193],\n",
       "             9258: [10015],\n",
       "             9259: [10017, 10018],\n",
       "             1068: [1194, 1195],\n",
       "             1069: [1197],\n",
       "             9262: [10019, 10020],\n",
       "             9266: [10022, 10024, 10026],\n",
       "             1075: [1200, 1203],\n",
       "             1077: [1205],\n",
       "             9270: [10029],\n",
       "             9272: [10030],\n",
       "             9273: [10032, 10033],\n",
       "             9275: [10034, 10035, 10036],\n",
       "             1084: [1208, 1209],\n",
       "             9277: [10041, 10042, 10043, 10044],\n",
       "             1086: [1210, 1212, 1213, 1214, 1215],\n",
       "             1087: [1217],\n",
       "             1088: [1218, 1219, 1220],\n",
       "             1089: [1222, 1226],\n",
       "             1090: [1228],\n",
       "             9283: [10048],\n",
       "             1092: [1233, 1235],\n",
       "             1094: [1236],\n",
       "             1095: [1237],\n",
       "             9288: [10049, 10050, 10051],\n",
       "             1098: [1238],\n",
       "             1099: [1240],\n",
       "             1100: [1242, 1243, 1244],\n",
       "             9293: [10055],\n",
       "             9296: [10056],\n",
       "             9299: [10059, 10060, 10061],\n",
       "             9301: [10062, 10063],\n",
       "             9302: [10064, 10065, 10066],\n",
       "             9304: [10068, 10069],\n",
       "             1113: [1253, 1256],\n",
       "             1115: [1257],\n",
       "             8378: [9128, 9129],\n",
       "             1120: [1258, 1259],\n",
       "             1121: [1260, 1262],\n",
       "             1122: [1263],\n",
       "             1123: [1265],\n",
       "             9319: [10087, 10088, 10089],\n",
       "             9320: [10090, 10091],\n",
       "             9322: [10093, 10096, 10097],\n",
       "             1131: [1271],\n",
       "             9325: [10099, 10100],\n",
       "             1134: [1274, 1276],\n",
       "             9327: [10103],\n",
       "             1138: [1277, 1279],\n",
       "             1139: [1280],\n",
       "             1140: [1281],\n",
       "             9334: [10109, 10110, 10111],\n",
       "             9338: [10114, 10115, 10117],\n",
       "             1147: [1284, 1285],\n",
       "             9340: [10120],\n",
       "             9341: [10121],\n",
       "             1151: [1288],\n",
       "             1153: [1289, 1290, 1291],\n",
       "             1558: [1742],\n",
       "             9350: [10127],\n",
       "             1159: [1293],\n",
       "             1160: [1294, 1295, 1296],\n",
       "             1161: [1298],\n",
       "             1164: [1299],\n",
       "             1165: [1302],\n",
       "             1166: [1305],\n",
       "             1167: [1307],\n",
       "             1169: [1310, 1311, 1314, 1315],\n",
       "             1171: [1320, 1321],\n",
       "             1172: [1324],\n",
       "             1176: [1327, 1328, 1329],\n",
       "             9371: [10135, 10136],\n",
       "             9373: [10138],\n",
       "             9374: [10139],\n",
       "             9376: [10141],\n",
       "             1185: [1330, 1332],\n",
       "             4294: [4809, 4810, 4811],\n",
       "             1190: [1333, 1335],\n",
       "             9383: [10142],\n",
       "             9385: [10143],\n",
       "             9388: [10146, 10148],\n",
       "             1197: [1343],\n",
       "             9390: [10153, 10154],\n",
       "             9393: [10158],\n",
       "             1203: [1344, 1345],\n",
       "             1204: [1346, 1347, 1348],\n",
       "             9398: [10160, 10161, 10162],\n",
       "             1209: [1350, 1351, 1352],\n",
       "             1210: [1353, 1354, 1355],\n",
       "             1211: [1357],\n",
       "             4298: [4814, 4815, 4816],\n",
       "             1214: [1365, 1366, 1367, 1368],\n",
       "             9407: [10174],\n",
       "             1217: [1372],\n",
       "             1218: [1373],\n",
       "             1219: [1377, 1378],\n",
       "             9413: [10180, 10181, 10184, 10185],\n",
       "             1223: [1379, 1380],\n",
       "             9416: [10188],\n",
       "             1226: [1381],\n",
       "             9419: [10192, 10193],\n",
       "             1570: [1749, 1750],\n",
       "             9423: [10194],\n",
       "             9424: [10197],\n",
       "             1237: [1388, 1389],\n",
       "             1239: [1391],\n",
       "             1241: [1393],\n",
       "             1242: [1395],\n",
       "             1244: [1397, 1399],\n",
       "             1246: [1400, 1401],\n",
       "             9440: [10206],\n",
       "             9441: [10207, 10208],\n",
       "             1251: [1403, 1404, 1405],\n",
       "             9444: [10212, 10214],\n",
       "             9446: [10215, 10216, 10218, 10219],\n",
       "             9448: [10222, 10223],\n",
       "             1258: [1408, 1409],\n",
       "             1261: [1412, 1413],\n",
       "             2941: [3295],\n",
       "             1264: [1414, 1415],\n",
       "             1265: [1418, 1419],\n",
       "             9459: [10233, 10237],\n",
       "             1268: [1421],\n",
       "             9463: [10239],\n",
       "             9464: [10240],\n",
       "             9465: [10243, 10246],\n",
       "             9466: [10250, 10251, 10253, 10254, 10255],\n",
       "             1275: [1426, 1428],\n",
       "             1277: [1429, 1431, 1433],\n",
       "             9470: [10256, 10258],\n",
       "             1279: [1435, 1436, 1437],\n",
       "             1280: [1439, 1440, 1441, 1443],\n",
       "             9473: [10264, 10265, 10266],\n",
       "             1282: [1444, 1445],\n",
       "             1283: [1446, 1447],\n",
       "             1284: [1449],\n",
       "             1286: [1452],\n",
       "             1288: [1453, 1455],\n",
       "             1292: [1456, 1457, 1458],\n",
       "             1295: [1460, 1461, 1462, 1463],\n",
       "             1296: [1466, 1467],\n",
       "             1298: [1468],\n",
       "             9492: [10275, 10277],\n",
       "             1582: [1763, 1764],\n",
       "             1302: [1469],\n",
       "             1305: [1473, 1474, 1475],\n",
       "             1307: [1477],\n",
       "             8410: [9160],\n",
       "             9502: [10290],\n",
       "             1314: [1485, 1486, 1488],\n",
       "             1315: [1489, 1490],\n",
       "             2950: [3308, 3309],\n",
       "             9510: [10294],\n",
       "             1319: [1492],\n",
       "             9512: [10295],\n",
       "             1321: [1494, 1495],\n",
       "             9514: [10296, 10297],\n",
       "             9515: [10299],\n",
       "             9516: [10300],\n",
       "             9517: [10305, 10307, 10309],\n",
       "             9518: [10312],\n",
       "             9519: [10314, 10315, 10317],\n",
       "             2952: [3314, 3315],\n",
       "             9522: [10319, 10321],\n",
       "             1333: [1496, 1498, 1499],\n",
       "             9526: [10322, 10323, 10324, 10325],\n",
       "             1338: [1500, 1501, 1502],\n",
       "             1339: [1503, 1504, 1505],\n",
       "             1342: [1507, 1509],\n",
       "             5685: [6236, 6243, 6244],\n",
       "             1344: [1511],\n",
       "             1346: [1513, 1515],\n",
       "             1349: [1516, 1518, 1519],\n",
       "             1351: [1520, 1521, 1522, 1523],\n",
       "             1353: [1524],\n",
       "             9546: [10329, 10330, 10331],\n",
       "             9547: [10332],\n",
       "             9550: [10334, 10335, 10336],\n",
       "             1359: [1526, 1527],\n",
       "             9552: [10337, 10338],\n",
       "             1362: [1532, 1534],\n",
       "             9558: [10342, 10343],\n",
       "             9559: [10344],\n",
       "             9560: [10346],\n",
       "             2959: [3320, 3321, 3322, 3323],\n",
       "             9564: [10351, 10352],\n",
       "             1373: [1543, 1544],\n",
       "             9567: [10353, 10355, 10356],\n",
       "             9569: [10360, 10361],\n",
       "             9571: [10362],\n",
       "             9572: [10363, 10364, 10365, 10366, 10367],\n",
       "             8422: [9171],\n",
       "             9576: [10371],\n",
       "             9577: [10374],\n",
       "             1386: [1553],\n",
       "             9579: [10379, 10380],\n",
       "             9580: [10381, 10384, 10385],\n",
       "             9582: [10386, 10388],\n",
       "             1391: [1557],\n",
       "             9584: [10389, 10390],\n",
       "             9585: [10393],\n",
       "             2963: [3325, 3326],\n",
       "             1396: [1563, 1565, 1566],\n",
       "             9589: [10395],\n",
       "             9590: [10396, 10397],\n",
       "             1399: [1572, 1573, 1574, 1575],\n",
       "             1401: [1576, 1577],\n",
       "             1403: [1578, 1579],\n",
       "             1404: [1580],\n",
       "             1405: [1581, 1582],\n",
       "             1406: [1584],\n",
       "             9599: [10402, 10403, 10404, 10405],\n",
       "             1408: [1585],\n",
       "             1409: [1586],\n",
       "             1410: [1589, 1590, 1591],\n",
       "             1411: [1594, 1595, 1596],\n",
       "             9605: [10409],\n",
       "             9606: [10410],\n",
       "             1415: [1601],\n",
       "             9608: [10415, 10416, 10417],\n",
       "             1417: [1605],\n",
       "             9611: [10418, 10419],\n",
       "             9612: [10420],\n",
       "             9613: [10422],\n",
       "             1425: [1606],\n",
       "             1426: [1609, 1610],\n",
       "             9622: [10427],\n",
       "             9624: [10429, 10430, 10431],\n",
       "             1434: [1617, 1618, 1619],\n",
       "             1435: [1620, 1621],\n",
       "             1436: [1622, 1623],\n",
       "             1438: [1624],\n",
       "             9631: [10434, 10435, 10436],\n",
       "             1440: [1625, 1626, 1627, 1628],\n",
       "             1442: [1629, 1632],\n",
       "             9635: [10443],\n",
       "             9636: [10445, 10447, 10448],\n",
       "             1446: [1639, 1640, 1641],\n",
       "             1447: [1642, 1643],\n",
       "             1448: [1644],\n",
       "             2972: [3338, 3339],\n",
       "             9642: [10453],\n",
       "             9644: [10457, 10459],\n",
       "             9645: [10461, 10462],\n",
       "             1455: [1645],\n",
       "             9648: [10466],\n",
       "             9650: [10471],\n",
       "             9652: [10472, 10473],\n",
       "             2974: [3343],\n",
       "             1465: [1650, 1652, 1654],\n",
       "             1466: [1655, 1657, 1658, 1659],\n",
       "             9659: [10479],\n",
       "             9663: [10480, 10481],\n",
       "             1474: [1660],\n",
       "             1611: [1801, 1802, 1803, 1804, 1805, 1806, 1807],\n",
       "             1476: [1662, 1663],\n",
       "             1477: [1666],\n",
       "             1478: [1668],\n",
       "             2977: [3347, 3348, 3349, 3350],\n",
       "             1481: [1669],\n",
       "             9674: [10486, 10487],\n",
       "             8439: [9176, 9177],\n",
       "             9680: [10489, 10490, 10491],\n",
       "             8440: [9178],\n",
       "             9682: [10494, 10496, 10497],\n",
       "             1491: [1672, 1673, 1674, 1675],\n",
       "             1493: [1676, 1677],\n",
       "             1494: [1680],\n",
       "             9687: [10500, 10501],\n",
       "             2980: [3354, 3356, 3357],\n",
       "             9690: [10504, 10505, 10506],\n",
       "             1503: [1681, 1682],\n",
       "             9696: [10512, 10513],\n",
       "             1507: [1688],\n",
       "             9701: [10514, 10515],\n",
       "             1510: [1689],\n",
       "             1511: [1690, 1691, 1692],\n",
       "             1512: [1693],\n",
       "             1513: [1695, 1696],\n",
       "             1515: [1699],\n",
       "             9708: [10517, 10518, 10519, 10521],\n",
       "             9709: [10523],\n",
       "             1518: [1703, 1704, 1705, 1706],\n",
       "             9713: [10524, 10525],\n",
       "             9714: [10527],\n",
       "             9598: [10399],\n",
       "             9718: [10532, 10533, 10534],\n",
       "             1527: [1709, 1710, 1711, 1712],\n",
       "             9720: [10536],\n",
       "             1535: [1716],\n",
       "             4352: [4879],\n",
       "             1538: [1719],\n",
       "             1539: [1721],\n",
       "             1541: [1722, 1723],\n",
       "             1542: [1724, 1725],\n",
       "             9735: [10549, 10550, 10551],\n",
       "             1545: [1726, 1727],\n",
       "             1547: [1728],\n",
       "             1549: [1731, 1733, 1734],\n",
       "             1555: [1738, 1739],\n",
       "             1556: [1740, 1741],\n",
       "             9750: [10554, 10555],\n",
       "             9751: [10557],\n",
       "             1560: [1743],\n",
       "             1561: [1744, 1745],\n",
       "             1565: [1746, 1747],\n",
       "             9759: [10563],\n",
       "             1568: [1748],\n",
       "             9761: [10567, 10568],\n",
       "             9762: [10569, 10570],\n",
       "             1571: [1751],\n",
       "             1572: [1752, 1754],\n",
       "             9766: [10580],\n",
       "             1575: [1756, 1757, 1758],\n",
       "             1577: [1760],\n",
       "             9770: [10582, 10583, 10584],\n",
       "             8455: [9199],\n",
       "             9772: [10585],\n",
       "             9774: [10586],\n",
       "             1629: [1821, 1822, 1823],\n",
       "             9776: [10587, 10588],\n",
       "             9778: [10591],\n",
       "             9779: [10593],\n",
       "             1588: [1766, 1767, 1768],\n",
       "             9781: [10595, 10596, 10597],\n",
       "             9785: [10598, 10599, 10601],\n",
       "             9786: [10603, 10606, 10611, 10618],\n",
       "             1595: [1769, 1770, 1771],\n",
       "             1596: [1773],\n",
       "             1597: [1780, 1785, 1790],\n",
       "             1598: [1791, 1792],\n",
       "             1599: [1793, 1794],\n",
       "             9796: [10633, 10634],\n",
       "             1606: [1795, 1796, 1797, 1798],\n",
       "             9799: [10640],\n",
       "             9800: [10641],\n",
       "             9802: [10644, 10645],\n",
       "             9803: [10646, 10647],\n",
       "             9804: [10651, 10653],\n",
       "             9809: [10654, 10655, 10656, 10658, 10659],\n",
       "             1618: [1808, 1809, 1810],\n",
       "             9812: [10660, 10661, 10662],\n",
       "             1621: [1812, 1813],\n",
       "             7419: [8101, 8103, 8105],\n",
       "             9817: [10666],\n",
       "             1626: [1816, 1817],\n",
       "             9819: [10667, 10668],\n",
       "             9820: [10670, 10671, 10674, 10675],\n",
       "             9821: [10676],\n",
       "             1630: [1826, 1827],\n",
       "             9824: [10680, 10681],\n",
       "             9825: [10682, 10683, 10684],\n",
       "             1636: [1831, 1833, 1834],\n",
       "             9830: [10685],\n",
       "             1639: [1840, 1841],\n",
       "             9832: [10688, 10689],\n",
       "             9833: [10690],\n",
       "             1642: [1843, 1844],\n",
       "             9602: [10407, 10408],\n",
       "             1647: [1846],\n",
       "             1649: [1847, 1849],\n",
       "             1652: [1851, 1852, 1853, 1855],\n",
       "             1654: [1858, 1859],\n",
       "             1641: [1842],\n",
       "             1656: [1867, 1868],\n",
       "             1657: [1870],\n",
       "             1658: [1871],\n",
       "             1659: [1874],\n",
       "             9853: [10707, 10708],\n",
       "             9855: [10710],\n",
       "             1665: [1877],\n",
       "             1666: [1879],\n",
       "             9859: [10712, 10713],\n",
       "             1668: [1880, 1881, 1882],\n",
       "             ...})"
      ]
     },
     "execution_count": 2,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "import first\n",
    "import nsfg\n",
    "\n",
    "live, _, _ = first.MakeFrames()\n",
    "live = live[live.prglngth >= 37]\n",
    "preg_map = nsfg.MakePregMap(live)\n",
    "preg_map"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "-6.0     0.001533\n",
       "-5.0     0.006897\n",
       "-4.0     0.015581\n",
       "-3.0     0.028608\n",
       "-2.0     0.084036\n",
       "-1.0     0.078672\n",
       " 0.0     0.652107\n",
       " 1.0     0.059770\n",
       " 2.0     0.048787\n",
       " 3.0     0.012261\n",
       " 4.0     0.007663\n",
       " 5.0     0.002299\n",
       " 6.0     0.000766\n",
       " 7.0     0.000511\n",
       " 10.0    0.000255\n",
       " 12.0    0.000255\n",
       "dtype: float64"
      ]
     },
     "execution_count": 3,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "data = np.concatenate([live.loc[indices[1:]].prglngth - live.loc[indices[0]].prglngth for indices in preg_map.values() if len(indices) > 1])\n",
    "pmf = pd.Series(data).value_counts().sort_index()\n",
    "pmf /= pmf.sum()\n",
    "pmf"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.axes._subplots.AxesSubplot at 0x7f7698010b90>"
      ]
     },
     "execution_count": 4,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAgAAAAF1CAYAAACef1IVAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAAPYQAAD2EBqD+naQAAIABJREFUeJzt3X2UXXV97/H3lwdJo3WqDiZQoYgPMa5elZliTX1sowaW\nV9parzBoRbBaKi7t1Nvae20virUstZCKGsFaDRQZjW0t1GuNhWofpBHNGHprA1oJRhBiRnSoJpGQ\nfO8fe49MDnNm5pzZZyYzv/drrbMy57f373x/OzNnn89+PJGZSJKkshyx0AOQJEnzzwAgSVKBDACS\nJBXIACBJUoEMAJIkFcgAIElSgQwAkiQVyAAgSVKBDACSJBXIACBJUoG6CgARcUFE7IiIvRGxJSJO\nnWbej0TEwYg4UP878fh/3Q9bkiTNRccBICLOBC4BLgROAW4GNkdEf5subwBWAsfV/z4GuAfY1M2A\nJUnS3EWnXwYUEVuAL2bmG+vnAXwLuCwz3zWL/r8C/CXw2Mz8VudDliRJc9XRHoCIOBoYBG6YaMsq\nQVwPrJnly5wHXO+HvyRJC+eoDufvB44EdrW07wJWzdQ5Io4DTgfOmmG+RwHrgNuBfR2OUZKkki0D\nTgI2Z+Z3283UaQCYq1cB3wOunWG+dcBHez4aSZKWrpcD17Sb2GkAGAMOACta2lcAd8+i/7nAVZl5\n/wzz3Q5w9dVXs3r16o4GODw8zPr16zvqM1fWXFo1F6quNa1pTWs2UXP79u284hWvgPqztJ2OAkBm\n7o+IrcBa4Dr48UmAa4HLpusbEc8DHgf8+SxK7QNYvXo1AwMDnQyRvr6+jvvMlTWXVs2FqmtNa1rT\nmg3XnPYQejeHAC4FNtZB4CZgGFgObASIiIuB4zPznJZ+r6a6emB7FzUlSVKDOg4Ambmpvub/Iqpd\n/9uAdZm5u55lJXDC5D4R8XDgV6nuCSBJkhZYVycBZuYGYEObaedO0XYv8LBuakmSpOYtue8CGBoa\nsqY1F2Vda1rTmtacz5od3wlwPkTEALB169atC3ICmCRJi9Xo6CiDg4MAg5k52m6+JbcHQJIkzcwA\nIElSgQwAkiQVyAAgSVKBDACSJBXIACBJUoEMAJIkFcgAIElSgQwAkiQVyAAgSVKBDACSJBXIACBJ\nUoEMAJIkFeiohR6AVIqdO3cyNjbWVd/+/n5OPPHEhkckqWQGAGke7Ny5k1WrVrNv356u+i9btpxb\nb91uCJDUGAOANA/GxsbqD/+rgdUd9t7Ovn2vYGxszAAgqTEGAGlerQYGFnoQkuRJgJIklcgAIElS\ngQwAkiQVyAAgSVKBDACSJBXIACBJUoEMAJIkFcgAIElSgQwAkiQVyAAgSVKBDACSJBXIACBJUoEM\nAJIkFcgAIElSgboKABFxQUTsiIi9EbElIk6dYf6HRMQ7IuL2iNgXEbdFxKu6GrEkSZqzozrtEBFn\nApcArwVuAoaBzRHxxMwca9PtE8CxwLnAN4DjcO+DJEkLpuMAQPWBf0VmXgUQEecDLwLOA97VOnNE\nnAY8Gzg5M79fN+/sbriSJKkJHW2FR8TRwCBww0RbZiZwPbCmTbcXA18G3hwRd0TErRHx7ohY1uWY\nJUnSHHW6B6AfOBLY1dK+C1jVps/JVHsA9gG/Ur/GB4BHAq/usL4kSWpAN4cAOnUEcBA4OzN/ABAR\nvwN8IiJel5k/mocxSJKkSToNAGPAAWBFS/sK4O42fe4C7pz48K9tBwJ4DNVJgVMaHh6mr6/vkLah\noSGGhoY6HLYkSUvPyMgIIyMjh7SNj4/Pqm9HASAz90fEVmAtcB1ARET9/LI23b4AvDQilmfmnrpt\nFdVegTumq7d+/XoGBgY6GaIkScWYaqN4dHSUwcHBGft2cynepcBrIuKVEfEk4HJgObARICIujogr\nJ81/DfBd4CMRsToinkN1tcCfu/tfkqSF0fE5AJm5KSL6gYuodv1vA9Zl5u56lpXACZPm/2FEvAB4\nL/AlqjDwceAP5zh2SZLUpa5OAszMDcCGNtPOnaLta8C6bmpJkqTmeTc+SZIKZACQJKlABgBJkgpk\nAJAkqUAGAEmSCmQAkCSpQAYASZIKZACQJKlABgBJkgpkAJAkqUAGAEmSCmQAkCSpQAYASZIKZACQ\nJKlABgBJkgpkAJAkqUAGAEmSCmQAkCSpQAYASZIKZACQJKlABgBJkgpkAJAkqUAGAEmSCmQAkCSp\nQAYASZIKZACQJKlABgBJkgpkAJAkqUAGAEmSCmQAkCSpQAYASZIKZACQJKlABgBJkgpkAJAkqUBd\nBYCIuCAidkTE3ojYEhGnTjPvcyPiYMvjQEQ8uvthS5Kkueg4AETEmcAlwIXAKcDNwOaI6J+mWwJP\nAFbWj+My8zudD1eSJDWhmz0Aw8AVmXlVZt4CnA/sAc6bod/uzPzOxKOLupIkqSEdBYCIOBoYBG6Y\naMvMBK4H1kzXFdgWEd+OiM9GxC90M1hJktSMTvcA9ANHArta2ndR7dqfyl3AbwK/BrwE+Bbw+Yh4\nWoe1JUlSQ47qdYHM/BrwtUlNWyLicVSHEs6Zru/w8DB9fX2HtA0NDTE0NNT4OCVJWmxGRkYYGRk5\npG18fHxWfTsNAGPAAWBFS/sK4O4OXucm4JkzzbR+/XoGBgY6eFlJksox1Ubx6Ogog4ODM/bt6BBA\nZu4HtgJrJ9oiIurnN3bwUk+jOjQgSZIWQDeHAC4FNkbEVqot+WFgObARICIuBo7PzHPq528EdgBf\nBZYBrwF+EXjBXAcvSZK603EAyMxN9TX/F1Ht+t8GrMvM3fUsK4ETJnV5CNV9A46nulzw34C1mflP\ncxm4JEnqXlcnAWbmBmBDm2nntjx/N/DubupIkqTe8LsAJEkqkAFAkqQCGQAkSSqQAUCSpAIZACRJ\nKpABQJKkAhkAJEkqkAFAkqQCGQAkSSqQAUCSpAIZACRJKpABQJKkAhkAJEkqkAFAkqQCGQAkSSqQ\nAUCSpAIZACRJKpABQJKkAhkAJEkqkAFAkqQCGQAkSSqQAUCSpAIZACRJKpABQJKkAhkAJEkqkAFA\nkqQCGQAkSSqQAUCSpAIZACRJKpABQJKkAhkAJEkqkAFAkqQCGQAkSSpQVwEgIi6IiB0RsTcitkTE\nqbPs98yI2B8Ro93UlSRJzeg4AETEmcAlwIXAKcDNwOaI6J+hXx9wJXB9F+OUJEkN6mYPwDBwRWZe\nlZm3AOcDe4DzZuh3OfBRYEsXNSVJUoM6CgARcTQwCNww0ZaZSbVVv2aafucCjwXe1t0wJUlSk47q\ncP5+4EhgV0v7LmDVVB0i4gnAHwPPysyDEdHxICVJUrM6DQAdiYgjqHb7X5iZ35honm3/4eFh+vr6\nDmkbGhpiaGiouUFKkrRIjYyMMDIyckjb+Pj4rPp2GgDGgAPAipb2FcDdU8z/k8DPAU+LiPfXbUcA\nERH3AS/MzM+3K7Z+/XoGBgY6HKIkSWWYaqN4dHSUwcHBGft2dA5AZu4HtgJrJ9qi2qe/Frhxii73\nAj8LPA14av24HLil/vmLndSXJEnN6OYQwKXAxojYCtxEdVXAcmAjQERcDByfmefUJwj+x+TOEfEd\nYF9mbp/LwCVJUvc6DgCZuam+5v8iql3/24B1mbm7nmUlcEJzQ5QkSU3r6iTAzNwAbGgz7dwZ+r4N\nLweUJGlB+V0AkiQVyAAgSVKBDACSJBXIACBJUoEMAJIkFcgAIElSgQwAkiQVyAAgSVKBDACSJBXI\nACBJUoEMAJIkFcgAIElSgQwAkiQVyAAgSVKBDACSJBXIACBJUoEMAJIkFcgAIElSgQwAkiQVyAAg\nSVKBDACSJBXIACBJUoEMAJIkFcgAIElSgQwAkiQVyAAgSVKBDACSJBXIACBJUoEMAJIkFcgAIElS\ngQwAkiQVyAAgSVKBDACSJBWoqwAQERdExI6I2BsRWyLi1GnmfWZE/EtEjEXEnojYHhG/3f2QJUnS\nXB3VaYeIOBO4BHgtcBMwDGyOiCdm5tgUXX4IvBf4t/rnZwEfjIgfZOaHuh65JEnqWjd7AIaBKzLz\nqsy8BTgf2AOcN9XMmbktMz+emdszc2dmXgNsBp7d9aglSdKcdBQAIuJoYBC4YaItMxO4Hlgzy9c4\npZ73853UliRJzen0EEA/cCSwq6V9F7Bquo4R8S3g2Lr/WzPzIx3WliRJDen4HIA5eBbwMOAZwDsj\n4j8z8+PzWF+SJNU6DQBjwAFgRUv7CuDu6Tpm5jfrH78aESuBtwLTBoDh4WH6+voOaRsaGmJoaKiD\nIUuStDSNjIwwMjJySNv4+Pis+nYUADJzf0RsBdYC1wFERNTPL+vgpY4EjplppvXr1zMwMNDJECVJ\nKsZUG8Wjo6MMDg7O2LebQwCXAhvrIDBxGeByYCNARFwMHJ+Z59TPXwfsBG6p+z8XeBPwp13UliRJ\nDeg4AGTmpojoBy6i2vW/DViXmbvrWVYCJ0zqcgRwMXAScD/wDeB3M/ODcxi3JEmag65OAszMDcCG\nNtPObXn+PuB93dSRJEm94XcBSJJUIAOAJEkFMgBIklQgA4AkSQUyAEiSVCADgCRJBTIASJJUIAOA\nJEkFMgBIklQgA4AkSQUyAEiSVCADgCRJBTIASJJUIAOAJEkFMgBIklQgA4AkSQUyAEiSVCADgCRJ\nBTIASJJUIAOAJEkFMgBIklQgA4AkSQUyAEiSVCADgCRJBTIASJJUIAOAJEkFMgBIklQgA4AkSQUy\nAEiSVCADgCRJBTIASJJUIAOAJEkFMgBIklSgrgJARFwQETsiYm9EbImIU6eZ91cj4rMR8Z2IGI+I\nGyPihd0PWZIkzVXHASAizgQuAS4ETgFuBjZHRH+bLs8BPgucDgwAnwP+NiKe2tWIJUnSnHWzB2AY\nuCIzr8rMW4DzgT3AeVPNnJnDmfknmbk1M7+RmW8Bvg68uOtRS5KkOekoAETE0cAgcMNEW2YmcD2w\nZpavEcBPAvd0UluSJDWn0z0A/cCRwK6W9l3Aylm+xu8CDwU2dVhbkiQ15Kj5LBYRZwN/CJyRmWPz\nWVuSJD2g0wAwBhwAVrS0rwDunq5jRJwFfBB4aWZ+bjbFhoeH6evrO6RtaGiIoaGhWQ9YkqSlamRk\nhJGRkUPaxsfHZ9W3owCQmfsjYiuwFrgOfnxMfy1wWbt+ETEEfAg4MzM/M9t669evZ2BgoJMhSpJU\njKk2ikdHRxkcHJyxbzeHAC4FNtZB4CaqqwKWAxsBIuJi4PjMPKd+fnY97Q3AlyJiYu/B3sy8t4v6\nkiRpjjoOAJm5qb7m/yKqXf/bgHWZubueZSVwwqQur6E6cfD99WPClbS5dFCSJPVWVycBZuYGYEOb\naee2PP/FbmpIkqTe8bsAJEkqkAFAkqQCGQAkSSqQAUCSpAIZACRJKpABQJKkAhkAJEkqkAFAkqQC\nGQAkSSqQAUCSpAIZACRJKpABQJKkAhkAJEkqkAFAkqQCGQAkSSqQAUCSpAIZACRJKpABQJKkAhkA\nJEkqkAFAkqQCGQAkSSqQAUCSpAIZACRJKpABQJKkAhkAJEkqkAFAkqQCGQAkSSqQAUCSpAIZACRJ\nKpABQJKkAhkAJEkqkAFAkqQCGQAkSSqQAUCSpAJ1FQAi4oKI2BEReyNiS0ScOs28KyPioxFxa0Qc\niIhLux+uJElqQscBICLOBC4BLgROAW4GNkdEf5suxwDfAd4ObOtynJIkqUHd7AEYBq7IzKsy8xbg\nfGAPcN5UM2fmNzNzODOvBu7tfqiSJKkpHQWAiDgaGARumGjLzASuB9Y0OzRJktQrne4B6AeOBHa1\ntO8CVjYyIkmS1HNHLfQApjM8PExfX98hbUNDQwwNDS3QiCRJOnyMjIwwMjJySNv4+Pis+nYaAMaA\nA8CKlvYVwN0dvtaM1q9fz8DAQNMvK0nSkjDVRvHo6CiDg4Mz9u3oEEBm7ge2Amsn2iIi6uc3dvJa\nkiRp4XRzCOBSYGNEbAVuoroqYDmwESAiLgaOz8xzJjpExFOBAB4GHFs/vy8zt89t+JIkqRsdB4DM\n3FRf838R1a7/bcC6zNxdz7ISOKGl21eArH8eAM4Gvgmc3M2gJUnS3HR1EmBmbgA2tJl27hRt3nJY\nkqTDyGF9FYDKsHPnTsbGxrrq29/fz4knntjwiCRp6TMAaEHt3LmTVatWs2/fnq76L1u2nFtv3W4I\nkKQOGQC0oMbGxuoP/6uB1R323s6+fa9gbGzMACBJHTIA6DCxmur80PnhYQdJpTMAqDgedpAkA4AK\n5GEHSTIAqGjze9hBkg4nXp8vSVKBDACSJBXIACBJUoEMAJIkFciTAKUlzPsdSGrHACAtUd7vQNJ0\nDADSEuX9DiRNxwAgLXne70DSg3kSoCRJBXIPgKRGeeKhtDgYACQ1xhMPpcXDACCpMZ54KC0eBgBJ\nPeCJh9LhzpMAJUkqkAFAkqQCGQAkSSqQAUCSpAIZACRJKpBXAegQ3sRFkspgANCPeRMXSSqHAUA/\n5k1cJKkcBgBNwZu4SNJS50mAkiQVyAAgSVKBPARwGPOMfElSryy5ADAyMsLQ0NCir3l4npE/Aszv\n/+3C1Fyoutaci+kC82c+8xlOO+20Kaf1KiwvlXWRNZduza4OAUTEBRGxIyL2RsSWiDh1hvmfFxFb\nI2JfRHwtIs7pbrgzGxkZ6dVLz2vNQ8/I3zrF4zlt2rcCV7Nv356u9x60N///twtTc6HqWrNbE4F5\ncHBwysdb3vKWttNWrVrNzp07Gx/TUlkXWXPp1ux4D0BEnAlcArwWuAkYBjZHxBMz80GfOBFxEvAp\nYANwNvB84EMR8e3M/Pvuhz6/ptu6GB8fZ3R0tG3fuW1htDsjv69Nu1SemS9hHQbWT9E+t8tXu10v\neIhOh4NuDgEMA1dk5lUAEXE+8CLgPOBdU8z/W8Btmfl79fNbI+JZ9essigAwm93xg4ODbad5gxxp\nvsxfYJ7LesF1gg4HHQWAiDgaGAT+eKItMzMirgfWtOn2DOD6lrbNTB3HD0vdb12AN8iRlib3Omix\n63QPQD9wJLCrpX0XsKpNn5Vt5n94RByTmT+aos8ygO3bt0/5grt37277Brjjjjv46Ec/OvXg+/s5\n9thj2wyzvQfGsaPNHP8FTD3WiT7tlmXmmp9u89p3AFMvpzXnWnO6utY83GrOrm4va3a6Xui+5l13\n3cVLXvJS7rtvX9t52u11eMhDlvHXf/2XHHfccR3Xne91rjXnVnPS39ay6WpHZs56oBFxHHAnsCYz\nvzip/Z3AczLzQXsBIuJW4MOZ+c5JbadTnRewfKoAEBFn036tL0mSZvbyzLym3cRO9wCMAQeAFS3t\nK4C72/S5u83897bZ+ofqEMHLgduB9lFXkiS1WgacRPVZ2lZHASAz90fEVmAtcB1ARET9/LI23f4V\nOL2l7YV1e7s63wXaphZJkjStG2eaoZv7AFwKvCYiXhkRTwIuB5YDGwEi4uKIuHLS/JcDJ0fEOyNi\nVUS8Dnhp/TqSJGkBdHwZYGZuioh+4CKqXfnbgHWZubueZSVwwqT5b4+IF1GdDvsGqrNxXp2ZrVcG\nSJKkedLRSYCSJGlp8NsAJUkqkAFAkqQCGQAkSSqQAUBaQBFxTEQcs9DjmA+lLKvLqcViUQeAiHhy\nRGyIiK9ExF314yt125Otac3DtOYLIuLTEfE9YA+wJyK+V7c9vxc1pxjDvKy8F3pZXc7G6xSxnKVY\ntFcB1LcT/htglOpuRxPfN7ACeAHVlxb9cmZOeycka1pznmueA3wI+Mspar6Q6h4Zr87Mv2iq5qTa\nL6D6hpo1wMPr5nupbsp1adOX5i7UsrqcLmdDdZ8MvL6uu7Juvruu+77M/I9FXzMzF+UDuBm4aJrp\nbwX+zZrWPMxqfg24YJrprwO+3mTN+nXPAfYDI8CrqO7OeXr98zXAfcCvL/ZldTldzobqng78iOqD\n961UX2v/W/XPX6C6Rf26RV+z6f+4+XoAe4FV00xfBey1pjUPs5r75rtm/boLsfKe92V1OV3OhuoW\nsUGymM8BuB140TTTXwR805rWPMxqfhV49TTTzwMa37UInAhMt6v0BuAxDddciGV1OSsu59w8kem/\nkXYEeMJir9nxrYAPI/8HuCYinkf1BzL5eNRa4DTgbGta8zCr+SbgUxFxWpuaJzN9KOnWxMr799pM\n78XKeyGW1eWsuJxzczvVstzaZnovN0jmreaiPQkQICJ+ger7BaY6YeI9mdn2Gwetac0FrHkS1bG9\nZ0xR8/LMvL0HNZ8HfAq4jWlW3pn5Tw3XPYl5XFaX0+VsqO7/oDrH4O/a1D0NODsz/2pR11zMAUDS\n7C1E8FgILqfL2VDdJb9BYgCQJKlAi/kkwGlFxB9HxIetac1FVvPKiPiH+ay5UEpZVpdTh6slGwCA\nnwZOsqY1F1nNb9P8yUUzWqCV953M87K6nD1VynIumQ2SxXwVwLQy8xxrWnMR1vxf812z9m3g4HwW\nzMz/PZ/1avO2nBERWVmI5byTMn6f876ctZ8GTljsNT0HQJpnEdFPdflS64k+NwIbM3P3Qo1NzYmI\n+4CnZub2hR6LNJVFGwAi4jHAvswcq58/Gzif6sYR3wTe34uzNFvGEMDzgMcDdwGbM3N/L2tOqn0b\n1W0hv96j138U8BTg5sy8p/7QejVwDPCJXq3U6t/r9zPzBy3tRwNrmr7cZ4axrAB+MzMvavA1T6W6\nj/oepr7UZznV7/XLTdWc5bhOAN6Wmec1/Lo/QfWdCvdky33MI2IZ8LLMvKrhmqupzhj/18y8JSKe\nBLyR6m/36sxsdJdxRFzaZtIbgauB7wJk5u80WbdlDA8FXsYD66KRzPxuwzUGgO9l5o76+a9z6Dr3\nfZn5sYZrvhfYlJn/3OTrqrKYA8AXgbdn5qci4peBv6a6XnQ71R2V/jvwksz8VIM1Pw0MZeZ4RDwS\n+DTwdGAMeBTVbSuf0+QWXES8oc2kS4F3UW05kpmXNVjz6cBnqb544/tUX4zzCeB+qvNGjgeelZmj\nDdY8DriW6sMiqa6Hfd1EEKg/jL+dmUc2VXMWY3oqMNpkzYjYQnXLz/Oz5c1XB8rLgadk5pqmas5y\nXL1Y1idS/R2dSPU7/RfgrMy8q57e+O+0vknNtcAPqMLUrwJXUf2fHwE8F3hhkyEgIg7Wr//9lknP\nBb4M/BDIzPylBmv+B9V78J46vP0T8AiqddDjqe6f/4yJD+uGat4MvCkzr4+I3wAuA/6Map27CvgN\n4I2Z2dhx6vr/NoFvAH8OXJmZdzf1+jPUfgjwK0y9p+7azLyvR3XnbyOoyfsKz+eD6g3+2PrnLcCb\nW6a/nmqF1mTNg8Cj6583UN2lamIMj6F6s3+gBzW/BexoeRwE7qh/vq3hmn9P9cb+SeB/1vX/bNL0\nDwOfbLjmlfXv8eeA59f/l18CHlFPXwEcbLjmU2Z4vAw40HDNvcCTppn+JHrzXQBnzPD47R4s6yep\nQnk/1YfSxA1dTpz0O2265o3AH9U/nwXcA7xj0vSLgc82XPP36+X6pZb2/cCTm/5d1q89eV10NdWX\nxfTVzx9Wv4evabjmHuBn6p9Hgde0TD8b+GoPlnMt8KfAbqov/7mWagPviF7839Z1H08VOvYCnwc+\nXj8+X7d9HXh8wzWPA24CDlBtbF0FPGzS9ObfL736D+z1gyptP6X+edfEz5OmPw74YcM1J7/pbgHO\naJm+luY/jC8HvgKsbmnv5crlnol6wNH1H+TTJ00fAO5ouOadLTWOAa6rl/2RPfqwOFgv28EpHhPt\nTdfcAbxymumvBG7vwe90umU92KNl3QX8t0nPA/gA1e7ik3v0Ox2fWDFTbfHvB06ZNP1ngbt78P97\nKtUtXP8EOLpum68A8A3gBS3TfwHY2XDNMWBw0u/2qS3THwfs6eFyHk0Vyj9D9QF5J/AOGv4grmv9\nPdVXhT98imkPr6dtbrjmvG8ELebLAP8RGKp//grVsfjJfpHqD6RpE7ttH0H1xpvsP6l2jzdXLPN8\n4CJgc0S8vsnXnsZDqFIuWZ3TsIfqzT9h4pBHk/qA7008ycwfAS+huj/254BHN1wPqqDzGuCxUzxO\nptrKaNqfAB+MiPdExBkR8fP144yIeA9V4HtXD+reRXVI7IipHlShrmk/QbWiBqp94Jn5W8DfUr1/\nn9iDmlC/RzPzINU32I1PmvZfVH9rzRbM/BLV4atjgS9HxM/ywLqiVyZefxnV73eyO+uxNOnvqO7I\nB9Xv76Ut019GtQ7siczcn5mbMvM0qvfnnwEvp/298+fimcAfZOa9U4zjXuAPgWc3XPP5wBsy88uZ\neX09hruAf6gPOUPDf1OL+TLA3wf+OSKOpzq2+I76BKuJ41FnUp2g0rSNEfEjqjT6WKrDABNW8uDj\ngHOWmZ+MiJuAqyLiRcC5Tddo8S2qN9jt9fOzOHQFcxyHBoIm3Ea12/3HJzVm5v31/bE/QbX7uGlb\ngeMzc8prlyPip6i2WhuTme+PiDFgmOqrTCeOfx+ox/OqzNzUZM3aVqoPqGvbDY2Gl5VqL9nPUb0n\nHyiU+frqdAeua7geVH+zT+CBcL4G2Dlp+ok8+MOyEVkdsz0nIs6iOsGz1+er3BAR91Ntka4C/n3S\ntJ+hPvmwQW8GvhAR/0i1dfqm+l79E+vcZ1Cdc9FzmbkTeGtEvI3qg7Np36e6B8i/t5l+Es2v6x+0\nERQRL6Fa/30OeEXD9RZvAMjM7RHx88AfUX1T1EOp0uD9VLtNzsrMv2m47JWTfr6W6iSjyX4N2NZw\nTQAy886IeD5V8PkKza+sJ/sYk7a4M/P/tkw/g+pYVZP+DngtcMgXXUwKAX9F81/7eTnV3007O+lB\n2MrMjwMfr0/q6a+bx7K3V5C8m+mX9T+p9po16ZNUe+n+onVCHQKOoPmQ/gEmffBmZusK/HSgpzeO\nycyPRcS/UAWuXt0Y520tz3/Q8vzFQKNnzmfmtyPiFKp10Iup1kFPp7o2/QvAM7P5q1e+SRWO240p\nqXbXN+1DVBtcb6f6yuHWq3X+AHhvwzXnfSNo0V4FMFl99vSjqVLotql22/Sw9jOBL9dp7aFUxzT3\n9bjmc6heLR8LAAABS0lEQVSO8V2Rmd+baf6Gak5ezuVUy/mjBl//KGD55N9dS82jgJ9ut7Xe4Dh+\nXLOXdSQd3iLizVSXcq7kgV3vQXUlwJ9mZqOH6iLincDTMnPdFNOOotoIenF9yK6ZmkshAEyIiHup\n/gNvs6Y1F0tNSYeviHgsky4DzAYvrWyp86CNoCmmN7oRtJhPApxKL3eLW9OakgqTmTsy81/rx8RN\nkE5o+r78mXn/DHuvjwMubLLmUgsAkiT12iOB+f6ukMZrLtqTANv4TR44WcOa1lwsNSUdRiLijBlm\nOXlJ1FxK5wBIkjRXk25BPN0hwcxmb2M97zU9BCBJ0qEW4sZZ817TACBJ0qEmbpzVTi9unDXvNZfa\nOQCSJM3VQtw4a95reg6AJEkF8hCAJEkFMgBIklQgA4AkSQUyAEiSVCADgCRJBTIASJJUIAOAJEkF\n+v9va5P8ueNR6AAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7f7698010290>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "pmf.plot.bar()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "0.0625"
      ]
     },
     "execution_count": 5,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "pmf.mean()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "0.00996168582375479"
      ]
     },
     "execution_count": 6,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "pmf.median()"
   ]
  }
 ],
 "metadata": {
  "anaconda-cloud": {},
  "kernelspec": {
   "display_name": "Python [default]",
   "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
}