tSNE_PCA_code.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 33,
   "metadata": {},
   "outputs": [],
   "source": [
    "#Hi! Thanks so much for helping us out with this!\n",
    "\n",
    "#my main question is: I'm not 100% sure of the arguments passed in the last two cells\n",
    "#should I do something differently? I'm not sure if this gives the best results\n",
    "\n",
    "#I used this site as a source: https://dataplatform.cloud.ibm.com/analytics/notebooks/v2/3bee909a-5b69-4a50-9cd3-fdcd98fa5bd1/view?access_token=178b412e00a89ec0bca53f6fa65e9be58e392b23f684cb7377418ead7b8af4c2\n",
    "#The code from above link was used in this story: https://pudding.cool/2018/11/dearabby/\n",
    "\n",
    "import pandas as pd\n",
    "import numpy as np\n",
    "import matplotlib.pyplot as plt\n",
    "\n",
    "import re, os, time, json, csv\n",
    "\n",
    "import nltk as nlp\n",
    "from nltk.stem import WordNetLemmatizer\n",
    "\n",
    "from sklearn.decomposition import PCA\n",
    "from sklearn.feature_extraction.text import TfidfVectorizer\n",
    "from sklearn.decomposition import TruncatedSVD\n",
    "from sklearn.manifold import TSNE\n",
    "\n",
    "from gensim.models import word2vec"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "#load master speeches dataframe\n",
    "df=pd.read_excel(r'C:\\Users\\Vivek Rao\\Desktop\\IDS\\Digital\\McRobbie NLP\\speeches\\master_speeches.xlsx')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>date</th>\n",
       "      <th>venue</th>\n",
       "      <th>title</th>\n",
       "      <th>contents</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "Empty DataFrame\n",
       "Columns: [date, venue, title, contents]\n",
       "Index: []"
      ]
     },
     "execution_count": 3,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "#find missing values in the contents column\n",
    "df[df.contents.isna()==True]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [],
   "source": [
    "#write function to clean sentences\n",
    "#will be used as tokenizer in next cell\n",
    "df_cleaned = df.copy()\n",
    "\n",
    "wordnet_lemmatizer = WordNetLemmatizer()\n",
    "\n",
    "#write function\n",
    "def clean_text(text):\n",
    "    #lowercase text first\n",
    "    text=text.lower()\n",
    "    text=text.replace('\\n', ' ')\n",
    "    text=text.replace('\\t', ' ')\n",
    "    tokens = nltk.word_tokenize(text)\n",
    "    text=[word for word in tokens if not word in set(nlp.corpus.stopwords.words('english'))]\n",
    "    stems = []\n",
    "    for item in text:\n",
    "        stems.append(wordnet_lemmatizer.lemmatize(item))\n",
    "    return stems"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "C:\\Users\\Vivek Rao\\Anaconda3\\lib\\site-packages\\sklearn\\feature_extraction\\text.py:300: UserWarning: Your stop_words may be inconsistent with your preprocessing. Tokenizing the stop words generated tokens ['le', 'u'] not in stop_words.\n",
      "  'stop_words.' % sorted(inconsistent))\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Time elapsed: 1027.7006571292877 seconds\n"
     ]
    }
   ],
   "source": [
    "t0=time.time()\n",
    "\n",
    "#first step: set up the model by loading in the sentence cleaner.\n",
    "#the following step does not run the model\n",
    "tfidf = TfidfVectorizer(stop_words='english', tokenizer=clean_text, min_df=0.025, max_df=.5)\n",
    "\n",
    "#next, run the model on speeches, i.e. the 'Contents' column of the DataFrame\n",
    "tfs = tfidf.fit_transform(df.contents)\n",
    "print('Time elapsed: {} seconds'.format(time.time()-t0))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {
    "scrolled": true
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[[ 0.27870188  0.13561265]\n",
      " [ 0.49859581 -0.28868892]\n",
      " [ 0.22744671 -0.14617388]\n",
      " ...\n",
      " [ 0.20220505  0.10243758]\n",
      " [ 0.24201907  0.05842697]\n",
      " [ 0.33800855 -0.07205232]] 501\n"
     ]
    },
    {
     "data": {
      "text/plain": [
       "['Dedication of the Elinor Ostrom Statue and Ostrom Commons',\n",
       " 'Economic Club of Indiana ',\n",
       " 'Report to the Trustees of Indiana University\\n',\n",
       " 'Recommendations to the Board of Trustees on the removal of David Starr Jordan’s name from campus structures and spaces\\n',\n",
       " 'Introduction to Bicentennial Strategic Plan Final Reports',\n",
       " 'Presentation of President’s Medal to Jon Vickers',\n",
       " 'NASEM Voting Webinar: Will our votes be secure and counted in 2020?\\n',\n",
       " 'Board of Trustees Business Meeting Retirement Remarks',\n",
       " 'Statement to Board of Trustees on FY2020 External Grant Funding and Private Philanthropy\\n',\n",
       " 'Remarks to the Board of Trustees Academic Affairs Committee',\n",
       " 'APAN50 Virtual Meeting and Conference – Hong Kong Keynote Address',\n",
       " 'Surprise Farewell Reception for William J. Lowe, Chancellor of IU Northwest',\n",
       " 'Board of Trustees President’s Report',\n",
       " \"Statement to the IU Board of Trustees on IU's Response to the COVID-19 Global Pandemic\\n\",\n",
       " 'Presentation of Richard G. Lugar Award to Ambassador Marie Yovanovitch ',\n",
       " \"Introduction of Senator Todd Young's Richard G. Lugar Lecture and moderator John Stehr\",\n",
       " 'Stephen L. and Connie J. Ferguson International Center Groundbreaking',\n",
       " 'Civil Rights Heritage Center Celebration',\n",
       " 'Luddy AI Center Groundbreaking',\n",
       " 'Big Red 200: Accelerating Research and Discovery at Indiana University',\n",
       " 'Founders’ Day Bicentennial Dinner',\n",
       " 'Indiana University 200th Anniversary Lunch',\n",
       " \"Celebrating the Class of 2019: Sharing in Indiana University's Essential Work\",\n",
       " '114th Annual Indiana Society of Chicago Dinner',\n",
       " 'Honoring Distinguished Alumnus Robert Gates',\n",
       " 'Celebrating the Hamilton Lugar School of Global and International Studies',\n",
       " 'George Taliaferro Statue Dedication Ceremony',\n",
       " 'Excerpt from \"Honoring Indiana University\\'s Traditions of Excellence\"',\n",
       " 'Sidney and Lois Eskenazi School of Art, Architecture + Design naming ceremony ',\n",
       " 'Celebrating the transformative generosity of Sid and Lois Eskenazi',\n",
       " 'Announcing a landmark gift to the School of Informatics, Computing and Engineering',\n",
       " 'Introduction of Michael Kirby\\'s keynote lecture at the \"Navigating the Backlash Against Global Law and Institutions\" workshop',\n",
       " 'Unveiling and dedication of the Elinor Ostrom historical marker',\n",
       " 'Announcement of the Enhancement to the IU Retirement and Savings Plan',\n",
       " \"A new social contract for Indiana University's third century\",\n",
       " 'National Service Archive Panel Discussion',\n",
       " 'Naming Ceremony for the Paul H. O’Neill School of Public and Environmental Affairs',\n",
       " 'The Installation of Susan Elrod as Chancellor of IU South Bend',\n",
       " 'The State of IU Athletics at the Bicentennial',\n",
       " 'Introduction of Sergei Eisenstein’s \"Battleship Potemkin\"',\n",
       " 'Pfau Course Dedication',\n",
       " 'Statement to Board of Trustees on FY2019 External Grant Funding',\n",
       " 'IU Bicentennial Tribute to Elinor Ostrom',\n",
       " '2019-20 and 2020-21 Tuition and Fee Recommendations',\n",
       " 'Announcement of the Selection of the new Big Ten Commissioner',\n",
       " 'Closing Ceremony for IUPUI’s 50th Anniversary Celebration',\n",
       " 'IU Bloomington Undergraduate Commencement ',\n",
       " 'IU Bloomington Graduate Commencement',\n",
       " 'A Celebration Honoring Alvin Rosenfeld',\n",
       " 'George Taliaferro Celebration of Life',\n",
       " 'Innovation Hall: Profiting by the Past, Providing for the Future',\n",
       " 'American Australian Business Council Annual CEO Summit',\n",
       " 'Announcement of the naming of the Sidney and Lois Eskenazi School of Art, Architecture + Design\\n',\n",
       " 'Patrick O’Meara: Accomplished Scholar, Teacher, and Internationalist',\n",
       " 'Groundbreaking for the Relocation and Renovation of the Arthur R. Metz Bicentennial Grand Carillon\\n',\n",
       " '\"Contending with Antisemitism in a Rapidly Changing Political Climate: An International Scholars Conference\"\\n',\n",
       " 'Celebrating Richard Lugar, Indiana’s Visionary Statesman',\n",
       " 'Conferral of Honorary Doctorate to Chi Pang-yuan',\n",
       " 'Dedication of the Republic Building',\n",
       " 'IUPUI 50th Anniversary Report to the Community',\n",
       " 'Wilkinson Hall Dedication',\n",
       " 'Winter Commencement',\n",
       " 'IU Kokomo Student Activities and Events Center Groundbreaking',\n",
       " 'BDEC2 Welcome',\n",
       " 'Celebration of J. Irwin Miller Architecture Program',\n",
       " 'Introduction of Benjamin Britten\\'s \"War Requiem\"',\n",
       " 'Office of International Services 75th Anniversary Luncheon',\n",
       " 'IU Hillel Memorial Service and Vigil for the Victims of the Tree of Life Synagogue Shooting\\n',\n",
       " 'Re-Dedication of Ernie Pyle Hall',\n",
       " 'Conferral of Honorary IU Doctorate to Kenneth \"Babyface\" Edmonds',\n",
       " 'Dedication of Ray E. Cramer Marching Hundred Hall',\n",
       " \"Honoring Two of America's Greatest Statesmen and Internationalists\",\n",
       " 'Office of International Services 75th Anniversary Celebration',\n",
       " 'Excellence Academy Dedication',\n",
       " 'Introduction of Michelangelo Antonioni’s \"Blow-Up\"',\n",
       " 'Reception to Celebrate 20th Anniversary of the IU International Networks Program',\n",
       " 'Indiana University Police Academy Graduation Ceremony',\n",
       " 'Stone Family Center for Health Sciences Dedication',\n",
       " 'David Zaret: Celebrating 41 Years of Outstanding Service to Indiana University',\n",
       " 'Air Force Research Lab Workshop Welcome at Science and Technology 2030 Initiative',\n",
       " 'IU Bloomington Undergraduate Commencement',\n",
       " 'Indiana University Bloomington Graduate Commencement',\n",
       " \"Special Event and Announcements in Celebration of IU's Master of Architecture Program\\n\",\n",
       " 'In Remembrance of Ken Gros Louis',\n",
       " 'Remembering ANZAC Day',\n",
       " 'The Inauguration of Ora Hirsch Pescovitz as the Seventh President of Oakland University\\n',\n",
       " 'Dedication of Luddy Hall, Remarks of Michael A. McRobbie\\n',\n",
       " 'Honoring Sir Salman Rushdie: Maker of Fiction, Teller of Truth',\n",
       " 'Introduction to \"Indiana in the World\" Panel \\n',\n",
       " 'Remembrance Event in Honor of Bruce Cole\\n\\n',\n",
       " 'Kelley School of Business, Prebys Career Services Center Dedication\\n\\n',\n",
       " 'James J. Fritts, D.D.S. Clinical Care Center dedication',\n",
       " 'IU Health Bloomington Regional Academic Health Center Groundbreaking',\n",
       " 'Monroe County Bicentennial Opening Celebration',\n",
       " 'Winter Commencement 2017: Diffusing light and truth at Indiana University',\n",
       " \"Introduction of Alain Resnais’ 'Last Year at Marienbad'\",\n",
       " 'Honoring Poland’s Greatest Living Composer: Krzysztof Penderecki',\n",
       " 'Wilkinson Hall Groundbreaking',\n",
       " 'Visions of Development Conference Welcome',\n",
       " 'The For All Bicentennial Campaign: \\nCelebrating Progress, Setting Bold Goals for Indiana University’s Future\\n',\n",
       " 'Induction as Honorary Member of I-Association',\n",
       " \"Deploying IU's Resources to Address a Deadly Public Health Crisis\",\n",
       " 'Living and Learning in the Heart of Campus: The Renovated Agnes E. Wells Quad Student Residences',\n",
       " 'Walking into the World of the Future in Northwest Indiana',\n",
       " 'From Bedford to Bloomington and Beyond: Building for the Future of South Central Indiana',\n",
       " 'From the Hills of Ancient Greece to the Hills of Southern Indiana',\n",
       " 'Reunion of the Wisconsin-Indiana-Purdue (WIP) Study Abroad Program in Madrid',\n",
       " '50th Anniversary of Madrid Study Abroad Program\\n',\n",
       " 'Announcing the Second Grand Challenges Research Initiative',\n",
       " 'The Class of 2017: Affirming Indiana University’s Commitment to the Continual Search for Truth',\n",
       " 'Celebrating the Graduate Class of 2017: Defenders of Truth in a Post-Truth Era',\n",
       " 'Celebration of the School of Art + Design and Rededication of Kirkwood Hall',\n",
       " 'Indiana’s Role in the World',\n",
       " 'Dedication of Paul H. O’Neill SPEA Graduate Center',\n",
       " 'Presentation of General Douglas MacArthur Award to IU Army ROTC',\n",
       " 'From Canberra to Bloomington: perspectives on higher education in Australia and the United States',\n",
       " 'Presentation of the President’s Medal for Excellence to John Lechleiter',\n",
       " 'Celebrating the life of David Baker',\n",
       " 'The university and civil society in the 21st century',\n",
       " 'Announcement of Don Brown’s gift to establish the Brown Center for Immunotherapy',\n",
       " 'The role of the world’s leading universities in the preservation of knowledge in the digital age',\n",
       " 'Opening remarks for the Ostrom Symposium on the study of the commons, governance, and collective decision',\n",
       " 'News conference to announce location of TASUS corporate headquarters in Bloomington Trades District',\n",
       " 'Kelley School of Business: Prebys Career Services addition groundbreaking',\n",
       " 'Introduction of Brian Schmidt’s lecture: “The Accelerating Universe”',\n",
       " 'Celebration and naming announcement for the Indiana University School of Nursing Champion Center for Cancer Control Research',\n",
       " 'Groundbreaking for the renovation of Riverside Hall',\n",
       " 'A life of scholarship and service: honoring Steven C. Beering',\n",
       " 'Dedication of the renovation and naming of Simon Skjodt Assembly Hall',\n",
       " 'Introduction of Peter Weir’s The Year of Living Dangerously',\n",
       " 'Sidney and Lois Eskenazi Museum of Art 75th anniversary gala',\n",
       " 'Introduction of Prime Minister Rui Maria de Araújo',\n",
       " 'Franklin Hall rededication',\n",
       " '2016 diversity reception at the 8th annual IUPUI Regatta',\n",
       " 'Groundbreaking for the IU School of Dentistry',\n",
       " 'The international university in the 21st century',\n",
       " 'Inauguration of Institute for Korean Studies',\n",
       " 'North Hall dedication ceremony',\n",
       " 'IU Athletics all-staff meeting',\n",
       " 'Engagement, collaboration, and innovation for the future of Columbus, South Central Indiana, and beyond',\n",
       " 'IU East Student Events and Activities Center dedication',\n",
       " 'Rededication of Ball Nurses’ sunken garden',\n",
       " 'Announcing the inaugural Grand Challenges research initiative',\n",
       " 'IU Natatorium grand re-opening',\n",
       " 'Celebrating a 40-year partnership',\n",
       " '50th anniversary celebration of IU Bologna study abroad program',\n",
       " 'Announcement of the partnership between the Uffizi Gallery and Indiana University',\n",
       " 'The Indiana University class of 2016: poised to enliven and strengthen our world',\n",
       " 'Profiting by the past: celebrating the renovation of the IU Kokomo main building',\n",
       " 'Indiana University’s Grand Challenges research program and its potential benefits to North Central Indiana, the State, and the Nation',\n",
       " 'Honoring Indiana University’s 2016 Founders Scholars and Meuller Scholars',\n",
       " 'IU Auditorium 75th anniversary gala',\n",
       " 'Universities and global sustainable economic development in the 21st century',\n",
       " 'Beyond the classroom walls: international engagement and online education',\n",
       " 'Planning for excellence in Indiana University’s third century: strategic iInitiatives and opportunities for IU, Central Indiana, and the State',\n",
       " 'Paul H. O’Neill SPEA Graduate Center groundbreaking',\n",
       " 'Local talent, global reach: engagement, collaboration, and innovation for the future of Wayne County, East-Central Indiana, Western Ohio and beyond',\n",
       " 'Bloomington economic development corporation annual meeting',\n",
       " 'Martin Luther King Jr. Day celebration leadership breakfast',\n",
       " 'Celebrating the Next Generation of Active and Engaged Global Citizens',\n",
       " 'Honoring Alice Rivlin: A Distinguished Career in Public Service',\n",
       " 'Regenstrief Institute Open House',\n",
       " 'Introducing Bernardo Bertolucci’s The Last Emperor',\n",
       " 'The Installation of Nasser Paydar as Chancellor of IUPUI',\n",
       " 'APLU Council of Presidents Panel on Student Debt Reduction',\n",
       " 'Welcome to Annual Meeting of the Association of Public and Land Grant Universities',\n",
       " 'IU Northwest Chancellor’s Medallion VIP Reception',\n",
       " 'IU Northwest Chancellor’s Medallion Dinner',\n",
       " 'Importance of Institutional Collaboration to Advance International Education',\n",
       " 'The Power of an International Education',\n",
       " 'The Data Dilemma: A Transatlantic Discussion on Privacy, Security, Innovation, Trade, and the Protection of Personal Data in the 21st Century',\n",
       " 'Honoring Twyla Tharp',\n",
       " 'Groundbreaking for Evansville Multi-Institutional Health Science Education and Research Campus',\n",
       " 'Official Opening of the Media Digitization and Preservation Initiative',\n",
       " 'Introduction of United States Secretary of State John Kerry',\n",
       " 'Dedication of the Global and International Studies Building',\n",
       " 'University Hall Dedication',\n",
       " \"Toward the Bicentennial: Transformational Progress for Indiana University's Third Century\",\n",
       " 'Farewell Reception for Bobby Schnabel',\n",
       " 'Luddy Hall Groundbreaking',\n",
       " 'For All: Indiana University Bicentennial Kickoff',\n",
       " 'Posthumous Presentation of Thomas Hart Benton Medallion to Michel duCille',\n",
       " 'Remarks at IU Health Announcement of Impending Retirement of Dan Evans as President and CEO and Appointment of Dennis Murphy as his Successor',\n",
       " 'IU Athletics All-Staff Meeting',\n",
       " 'Presentation of the IU Lilly Family School of Philanthropy Stead Medal in International Philanthropy Award Ceremony',\n",
       " 'Indiana Black Expo Corporate Luncheon',\n",
       " 'National Summit on Collegiate Financial Wellness',\n",
       " 'Reception in Honor of Vice President for Research Jorge José',\n",
       " 'An Evening of Celebration with IUPUI Chancellor Charles Bantz and Professor Sandra Petronio',\n",
       " 'U.S.-China Presidents Roundtable',\n",
       " 'United in Common Purpose: Breaking Ground on a Partnership for Excellence',\n",
       " 'Announcement of Mark Cuban Center for Sports Media and Technology',\n",
       " 'Tribute to Chancellor Charles Bantz',\n",
       " 'Leading for Change: Innovation and Research',\n",
       " 'Cream and Crimson Awards Luncheon',\n",
       " 'The Value of Higher Education: Benefitting Graduates, Advancing the Common Good',\n",
       " 'The Fruits of Graduate Education in the 21st Century',\n",
       " 'North Hall Groundbreaking',\n",
       " 'Orson Welles: A Centennial Celebration and Symposium',\n",
       " 'Remembering ANZAC Day: 100 Years Later',\n",
       " 'Australian Intervention in the Great War: The Definition of a Nation',\n",
       " 'A Legacy of Leadership',\n",
       " 'News Conference to Announce New Location of IU Health Bloomington Hospital and Establishment of Academic Health Center',\n",
       " 'Celebrating Gathering Spaces that Reflect the Values of Indiana University',\n",
       " 'Groundbreaking for the Renovation of Assembly Hall',\n",
       " 'Australian National University 2015 Alumnus of the Year Award',\n",
       " 'International Conference on China in the Middle East',\n",
       " 'World War I: 100 Years Removed',\n",
       " 'Celebrating the Legacy of David Starr Jordan, Eminent Educator, Philosopher, and Scientist',\n",
       " 'Acceptance of the U.S. Midwest Award for Outstanding Achievement',\n",
       " 'The Transformative Power of the Arts',\n",
       " 'The Installation of Ray Wallace as Chancellor of IU Southeast',\n",
       " 'The IU East Student Activities and Events Center',\n",
       " \"Strengthening Indiana University's Longstanding Connections with India\",\n",
       " 'The Preservation of Knowledge and Culture in India, Indiana, and Beyond',\n",
       " 'Honoring George Ortman',\n",
       " 'Building on Time-Honored Strengths to Meet the Emerging Challenges of the Future',\n",
       " 'Toward Enduring Excellence: The Bicentennial Strategic Plan for Indiana University',\n",
       " 'Improving the Health of the People of Indiana and Beyond \\nThrough Neurosciences Research',\n",
       " 'Hodge Hall: Improved Equipment in the Service of Learning',\n",
       " 'Honoring \\x91the American Olivier:\\x92 IU Alumnus Kevin Kline',\n",
       " 'Acceptance of the IU Hillel Chaver Award',\n",
       " 'A Celebration of Leadership',\n",
       " 'Announcement of a Gift from Irene Lilly McCutchen to the IU Lilly Family School of Philanthropy',\n",
       " 'Dedication of the Renovation of the Rotary Building',\n",
       " 'Center for Applied Cybersecurity Research Reception',\n",
       " 'Education for a World of Change',\n",
       " 'Celebrating a Legacy of Leadership in Education, Research, and Practice',\n",
       " 'Honoring A Medical Pioneer: Distinguished Professor of Medicine Harvey Feigenbaum',\n",
       " 'Singapore Alumni Reception and Dinner',\n",
       " 'Opening Ceremony for the IU China Gateway Office',\n",
       " 'Tokyo Alumni Reception and Presentation of Benton Medallion to Tsuyoshi Tsutsumi',\n",
       " 'Democracy, Education, and the Freedom to Know',\n",
       " 'IU Honors His Excellency Michael D. Higgins, President of Ireland and Former U.S. Secretary of the Treasury Paul O’Neill',\n",
       " 'Graduate Education and the Renewal of Democracy in the 21st Century',\n",
       " 'Rotary Club of Evansville',\n",
       " 'Toleration, Understanding and the Duty of a University to Defend its Fundamental Values',\n",
       " 'An Unparalleled Artist of Stage and Screen: Celebrating Meryl Streep',\n",
       " 'Decadence, Detail, and Lavish Beauty: the Films of Luchino Visconti',\n",
       " 'Recommendation to the Board of Trustees on the Location of the Evansville Health Sciences Campus',\n",
       " 'Remarks at the opening of Internet2/NKN Partnership Discussion',\n",
       " 'Profiting by the Past: Celebrating the Renovations of Jordan and Owen Halls',\n",
       " 'Celebrating the Refabricated George Ortman Banners',\n",
       " 'New Space for Two of Indiana University’s Pioneering Schools',\n",
       " 'Celebrating the Vision and Energy of the Rising Generation',\n",
       " 'Outstanding Leadership and Service in Higher Education',\n",
       " 'Celebrating the Transformative Generosity of Cindy Simon Skjodt',\n",
       " 'The Rise and Fall of the Third Reich on Film',\n",
       " 'Education, Knowledge, and Innovation in Science and Engineering',\n",
       " 'Celebrating Jane Fortune',\n",
       " 'Maintaining a Noble Course in the Musical Arts',\n",
       " 'Staff Appreciation Celebration',\n",
       " 'Building Friendships, Brightening Futures',\n",
       " 'The Installation of Kathryn Cruz-Uribe as Chancellor IU East',\n",
       " 'Psychological and Brain Sciences’ 125th Anniversary Celebration: 1888-2013',\n",
       " 'The Installation of Terry Allison as Chancellor of IU South Bend',\n",
       " 'Violette Verdy: Celebrating One of the Great Ballerinas of the 20th Century',\n",
       " '2013 Indiana University Intercollegiate Athletics Hall of Fame Induction',\n",
       " 'Looking to the Future: Preparing for Indiana University’s Bicentenary',\n",
       " 'A Great History and a Brilliant Future: IU’s Distinguished and Titled Faculty',\n",
       " 'Celebrating the Impact of IUPUI\\x92s Generous Supporters',\n",
       " 'Inauguration of the School of Informatics and Computing',\n",
       " 'Orphans Midwest Film Symposium',\n",
       " 'King of Cellists: Celebrating the Life of Janos Starker',\n",
       " 'Improving Health and Well-Being in North Central Indiana',\n",
       " 'Honoring Ed Marshall’s Legacy of Outstanding Leadership',\n",
       " 'Excerpt from Welcome School of Medicine Dean Jay Hess',\n",
       " 'Celebrating the GlobalNOC’s 15th Anniversary',\n",
       " 'Anchors of Progress and Freedom: Honoring the USS Indiana and her Crew',\n",
       " 'Hine Hall: A Monument to a Visionary Leader',\n",
       " 'The Importance of Political Networks',\n",
       " 'Dinner in Celebration of Dean Craig Brater',\n",
       " 'Growth and Change in Allen County',\n",
       " 'Tribute to Keith Cash',\n",
       " 'Honoring Jeanne Madison: Presentation of the E. Ross Bartley Award',\n",
       " 'Retirement Reception for Chancellor Sandra Patterson-Randles',\n",
       " 'Rededication of the Renovated Alumni Hall',\n",
       " '2013-14 and 2014-15 Tuition and Fee Recommendations',\n",
       " 'Excerpt from Excellence, Perseverance, and Dignity  Through Adversity: Honoring Ju Tong-Kui',\n",
       " 'The Promises and Challenges of an Era of Change',\n",
       " 'Indiana University Honors David Brooks, Alecia DeCoudreaux, and Vi Simpson',\n",
       " 'The Indiana University Class of 2013: Poised to Enliven and Strengthen Our World',\n",
       " 'Widening Horizons: Celebrating the School of Global and International Studies',\n",
       " 'Coming Home: Celebrating Bart Kaufman Field',\n",
       " 'Big Red II: Accelerating Research and Discovery at Indiana University',\n",
       " 'Menahem Pressler and Friends',\n",
       " 'Cultivating the Talents of the Sons and Daughters of Indiana: The IU South Bend Education and Arts Building',\n",
       " 'Presidents Hall Celebration',\n",
       " 'Pioneering Philanthropy Education, Research, and Training',\n",
       " 'Preparing to Address Society\\x92s Paramount Challenges: 2013 Founders Scholars',\n",
       " 'A Sound and Lasting International Partnership',\n",
       " 'Higher Education at the Vanguard of a New International Era',\n",
       " 'American Higher Education in the 21st Century',\n",
       " 'Playing to the Camera: Great Concert Films of the Classic Rock Era',\n",
       " 'Cooperation and Collaboration: A Vision for Progress in Northeast Indiana',\n",
       " 'Strengthening Indiana University’s Ties  with Two of America’s Greatest Statesmen and Internationalists',\n",
       " 'High Performance Networks in Higher Education: Past, Present, and Future',\n",
       " 'The Creation of the Future at Indiana University',\n",
       " 'Farewell Reception for Neil Theobald and Sheona Mackenzie',\n",
       " 'Educating Generations of Future Leaders',\n",
       " 'The Mind and Soul of the University: Celebrating the Lives of \\nElinor and Vincent Ostrom',\n",
       " 'Transforming Indiana University for the 21st Century',\n",
       " 'Promoting Improved Public Health Through Education, Research, and Service',\n",
       " 'Transforming Public Health in the 21st Century',\n",
       " 'Naming Ceremony for the IU South Bend Judd Leighton School of Business and Economics',\n",
       " 'An ‘Impassioned and Adventurous Career’: Celebrating Filmmaker Werner Herzog',\n",
       " 'IU Health Neuroscience Center Dedication',\n",
       " 'Freshman Induction',\n",
       " 'Unlocking the Secrets of the Brain: Breaking Ground for the Neurosciences Research Building',\n",
       " 'Affirming the American Dream',\n",
       " 'Dottie Frapwell: Celebrating 35 Years of Extraordinary Service to Indiana University',\n",
       " 'Excerpt from Honoring Professor H.A.R. Alex Tilaar',\n",
       " 'Global Impact: Indiana University and Institution-Building',\n",
       " 'The Diffusion of Light and Education: Serving the Public Good at Indiana University',\n",
       " 'The Public Mission of the Research University: Serving the Public Good in the 21st Century',\n",
       " 'A Founder of American Soul Music: Indiana University Honors Booker T. Jones',\n",
       " 'Innovation, Transformation, Community',\n",
       " 'Opening of IUPUI Campus Center Theatre and The History of Cardenio',\n",
       " 'Simic Commons Dedication',\n",
       " 'The Legacy of Herman B Wells',\n",
       " 'Honorary Degree Conferral: Yi Gang',\n",
       " 'Ambassador Nirupama Rao: Poet, Diplomat, and Trailblazer',\n",
       " 'Presentation of Distinguished Service Medal to Professor Michael Schwartzkopf',\n",
       " 'James R. Hodge Hall Naming Ceremony',\n",
       " 'Groundbreaking Ceremony for the Science and Engineering Laboratory Building',\n",
       " 'David Baker: Celebrating a Jazz Legend',\n",
       " 'A Legacy of Academic Leadership: Celebrating the Service of Provost Karen Hanson',\n",
       " 'The Lilly Endowment:  Astonishing Generosity, Inestimable Impact',\n",
       " 'The Value of Knowledge in an Age of Change',\n",
       " 'Dedicated to Justice, Integrity, and Civility: Indiana Chief Justice Randall Shepard',\n",
       " 'Excerpt from Honoring a Remarkable Leader: Former State Representative Bill Cochran',\n",
       " 'Extending the Tradition of Leadership and Vision: The Indiana University \\nRobert H. McKinney School of Law',\n",
       " 'Press Conference and Restitution Ceremony',\n",
       " 'A Beacon to Europe and the World: Celebrating Ten Years of Progress at South East European University',\n",
       " 'The Past, Present, and Future of Software Serving the Academy',\n",
       " 'A Natural Teacher and IU’s Ambassador to the World',\n",
       " 'Andrei Tarkovsky and the Nuance of the Human Condition',\n",
       " 'The Fabric of a Great American City: Art and Design in Columbus, Indiana',\n",
       " 'A Great Person with a Great Mind and a Gentle Soul: Remembering IU President Emeritus John Ryan',\n",
       " 'A Century of Partnership and Progress: The IU School of Social Work',\n",
       " 'Making a World-Class Education More Affordable',\n",
       " 'The Cyberinfrastructure Building: Continuing the Journey',\n",
       " 'Economic Engagement and the Research Enterprise: \\nThe Role of American Higher Education',\n",
       " 'Climbing Our Way to the Top: Dedication of the HPER Courtyard Project',\n",
       " 'Excellence in an Evolving Environment',\n",
       " 'Defining the 21st Century Public University',\n",
       " 'Not for the Glory but for the Story:  \\nOne Hundred Years of Change at the IU School of Journalism',\n",
       " 'A Collective Act of Memory: Remembering 9/11',\n",
       " 'Getting in the Game:  IU Student-Athletes and Academic Excellence',\n",
       " 'Confederation of Indian Industry',\n",
       " 'Strengthening Global Partnerships and Honoring Distinguished Alumni: IU Delegation Visits OP Jindal Global University',\n",
       " 'Proud Traditions, Limitless Opportunity, and Liberal Education at Indiana University',\n",
       " 'A New Language of Hope: Dedicating the Glick Eye Institute',\n",
       " 'IU Health Neuroscience Center (Phase I) Topping Out Ceremony',\n",
       " 'Paying Tribute to an Indiana Treasure:  Honoring the Lilly Endowment',\n",
       " 'A Life of Integrity, Scholarship, and Service: Honoring Bruce Jaffee',\n",
       " 'Honoring Marilyn Glick: A Visionary in the Community',\n",
       " 'Expanding for Campus and Community: IU South Bend Education and Arts Building',\n",
       " 'Mini University: Building a Community of Scholars',\n",
       " 'A Longstanding Commitment to the Arts and Humanities: The New West Wing of the Lee Norvelle Theatre and Drama Center',\n",
       " 'Opportunity for Columbus and Beyond: The Dedication of the Advanced Manufacturing Center of Excellence',\n",
       " 'Bridging the Gap Between Theory and Practice: 2011 Public Economic Theory Conference',\n",
       " 'A Visionary and a Humanitarian: Celebrating the Life of Bill Cook',\n",
       " '2011-12 and 2012-13 Tuition and Fee Recommendations',\n",
       " 'Honorary Degree Luncheon Honoring Bill Hunt',\n",
       " 'Inheriting the Power of a Liberal Education and Investing in the Nation’s Future',\n",
       " 'Possibilities and Progress Through Higher Education: William G. Bowen, Anthony J. Adams, and Patricia R. Miller',\n",
       " 'Investing in Excellence and Imagining the Next Horizon',\n",
       " 'Building on a Legacy of Teaching and Research:\\nThe Rededication of Integrated Science and Accelerator Technology Hall',\n",
       " 'Building a Scholarly Community',\n",
       " 'To Swell the Soul with Music',\n",
       " 'An Unmatched Choice: Ambassador Kim Beazley',\n",
       " 'Strengthening the Foundation of Tolerance at Indiana University',\n",
       " 'A History of Artistic Collaboration: \\nStanford, Indiana University, and Domenico Tiepolo',\n",
       " 'A Model of Collaboration in the Search for a Cure: The Dedication of Harper Hall',\n",
       " 'The Difference the Big Screen Makes: The Dedication of the IU Cinema',\n",
       " 'Printing the Legend and Honoring Peter Bogdanovich',\n",
       " 'IU Health Dedication Ceremony and News Conference',\n",
       " 'A Legacy of Leadership and Achievement: Kappa Alpha Psi Centennial Celebration',\n",
       " '“Dreaming of New Worlds and Opening Doors of Innovation”',\n",
       " '“Worldwide Service, Commitment, and Compassion: \\nHonoring Her Royal Highness Princess Maha Chakri Sirindhorn and Vice Admiral Adam M. Robinson”',\n",
       " '\"IU Center for Art and Design-Columbus Press Conference\"',\n",
       " '\"Women of Distinction at Indiana University: The Colloquium for Women of Indiana University and Meryl Streep\"',\n",
       " '“Carrying and Transmitting the Essence of Our Humanity: The Lilly Library Celebrates Fifty Years”',\n",
       " '\"Shared Dreams and Aspirations: The Installation of William Lowe as Chancellor of IU Northwest\"',\n",
       " '\"Residential Housing in this New Era of Imagination\"',\n",
       " '\"The Global Power of Collaboration:  \\nAnnouncing the IU Center for Global Health\"',\n",
       " '\"Imagine What the Future Holds: The Installation of Michael Harris as Chancellor of IU Kokomo\"',\n",
       " '\"Impact IUPUI:  The Celebration Continues\"',\n",
       " '\"Galvanizing Support and Celebrating Matching the Promise\"',\n",
       " '“Principles of Excellence”',\n",
       " '\"Introducing a Distinguished Public Servant and Hall of Fame Wrestler: NEH Chairman James A. Leach\"',\n",
       " '\"Recommendation for Faculty and Staff Salary Increase\"',\n",
       " '\"Statement for Academic Affairs Committee Regarding Research Numbers\"',\n",
       " '\"International Computing Networks at IU Receive NSF Awards\"',\n",
       " '\"Let Your Minds Carry You Towards the Future\"',\n",
       " '“‘It’s All About Relationships’: Interim Chancellor Stuart Green”',\n",
       " \"“The People Say 'We Did It Ourselves': \\nCelebrating the Chancellorship of Bruce W. Bergland”\",\n",
       " '“National Institutes of Health Press Conference”',\n",
       " '\"Continuing a Partnership Toward Better Health\"',\n",
       " '\"The Currency of Knowledge at Indiana University and Beyond\"',\n",
       " '\"A Legacy of Learning\"',\n",
       " '\"A Musical Innovator, Pioneer, and Global Humanitarian: Indiana University Honors Quincy Jones\"',\n",
       " '\"Breaking Ground and Continuing Leadership in Information Technology at Indiana University\"',\n",
       " '\"Celebrating Cook Hall and the People Our Players Become\"',\n",
       " '\"Ambassador Nigel Sheinwald: Serving the United Kingdom and the Greater Good\"',\n",
       " '\"The Best of What We Know: 2010 Founders Scholars\"',\n",
       " '\"IU  Bloomington Libraries:\\xa0 Leading the Way  to the Future\"',\n",
       " '\"Indiana University: A Vital Partner in the Regional Economy\"',\n",
       " '\"Sustaining Educational and Research Excellence in Times of Challenge: The Central Importance of Indiana University\\'s Academic and Facility Core\"',\n",
       " '\"Honoring Excellence at Indiana University: Professor Elinor Ostrom Delivers Nobel Lecture\"',\n",
       " '\"Clarian Health-Bloomington Hospital Press Conference\"',\n",
       " '\"Vital Skills in a Changing World\"',\n",
       " '\"Incubating Excellence at Indiana University\"',\n",
       " '\"Securing IU\\'s Information Technology Future\"',\n",
       " '\"A Laboratory of Leadership: Celebrating 100 Years of the Indiana Memorial Union Board\"',\n",
       " '\"Community Partnership for Sound Minds\\nand Sound Bodies\"',\n",
       " '\"A Transformative Leader of Wisdom, Dignity, and Courage:\\nIn Honor and Memory of Myles Brand\"',\n",
       " '\"A Tribute to Leadership\"',\n",
       " '\"The Unfolding of Great Possibility:\\nDedicating a New Space for Science at Indiana University\"',\n",
       " '\"Building on a Glorious IU Tradition: Breaking New Ground in the Arts and Humanities\"',\n",
       " '\"Honoring a Great Benefactor: Tom Raper\\'s Support of IU East\"',\n",
       " '\"The Natural History of Creation: An Evolutionary Snapshot\"',\n",
       " '\"A Voice for Democracy: The Center on Congress at Indiana University\"',\n",
       " '\"Celebrating Unsung Heroes of the Clinic and the Lab\"',\n",
       " '\"The Spirit of Generosity: The Maurer School of Law\"',\n",
       " '\"Liberal Education and the Beginning of Success at Indiana University\"',\n",
       " '\"Research for a Thousand Years: Sustainability, Clean Energy, and Our Shared Future\"',\n",
       " '\"Indiana University 2009-10 Budget Proposal to the Trustees of Indiana University\"',\n",
       " '\"A Beautiful Vision for IU East: The Leadership of Nasser Paydar\"',\n",
       " '\"The Lessons of History: Guide towards a Future of Renewal and Innovation\"',\n",
       " '\"A Humanitarian for the Times: Indiana University Honors Justice Michael Kirby\"',\n",
       " '\"Strategic Planning and Innovation: International Education and Research at Indiana University\"',\n",
       " '\"Celebrating the Spirit of Ed Hutton and a New Intellectual Home\"',\n",
       " '\"A Man of Action and Integrity: IU President Emeritus Myles Brand\"',\n",
       " '\"The Wisdom and Courage to Change the World: 2009 Founders Scholars\"',\n",
       " '\"The Fire of the Imagination: Master Planning and Student Life at Indiana University\"',\n",
       " '\"Political Art and Artful Politics: ArtsWeek 2009\"',\n",
       " '\"The Evolving Partnership between Ivy Tech Community College and Indiana University\"',\n",
       " '\"Raising the Bar for Matching the Promise: Heightened Goals in Challenging Times\"',\n",
       " '\"IUPUI’s Dynamic Culture of Collaboration and Partnership\"',\n",
       " '\"\\'A Generation of Change\\': Herbert Presidential Scholars\"',\n",
       " '\"Atkinson and Seffrin: Exemplars of Excellence at Indiana University\"',\n",
       " '\"The Michael Maurer School of Law at Indiana University\"',\n",
       " '\"A Presidential Vision for an Internationalized Institution\"',\n",
       " '\"Cultivating a Spirit of Collaboration\"',\n",
       " '\"A Century of Progress: IU School of Education\"',\n",
       " '\"Searching for Truth on the Global Horizon: Academic Freedom and International Education\"',\n",
       " '\"Announcement of New IU Athletic Director\"',\n",
       " '\"A Great History and an Auspicious Beginning: Group 43 (1968) to Group 43 (2007)\"',\n",
       " '\"A Campus Home at Indiana University Southeast\"',\n",
       " '\"Building Excellence at Indiana University: A Year of Progress\"',\n",
       " '\"A New Era of Visual Health: Dedicating the School of Optometry Eye Care Center\"',\n",
       " '\"The Glowing Promise of the Future: Indiana University and the Eugene and Marilyn Glick Eye Institute\"',\n",
       " '\"Uniting Generations and Cultures: A Tribute to Hoagy Carmichael\"',\n",
       " '\"Finding a Place at Indiana University\"',\n",
       " '\"In the Service of Mankind: The IU Melvin and Bren Simon Cancer Center\"',\n",
       " '\"A Farewell Tribute to Ruth J. Person\"',\n",
       " '\"Cultivating Achievement, Success, and Excellence at Indiana University\"',\n",
       " '\"The Life Sciences of Today and Tomorrow: \\'Genes and Your Health\\'\"',\n",
       " '\"Indiana University: An Era of Progress\"',\n",
       " '\"Back Home Again in Indiana: Celebrating the Career of Curt Simic\"',\n",
       " '\"Action, Passion, and Hope: A Foundation for Our Shared Future\"',\n",
       " '\"The Iron of Courage: President Ellen Johnson Sirleaf\"',\n",
       " '\"Indiana University and Traditions of Leadership\"',\n",
       " '\"The Transcendent Power of Liberal Education\"',\n",
       " '\"A Figure of Excellence: Peter G. Hall\"',\n",
       " '\"From Pre-School to Med School and Beyond: Comprehensive Engagement in the Life of the State\"',\n",
       " '\"A History of Partnership: IU and Ivy Tech\"',\n",
       " '\"Partnerships in Education and Industry: China, Indiana, and Indiana University\"',\n",
       " '\"Investing in Community and Building Our Shared Future\"',\n",
       " '\"Looking Towards a Brilliant Future\"',\n",
       " '\"A United and Powerful Voice for Indiana University\"',\n",
       " '\"Limitless Vision at Indiana University\"',\n",
       " '\"Social Justice, Education, and the Information Revolution\"',\n",
       " '\"The Indiana Edge: Excellence, Innovation, and Economic Development\"',\n",
       " '\"Architects of a Global Future: Knowledge and the Cardinal Virtues\"',\n",
       " '“Gifts to the Future from the Lilly Endowment”',\n",
       " '\"The Reciprocal Gifts of Education: Our Hope for the Future\"',\n",
       " '\"American Universities and Businesses: New Opportunities and Challenges\"',\n",
       " '\"Commitment to Community: Anthony Thompson and Indiana University\"',\n",
       " '\"Rising to the Challenge and Extending IU\\'s Commitment to Diversity and Equity\"',\n",
       " '“Endurance, Excellence, and the Energy of Change at Indiana University”',\n",
       " '“The Wonders and Possibilities of Interdisciplinary Science”',\n",
       " '\"Synergy, Community, and Transactional Experience\"',\n",
       " '\"Allied in the Pursuit of Academic Excellence\"',\n",
       " '“A Transformative Journey: Freshman Induction 2007”',\n",
       " '\"Building on Strength at the Purdue School of Engineering and Technology\"',\n",
       " '“Diversity and Equity at Indiana University”',\n",
       " '\"Traditions and Persistence\"',\n",
       " '“Increasing Diversity at IU”',\n",
       " '“An Education for the Future: Presenting the Class of 2007”',\n",
       " '“Cooperation and Collaboration:  A Vision for Progress”',\n",
       " '“The Life of a Senator: Birch Bayh”',\n",
       " '“To Forge a Leader: Experience, Community, and Character”',\n",
       " '\"Presidential Announcement, Board of Trustees Meeting\"']"
      ]
     },
     "execution_count": 14,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "#reduce number of dimensions of text to make analysis easier\n",
    "#n_components is 3 because we want the visualization to be three-dimensional\n",
    "tfs_reduced = TruncatedSVD(n_components=2).fit_transform(tfs)\n",
    "print(tfs_reduced, len(tfs_reduced))\n",
    "\n",
    "#convert speech titles to a list to append to each x-y-z coordinate,\n",
    "#so users know what speech they're hovering over\n",
    "#to see a demo of the visualization, check out tsne3d.html\n",
    "speech_title = df.title.to_list()\n",
    "speech = df.contents.to_list()\n",
    "dict={}\n",
    "for speech_title, speech in zip(speech_title, speech):\n",
    "    dict[speech_title]=speech\n",
    "    \n",
    "#print for my sanity\n",
    "list(dict.keys())"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 32,
   "metadata": {
    "scrolled": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[t-SNE] Computing 70 nearest neighbors...\n",
      "[t-SNE] Indexed 501 samples in 0.001s...\n",
      "[t-SNE] Computed neighbors for 501 samples in 0.008s...\n",
      "[t-SNE] Computed conditional probabilities for sample 501 / 501\n",
      "[t-SNE] Mean sigma: 0.022359\n",
      "[t-SNE] Computed conditional probabilities in 0.021s\n",
      "[t-SNE] Iteration 50: error = 59.6349945, gradient norm = 0.5257760 (50 iterations in 0.140s)\n",
      "[t-SNE] Iteration 100: error = 60.3977928, gradient norm = 0.4967720 (50 iterations in 0.147s)\n",
      "[t-SNE] Iteration 150: error = 59.1676216, gradient norm = 0.4961248 (50 iterations in 0.171s)\n",
      "[t-SNE] Iteration 200: error = 59.4349670, gradient norm = 0.5022141 (50 iterations in 0.220s)\n",
      "[t-SNE] Iteration 250: error = 58.9213257, gradient norm = 0.5147982 (50 iterations in 0.187s)\n",
      "[t-SNE] KL divergence after 250 iterations with early exaggeration: 58.921326\n",
      "[t-SNE] Iteration 300: error = 0.4658591, gradient norm = 0.0009254 (50 iterations in 0.183s)\n",
      "[t-SNE] Iteration 350: error = 0.4131903, gradient norm = 0.0012081 (50 iterations in 0.151s)\n",
      "[t-SNE] Iteration 400: error = 0.4043961, gradient norm = 0.0002171 (50 iterations in 0.137s)\n",
      "[t-SNE] Iteration 450: error = 0.4015154, gradient norm = 0.0001492 (50 iterations in 0.180s)\n",
      "[t-SNE] Iteration 500: error = 0.3999741, gradient norm = 0.0001287 (50 iterations in 0.178s)\n",
      "[t-SNE] Iteration 550: error = 0.3882450, gradient norm = 0.0008334 (50 iterations in 0.146s)\n",
      "[t-SNE] Iteration 600: error = 0.3822435, gradient norm = 0.0001206 (50 iterations in 0.133s)\n",
      "[t-SNE] Iteration 650: error = 0.3819966, gradient norm = 0.0000973 (50 iterations in 0.137s)\n",
      "[t-SNE] Iteration 700: error = 0.3818411, gradient norm = 0.0000909 (50 iterations in 0.131s)\n",
      "[t-SNE] Iteration 750: error = 0.3813442, gradient norm = 0.0000580 (50 iterations in 0.138s)\n",
      "[t-SNE] Iteration 800: error = 0.3811342, gradient norm = 0.0000885 (50 iterations in 0.148s)\n",
      "[t-SNE] Iteration 850: error = 0.3808340, gradient norm = 0.0000896 (50 iterations in 0.127s)\n",
      "[t-SNE] Iteration 900: error = 0.3807232, gradient norm = 0.0001005 (50 iterations in 0.127s)\n",
      "[t-SNE] Iteration 950: error = 0.3805510, gradient norm = 0.0000839 (50 iterations in 0.127s)\n",
      "[t-SNE] Iteration 1000: error = 0.3805502, gradient norm = 0.0001039 (50 iterations in 0.142s)\n",
      "[t-SNE] Iteration 1050: error = 0.3805763, gradient norm = 0.0000854 (50 iterations in 0.134s)\n",
      "[t-SNE] Iteration 1100: error = 0.3803219, gradient norm = 0.0000916 (50 iterations in 0.132s)\n",
      "[t-SNE] Iteration 1150: error = 0.3804070, gradient norm = 0.0000877 (50 iterations in 0.135s)\n",
      "[t-SNE] Iteration 1200: error = 0.3801246, gradient norm = 0.0000633 (50 iterations in 0.142s)\n",
      "[t-SNE] Iteration 1250: error = 0.3800744, gradient norm = 0.0000689 (50 iterations in 0.136s)\n",
      "[t-SNE] Iteration 1300: error = 0.3799809, gradient norm = 0.0000804 (50 iterations in 0.136s)\n",
      "[t-SNE] Iteration 1350: error = 0.3799284, gradient norm = 0.0000523 (50 iterations in 0.136s)\n",
      "[t-SNE] Iteration 1400: error = 0.3797877, gradient norm = 0.0001074 (50 iterations in 0.124s)\n",
      "[t-SNE] Iteration 1450: error = 0.3797837, gradient norm = 0.0001252 (50 iterations in 0.112s)\n",
      "[t-SNE] Iteration 1500: error = 0.3797694, gradient norm = 0.0001000 (50 iterations in 0.101s)\n",
      "[t-SNE] Iteration 1550: error = 0.3796790, gradient norm = 0.0000928 (50 iterations in 0.105s)\n",
      "[t-SNE] Iteration 1600: error = 0.3797085, gradient norm = 0.0001065 (50 iterations in 0.101s)\n",
      "[t-SNE] Iteration 1650: error = 0.3796883, gradient norm = 0.0000792 (50 iterations in 0.105s)\n",
      "[t-SNE] Iteration 1700: error = 0.3796811, gradient norm = 0.0001090 (50 iterations in 0.102s)\n",
      "[t-SNE] Iteration 1750: error = 0.3797168, gradient norm = 0.0000698 (50 iterations in 0.104s)\n",
      "[t-SNE] Iteration 1800: error = 0.3796596, gradient norm = 0.0001074 (50 iterations in 0.103s)\n",
      "[t-SNE] Iteration 1850: error = 0.3795515, gradient norm = 0.0000628 (50 iterations in 0.101s)\n",
      "[t-SNE] Iteration 1900: error = 0.3795482, gradient norm = 0.0000851 (50 iterations in 0.106s)\n",
      "[t-SNE] Iteration 1950: error = 0.3795261, gradient norm = 0.0001018 (50 iterations in 0.101s)\n",
      "[t-SNE] Iteration 2000: error = 0.3795269, gradient norm = 0.0001133 (50 iterations in 0.140s)\n",
      "[t-SNE] Iteration 2050: error = 0.3793516, gradient norm = 0.0001238 (50 iterations in 0.142s)\n",
      "[t-SNE] Iteration 2100: error = 0.3794540, gradient norm = 0.0001105 (50 iterations in 0.125s)\n",
      "[t-SNE] Iteration 2150: error = 0.3794883, gradient norm = 0.0000986 (50 iterations in 0.154s)\n",
      "[t-SNE] Iteration 2200: error = 0.3793907, gradient norm = 0.0001382 (50 iterations in 0.141s)\n",
      "[t-SNE] Iteration 2250: error = 0.3793393, gradient norm = 0.0000913 (50 iterations in 0.130s)\n",
      "[t-SNE] Iteration 2300: error = 0.3792703, gradient norm = 0.0001296 (50 iterations in 0.125s)\n",
      "[t-SNE] Iteration 2350: error = 0.3794485, gradient norm = 0.0001128 (50 iterations in 0.137s)\n",
      "[t-SNE] Iteration 2400: error = 0.3794475, gradient norm = 0.0001091 (50 iterations in 0.130s)\n",
      "[t-SNE] Iteration 2450: error = 0.3794048, gradient norm = 0.0000903 (50 iterations in 0.141s)\n",
      "[t-SNE] Iteration 2500: error = 0.3793595, gradient norm = 0.0000784 (50 iterations in 0.134s)\n",
      "[t-SNE] Iteration 2550: error = 0.3793159, gradient norm = 0.0000874 (50 iterations in 0.139s)\n",
      "[t-SNE] Iteration 2600: error = 0.3793717, gradient norm = 0.0000864 (50 iterations in 0.128s)\n",
      "[t-SNE] Iteration 2650: error = 0.3792337, gradient norm = 0.0001001 (50 iterations in 0.134s)\n",
      "[t-SNE] Iteration 2700: error = 0.3792400, gradient norm = 0.0001030 (50 iterations in 0.148s)\n",
      "[t-SNE] Iteration 2750: error = 0.3793135, gradient norm = 0.0001012 (50 iterations in 0.136s)\n",
      "[t-SNE] Iteration 2800: error = 0.3793025, gradient norm = 0.0000878 (50 iterations in 0.138s)\n",
      "[t-SNE] Iteration 2850: error = 0.3792403, gradient norm = 0.0000495 (50 iterations in 0.169s)\n",
      "[t-SNE] Iteration 2900: error = 0.3792832, gradient norm = 0.0001296 (50 iterations in 0.141s)\n",
      "[t-SNE] Iteration 2950: error = 0.3792243, gradient norm = 0.0000750 (50 iterations in 0.114s)\n",
      "[t-SNE] Iteration 3000: error = 0.3791953, gradient norm = 0.0000815 (50 iterations in 0.106s)\n",
      "[t-SNE] Iteration 3050: error = 0.3791305, gradient norm = 0.0001372 (50 iterations in 0.104s)\n",
      "[t-SNE] Iteration 3100: error = 0.3791908, gradient norm = 0.0001125 (50 iterations in 0.102s)\n",
      "[t-SNE] Iteration 3150: error = 0.3791806, gradient norm = 0.0001214 (50 iterations in 0.108s)\n",
      "[t-SNE] Iteration 3200: error = 0.3791544, gradient norm = 0.0000755 (50 iterations in 0.101s)\n",
      "[t-SNE] Iteration 3250: error = 0.3791632, gradient norm = 0.0001014 (50 iterations in 0.101s)\n",
      "[t-SNE] Iteration 3300: error = 0.3792209, gradient norm = 0.0000779 (50 iterations in 0.108s)\n",
      "[t-SNE] Iteration 3350: error = 0.3791878, gradient norm = 0.0000739 (50 iterations in 0.101s)\n",
      "[t-SNE] Iteration 3400: error = 0.3791141, gradient norm = 0.0000850 (50 iterations in 0.100s)\n",
      "[t-SNE] Iteration 3450: error = 0.3791756, gradient norm = 0.0000860 (50 iterations in 0.106s)\n",
      "[t-SNE] Iteration 3500: error = 0.3790947, gradient norm = 0.0000819 (50 iterations in 0.150s)\n",
      "[t-SNE] Iteration 3550: error = 0.3790985, gradient norm = 0.0000924 (50 iterations in 0.139s)\n",
      "[t-SNE] Iteration 3600: error = 0.3791698, gradient norm = 0.0000802 (50 iterations in 0.137s)\n",
      "[t-SNE] Iteration 3650: error = 0.3790851, gradient norm = 0.0000705 (50 iterations in 0.121s)\n",
      "[t-SNE] Iteration 3700: error = 0.3791291, gradient norm = 0.0001094 (50 iterations in 0.129s)\n",
      "[t-SNE] Iteration 3750: error = 0.3791066, gradient norm = 0.0001044 (50 iterations in 0.132s)\n",
      "[t-SNE] Iteration 3800: error = 0.3791301, gradient norm = 0.0001102 (50 iterations in 0.134s)\n",
      "[t-SNE] Iteration 3850: error = 0.3790902, gradient norm = 0.0000898 (50 iterations in 0.129s)\n",
      "[t-SNE] Iteration 3900: error = 0.3790682, gradient norm = 0.0000945 (50 iterations in 0.131s)\n",
      "[t-SNE] Iteration 3950: error = 0.3790022, gradient norm = 0.0001110 (50 iterations in 0.130s)\n",
      "[t-SNE] Iteration 4000: error = 0.3790168, gradient norm = 0.0001039 (50 iterations in 0.144s)\n",
      "[t-SNE] Iteration 4050: error = 0.3790686, gradient norm = 0.0001000 (50 iterations in 0.135s)\n",
      "[t-SNE] Iteration 4100: error = 0.3790564, gradient norm = 0.0000887 (50 iterations in 0.137s)\n",
      "[t-SNE] Iteration 4150: error = 0.3790047, gradient norm = 0.0001007 (50 iterations in 0.132s)\n",
      "[t-SNE] Iteration 4200: error = 0.3790505, gradient norm = 0.0000813 (50 iterations in 0.133s)\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[t-SNE] Iteration 4250: error = 0.3790710, gradient norm = 0.0000971 (50 iterations in 0.144s)\n",
      "[t-SNE] Iteration 4300: error = 0.3790406, gradient norm = 0.0001038 (50 iterations in 0.141s)\n",
      "[t-SNE] Iteration 4300: did not make any progress during the last 300 episodes. Finished.\n",
      "[t-SNE] KL divergence after 4300 iterations: 0.379041\n"
     ]
    }
   ],
   "source": [
    "#run above dictionary through the TSNE model\n",
    "#n_components is 3 because we want this to be three-dimensional\n",
    "#function documentation can be found here: https://scikit-learn.org/stable/modules/generated/sklearn.manifold.TSNE.html\n",
    "model = TSNE(n_components=2, perplexity=23, n_iter=25000, verbose=3).fit_transform(tfs_reduced) #, method='exact').fit_transform(tfs_reduced)\n",
    "\n",
    "# save to json\n",
    "x_axis=model[:,0]\n",
    "y_axis=model[:,1]\n",
    "# z_axis=model[:,2]\n",
    "\n",
    "x_norm = (x_axis-np.min(x_axis)) / (np.max(x_axis) - np.min(x_axis))\n",
    "y_norm = (y_axis-np.min(y_axis)) / (np.max(y_axis) - np.min(y_axis))\n",
    "# z_norm = (z_axis-np.min(z_axis)) / (np.max(z_axis) - np.min(z_axis))\n",
    "\n",
    "data = {\"x\":x_norm.tolist(), \"y\":y_norm.tolist(), \"names\":list(dict.keys())}\n",
    "\n",
    "with open(r'C:\\Users\\Vivek Rao\\Desktop\\IDS\\Digital\\McRobbie NLP\\models\\perp23_ncomp2_verb3.json', 'w') as outfile:\n",
    "    json.dump(data, outfile)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 40,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Some helper functions for plotting annotated t-SNE visualizations\n",
    "\n",
    "# TODO: adjust_text not available in kernels\n",
    "try:\n",
    "    from adjustText import adjust_text\n",
    "except ImportError:\n",
    "    def adjust_text(*args, **kwargs):\n",
    "        pass\n",
    "\n",
    "def adjust_text(*args, **kwargs):\n",
    "    pass\n",
    "\n",
    "def plot_bg(bg_alpha=.01, figsize=(13, 9), emb_2d=None):\n",
    "    \"\"\"Create and return a plot of all our movie embeddings with very low opacity.\n",
    "    (Intended to be used as a basis for further - more prominent - plotting of a \n",
    "    subset of movies. Having the overall shape of the map space in the background is\n",
    "    useful for context.)\n",
    "    \"\"\"\n",
    "    if emb_2d is None:\n",
    "        emb_2d = model\n",
    "    fig, ax = plt.subplots(figsize=figsize)\n",
    "    X = emb_2d[:, 0]\n",
    "    Y = emb_2d[:, 1]\n",
    "    ax.scatter(X, Y, alpha=bg_alpha)\n",
    "    return ax\n",
    "\n",
    "def annotate_sample(n, n_ratings_thresh=0):\n",
    "    \"\"\"Plot our embeddings with a random sample of n movies annotated.\n",
    "    Only selects movies where the number of ratings is at least n_ratings_thresh.\n",
    "    \"\"\"\n",
    "    sample = mainstream_movies[mainstream_movies.n_ratings >= n_ratings_thresh].sample(\n",
    "        n, random_state=1)\n",
    "    plot_with_annotations(sample.index)\n",
    "\n",
    "def plot_by_title_pattern(pattern, **kwargs):\n",
    "    \"\"\"Plot all movies whose titles match the given regex pattern.\n",
    "    \"\"\"\n",
    "    match = df[df.title.str.contains(pattern)]\n",
    "    return plot_with_annotations(match.index, **kwargs)\n",
    "\n",
    "def add_annotations(ax, label_indices, emb_2d=None, **kwargs):\n",
    "    if emb_2d is None:\n",
    "        emb_2d = embs\n",
    "    X = emb_2d[label_indices, 0]\n",
    "    Y = emb_2d[label_indices, 1]\n",
    "    ax.scatter(X, Y, **kwargs)\n",
    "\n",
    "def plot_with_annotations(label_indices, text=True, labels=None, alpha=1, **kwargs):\n",
    "    ax = plot_bg(**kwargs)\n",
    "    Xlabeled = model[label_indices, 0]\n",
    "    Ylabeled = model[label_indices, 1]\n",
    "    if labels is not None:\n",
    "        for x, y, label in zip(Xlabeled, Ylabeled, labels):\n",
    "            ax.scatter(x, y, alpha=alpha, label=label, marker='1',\n",
    "                       s=90,\n",
    "                      )\n",
    "        fig.legend()\n",
    "    else:\n",
    "        ax.scatter(Xlabeled, Ylabeled, alpha=alpha, color='green')\n",
    "    \n",
    "    if text:\n",
    "        # TODO: Add abbreviated title column\n",
    "        titles = mainstream_movies.loc[label_indices, 'speech_title'].values\n",
    "        texts = []\n",
    "        for label, x, y in zip(titles, Xlabeled, Ylabeled):\n",
    "            t = ax.annotate(label, xy=(x, y))\n",
    "            texts.append(t)\n",
    "        adjust_text(texts, \n",
    "                    #expand_text=(1.01, 1.05),\n",
    "                    arrowprops=dict(arrowstyle='->', color='red'),\n",
    "                   )\n",
    "    return ax\n",
    "\n",
    "FS = (13, 9)\n",
    "def plot_region(x0, x1, y0, y1, text=True):\n",
    "    \"\"\"Plot the region of the mapping space bounded by the given x and y limits.\n",
    "    \"\"\"\n",
    "    fig, ax = plt.subplots(figsize=FS)\n",
    "    pts = df[\n",
    "        (df.x >= x0) & (df.x <= x1)\n",
    "        & (df.y >= y0) & (df.y <= y1)\n",
    "    ]\n",
    "    ax.scatter(pts.x, pts.y, alpha=.6)\n",
    "    ax.set_xlim(x0, x1)\n",
    "    ax.set_ylim(y0, y1)\n",
    "    if text:\n",
    "        texts = []\n",
    "        for label, x, y in zip(pts.title.values, pts.x.values, pts.y.values):\n",
    "            t = ax.annotate(label, xy=(x, y))\n",
    "            texts.append(t)\n",
    "        adjust_text(texts, expand_text=(1.01, 1.05))\n",
    "    return ax\n",
    "\n",
    "def plot_region_around(title, margin=5, **kwargs):\n",
    "    \"\"\"Plot the region of the mapping space in the neighbourhood of the the movie with\n",
    "    the given title. The margin parameter controls the size of the neighbourhood around\n",
    "    the movie.\n",
    "    \"\"\"\n",
    "    xmargin = ymargin = margin\n",
    "    match = df[df.title == title]\n",
    "    assert len(match) == 1\n",
    "    row = match.iloc[0]\n",
    "    return plot_region(row.x-xmargin, row.x+xmargin, row.y-ymargin, row.y+ymargin, **kwargs)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 38,
   "metadata": {
    "scrolled": true
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.collections.PathCollection at 0x21a7f36ea48>"
      ]
     },
     "execution_count": 38,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAlgAAAHSCAYAAAAwpbX/AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nO3dfYwl2Xnf99+ZvrffezyzszMkd4fr5SprQpRAU06DYZD8oUg0RBmO1nYsgEISETKBhQAJsRD/QdIL2HAMAjIE2IET28gCEkwDtGnBcsKFZEUiCQpKgFDUMKFl0iuay6XI3R2aO5yd4U6/THff7pM/+t6e23fq5ZyqU1Wn6nw/AMGd7tv31q3Xp57znKeMtVYAAAAI50LXCwAAADA0BFgAAACBEWABAAAERoAFAAAQGAEWAABAYARYAAAAgY26XoB5jz76qH3yySe7XgwAAIBSX/rSl75nrb2a9buoAqwnn3xSN27c6HoxAAAAShljvpX3O4YIAQAAAiPAAgAACIwACwAAIDACLAAAgMAIsAAAAAIjwAIAAAiMAAsAACAwAiwAAIDACLAAAAACI8ACAAAIjAALAAAgMAIsAACAwAiwAAAAAiPAAgAACIwACwAAILBR1wsAAAD66f7RsW7vHOhgcqKV0QVd2VzR6nip68WKAhksAADg7f7RsV67s6cTK60vL+nESq/d2dP9o+OuFy0KZLAAAIC32zsHWh4taXl0mqtZHpmznz9+eV1S2hkuAiwAAFAoK1A6mJxoffl8sDReMto7PD77m9fu7Gl5tKT15SUdHVu9dmdPj19eTyLIYogQAADkyhsKlLU6OrbnXnt0bLUyzWjNZ7iMMVoeXdDyaEm3dw66+BqtI8ACAAC58gIlGaPDybEOJyey1upwcqLDybGubK5Ikg4mJxovmXPvNV4yOpicdPE1WkeABQAAcuUFSpL0+OV1XTDS3uGxLhidG/5bGV0ozHANXRrfEgAAVFIUKK2Ol/T45XU9dmlNknTz7v7ZTMIrmyuFGa6hI8ACAACSHtRbvXxrxzlQyq3RUnGGa+iYRQgAAApn/T1+eV23dw60d3isldGFc4FSWbuGVFs2kMECAACFs/6KhgJditlTbEpKgAUAAEoDpaKhwLJi9hRbNhBgAQCA0ll/ue0arC0tZk+xZQMBFgAAKC1mLwqSjE6zWd+8taOj45OHitlTbNkw3G8GAACczeqsfPpa7dyf6M7eocajJT356IauP7Iha+1D751iywZmEQIAAEkPgqwsVzZXzmquxktGR8dWt+7d19Wt1dwZhPMzB42ko+MTHR3roZmIQ0SABQBA4lxaKMyCr/l2DZc2VrS1ej6UmD3wOavtw+HkePCB1UztIUJjzKox5ovGmH9rjPmqMebvTH/+DmPMHxhjvm6M+ZfGmOX6iwsAAELyaaEwC7Keurqpxy+v6+LqKLe2KsWZg/NC1GAdSPoxa+2flfQeSR8wxrxP0t+T9A+stU9LuiPpwwE+CwAABFQnECqqrUpx5uC82gGWPbUz/ed4+j8r6cck/avpzz8h6S/V/SwAABBWUSCU9eiceUWF8SnOHJwX5FsaY5aMMV+W9Lqkz0j6hqS71trJ9CWvSno8xGcBAIBw8gIhWes0dLg4bDirr0px5uC8IAGWtfbYWvseSdclvVfSD2a9LOtvjTHPGmNuGGNu3Lp1K8TiAAAAR3mBkIwpHDqsk91KQdA8nbX2rqTfk/Q+SZeMMbOpBdcl3cz5m+ettdvW2u2rV6+GXBwAAFAiLxCSVDp0WDW7lYIQswivGmMuTf97TdL7Jb0o6fOS/ur0ZR+S9Om6nwUAAMLLCoSKaqhSnyHoIkQG622SPm+M+SNJfyjpM9ba35T0EUn/ozHmJUlXJP1qgM8CAAAtYIZgPbUbjVpr/0jSj2T8/GWd1mMBAICeyWosupjdmnVul9KaIeiCTu4AACBT3qNzsh6bM+vSjlMEWAAAJMbl0ThFirJbOEWABQBAQrKeEfjanT3vAKnowdAgwAIAoPd8MlLzMwAlndVR3d45IGAKiGo0AAB6zOdhzVK9R+PAHQEWAAA95tuTqu6jceCGAAsAgB7z7UlV9dE48EOABQBAjxV1XM9S5dE48EeABQBAjxV1XM+SVxDvG6ihGGsNAIAey8tIZc0iLCqIXwzU3tw/0itv7OrN/SNqsSqgTQMAAD3n2pOqrEXDrHnonb0j3d090LWtVW2ujir3ykoZGSwAABJRVhA/C9Quro50/ZENba2NKXiviAALAIBEuNZZ+c5MxMMYIgQADEbdZ+wNXdFDmufX3Z29A52cLGtrbXz2txS8+2FNAQAGwbejeYqKWjTMr7uLq8t65Y1d3ds/cpqZiIeRwQIADEKVZ+ylmPHKKoifPfx5tu4uro11/ZENff/+oS5cMFoZXaDA3RMBFgBgEA4mJ1pfPh8AjJeM9g6zM1izjNfyaEnry0tJz5TLWndbqyMtXTB66upmR0vVbwwRAgAGwbdRpu8z/IaMJqPhseYAAIPg29GcmXIP+K47lGOIEAAwCLPaots7B9o7PH6obmix3ko6zdLMarVm/04xa1O27uCPAAsAMBh5Hc2z6q3uH060b4wuro4falmQItdu8HBDgAUAGLysGYZba8s6mhyftSwYetYmxRmTXSLAAgD0lmvQkDfD8OjYJJG1YcZk+wiwAAC95BM0zGbJtVVv5RL4tZlRqtIjDPWkV8kHABgEnzYLbc6Sc+ko33bXeWZMto8ACwDQSz5BQ94jYprIGLkEfm334KLPVfsYIgQA9NLK6IJ27k+0ezg5G2bbWB5pbTk7aGprlpxLR3nfrvN1FT3kGc0gdAUA9NLGykivvLGrg6MTrY4u6ODoRK+8sauNlW5zBy7ZorYzSm1m8HCKDBYAoJd2Dya6/siGdg8muj850fLogi5vnP770vpyZ8uVlS168/6R1kYX9PKtndNM28robDiwrYwSfa7aRYAFAOilg8mJtlZHurg2PvuZtdZrmK2JmXyLXdElyVir8WjpLJi6vXOgK5sr2j2YJNGDK0UEWACAXqrbeqHJ3lDz2aLX7uxpvLT8UIuE3YMJGaUBowYLANBLdVsvVJnJNwvKXr6149xWgRYJaSLAAgD0Ut3Cbd/Ap2rvKlokpIkhQgBAb9Up3PYdYqzaDb2sRQLPCBwmwmcAQJJmQ4z39o/0nbt7+sZ37+nVgjYPVYf6ijJtbXd0R3vIYAEAkrQ6XtKVzRX9+5vfl9VpgLO5OtbtnQOtjpeCPs8wL9PGMwKHiwALAJCs3YOJ3v7IxlmAI0mHk5PMAKeJbuh1O7ozvBgvhggBAMlq43mGRTMP6xTAM7wYNwIsAECyfAOcWZD11NXNc32u8to2lAVBdVpNtP3AaPghwAIAJKtOgOOSQSoLguq0mqC/VtyowQIAJGvxsTY+j6xxKVB3qbGq2mqibid7NIsACwCQtKoBjkvw1GQQ1ETRPcIhwAKAgJjVNSxF29MleGoyCKqTfUPzCLAAIJBQDw8mSItD2fZ0CZ6aDoKKsm/sR90iwAKAQMpqclwueKGCtCHpKlAo256uwVOdx/lUxX7UPSrhACCQolldrj2LmHp/Xpe9nlxm6S22bWgyeCnqp7WI/ah7BFgAEEhRTyXXCx5T789rKlBwCVbqNAENuRzzr3MNNNmPukeABQCBFPVUcr3gtXFR75MmAgXXYKVOj6yQyyH5B5rsR91jTQNABVmZh6Kmka4XvKYv6jHpKovkGqzUaQIacjkk/0Azpf0oVgRYAOCpKPOQV5PjesFr+qIeiy6zSFWeP9hEjZXPclR9pM/Q96OYMYsQADy5dPBeNLvg3by7r1ff2JcxRte2soOELmadtc11Hfq2OXCZcRhLB3Sf5ajSTyuF/ShmtfcmY8zbjTGfN8a8aIz5qjHmr09//ogx5jPGmK9P//9y/cUFgO7VqQuy1ur6Ixt68tENjUdLQWfE+cwya4LP5zeRRYqltsqVz3KQkeqfEOH6RNLfsNb+oKT3SfoFY8y7JH1U0uestU9L+tz03wDQe1XrgpqcOt9lO4Mqn5+1DnfuT3Rn76BygBhLbZUr3+VosyUE6qs9RGit/Y6k70z/+54x5kVJj0t6RtKPTl/2CUm/J+kjdT8PALpW9fEnLs+uq6rKsGVIvp+/uA537k/0yhu7uv7IRmljzLxhQJ/1G8vwWSzLgfCCDjgbY56U9COS/kDSW6bB1ywIuxbyswCgK1UzIHVnxBUNwXXd98j38xfX4ffvH+r6Ixu6uDYuzD4VZcpoTYCYBNvrjDGbkn5D0i9Za9/0+LtnjTE3jDE3bt26FWpxAKBRVYZr6tT+lA3BdR1cVPn8+XV4eX1FW6vnB1WyArSiYcBYaqsAKVCAZYwZ6zS4+qS19l9Pf/xdY8zbpr9/m6TXs/7WWvu8tXbbWrt99erVEIsDAFHKy3xJKi0OL6sv6jq4qPv5rgFaUaYsltoqQAozi9BI+lVJL1pr//7cr16Q9KHpf39I0qfrfhYA9N1i5kuSU3F42RBc18FF3c93DdDKAjEKwRGLEH2w/gtJ/72kf2eM+fL0Z39T0i9L+nVjzIclfVvSTwf4LAAYFNficJeeSV0XTNf5fNd+V1UnGABtCzGL8P+WZHJ+/eN13x8Ahsxl5tv9o2MdHB3r22/saWNlpKtbK1q6cOFcYOHSYLNtvsvkEqD5Nh4FusLUCgDoUNmQ16y4fTxa0pOPbkhG+pPv7ero+OQssOi6B1YW12Wq0hyVYUD0AQEWAHSorPZofghxbXmkJx7Z0A9c29LK6MJZYNFkA9OqXJYpxsAQCIVnEQJAh8qGvFyGEJtsYJqnbPjPZZmabI4a45Ap0kKABQAdK6o9cilub/vhxbPM0/JoKbfrussyVQ0My4Inl+UDmsYQIQBEzKV9Qds9sFyG/1yWqUpzUpdhxRiHTJEeAiwAiJhLf6m2e2C5PBZnfpnu7B3p9Xv7Opyc6PbOwVkwVCUwdAmeun5sECAxRAgAUSga9vJpX9AG1yHJ1fGSrmyu6P7RsTa31s76Vs0P1/m2XHAZVmx7yBTIwt4GAB3r22w6n8xTWcbJt+WCy7Bi148NAiQCLADoXN9qhnyGJEMP17kET0XLV6XvFlAFQ4QA0LEu2iyUKZup5zokGXq4znVYMWv5mF2INpHBAoCOVZlN16SQQ5ZNDNdV7eTet0wh+o0ACwA6FlvNUMhApO0ZjkWYXYg2MUQIAB3r+gHGi8OBb96f6PL6+Nxr6gxZtjnDsQizC9EmAiwAiEBXQUhWXdLd3QONLhhdXHsQZA0hELmyuaLX7uxJ0lnLiMPJcRTBH4an30cLAKCWrOHAq1urunXvfjRDlqHENFyJ4SODBQAJy5rBuLk60tHx8lkg0vaQZZNiGa7E8BFgAUDC8uqSLq6NzwKRspYNAB7GECEAJKxsBmPfuswDsSDAAoCEldUl0TsKqIYhQgBIXFFdUpNd5hl6xJCRwQKAhgzhuXdNdZln6BFDR4AFAA0YSgDRVJd5hh4xdARYANCAoQQQTfWO4rE1GDpqsACgAU3WLrWtid5RPLYGQ8eeDAANaKp2aShie8A1EBpHOgA0gACiGI+twdAxRAgADZgFELd3Dgb3uJlQeGwNhowACwAaQgCBoaBnmT8CLABAKS6w6Zq1HFkeLWl9eUlHx1av3dkjI1uCGiwAqGgIjURdDKWnF6oZSsuRthFgAUAFKQUdXGDTRs+yagiwAKCClIIOLrBpo+VINawdAKggpaCDC2zaaDlSDUcHAFSQUtDBBTYti7WFkuhZVgGzCAGggiubK2cXn/GS0dGx1eHkeJBtGWLt6cXMxvCKZgwOcd9u0vButQCgBal1Ip9936eubkbxPVOaZNCmlGoLm0YGC8A5ZAXcpdxItGw/aXo/mg8EJJ09NPr2zsFD24R92t2QHlLeNTJYAM6QFYCLsv2kjf3IdZIB+7SfkLWFqfSJy0OABQyY7wmO4QG4KNtP2tiPXAMB9mk/oSY0ENgSYAHRCXXXV+UE19fWA03dKad+B56nbD9pYz9yDQT6uk93JVRtIYEtNVhAVEI+88ulRmWxNkU6zQLMXjv7d8ytB3zXmWs9Ds9fyzfLHuXtJ2W/D8F1ZmMbyzI0IWoLm6zl6ktNHXsYEBHfu76iDEvZnXtWhuv+4URv3j/qVb8jn3Xmk9XjDjxfWfaorb5ZLjMb6eHVjod6Z1nbSJ+4Pg09EmABEfEZzig70ZTVqGQFEFtry1obXehV6wGfdeYTNDG0lK9sGCnr91c2V3R756D14dbU2ml0IetctD850b39w+CBbZ9ufBgiBCLiM5xRNgRY1ggzL4V/dGx61XrAZ535DFswtFSsbBhp/vddD7em3E6jKp9huKxz0cXVsY6OT84C21DNafvURoIzBRARn+GMsgxL2Z37UB714rPOfL4zQ0vh9CnrAP9huLxzkaTgzWn7dN6Kb4mAhPkMZ7icaIpqVIYSQPisM5/vzNBSOAy39otvQNx00DNf33XQ0NBjExgiBDqUl4Z3Gc6o+yy8WJ8vV4XrOvP9zgwtnVe1ezvDrf3iOwzX5HM5s4aX943R0eRYR8cm6vMWARbQkbp1KSECpBQDiBS/cwhl+2vR71N6MPYQ+AbEIc5FecF5Xn3XLKMcMwIsoCM+z1LLUydY6EsvmbpS+Z5V1C1knv18dnEt+n3WBVg6re1h28SlSkBc91yUF5z3qah9EQEW0BHfE0fIQKFK9iyGQMV3GbqevRYz33VTtr+W/X7xAsy2qa7pY7Ht8oGi4LwsmxbDeSkPA+BAR3wKQ0M316va0LTL5n5VlqHK7LVUHo8TupDZt9CZmYXV1D0WXfdvlyauoRRNgiiamBLDealIkADLGPNrxpjXjTFfmfvZI8aYzxhjvj79/8shPgsYCp8ZbS4XI5/AwHdWVwwXwyrL4Ps9Yz9hh+S7bkJ3b2dmYTV1jsUm9u8QNyRFwXnRbN4YzktFQmWw/qmkDyz87KOSPmetfVrS56b/BjDl0wagymNvik52vtmGKh3mQ2eAqlyQyark8103Vbq3F2U++tTPKCZ1AtPQ+3eogK0sOM/LpsUepAepwbLW/r4x5smFHz8j6Uen//0JSb8n6SMhPg8YCtfC0LI6BN+Ced8iVtdZRaHqarLqKqpM9ff9nn0uqPXVRCGzT6EzMwurqVOTVHX/9pnhJ/lN1JHKa7762v6jyaV4i7X2O5I0/f9rDX4WMGhld3i+d3K+2QbX4Z8Qd8h5d8UbKyPvxqghsio79ye6s9f+M/TK1M0Udt1ItevP76s6NUlVsoZF7xkyg5SXpSr6/NibJRtrbfmrXN7oNIP1m9baH57++6619tLc7+9Yax+qwzLGPCvpWUl64okn/tNvfetbQZYHqCrWWSlFyzU7AS3PnSgPJydBe8W4rJeXb+1ofXlJxjw46VprtXd4rKeubjp9TtF3mT0wuGgZ6my/+QzceMlo5/5Er7yxq+uPbGhrdXQuy9LlPrG4nC7LVWW9xHospC5vu5SdB6rsN0XvKanx847Ld+pyHzXGfMlau531uyYzWN81xrxtugBvk/R61oustc9ba7ettdtXr15tcHGAcjEXOXf92BuXWUUh6mqK7orLlqHu9lvMqnz//qGuP7Khi2vjqGqy2pgFGvOxkLqqNUlVsoZVZ/iFEnudVZEmA6wXJH1o+t8fkvTpBj8LCKKvRc4hh1tch56yXhfihFsnSAux/eYvXpfXV7S1er5UNYaTexuzQPt6LKSsiYkDVWf4hVL0+bHfBIRq0/AvJP0/kt5pjHnVGPNhSb8s6c8bY74u6c9P/w1Erc93SyH61riesPJeJ6n2CTcrSLu3f6iDyUlp0Oey/Xxql2Kd6dbkLNCqf5NK/7CYld3gVAlIqs7wa+M7xX4TEGoW4c/k/OrHQ7w/4KrueHyMs1LarDFwnRXk8liUqhZnFMlaWWM0XrpwVjeSNzPRZYaVzyzHWGe6NTULtOrf0JU9DmWz8arM+uv6ofCr4yVd2VzRN27d0879iTZXR/qBq1taHS9FP+uXR+VgMEKc5MsuXE0EO0Xv2fSFa/Gz37w/0eX18bnXZJ2wmj6xzU/3f+3OnsaOF4Wy7edygcnaHrsHk04uLnl8L3pVAkWfvwk1XR/1FbXKqHrcdvmA9NnxeG1rTY9fOt0Pb+8cJN+mAWhVyPqbrCGuJrsg571nkynwrM++u3uge/cn516XdcJqc+jMZ6iqrCakSsPW2zsHurK50sojQ3z4DM1UqZUJ2QgXzXIdng1x3LY9FFx0Doy9TQMBFgYj1Ek+78LVRLBT9p5NXriyPvvq1qpu3btf2GPn5Vs7Ojg61pv3j1o5sfk+s7Eow1j2XrHXdNThWyvjk62NtVYtBT43fnUDki6Kyl1mFcfaS429H4PR9Em+iWCn7D2b/E5Zn725OtLl9eWHTliSzp1Yx6MlGWt1dHzS+InN9aLgcvIP3bB1qHwvpLFnEobM56agbkDSxQ1I2TmwzYdS+6IGC721eIe9sTI6O9CbKEhuYry/7D2bLLLO++yLa+OH3n9WBzZfY7O1thy0oWAe13ojlzqgsveKvaajLb41VV0XQqfMt66qrJ6qiUft1BHrRBMXaZ01MBhFtTJNpYubuEt3nQLdxHfy+T5dZ3Zc7lJdl7Hrhq19UGV7x5xJGLKQWe4mHrVTV+zDgEXIYKGX8u6wdw8mjd3ZNHGXXvUhp21/nxgyO671VXWWkUzMqRi2N9yEzPCUZS67yiZ1OYuxDgIs9FLop8K7auJAz3vPNnoLuX6frtP0Lusi1DL29WQeUtfbG+5C3hSUnVf7cgPS9fMJZ7gdQS+Ffip8jGKa0dZ1mt5lXXS9jENSZ13S0b19oYZny86rsQQuRWI6z5PBQq/MDvA37090d/dAV7dWtbk6crrD7lszxNi6FHeZ2XFdF2Sfwqh6IaWjezhNBzNZ71+UuWy76XHV7xvTeZ4MFnpj/s7k8vpYj26t6vV793Vn97D2U+FjRG+hB1gX7amTAYgp69pnvtvAN2tY5VmibTc9rpp1iuk8TwYLvbF4Z3JxbazV8VJmq4Csu6G2CndD3YlRB/MA66I9dTIAsWVd+8pnG1TJLFV5lmiT2zZk1immCRrc/qE3XO9M8u6GNlZGjU/BD3knRk3RA6yL9tTJAJBpDMNnG1TJLFXZxm03Pa6adYqp1Qp7PXrD9QDPO+HMWjg0eZEOnUant9ADrIt21LmQxnRx6zOfbdBWsNTktg0ZvMV0M8YQIXrDdZioKJXddBG0axq9D7Nx+oz1W13V4djZOj+cnOju/qHWl8e6uDoiGK7AZxtUGRKrso1Dtmho+ikcsUx2IYOF3nC9M+lymMLls2OaRjxErN96qmQAzk1A2VjWta01LS8ZAtuKfLaBS2ZpsQheyi9md1muOlnkJp/CEVuLEGOtLX9VS7a3t+2NGze6Xgz03HzR5+LdUNMne5fPnp1clueCrsPJSSvP9UsB67d9rPNuFWVsuzwfZmlqX+nqexpjvmSt3c76HUOE6L2sk8sslX1n91B7RxOtL4/P7pKaPNhc0ujMtPLnM+RXtn5THT5s8nuzT7cr75yXJaa+UFJz+0ps31NiiBA9V9TP5crmipZHF3Rta02X18etDRXNgqzHLq1Jkm7e3e/8gal95jvkV7R+Ux0+bPp7s0+3x3dbxtQXSmpuX4nte0oEWOi5oll7XTY9LDoJMtPKj+92LFq/qTbCbPp7s0+3x3dbxhb8NrWvxPY9JQIs9FzRXUtbdzRZhZVFJ8GYphH3ge92LFq/Md7ltqHp780+3R7fbdlW8OtaYN7UvhJjkE8NFnqtbIpy0x1987ooHx5bXV4fn3tt1lPpUa7KNPS89RtTl+c2hfreRXVc7NPt8N2WIdsr5PHtJh9iXymqvW3qe/oa9lkFg1d019LGHU1epmrv8Mg7XR3bFONYFG1H33UW411uG8q+t8t6TLV+LTZV9uH59gqzofKQ55m2h97LnqWYV//aNgIs9FpRurmNYYu8dP362O+xPFy88uVtR0ne6yzVoayi7+2676VavxabOvtwU+eZtofei/bFmM6lDBGi94rSzU0PW+Sl6y+ujc/uFF3S1TFOMY5J1nacDUnkrbO84axUh7LyvrfrvkcrhnhU3YebOs+0PfRetC/GdC4lgwXUUJSuz+p6nDcUk2rxdR1F6yymu9jYue57Mc7Sgp+mzjNtD70X7YsxnUs5MoAafNL1RRd9Ll7+itZZ6sNZPrVprvteqvVrfZW1DzR1nml76L1oX4zpXMrZG6jJ9flcRRd9Ll7+itZZ2V3skCcU+GbvXPe9VOvX+ihvH9hY8asN9RHiOYW+n5W1L8Z0LiXAAlpSdNHn4uWvaJ2l3M3dN3vns++1eRFFdXn7wO7BZDDnmbx9MaZzKUXuQEvKCkFTKL4O/Ty8vHV2ZXPlbNr24oNfYyqCbUKVYvQU9r2UFO0DKWzrWL4jGSygJTGlrrvgmjkKMXyXcjf3mGpQ0A2ffWDIw+Vd44gDWhJT6roLLkNXIYbvZu9x8+6+JOmxS2vn1vPQA5DUA3m47wNDHy7v2jDOKEBPNN1ROWYumaO6s/9cLhhDCECKsg7zgfydvSO9fm9fh5OTsyaMGL5Zsffr9/b1tf/4pl6/t585HJ/6bNumEWABHUjxztElc1R3+M7lgtH3TKLLvjO7wC4vGV3bWtPljeUk9rFUlA3rzWodr22t6Z1vvahrW2uZAfbQh8u7RoAFNKDsBJjinaNL5qju8J3rBaPPs+Fc950U97EUuATYrtt+6MPlXWMtAoG5nABTvHN0yRzVHb5L4YLhuu+kuI+lwCV4ct32QxgujxltGhCt0FP62+LSBqDtZ3fFomz69Oz3Wc9wdNkfitozDIXrvpPqPjZ0Lm04XLd90fGG+jjSEKU+1yi53D1y55g/jFr0DMey/aHv9VUuXPcd9rFhcsnS+swi7ONNbF8QYCFKfa4fcTkB1g0E+t67xjeA9tkf+lxf5cJ130kh2EyRS/Dksu37fBPbFwwRIkpVulHHwnWYajbTa3YHOXsmYdkFcHZiXB4taX15SUfHVq/d2evVxdO3m3qf9wzP1McAACAASURBVIdQqmQbYulojXBch/XKtv3Qn2gQAwIsRMm1hiDGFLfrCbBqoDSEE6NvwJR6PZHPvhLjMYGwQgTO3LQ0jwALUXLJAoXO5IS8MLmcAKsGSkM4MfoGTCkUrxdx3VeGkN3EeXXOS0V/m/pNSxtYk4iSSw1ByDqtLuoRqk6jH0IrAt8C7Bjqibqse3PdV/pcu4iH1Tkvlf0tkyCaRwYL0SrLAoXM5HQx7Fb1DnII2Zwq7Ri6rCfqOjPkuq8MIbuJB+qcl4r+dlb7eTg50d39Q60vj3VxdUSmMzACLPRWyBR3FxemqoHSUHrXZAVMXQcyebque3PdVxj2GRaX81LeDUne397ZO9L9o2Mtj5Z0eWNZm8fjs8zV7Bijji8Mjjr0VsgUdxfDbnWGvYbaiiDWIa6uu6K77isM+wxL2XmpaBgw72/3Do8KjzHaN4RDBgu9FTKT45ohCH1nxzT682Id4oohM+Syrwwlu4lTZeelsmHArL9dH48ybxZmx1jX2dohIYOFXpv1kloZXTjrJVXlTovGfM1yLRCPtYC/T5mhoWY3U1R2XirKrOb97cW1ceEx1nW2dkjIYKHXQtbsNNGYj1oGv20UawE/mSF0pei8VJZZzfrbsmMshmztULDG0Gtt1uz43tmR8TpVto3ms1uzoY0YH+/SVmao749BQnuqZFbLsmJ9ytbGjgALveYT9NS9cPkOX8VasN22om2UFYTOgqwUh7hcgnICMMw00R8uhp5zQ9F4gGWM+YAx5mvGmJeMMR9t+vOQFtegJ0Q2yffOjlqGU0XbiCD0PNdsX+pZ0aGqEjwvZlYlFb6Hyz5EHV8YjQZYxpglSf9I0k9KepeknzHGvKvJz0RaXIOeEBdy3zu7WAu221a0jQhCzytbHwSkwxUieHZ5D/ah9jR9pn+vpJestS9baw8lfUrSMw1/JhLiGvSEupD73NlRy3CqaBsRhJ5Xtj4ISIcrRODj8h7sQ+1pehbh45Jemfv3q5L+s4Y/E4lx6Q/UxcwYZp49kLeNYp012BVmeKUrRA84l/dgH2pP02vUZPzs3O2ZMeZZY8wNY8yNW7duNbw4SFVX2SRqGYpRUHseM7zSFSKb6/Ie7EPtaTqD9aqkt8/9+7qkm/MvsNY+L+l5Sdre3j6/ZwAZfHpLzb/WSDo6PtHRsZLOJsWGbvbnFa0PsqLDFSKb6/IedfYh+vr5aTqD9YeSnjbGvMMYsyzpg5JeaPgzMWA+haCLrx2PlmSt1WOX1rgoobfIig5TiGxukxlhZrD6azSDZa2dGGN+UdLvSFqS9GvW2q82+ZkYNp9u6jxTC0NFJmGYQmRzy96j6tMvOJ/6a7yqzVr7b6y1f8Za+wPW2o83/XkYNp8ZMMyWwRCRSUAdVWcrcj71x7QB9IpPISgtADBE9DFCHVUDJc6n/lgz6JWiGTCLXZA3VkbMlsHgkElAHVUDJWYf+iPAQq/kFXFKyn2mHS0A0Ed5j00hk4A6XAKlrH2Plir+mm7TAASXVcQ5K9pcLMDcPZhQgIneKSpEpjkr6ihr01BWBM9+5o4AC4MQogsyEIuyGVv0wkIdRYESswXDIcDCIPD4BwxJ2Q0DmYSwaHvxADer4XD1wSBQgIkhoc6qPbS9OI99LxzWGAbBpQAzr2gYiA03DO2h7cV57HvhEGBhEMpS/Nylok+YsdUe2l6cx74XDjVY6D2XRz9QuIm+oc6qPpfaKuo3H8a+F0a6exAGwyXFz10qkBbXrDVDYmgKGSz0nsusF+5SgbS4Zq2L+kL5zi4c2mzEoX2ftnF1Qe+5zHrhLhV90MREjFQnd7hmrfOCCN+6zaHVeQ7t+3SBAAu95xI8UbiJ2DVxQUv5Iuly41W0fnxnFw5tNuLQvk8XCLDQe67B0+x1T13dJLhCdJq4oKV8kXS58SpaP751m0Or8xza9+kCNVgYBGa9oO+a6KCdclfu1fGSrmyu6Bu37mnn/kSbqyP9wNWtczdWRevHt25zaHWeQ/s+XWBNAUAEmuignXJX7tkw37WtNb3zrRd1bWtNt3cOzg2PFq0f37rNodV5Du37dGH4RxkA9EATF7SUL5Iuw6Nl68fotCbrm7d2dHR8UlhaMLQ6z9n3OZoc60++t6tX39iVMab8D3GGIUIAiEBRu4A8ZdPoq7znULgMj+atH0lnzYuffHRDR8dWh5PyYdUhlipYnQaK4yWT2cQZ+QiwACASPhdolycY+L7nkLjWEGWtn9l6Tf3JDzwBox6GCCEp3V45QF+lPEPQRZ3hUWbQnWI91EMGC853wl2iozBw3hBnCNY5zrP+turwKDPoTrEe6mEtIfo7YZ9miWTikIqhzRCs0xQ1728lVep955P9GvI5J+VJEiH080hEULGngV0DwJS7ViM9Q7v41bnRC3GTOB8o3d450JXNldIZgUM/5wxtZmTbGCJE9Glg16EQCjKRkqHNEKwz5Fl3uDSrTGJ23ihanymcc1KdJBECAVZisuoUrmyunKXTZ1NxDyfH0RxUrgHgEGtSgCJDuvjVudGre5PoEihlnTs556BIHCkKtKKsTiHWNLDrUMjQalKAlNQZ8qw7XFpWJpF77rSWcw5ysRckpKhOIeYHIbvWAQytJgVISZ16n7q1QmU3Z3nnThnDOQe5GCJMiGs6O8aWCC5DIUOrSQEWxXhshlRnyLPO35aVSeSdO4+OxTkHuQiwEuJSp9CHnlhFhlSTAszr+7EZs7Kbs6JzJ+cc5CHASohLMXsKs2JiV5alGHoWA9k4NptVFCjFPhEIcaIGa4DyGt+51CnE3hMrJk00GCzrqzP0vjvIx7HZHfpBoQoyWANTNoxQls6OvSdWLJoarinLUpDFGK6yzCTHZhhVM8AMBcIXR+bA1O1o7DoTb8iPh3DR1OOFyrIUZDGGySUzmXVs3ts/1MHkJNnj0BcZYLSJAGtg6l6AXVLhRSepVAKvpgKdsuni9PoaJpeAffHYPJocyxqj8dIFggVHsT93FcPCWXlgQlyAy3pi5Z2kbt7dT+busKlApyyDSK+vYXIN2OePzZXxki6ujgkWPJABRpsIsAamjQtw3knqu9/fD3p3GHM2rKn1XJZBpNh2mKoE7AQL/ureGMV8TkJ8CLB6qOggb+MCnHeSMsYEO+HHXiuRNVxjjDnL4jW5nDF33Uc1VQJ2hov9la3nonNr7OckxIcjsWdcDvKmL8B5J6lrWyvBTvh9qJWYrefHLq3JSkFqYTiJD1foGyOGi7NVXc9lx14fzkmICwFWz8RwkOedpB67vB7shN+n4Y+Q2ySG7Yvw6twY1elrl5o667ns2OvTOQlxoA9WB1z6sOS9xvV5gk3L6wkT6rlcMfb8aWObxLJ9EVbV/mV1+9qlpk6fuLJjL8ZzEuLGntGQvLtOlzusotfEXncRanjSZ/ijjcLTtrZJ7NsX1VTNfpDR9FMny1R27OWdkzZWRhS+IxNn7QYUXYxdTphFr0ml7sJ1+KOtmqW2tkkq2zc1VQNnhqX81LlBKTv2ss5JVzZXdHvngJpJZGKIsAFFaWqXIaCi15Q99X1IXIY/fIYE6jwkua1t0ub25aHS7an6sGCGpfzUeSizy7G3eE6aDd/y6Cpk4ShtQNFdp8sdVtlrmKb/gOsdft1MV5vbpI3ty0Ol21W1IJ2Mpp+2C//JMKIIGawGFN11utxh1bkLS43rHX7dhyRX2Sa+GaA2M0Y8VLp9VQrSU8pYh1K18L/KA9zJMKIIe0EDiu46Xe6wmH7tzvUOv41nNM7zzQBVfX3VwloeKt0fZKzbUWVCARlGFCGDVUFZpqHsrtPlDovp125c7/BD3Gn6bBPfDJBvLZnvnfaisvXBnTlSU6VFChlGFCHAmnIdnnG9uBEgNct3OK3tYVffk7XP60MM35WtD4apH0bR/7C53FTk7QMpHxfIx+2o/IZn6EvTvSoF2G0Pu/pOF/d5fYjhOx4q7cd1n+NhwP3l+pxCJn7AFRks+WUE6LTdvaoZnDbvNH0zQD6vDzV8t7g+ZheQ+btz7sxPuexzIYZu0Z2y4T4mfsBXrQyWMeanjTFfNcacGGO2F373MWPMS8aYrxljfqLeYjbLJyNAp+3u9aEA2zcD5PP6OoW1dZ4wkDKXfY7sdv8VTSjow3kHcambwfqKpL8i6X+b/6Ex5l2SPijphyQ9Jumzxpg/Y62N8mztkxHosjaFGpBToQuw667XvL/3zZi5vr5qYW1RhoW782Iu+xzZ7WFj4gd81dozrLUvWmu/lvGrZyR9ylp7YK39pqSXJL23zmc1yScj0FVtChmGB0JOja67XrvaLlWm7hdlWLg7L+ayz5Hdjl+dGjlaMsBXUzVYj0v6wty/X53+LEqzi9XNO3t67c6BrLV6y59aK329qxCZJzIMD4ScGl13vfZpuxRlWNq+O485G1s0U6xon2PmZdxca+TKMtK0ZICr0gDLGPNZSW/N+NVz1tpP5/1Zxs9sxs9kjHlW0rOS9MQTT5QtTqOsTrNRs5NjiALVUIWvDD+cF6pgve567dN2qfuEgVBiLgYvWzaXSRRcgOMUYqLC/D4Q800C4lAaYFlr31/hfV+V9Pa5f1+XdDPn/Z+X9LwkbW9vZwZhbaiaiSg7yEJlOBj/b0bd9dqn7VIURNUNDnwuNjFn/eouGz2R4uVyM+S6/WO+SUA8mroKvCDpg8aYFWPMOyQ9LemLDX1WEGU1KFlj9y71N6FqWxj/b0bd9dqn7eLa+8r3kSy+dWhd13sV1eF0vWxojkuNnOv2Z8YoXNRt0/CXjTGvSvrPJf2WMeZ3JMla+1VJvy7p30v6PyX9QqwzCGeKDr68C8jNu/ulB1mowlcaPzaj7nrt23Zp4rl2vhebLovBy4JBCtWHK+REBQJxuDDWdjYq95Dt7W1748aNxt6/aBhjPuW7OHxye+dAJ1ZnaWNJOpyc6NU3dvWOq5sy5sGBZq3V3uGxnrq6Wfq+sV6EUU9qtRkv39rR+vJS4XEwr8tjYhZcLR7Ls8CY43XYyo5N1+1fth8hHcaYL1lrt7N+l8xtWdmda1EmIu9uxRhTerfTtwwH6ulzO42qU9h9sz5dHhNlmQeO12Ery+C6bv8+lQagO8lksOrcceT97dHkWFbibhdn+npnWydzU/dv28z29XX7ID6pZaqRrSiDlcyzCOtMpy+bwl4084qDMC19atswr+4s2sPJie7uH2p9eayLq6Nzx0HeMdDFTCx6VaGKot5oKeK65iaZIcI6xatFaeOilHOfh4uQr2gora9F0lWKduf378sby7q2tablJZNZ25h1DHQxE4shQPjiPH4e68NdMhksnzvXoruV2e9u3t3vdb8fVFOWdelrhqRKPy+X/bvoNSGyfS5Fy2QeUAfn8fPK1gfZrQfivq0OyPXOtSg676rfT53nZyGssqxL7BmSvH2pStGuy/5d9Jq62b6y45E7bYRAS4bzitYHx9x5yWSwJLcuy0XR+em/3e9kQnT5pmNwXFyyLrFmSFweA+PTyd1l/27y8Txld9I+mQfuupGnT09rCKXoeChaH2T7zhvuHlJRUXTueycTYiovHYPj0tcaK8k9++bahNRl/y56Td1sX9nx6Hq8cteNIqm1ZCg7Hq5sruje/qG+/cau/uR7O/r2G7u6t3+oK5srZPsWxH9VaFnRBbSpfj88uqM/YjvZ+gwfh96XXPbvph7PI5UHu67HKzcxKBL7sH9oLseDNUaaHVp2+m/1+wa0CUkNEbooG7bwHdIoGy4qG7ZJMT0ds7oPRQ7Jd/i4iX3JZTi0qSHTsmPVdQiyr6010J5Yh/2bUHY83N450MXVsR6du6k8nJzo9s5Bbyf5NIWr9AKXlgwh72TK7hZiy5igmef5VVG27yxmtzZWRoPal1yzY2XHK3fdiFnbk5zKjoeiTHhq2b4ynEEyFF1AZ7977NKaJOnm3f1aOz2P7kBVvrN5ZneYQ9qXXB59cmVzRSujCzqY3mUvHqvcxCBWXdQHlh0P3JC4Y41M+dwlhNzpXXbWWDImiEvRvpOX3do9mCS1L7kcq9zEIFYxNuMtCsCYMHIeNVjyr2UJ2WiNMet+imFaf9G+c/PuPnVFcm8SmVKNTYxiOJ5i1EYz3ix5x0PZo7Fm11HaNJwigyX/u4SQjdZivXumuWn+OojlLq1o3yGNf4pZuPGL5XiKUdZxvHN/ojt7B62Ptrg8Govj7by0zrY5fHeKKkMzRSnd2IYAOeHF9wy9PHn7DnVFpwg04xfT8dS2shvZxeP43v6RXnljVxdXl53OzSHXrct7cbydl+a3XuC7UxRdvIYQwad8wpspWgd92MaxZkbb5htokrltXx+OpyZUqQ/8/v1DXX9kQxfXxrVHW3y5vBc3ducRYMl/pxj60EyqJ7x5TT5Dry2xZUa74BNokrntRl+Op9Bcb2Tnj+PL6yvaWj1fOl11tMWXz4Ss1G/sZihyV3bzyCubK4WFgXlFgEMoWqe5abPP0EO7XAvYeY5aN1I9nqoUsPuem+us28Xi+I2V0VnwV/ReTBh5IJ0rZon5u4RZcFXlTnYIETxp3mafoYd2uQ77kbntRqrHU5XskkuLhPn9vOq6baOPXgrD8cZaW/6qlmxvb9sbN250vRhnO9by3I5+ODk526FSwLRp1sEQzLdgWbzrXtyWHPfI0tR5wGffLFseSZXeK0/Tx0LV7x4jY8yXrLXbWb9jiDADzyZLM82bdeJKbR0Mjc+wX6pDVcjn2yPRR95zTaXTAMenPCV0/6mmr4GpDMczRJgh1aLLlFHgPEw+w37zwyl39o70+r39s4fYsh+kqekZ1YsTUSRVOg+FHt72vQb6DvelMhxPxJCBGqT00JpimHwvFLNnFy4vGV3bWtPljWWC7YS1HQhUPQ+FTgr4XAOr3JymksQY1rcJJNWiy5SlckeVmio3SwTbmGk7EKh6HgqdFPC5BlY5XlJJYlCDlSPFGqQsqRR605pimPLqXIr2YWowMROiLs/nHFr1PFRlPy/jeg2scrw0sbwxIsBCriYLPGNDgfNw+d4sEWwPn2vQUzcQ8D2H1jkPdZUUqBsUDhkBFnJVnenRx6yX64m0j98tVVW3VarBdpv7dpfHkW/QUycQ8D2H9jGzU3a8pHzO5JYMuarUA/R5Nl7Zo2X6/N1SU2dbpViDGWrfdplN1vVx1GaNXZVzaIyPuCrarkXHS9fbumtksJCrSup3yP1NhvzdhsZnW+XdYae0TUPs266Zoa6PozZr7FzOobFneFy2a97x0vW27hoZLOSqMtNjyLPxhvzdhsZ1W6V+hz0TYt92zQx1fRy1OTOw7BwaYv9r+pEzdTJ+XW/rrhFgIVeVoZIh9zeJ+bul8FwvH67bipYMp0Ls264XU9fPqrNPF/1tmy0Cys6hdfe/Nm4Q6gRJMZ8z25DGt0RlvvUAQ+5vEut3IwvzMNdtlfod9kyIfdv1YuryWXX26bK/bbvGrugcWnf/a+MGoU6QFOs5sy0EWAhqyAXCsX43sjAPc91Wqd9hzzPG6NU3dvUn39vVUcaDd8sySq4X01m3/Nfv7etr//FNvX5v/6G6ozr7tMvfxlJIXnf/a+MGoU6QFOs5sy0UuSO4IRcIx/jdaIyZzWVbpdqSYd58EfM7rm6erYO815QVOru2Orm2tabHL52u89s7B1odL529ts4+3afjoe7+10bPtqLt6lKgH+M5sy0EWKgk9pkvKaExZnV97DsUmstML9fZYPMX07xzhMt71dmn+3Q81N3/qgRoVc7dWUFSSo2oqyLAgjcOrLiQhamnyzvsKhe70Dc3Lhkf36xQ0TnC5b3q7NN9Ox6q7H/z+4CRdHR8oqNjlQZoIc/dqbdgcBFfSI/W+c7WKatxYEZbu1Kvc+irKoXcTUxocKkD8q0VKjpHuLxXnX166MfD4j4wHi3JWqvHLq2Vfs+Q9ZpMEClHBitxVe5oiu5AyW51w/cuuK0hXoaS81XJALj8je86d8n4+GaFis4Rj11ac3qvOpnFIdf91MkchaxP69NQbFdYE4lzuaNZzEhJyr0DZUZb/Npq60D7iGJVMgBlf1NlnbtkfHyzQkVZKpdHq5D9zhdLX6rUWzC4IIOVuLI7mqyM1P3DifaN0cXV8UN3oDfv7vdmBk+qiu6Ar2yuBMs4UaNRrEoGoOxvqq5zl4yPT1aoLOMVsmjaNWM3lGxqncxRyPo0JoiUI4OVuLI7mqyM1NbastZGFzLvQOkrFL+8O+A394+CZpyo0SiWlwHYWBlV7kIeyzqvkqWqkv12zdgNKZsaU1+qxX5ikshAziGDlbiyO5q8DNfRscm86+nbDJ4U5d0B7x1NtLm6Vpr9cM0EUKNRLCsDMMsg5mVxyrIGXa1z1wdm151duMg1YzekbKpP5qjNB5lTf/swznSJK7uj8c1IDX0GzxDk3QGvL49Lsx8+mQBqNMotZgB2Dya5WZzZur95d1+SMmeNdbHOffaJurMLF7lm7GLJ7IXi0om+7awd9bcPI8BKXFk2osoJO5bHUCBbXhB8cXVUeoHzOYkSbPurM3w7O5YPj61ev7evO7uHraxzn32iKNCpcq5xDcpSLF1oO+AZWhAbwnD3LpzJq3lwucPhIjlMWUGwywXO9ySaerDtOysuLxDYO8rPbM1/zomVLq+PdW3rdKi3jUJun32i6uzCPK5BWYrZ1LYDnhSD2DLpfvNEFAVRrnc4qV8kU+FygeMk6q7KEE3V4dsuh2d89omyQMf3XOMalKV4o9j2sZpiEFuGIveBKyru7NNDUdGOsgJYJjG4q1JYnVfAfHvnoLB4vctj2WefaGJqv2vR9pCbj2Zp+1ilbcPDCLAGrujEWzbjKLX+MjGKad2e1fhMTnR3/1Dry2NdXB0lfxLNEzLoKbtYdjlj0/fCmlqg05UuAh627Xnk9QeuKE1clNJNsb9MbGJat+dqfDaWT2t8lgzBdIEqQzR521xS4RCXy7HcZG8iygji4jLjFM0jwBq4ohNvUV2Ca00HU3ObE9O6jWlZ+qJKTUrRei4KYvKOZUm1gnQeXdM/Md2Ypa7WEKEx5lck/deSDiV9Q9LPWWvvTn/3MUkflnQs6X+w1v5OzWVFgbyhpLI0cV5K13V4o+zBz7EMb/VRm3U1ZduKej1/VYZoXNZz2bE+b9b4sUqDTRpH9tOQmqr2Xd0M1mck/bC19t2S/oOkj0mSMeZdkj4o6YckfUDSPzbGcEQ2pOyOpUr6vm5/GVnLXVRNbc0CcrnjrbssXWVCyj636eXyPfbK1rNvdqLOVH2ylv1EP6p41DpTW2t/11o7mf7zC5KuT//7GUmfstYeWGu/KeklSe+t81nI18SJsG5/GRnDybmmtqY9u+w/dZYl1JCFbzBU9rkxDqWUrWffY71OYMyFup9opRKPkGv8r0n67el/Py7plbnfvTr9GRrQxImwbn+Z2TKEXKbUtNW7x2X/qbMsIW4AqgRDZZ8bY4ambD37Hut1AmMu1P1EP6p4lNZgGWM+K+mtGb96zlr76elrnpM0kfTJ2Z9lvN5m/EzGmGclPStJTzzxhMMiY1FTU7TrTLnlQb9htDHt2XVbVV2WEPVbVepKyj431rqyovXse1zVmapPz7N+oh9VPEoDLGvt+4t+b4z5kKS/KOnHrbWzIOpVSW+fe9l1STdz3v95Sc9L0vb2dmYQhmJVToShCtDzCmGvbK6cZQKGfnJuq5i/qc9x2X/qfHaIYLtKMFT2uX28CahyrFcNjLlQ9xf9qOJQ60xijPmApI9I+ilr7d7cr16Q9EFjzIox5h2Snpb0xTqfhXy+wzcha0/yhll2DyZJPJqirTqeJj+nbP+p+9khejRVGa4qGyrpundUFW0/8oX+VkB15kHSqcIfG/OSpBVJt6c/+oK19uenv3tOp3VZE0m/ZK397ex3eWB7e9veuHGj8vLAzexiuTx3cTqcnJyrn3L18q0drS8vyZgHWQBrrfYOj/XU1c2znw21ZUPIdRnD5zT12Vnbf/bey6Olh7Ixi/vGfKa07LVlnzv/+rrLBSBtxpgvWWu3s35Xqw+WtfY/KfjdxyV9vM77oxkha09chlmG3E8ndB1PXkDQZb1QiM+u26Op6nBV2VBJ6N5RADDDswgTFLL2xKUmZMiN70I9z3H22rxAtMt6oaY+2zdwC11XEmMwC2A44q3mRGNcp/G61qEYnd71f/PWjo6OTx7KLAy5n06I5znOFLUN6HLqdVOf3WUbgKJtQ3sCACFwxkiQS6GsS3Awe814tKQnH93Q9Uc2lFXTN+QLVojnOc4UBaJtFze7fsc6ugwaYw1mAQwHQ4SJKhtucRnWcx36G3o/nbrPc5wpG4rrcup1E5/dZRuAom1DewIAIRBgIZNLcOAaQKR6wfKtXRp6IJqlq6Ax5mAWwDAQYCGTS3DgE0CkcsGaL5yWpPuHE22tLTsFTKkGol1IMZgF0C4CLGRyuQBxkTovaxbgvjE6mhzr6NgUBkyLM9oeu7Q2yMAqln5oBLMAmkaA1RNtX5hcLkBcpM7Lqkm7uDoubchZp09YLAGLi9j6oaWSVQXQDQKsHujqwuRyAeIi9UDV/klV+4TFFrCUGXI/NABY1P958gnwne6PblRtR1G1T1jf9ouu+qHF+lxBAMNGBqsH2u4s7Trs1KfhqTZUrUlznSywuL7f3D/S5Y3lc6+JueN4F93o+5blAzAcZLB6oM1Gna7dx327lKegakNOl8aWWev7zt6hdu5Pzr1XzA1cm2zgmZelCpHlIwMGoIo4z8Q4p83O0q4XpL4NT7VlFmQ9dXXTOUviEphlre+rW6u6de9+bzqON9URvijYrzssyY0EgKoYIuyBNmfruQ5H8kDcdmWt763VkSYnK2cBSx9mcTYxKaKoeL7usCSF+QCqIsDqidAXprz6KdcLUhf1NEPlUieUt74vro6iutB3UZdXFOw/dmmtVq82biQAVMXVMEFFwx6uw5E8uWjVpAAADyJJREFUENeNS/2Oy3BrH9Z3V8NpRTWKdYclh/ygcgDN4iyRoKILuusFqal6miFxDThc6oT6sL67qssrCz6r1MW5vjcA5GGIMEFlwx6uw5E0GS3mWr/jOtwa+/ruajityRpFnlYAoCoCrASFqJ+iB1Y514Ajhmc6htieXdblNRl8xh7YAogTQ4Q9E6InT91hjyq1Nin2EnKt3+l6+C9U7RTDaQDwAAFWj4S6ENa9oPvW2qTaS8gn4PCpEwodrFbdnouf33WgCAAxIcDqkZBFxHUKf32bN6balLSJgKOJYNVne5Z9fp39CgCGhBqsHomlJ49vrU0sy92Uovoll/odn/qnJhpf+mxPGm8CgBsyWD0SS08e31qbWJa7CXUzSr5/X/fRL1l8tmcTnw8AQ9T/K9yALda6bKyMoigi9h36GnLxc93hT9+/byJY9dmeQw6WASAkhggjlfX4lNs7B7qyuaLdg0nnPXl8pq4PuZdQ3eFP379vqqWD6/aMoaUEAPQBAVak8mpddg8mUV/M8uqJhtpLqG7vJ9+/7zpYLft8+qMBwCkCrEj1sTDc5aHF868dwoW4bkanyt/XDVbrrvu8z/fZ/gAwdBRORKqPtS6u9URD6otVtxVD272jmlz3qbbjAIAsZLAi1cdaF9es29Cm+tfNKLU5fNrkuu9j1hUAmhJvOiRxrpmNmB5B45p1Y6p/d5pc933MugJAU8hgRawssxFbzYtr1q3LhwJn6Vs9WJ3lbXLd9zHrCgBN4dYyIkXZqKzfxVbz4pp1i6kvVt/qweoub5PrnmcRAsADZLAiUZSNkpT5u8PJiS5vLJ97n65rXlzqibpuNTCvb/VgdZfXdd1XzZINtR0HAPgiwPLkeuHxvUAVXThP//3w7+7uH2rzeBzNUJuPqhfi0MN5fSvMDrG8dYee+zakCgBdiP9KHBHX4ZkqwzhFxcd5v1tfHkcz1NaGJobz+laYHWp5i4aji4ae+zakCgBdifMqEinXmqcqtVFFF868311cHSVV89JEzVlM9WAuQixvWZBUFOzHVvcHALEiwPLgOsW9ylT4ogtn0e9mwz1PXd0cdHAlNdNioG+F2SGWtyxIKgr2abEBAG6owfLgOsW9ylT4+eLjO7uH2juaaH15fPaA56LfxRoMhOa6Xn1rhPpWmJ21vD7fuayOq6jdwu2dg6habABArDgrenAdnqk6jLM6XtKVzRUtjy7o2taaLq+Pz4ZvZu+b9btU6l9c1uuQa4Ty6qZ8v3NZHVdRlqxvQ6oA0BVjrS1/VUu2t7ftjRs3ul6MQnmZgsWfb6yMtHsw8Z5pNbtQLs9lBA4nJ7owTRjk/a5PGZg6yjI1Reuvz+tofmZfVlbJ5zsXvZfLPuqbIWTWIYChMsZ8yVq7nfU7hgg95Q3PLE5rn/Ul8r2QlA3f9KmlQBNmWZTZBXtxmLRvbRdcFbXx8P3OdfuQ+Qypxva0AQBoCwFWACGbVZbVGaVe/1J2wY7tMTyhFAVRdWr+mtZ0I1eyYwBi1e+rTiRCzqyqOpswFWUz4Ia6jorqpmL+zqGOjaz6syHX2wHoPwKsAFyaPxY1dpxXVGDct5YCTSi7YA91Hbm06ojxO4dojJoXSN2cZjLpyQUgRgwRBlA0rV3yr0MpGr6Z/W42NHLz7v5gh0ayhn9chsP61nbBRVndVAzfOWt7lR0bLvKGGV+7s6cnH90499oh1NsBGAYyWDXMAqebd/dlJB0dn2RmEEJ3v/Z9ZE9Z1ixGed9xY2UU7XBY02JuKpu3vSTVzq7lZS2ttb16zBGAtJDBqigrK5U31b1sllfIB0P7Zs1iLRLO+467BxPnGXCxfrchKtsn62TX8rKWb/lTazqcnB5DVbNjANAUbvUq8slKFdWhhH4wtM/yVfnstrJiRd/RJZMzpALoPmQim3yETl792WOX1qKtPQMAAqyKfC4oRQXKoR8M7bN8vp/dZtBStzh6KA8l7kugGKKYfZHLEHzMw6YA0kaAVZHPBaVollfoB0P7LJ/vZ7cZtNRtPTCUhxL3JVAM3SpiMbAcj5ZkrT3LWhFIAYhdrQDLGPN3jTF/ZIz5sjHmd40xj01/bowx/9AY89L0938uzOLGw/eCknenXeXO32VafqggbF6bQYvLdywaOmsio9KFvgSKWdtrlqGtMrTZl8ASAPLUvdr8irX23dba90j6TUl/a/rzn5T09PR/z0r6JzU/Jzqheg/VeTD0LGDLupAtLt/R8YmMpJt3989e4/vZbQctRcM/ZUNnMTff9NGnQDFrn6w6tNmXwBIA8tQ6S1tr35z754ak2ZXgGUn/zJ76gqRLxpi31fmsGIWo/6gbqBUFGrP3fuzSmqy1Gk9nFFadQh9T0FKW4Yi5+aaPmNa5j7oZqD4FlgCQpXabBmPMxyX9rKTvS/qvpj9+XNIrcy97dfqz79T9vCGq0yTSpWVDqCn0dR8SHJLLA45jaL5ZV0zr3EfR9nFpn3Flc0Uvv35Pu0fHOjmxunDBaGO8pKeubbX5NQCgstIAyxjzWUlvzfjVc9baT1trn5P0nDHmY5J+UdLflmQyXm8zfiZjzLM6HUbUE0884brcmHIJNFxe4yqWoGWoD3XOEss6z5IXLOVtH0nOTzWwxjw4a9jpvwGgJ0oDLGvt+x3f659L+i2dBlivSnr73O+uS7qZ8/7PS3pekra3tzODMORzCTSGGIyEeAQL6ilqZpu3fYxUmnGd/fvi6liPzg2FHk5OHnodAMSq7izCp+f++VOS/nj63y9I+tnpbML3Sfq+tZbhwQa41Oj0tY6nyFBqrGJQtZFpUZ1V3vaRMU7F6xS5A+i7ujVYv2yMeaekE0nfkvTz05//G0l/QdJLkvYk/VzNz0EOlxqdvtbxlIl56KwvfB9EPq9s6Dlr+7hmU4eYdQWQlloBlrX2v8n5uZX0C3XeG+5cAg2CEWRxmSSRp0oQ5Dq0yxAwgL7jdhBIWJ2huCpDz65DuwwBA+i72m0aAPRXnaG4qkPPrtlUsq4A+owAC0hYlaE4lz5WAJA6hgiBhPkOxZU9oijrtVWeRQgAfUcGC0icz1Cca1F8ndmJADAEZLAAOHMtiq/7LEIA6DsCLADOXB/CTKNQAKkjwALgzLU1g2sgBgBDxdkOgDPXovghPp4JAHxQ5A7AS1FR/HwLByPp6PhER8cazOOZAMAVARaAILJmDs56ahFYAUgNQ4QAgmDmIAA8QIAFIAhmDgLAAwRYAIJg5iAAPMCZD0AQzBwEgAcIsAAE4ftcQwAYMmYRAgjG57mGADBkZLAAAAACI8ACAAAIjAALAAAgMAIsAACAwAiwAAAAAiPAAgAACIwACwAAIDACLAAAgMAIsAAAAAIjwAIAAAiMAAsAACAwAiwAAIDACLAAAAACI8ACAAAIjAALAAAgMGOt7XoZzhhjbkn6VtfLEZlHJX2v64WIGOunGOunGOunGOsnH+umWCrr509ba69m/SKqAAsPM8bcsNZud70csWL9FGP9FGP9FGP95GPdFGP9MEQIAAAQHAEWAABAYARY8Xu+6wWIHOunGOunGOunGOsnH+umWPLrhxosAACAwMhgAQAABEaAFSFjzN81xvyRMebLxpjfNcY8Nv25Mcb8Q2PMS9Pf/7mul7ULxphfMcb88XQd/O/GmEtzv/vYdP18zRjzE10uZ1eMMT9tjPmqMebEGLO98Lvk148kGWM+MF0HLxljPtr18nTNGPNrxpjXjTFfmfvZI8aYzxhjvj79/8tdLmOXjDFvN8Z83hjz4vTY+uvTn7OOJBljVo0xXzTG/Nvp+vk705+/wxjzB9P18y+NMctdL2ubCLDi9CvW2ndba98j6Tcl/a3pz39S0tPT/z0r6Z90tHxd+4ykH7bWvlvSf5D0MUkyxrxL0gcl/ZCkD0j6x8aYpc6WsjtfkfRXJP3+/A9ZP6em3/kf6fR4epekn5mum5T9U53uE/M+Kulz1tqnJX1u+u9UTST9DWvtD0p6n6RfmO4zrKNTB5J+zFr7ZyW9R9IHjDHvk/T3JP2D6fq5I+nDHS5j6wiwImStfXPunxuSZoVyz0j6Z/bUFyRdMsa8rfUF7Ji19nettZPpP78g6fr0v5+R9Clr7YG19puSXpL03i6WsUvW2hettV/L+BXr59R7Jb1krX3ZWnso6VM6XTfJstb+vqQ3Fn78jKRPTP/7E5L+UqsLFRFr7Xestf/v9L/vSXpR0uNiHUmSpteknek/x9P/WUk/JulfTX+e3PohwIqUMebjxphXJP23epDBelzSK3Mve3X6s5T9NUm/Pf1v1k8x1s8p1oObt1hrvyOdBhiSrnW8PFEwxjwp6Uck/YFYR2eMMUvGmC9Lel2nowzfkHR37mY4ueOMAKsjxpjPGmO+kvG/ZyTJWvuctfbtkj4p6Rdnf5bxVoOcBlq2fqaveU6nqftPzn6U8VbJrp+sP8v42SDXTwnWAyoxxmxK+g1Jv7Qw0pA8a+3xtKzluk6zxD+Y9bJ2l6pbo64XIFXW2vc7vvSfS/otSX9bp3cAb5/73XVJNwMvWhTK1o8x5kOS/qKkH7cPeo2wfools35KsB7cfNcY8zZr7XempQivd71AXTLGjHUaXH3SWvuvpz9mHS2w1t41xvyeTmvVLhljRtMsVnLHGRmsCBljnp77509J+uPpf78g6WenswnfJ+n7s/R0SowxH5D0EUk/Za3dm/vVC5I+aIxZMca8Q6eTAb7YxTJGivVz6g8lPT2d4bSs08L/Fzpephi9IOlD0//+kKRPd7gsnTLGGEm/KulFa+3fn/sV60iSMebqbDa3MWZN0vt1Wqf2eUl/dfqy5NYPjUYjZIz5DUnvlHQi6VuSft5a+9r0IP9fdTrbZ0/Sz1lrb3S3pN0wxrwkaUXS7emPvmCt/fnp757TaV3WRKdp/N/OfpfhMsb8ZUn/i6Srku5K+rK19iemv0t+/UiSMeYvSPqfJS1J+jVr7cc7XqROGWP+haQflfSopO/qNGP+f0j6dUlPSPq2pJ+21i4WwifBGPNfSvq/JP07nZ6XJelv6rQOK/l1ZIx5t06L2Jd0mrj5dWvt/2SMeUqnk0gekfT/SfrvrLUH3S1puwiwAAAAAmOIEAAAIDACLAAAgMAIsAAAAAIjwAIAAAiMAAsAACAwAiwAAIDACLAAAAACI8ACAAAI7P8Hn5kEYRb948sAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 720x576 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "FS = (10, 8)\n",
    "fig, ax = plt.subplots(figsize=FS)\n",
    "# Make points translucent so we can visually identify regions with a high density of overlapping points\n",
    "ax.scatter(x_axis, y_axis, alpha=.1)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 56,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAA28AAAIICAYAAAD5SdaqAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nO3dbaxl133f99+6D4dzOeKtxIh2Zc1cSm3ZwLLhsYMLWUX7Ih0LsGIwppNGgOSyIzgCyAHkJlHzwlKGqMFmBnDgINO6cUoSsBENRrYixGkk0HZtSVTg9oUsD1ONI0WWRdvQzESCxZSW75i8w3PPvasvzjkz9+Gce/Z62nuttb8fQBDn3Iez79lr773+a/3XfxlrrQAAAAAAeVvq+gAAAAAAAIsRvAEAAABAAQjeAAAAAKAABG8AAAAAUACCNwAAAAAoAMEbAAAAABRgpesD2O/Nb36zfdvb3tb1YQAAAABAJ1588cX/aK19aNbXsgre3va2t+natWtdHwYAAAAAdMIY8415XyNtEgAAAAAKQPAGAAAAAAUgeAMAAACAAhC8AQAAAEABCN4AAAAAoAAEbwAAAABQAII3AAAAACgAwRsAAAAAFIDgDQAAAAAKQPAGAAAAAAUgeAMAAACAAhC8AQAAAEABCN4AAAAAoAAEbwAAAABQAII3AAAAACgAwRsAAAAAFIDgDQAAAAAKsNL1AQAAAADIz87unraHI432pJUlaW2wotVl5n66xKcPAAAA4ICd3T3dvrMjK6PBypKsjG7f2dHO7l7Xh9ZrzLwBAAAAOGB7ONLy0pKWl4wkTf5/SdvDkVbXBke+n1m6dvCJAgAAAD2ys7unre2hXnl1qK3t4czZtNGe7gZuU8tLRqMZE2/M0rWH4A0AAADoiaaB1sqStLtnD7y2u2e1MiN6mDVLt7w0nqVDXARvAAAAQE80DbTWBiva3du7G8Dt7lnt7u1pbXB01ZXLLB3CELwBAAAAPdE00FpdXtIDJ1ZlZDUc7cnI6oETqzPXsbnM0iEMBUsAAACAnpgGWvsDuHmB1ury0sziJIetDVZ0+86OpPGM3nSW7oETqxGPHBIzbwAAAEBvuKRDTi0qcOIyS4cwzLwBAAAAFTmubP800NoejjQcjWfcjgu0pgVOlpeWNFgZz6rdvrNz5GeaztIhDMEbAAAAUIkmwZZLoOW639v+42Dft/j4BAEAAIBKuJTtj73f2/7fy75vaRC8AQAAAJVoGmyl2O9tin3f0iF4AwAAACrRNNhKsd/bFPu+pUPwBgAAAFSiabCVYr+3KfZ9S4ePEAAAAKhE02CrSYA1XRN3+854Nu6BEytaXxssLDziM1uHZvgEAQAAgIo0qSa5aGPtplsEzHt/l+0I0BzBGwAAANAziwIs1y0CZm0NsM6+b9ERvAEAAAAF891T7bgZutGeNFg5uiZuOLJHvjdklg5u+DQBAACAQqXaU82l6AhbA7QnOHgzxpwwxnzRGHPdGPMVY8zTk9ffboz5XWPM140x/8IYw7wpAAAAEFGqwMml6AhbA7Qnxszb65LOWmvPSPpBSe8xxrxL0j+SdNla+4ikP5P0wQjvBQAAAGDCJXCaVo985dWhtraHx87OuWwRwNYA7Qn+SO3YX0z+uTr5n5V0VtK/nLz+MUk/EfpeAAAAAO5pGjj5pFeuLi9pfW2gB08Ojt0igK0B2hMlHjbGLBtjviTp25I+I+mPJH3HWjudr70l6a1zfvYJY8w1Y8y1l19+OcbhAAAAAL3QNHBKuS7NZyNv+InyiVprd621PyjplKR3SvreWd8252efs9ZuWms3H3rooRiHAwAAAPRC08CpaXqlS2rl4eNoMkuHMFHnMq213zHG/BtJ75L0RmPMymT27ZSkb8Z8LwAAAADNNuWeplfuD+AOp1dS8j9/MapNPmSMeePkv9ckvVvSVyV9XtLfmnzbByR9KvS9AAAAALhrkl5Jyf/8xQih3yLp88aY35f0e5I+Y619XtLPSPqfjDEvSfpLkn4pwnsBAAAAcNQkvZKS//kLTpu01v6+pB+a8fofa7z+DQAAAEDHFqVXNkmtRLc4FQAAAAAo+V8AzgQAAABQgJ3dPW0PRxrtjWfJ1gYrUQuJTFMrt4cjDUfjGTeKleSF4A0AAADIXFuVIJtUrkR3CKMBAACAzFEJEhIzbwAAAEAysVIdR3vSYOVoJcjhyMY6VBSA4A0AAABIIGaqY9NKkKnXxaFbnEkAAAAggZipjk0qQU6DRSujwcqSrIxu39nRzi4btdWC4A0AAABIIOam10022WZdXP1ImwQAAAASiL3p9aJKkKyLqx/BGwAAAJDA2mBFt+/sSBrPhk1THR84ser8u5qsZYsdLCI/nEoAAAAggSapjk00Xcs2a13cnZ2Rdnb39MqrQ21tD1n/Vjhm3gAAAIBEYmx6PWstmzRey7b/d0+Dxe3hSMORlbXjQG1lefnuzF+Kjb3RHoI3AAAAIGMua9n2B4tb20OtLJuFQR/KQcgNAAAAZGy6lm2/JmvZYla7RB4I3gAAAICMNdnjbRbfoA/5Im0SAACgR5pULUReDq9lW1nSkXVrs85rzGqXyANXKgAAQE80rVqI/KwuL2l9baAHTw60vjY4ErjNOq+SolS7RD6YeQMAAOiJplULm2AGLx/Hndf1tQHFSSrCFQYAANATsQpYMIOXFwqT9AfBGwAAQE/EKmAxa6ZneWk804P2UZikPzilAAAAPeFbtfAwZnryEuu8In8EbwAAAD0xrVoYWsCCmZ68xDqvyB/hOAAAQI+sLi81KmBxXEESStDnp+l5RdkIxwEAAHDAooIkzPQA3WDmDQAAAAc02VKAmZ702I4Bh3H2AQAAKnTl+lVtXH5YS08vaePyw7py/Wrjn6UgSffYjgGzMPMGAABQmSvXr+r8809qe/SaJOnm1g2df/5JSdK5M48v/PlpQZL9AVxuBUl8Z6VKmc2KuaE66pFfSwUAAECQp164cDdwm9oevaanXrjQ6OdzLz3vOytV0mwWs5+YheANAACgMre2bjq9fljuBUl8NwkvaXNxtmPALJx+AACAypxaP+30+iyry0taXxvowZMDra8NsgncJP9ZqZJms3Kf/UQ38rkKAQAAEMXFs5e0trJ24LW1lTVdPHupoyOKy3dWqqTZrNxnP9ENzj4AAEBl3v/9P6l//O5f1FsfOC0jo7c+cFr/+N2/qPd//092fWhR+M5KlTablfPsJ7phrLWLv6slm5ub9tq1a10fBgAAQNG2toeyMkeqRRpZrVdSqXBe1chF1SRLqTaJ/jLGvGit3Zz1tTyHGQAAAOBttCcNVo6u7RqO4g3adx0EzdokfFpNcnlpSYMVo909q9t3dg6kG7K5OErGMAMAAEBlUq/tyrXkfknVJAEfBG8AAACVSb22K9cgqaRqkoAP0iYBAAAqM61UuD0caTgaz7jFrFQYMy0zZvrldMbx8Fq/HKtJAj4I3gAAACqUcm1XrCCpyRo1F2uDFd2+syNpPCs4nXF84MSq8+8CcsQ4BAAAAJzESsuMnX7pujfazu6etraHeuXVoba2h52v2QMWYeYNAAAATvanZb76+q5Gu7taXV4eB10OaY8pqmI2nXGMPesHtIGWCQAAAGery0uT9WlGbzgx0P33rThXnUxdFfM4uRZdAY5D8AYAAAAvoQFQ6qqYx0ldmZKUTKRA8AYAAAAvoQGQ6xq1Ra5cv6qNyw9r6eklbVx+WFeuX537vSln/XLdBw/lY80bAAAAvMSoOjlvjZrrFgJXrl/V+eef1PboNUnSza0bOv/8k5Kkc2ceP/L9KStTzpqRlMYzkqkqgKIfmHkDAACAl1Rpjz4zV0+9cOFu4Da1PXpNT71wYeb3x57124/NwpEKM28AAADwkmozcJ+Zq1tbN51enx5/ipkwNgtHKjQhAAAAeFtdXtL62kAPnhxofW3Q2czVqfXTTq+n1GUhFtSN4A0AAABJuRQSkfyKiVw8e0lrK/cfeG1t5X5dPHvJ+7h9pUzJRL8R/gMAUBjXQg5Al1wLiUh+xUSmv+upFy7o1tZNnVo/rYtnL819j9RSpWSi34y1/jvYx7a5uWmvXbvW9WEAAJCtaSGH6Xqg/Z3amAHcletXs+kEo2wblx/Wza0bR14/vb6hGx/+xtyf68sgRV/+TjRnjHnRWrs562vMvAEAUBCXQg6+nUKfmRK0o8SOvk8hEakfM1f7B2MGK+PBmNt3dkixxFy0CgAACtK0kEPIJsGuJdfRjlI3fs6pkEhsO7t72toe6pVXh9raHjqfi1mDMctL48EYYBaCNwAACtK0kENIp9B3pgRpldrRz6mQyH6hgVeMYJr94OCK4A0AgII0LUEe0imseaakZDl09H0CnnNnHtczjz6r0+sbMjI6vb6hZx59ttMU3BiBV4xg2qeqJvqNNW8AABSk6abIIZsEXzx76cCaNymPmZJa+K5b63rj55D1WefOPJ7VekmfTcAPG+1Jg5WjwfRw1LwYoE9VTfQbcT0AABlpMrPRZFPkkE2Cc5wpqUXIjE/XGz+XmrY5S4xZzBizZuwHB1dsFQAAQCZibwNQYmXC2m1tD2VljsyeGVmtN5jx6fKcvvLqUIMZkclwtKcHT5ZVFTL0PEjtbduB/mGrAAAAChAjlWu/1eUlabByt7O/PRxJHQZwBJPhqXah5fNDzkHXaZsxxUhXbJrCDMQU3LqMMaeNMZ83xnzVGPMVY8zfnbz+oDHmM8aYr0/+/03hhwsAQL1iF6TIqbR8TsfiK7Q6odRtgYrQc9B12mZMsdIVm6QwAzHFaGEjSX/fWvu9kt4l6UPGmHdI+oikz1lrH5H0ucm/AQDAHLE79jmtUcrpWHzECj4XBUAxAsR5Qs9BbeuzCLxQouBWaq39lrX2307++7akr0p6q6THJH1s8m0fk/QToe8FAEDNYs9s5FBaPsdj8REr+DwuAEo9OxnjHBDwAN2KOs9tjHmbpB+S9LuSvtta+y1pHOAZY74r5nsBAFCb2GtoUq9Rclk/Vfp6qRhl4afmrVsLWfPY5FyUfg4ARNwqwBjzBkm/JunvWWu3HH7uCWPMNWPMtZdffjnW4QAAUKSYMxsp1yi5zhKVvl6qjbVqvjNjTc9F6ecAQKTgzRizqnHg9nFr7b+avPynxpi3TL7+FknfnvWz1trnrLWb1trNhx56KMbhAAAApV2j5JpGWPp6qTYCH98Asem5KP0cAIiQNmmMMZJ+SdJXrbX/ZN+XPi3pA5J+bvL/nwp9LwAA4GZRaXnf0vE+aYShZe671EZZeN/y9S7nouRzACDOmrf/WtL/IOnfGWO+NHntH2gctH3SGPNBSTckvTfCewEAgEj2bzI8WBkHC7fv7DQKSvq4fip14OMbIPbxXAB9FRy8WWv/H0lmzpd/JPT3AwCANEIKZLjOErFBdzM+AWKMDacBlIG7JgAAPRVSIGN7ONLOrtVf3BnqtddHx66fqmGD7uOk3JutCdayAf1BeSEAAHrKJ91uf6rlyftWtLu3fLdwx7xgIWSGL3chqacxZyNZywb0A0MyAAD0lE8FRZ/NqkvfoPs4vpt31z4bCSANZt4AAOgpnwIZPlUmSyqo4Tob5rt5dy6zkaxFBMpC8AYAQI+5ptv5BGKlFNTwSYH0DUx9g75Fx+8SiIWkfALoBlcmAABozCfVspSCGj4pkL6bd/tuyD2PTxqmb8ongO4w8wYAABrz3YushIIavhuPH/48Tqwua3s40u1jZsBiz0b6pGGmmP0DkBbBGwAAcFJCIObDNwVy/+fRNBXRNwiep/a1iADGCN4AAACFKxRnNsxlBixmEFzzWkQA9/TrrgwAAI6gbP1YjLV5XW2LUPNaRAD3MPMGAEDP5VK2Pgehs2FdpSKmXIvIrCyQD4I3AAB6rubCFW0HHl2mIqZYi8h2AkBeuOoAAOi52GXrc9FFOmhtqYhsJwDkhZk3AAB6rtbCFV2lg9ZUjbPmWVmgRGUOAwEAgGhqmy2a6qp4SE1qnZUFSsXMGwAAKHa26Lg1bexjFq7WWVmgVNy+AABAkRatafMpn4+Dap2VBUrF3QsAABRp0Zo23/L5OKjUWVmgRgRvAACgSE2KaRB4AKgJwRsAACiSz5o2NpwGUDLuVgAAoEiua9q62PcNAGIieAMAAEVyLabBhtMASkfaJAAAKJbLmracNpwmfROAD+4SAACgF3LZcJr0TQC+CN4AAKjAzu6etraHeuXVoba2hwQCM+Sy7xvpmwB8EbwBAFA4ZnKayWXD6dGeDlTIlMb/HnG6ACzAmjcAAAq3aLNq3JPDvm8+WxwAgMTMGwAAxWMmpyy5pG8CKA/BGwAAhculEAeaySV9E0B5GOIBAKBwa4MV3b6zI2mcOjmdyXngxGrXh4Y5ckjfBFAehngAACgcMzkA0A/MvAEAUAFmcgCgfgRvAACgd3Z297Q9HGm0N14zuDZYYaYSvUH7LxdnCQCADLHpdjrsi4c+o/2XjeANAIDM0LlKa9a+eMtL433xgNrR/stG8AYAQGboXKXFvnjoM9p/2QjeAADIDJ2rtNgXD31G+y8bpwkAgMzQuUprbbCi3b29u5/xdF+8tQF13FA/2n/ZeAwAAJAZOldpsS8e+mBe0SPaf9l4CgAAkJlp52p7ONJwNJ5xo3MVV8374lEGHtOiR8tLSxqsGO3uWd2+s3P3PlJz+68dwRsAABmicwUfizrt6IdZRY+kcdEj7itl4yoGAACoBJVKIVH0qGbMvAFAhUibAty5Xjc5XmejPWmwcrTTPhzZOT9xVI5/F9xMix7tD+AoelQHTiEAVIYNngF3rtdNrtdZaKXSXP8uuGmz6NG8wihIg+ANADIQ8+FH2lSe6ODkzfW6yfU6C+205/p3wU1bFSUJ9ttH2iQAzNBm2lDsAgMx0qZqEnouY7QFikjkz/W6yfU6C61UmuvfBXdtFD2iMEr7CN4A4JC2O9q+D795QQVrHe4JPZex2gIdnPy5Xjc5X2chnfac/y7kJ5dgv0/rNOv8qwAgQKy0oaZpcj5VwY5LVWGD53tCz2WstkDlt/y5Xje1Xme1/l0IM+95FrrGMtax9Sl1k+ANAA6J0dF2eZj4PPyOCyraWutQgtBzGSvoyqGDg+O5XjdNvr/EdY7cP3BY7oOFfVunyTAKABwSI23IJU1ubbCi23d2JI2/f/rwe+DE6tzfvyhVhQ2ex0LPZawUMp9zjPa5XjfHfX/J6xy5f9QjRjrhcc+z9bVB0BrLGHJJ3WxL3ncPAOhAjJFElxkbn5FuZnKaCT2XsUaVmc3on77NBiA/sdIJFz3PVpeXtL420IMnB1pfG7R+X+vb85CZNwA4JLRam+Q+Y+M60s1MTjOh5zJGW9j/u5jN6I++zQYgP7EKJeVexKZvz0OCNwCYIbSjnfphEjOoqF3ouSToqlfKCnW5d3hRv1gDCDkFR/Ou2T49DwneAPRWyo5bGw8TggrAn8+aNJd7Rk4dXvRTrAGENp5nTa6tRddsX56HBG8AeqmNYgJ9epjUok97BfWda0qZ6z2jb7MByE/MAYSUz7Om1xb7ZY4RvAHopZIfAgQYaZRcHbBvYlwDrillPveMeR1ermG0oZQBhKbXFutIxwjeAPRSrIdA252wmAFGDR3ImH9DyQF9n8S6BlxTymLeMxgkqE+u99MSMkCaXluu12yu5yRU+X8BAHiIUVo4VhlmF7HKj3dx7LHF/htibchd4sbMJYl1DbhuAxGrHDlbCNQn9f209ntK02vL5Zqt4Rk3D8EbgF6KsX+XaycsxgM4VoBRQwcy9t9QakDfN7GuAde992Lt+Rfr+JGPlPfTPtxTml5bLtdsDc+4eaIEb8aYXzbGfNsY8+V9rz1ojPmMMebrk/9/U4z3AoAYYmya7NIJi/UAjjX6H9qBzGEkOHYnuIuAHu5ibsjrsrlwrI3W+7ahcB+kDMi7vKe0dZ93ubaaXrM1D5LEulX8c0nvOfTaRyR9zlr7iKTPTf4NANlw6bjN4tIJ6yrVK8axH5bLSHDsTnDbAT38xLoGfITeM6Rujx9ppAzIY99TmgZkbd/nY1xb+9U8SBLlT7DW/o6kVw69/Jikj03++2OSfiLGewFALlw6YV2lesU49sPaGgle1MlI0QluM6CHn1jXQFdKP34c5XovcpnRinlPcQnIcs0iaPrZ1TxIkvIv+G5r7bckyVr7LWPMd836JmPME5KekKSNjY2EhwOgTbVWedrPpQxzrM1Sp+8bWj0spIR0G+Wam1Tky7EM9qJ9lfpwXcT6G4/7PSVU0DtO6cePg1zuRa7VRmPu1eZSUTfHsvwun12Oz4dYOg8/rbXPSXpOkjY3N/u1UQNQqT6Vwm7aCYv5AI7FtwMZMxCdp2knI6QTnCKQOq7D0IfrItbfGPuz6kPQjG41vRe5bkkSMwhxCcjauM+78vnsahwkSXkK/tQY8xZJmvz/txO+F4CM5Jpu0aWaUqXaSEdJvXYs5XqOeamXfbguYv2NMT+rXNZoApLfvS3WejCXFMwc0w5dP7scCmulkLLX8GlJH5j89wckfSrhewHICEUbZou9IPs4vg+tJj/XRiCaeu1YF4FUH66LWH9jzM+qD0EzyuFzb4sVhLgEZDkOOLp8djUP2kQJn40xvyrpr0p6szHmlqSflfRzkj5pjPmgpBuS3hvjvQDkL8d0iz7xTTlzXU+QMh0l9dox1/UcMdLu+nBdxPobY35WMdbukHaJWFxT6GOmELumYOaWdujy2bmmWJYkSvBmrX3/nC/9SIzfD6AsOa7v6hPfh1ZOD7vUa8dcgoNYnac+XBex/saYn1VoINiHtYpoj2sAFfu+nFtA5mJ1eUknVpf1Z6++rtdHVvetGL3p5H0zP7scC67E0nnBEgDtamME2bXyVikj2qUcq+9DK7eH3bxORozOTBsjuLPaS63Vz6ZiFVeIWaQhNBDMaVADdXAJoHK7L3dpZ3dPd3Z29YYTA/0nk2v5zs7u+PM8dG+oOdOB4A3okTZHkJs8nHIY0W4akOVwrLOOadax+z60SnnYxejMuAQHPu93XHtZr7zDH2tkP+bvCQkE6TyjS23cl0sZmHQZSKk50yG/MwMgmdwW7nd9PCVvWHrcsftWCcuxutgsIcVM9i/83x6OtDZYWVhAxuf9cmsvfRbaMWXjdYQILTaS+r5cUmEPl0JGORZciaX8vwBAY7lVu+v6eFw62F0f62HHHXvTh9bhToWkIh52vp0Z306Kz/vl1l76KkbHtJRBDeQnRvtLHYSUNNDkOpDSZoXnNnHnASpz3ChzbmlxXR9PyRuWLjr2RSlnJaf1+abB+a5d8nm/3NpLX8VYrxZz/R36JdZ6SdcUYpfZ5pLSgmtOhXTBnQeoyKJRvtxGkLs+npI3LA1N5SpptHUWnxHVkNkw1/fLrb30VawZ0FpH8JFWFzPwrrN9JaUF15wK6aJffy1QuUUd8txufF0fT8kbloYGB31M62uzk5Jbe+mrkjqmqE8X7c91YK60gSYGUkibBKrSJP0htz1eujyekjcsDU3lqiGtz7UQRdspNzm1l74izQpd6qL9uaZBkhZcHoI3oCIpOuSllBD21bSDnePnEBIclN6p9dm6gU5K/3DO0aUu2p9PP6CPA005PtObIngDKhK7Q57j3mZdqPFzKL1TG1J8pG+dlL5Lcc5L7vihXW3fc1z6AX1tx6U/0/M/QgALTUu+374zkrVWo93dKOtsSi9qEUutn0PJawf6uGYPeShpXyz0j8tWMX1tx6U/05l5Awp3dATJTEbZwkfQSiohnBKfQ35qWLOH9FLMLMQq/44y5DY71eR4msz25daO2/ycS3+m85gDCpdyBIlKbWN8DvkprUIa2pdqZoFZ3/6I0YammTGvvDrU1vYwqP3FbNM5teO2ZwFLf6YXcpgA5gm9AR/3YMm9gxzzoXic3D+HPqIUPxZJNbBVescPzYW2odhBScw2nVM7bjuNsfRnOrcaoHAhN+BFD5acO8htjtTl/Dn0Wclr9pBeqpmF0jt+aC60DcUOSmK26ZzacduzgKU/07nTAIULqTDZJOc91+p8befr5/o5AJgt1brI0iu1ornQNhR7bVXMNp1TO+5iDXPJz3SCN6BwITfgkhft+hx7bgvP0T3aRL1SbJ1CW+mX0DYUOyiJ3abbDGCOu35K33e0bdx1gAr4po/llPPuyvXY+1wWGbPRJuoWMzWKttJPoW3IJzXxuLXcpab7lbxEI0fMvAE9VvJol+ux51YWGd2jTdQv1swCbaW/QtqQa2ZMk82jS0z3y2WJRi2z5wRvQI/llPPuyvXYS04RRRpN20QtD/zS5PS5c/+AL5egpNZBghyunyaBcSkI3oDKLeoA7X+wTL/3dgadpSZcHops6lyPWJ36Jm2ipgd+SXL73Ll/YJEY96UcgpwUcrh+agqMue0AFXNZp1Hymo4m+73lVBYZ/mK20yZtou39hzCW2+fO/QPHiXVfKnkd+nFyuH5y2pQ8VOHNAcBxXDpAuXWWmmr60GRBdB1ittMmbaKmB35JcvvcuX/gOLHuSzkEOSnkcP3UFBiX3RoAHMslBSPHdI0maSguqRAlLvTGQbHb6aI2kUO6Tx/l+Llz/8A8se5LOa9DD00L7fr6KblA22HdtwYAybiMNOU2KtV0Ri23EXqk1XY7rXUkPHdtfu5N0q6B48S8L/lu/ZNSScsq5l3POcz+xVLeEQNozKUDlFsntWkaSuzOPB25vLXdTmt64JfE5XMPuWZL6pQiX6nvS10/l0pZVtFkP7lpYLw2WNH2cFTks56nD1Axlw5Qbp3UpjNqMR+adOTy11an/vB75jYS3gdNPvfQa7aUTinylvL5mcNzqZQMl6bXcw6faQjyPoDKueSZd52Tvl/TNS8x1wjUVEq4Zk3aaW6l5pFG6DWb41pflCnV8zOH51KOa1BnaXo95/CZhiB4A5All8XFsR6adOTq0fThnNNG0HAXes2W0ilFf+XwXCql2EfT6zmHzzQEtycAWYqZhtI0fS63oi3w1yTNp/TUGYRfs7mt9UX/LHo+5fBcym1ZxTxNr+ccPtMQhRwmgD6KsdbIpYNOR64eTR7OrHfqVow1iaHXbCmdUtSpyZ2bme4AACAASURBVPMpl+dSCWt/m17PuXymvso4SgDw5LoPXK577MBNkzQfn9QZ0izjiLUmMcY1m9NaX/RLk+cTzyU3Ta7n0j9TgjcAVXPtoNORq0OTh7PreieKoMQTs2AA1yxK1fT5RBuPr+TPlOANQNUoSNCtLmeqFj2cXRfhl16hLCelFwwAYuD5BB80DwBVKz23vWQhBUHa2JTWdb1TKXsdlaD0ggFADDGeT11v4I32cZsEUDUKEnTHtyBIG1Ugpx2e23fGx/LAiZWFi/AJOOJhUAUIfz5RMbefeOQAqN7hKlm/+uVf0cblh7X09JI2Lj+sK9evdn2IVfKdqUpdBdK3w0PAEc+sTuuJ1WVtD0fMIKBXQqo4UjG3n3jiAOiVK9ev6vzzT2p79Jok6ebWDZ1//klJ0rkzj3d5aNXxXc+Rej2U79q10iuUtcFljeP+NYkUg0HpfNf3hqwLZu1oP3FHBNArT71w4W7gNrU9ek1PvXChoyOql+9MVer0xJC1ayXsddSVkBQuZhBQMt+2H5r2SCp3P3F6ARTNdbH2ra2bTq/Dn+96jtTpiXR40ggJwCgGg5L5tv3QQQtSufuJRxWAYvmMWp5aP+30OsL4zFSlLjJDhyeNkACMgBol8237oYMWFOTqJ84uUABKAc/mM2p58ewlra3cf+C1tZX7dfHspaTHCjeLgr6Qa4IOTxohARgBNUrm2/ZjDFqQyt0/nGEgc5QCns9n1PLcmcf1zKPP6vT6hoyMTq9v6JlHn6VYSctCgq8Y1wQdnvhCAjACapTMt+2HXDMM6vaXsTafijSbm5v22rVrXR8GkJWt7aGszJGKfUZW68dUxuuDHD+bkMphfbG/suDykrnbYWnaWc/xvGOM9o++arPaZOg9FPkzxrxord2c9TXyEYDMUQp4vrXBim7f2ZF09AHWROyOJuXOm/Et1T/FNZGv/eX/gT7xbfs+Pxd6D0XZ6E0AmWMh/3whqVYp0lEpd95M6CJ9rok8hKZtkfYF+KE6a78x8wZkLmR2qQ8pTL6jnSlGLpkRasZ38+6p0BlXhAudZWaWGrjH9Vkdeg9F2TjNQOZ8Z5faLnRS2ih6ipFLZoSaCa0sWEtxi9Kumf1CZ5mZpUaJUlyzPs9qqrP2W1lPOqCnfCrjtdk5KrEiZopAiwdqMzGCr9KrRZZ4zewXOvhB2hdKk+qa9XlW1zKABT/0KIBKtZnCV+Li6RSpd9MH6vZwpOFoHAjyQJ3tuHTXPqT7lnjN7BeatkXaF0qT6pp1eVb34d6IxTjjQKXaTOErcRQ91chl6TNCXSt9RqqpEq+Z/UJnmZmlRml8rtkmaZZNn9U+98aSU7MxH70KoFJtdo5KXetFoJWfvqyFKvWamQod/CDtC6VxvWabBltNn9Wu98a+DIT1EUNcQKXaTOHLbb81lKsvFTtrqJgZuqcbe8KhJK7XbNM0y6bPatd7Y+mp2ZiP3hFQsenM0gMnxuM0t++MkqRO5LbfGspV+oxUU8w8AWVxvWZd0iybZIG43htLT83GfMy8AZVraz+lrvdbY/YuHyHnooYZqaaYeQLK4nLNxi7K43pvpChQvTiFQOVyX0MUY3SQ2bt8hJ4LZqQA1CD2unPXeyNFgerFGQQqF7qGKPWMVozRQXL78+FzLma1sXXOW6eYyQbCpFh37jLzx9Y19SJ4AyoXEhy1kXIZI02uL0UuSuB6LtpK60VzPueEYA84al6w5Xu9uP4cqdl1Sn5nNca8xxjzNWPMS8aYj6R+PwAHhaROtJFyGSNNri9FLkrgei5yT+vtI0qSo6/a2BfN93rhOsNU0q6NMWZZ0i9K+muS3iHp/caYd6R8TwAHhQRHbVWrCt1vjdz+fLieCyqi5cf1nBCAowZtBUe+1wvXGaZSj0u/U9JL1to/ttYOJX1C0mOJ3xPAIb7BUSkzWhS5yIfruSiljfUJJcnRR20FR77XC9cZplIPS79V0s19/74l6YcTvyeASEoq205ufz5czkVJbawvKEmOPmpr7bTv9cJ1hqnUp9zMeO3AVWCMecIYc80Yc+3ll19OfDgAXDCjhdRoY/mhJDn6qK0sAN/rhesMU6nP+C1Jp/f9+5Skb+7/Bmvtc5Kek6TNzU1KwwEdmVfFihktpEYbyw8lydE3bWUB+F4vqa4zKsWWJ3Xw9nuSHjHGvF3Sf5D0Pkk/mfg9gd4KKT9MuXYAvgjAUbo2ByF8r5fY1xnP/jIlDd6stSNjzE9L+i1Jy5J+2Vr7lZTvCfRVyE2YTa4BAH2X+yBE7Fkynv1lSp4oa639DUm/kfp9gL4LuQmzyTWAUKRfAemkmCXj2V8m7qpAJULKCFOuHUAINhAG0kqxlQHP/jJxeoBKhNyEqWIFIAQbCANppdjnjWd/mQjegEo0uQnv7O5pa3uoV14damt7eHdUnHLtAEKwgTCQVopZMp79ZSK0BiqxqFLWonz53BdqA2hf03VsbCAMpOWzlUGT65dnf3m4rQIVWV1e0vraQA+eHGh9bXDgJk1aEwAXLuvYSL8C0nKdJWMdar24qwI9QVUpAC5cKtiyUTeQnsssGdsA1IvgDegJ0poAuHAd8CH9qp/YIiJPDNjWi6sL6AnSmgC4oIw4FiE1L19cv/XiFAI94VNVal51SgD1Y8AHi7CWOl9cv/UieAN65LiCJocxogr0G2XEsQhbROSL67dehN9AD/isSWCxMwDWseE4rKXOG9dvnbi8gMr5zqAxogoA/eCbIk9qHtA+gjegcr5rEljsDAD1C0mRJzUPaB9DI0DlfMsFrw1WdPvOjqRx4DcdUX3gxGrCowUAtCk0Rb7t1Dy2JvDD51YPgjegcr5rEth0F0AX6GS2y3eAz/U8xTiv01nC5aUlDVbGg4q37+zwbFqAz60unDGgciFrElyqUwKoQ5dbhFDltn0+KfKu5ynWeWVrAj98bnWhJwZUjjUJAJrqOniik9k+nwE+1/MU67xSSMsPn1tdSJsEeoBywQCa6HqLEN8UPvjzSZF3PU+xzitbE/jhc6sLpw0AAEjqfoSeKrdlcD1Psc4rWxP44XOrC7dDAAAgqfvgiU5m+3xSZV3PU6zzyjIAP3xudeFuCAAAJHW/RQhVbtvnkyrb9DztrzBprdVod1e7e0tB55VlAH5Wl5ekwcrd87E9HElUci0SwRsAAJAUN3jyLQ1P57xdvuvRFp2no+XpzWQggIChC2wXUA+CNwAAcFeM4ImOYjlSFbPouvgNDuJ81IM7KKrQ5b5EAICDKPlfjlTrDLsufoODOB/1YOYNxWOEFwDyQsl/P76ppiG/O9U6Q8rT54XzUQ+CNxSvj6kAIQ/4lJ0DAJDoKPpIORC56HenWGfYdfEbHMT5qAe3URSvb6kAPmWdY/wsADRFyX93KVNNu0hjjVGeniUR8bBdQD24i6J4fRvhDZlp7OMsJYD2UfLfXcpU09RprMelZPo+W1gSER+VXOtA60fx+jbCGzLT2LdZSgDdWV1e0vraQA+eHGh9bUCHe4GUG6Sn/N2pMjooegPMxp0URTgudaJvqQAhD+GUD3AAgL+UA5Epf7dvkLUoJZLBRmA2umzIXpNRvT6N8IY8hPs2SwkApUg5EJnyd/sEWU2e6ww2ArPRY0P2WKd1UMhaEtahAEC+Uq5JSvW7fdadN3muUx0RmI3gDdnzWWhdezn8kIcwC5YBALH4BFlNnusMNgKzEbwhe66jelSoAoB21D5QhsV8gqymz3UGG4GjuMMie67rtKhQBQDpsW8kplzXnbP+GvDHVYLsuY7qpd7PBv3gOqPADAT6hvXI8EVKJOCP4A2dcensuqRO9G3T7hp1HQi5pt6SqovS+VxzDJQhBCmRgB96FehEynQb0jHKlkMqlmvqLam6KJnvNUcp9/5YtCcbgPZwi0UnUnZ2Q/ez4SHVrRwCIdd9i9hMFiXzveaOGyjjPlqPHAbUANxD8IZOpO7s+m7a3eQhdeX6VW1cflhLTy9p4/LDunL9apyDhqQ8AiHXGQVmIFAy32tu3kCZJDr7FclhQA3APeSRoRO5rktbtAD/yvWrOv/8k9oevSZJurl1Q+eff1KSdO7M410ddlVyaBuu+xaxmSxKFnLNzVq3tLU9pJBJRVjbCOSFcWF0Itd1aYtGoJ964cLdwG1qe/SannrhQtLj6lMKUg5twzX1NjRVF+hS7Gsuh9lzxJMys6BPzzYgFnoWiML1BpxrZ3fRQ+rW1s2ZPzfv9Rj6tt5gVts4sbqs7eGotQe8T+U931RdoGux78ekEdcl1YBa355tQCykTSKYb5n0HMsEL0p/O7V+Wje3bhz5uVPrp5MdUx/3UtrfNtouw0/Zf9TAdQAi5v2YNOK6uO7J1rTt9fHZBsRATwTBalrMvGgE+uLZS1pbuf/Az6yt3K+LZy8lO6a+pyC13b5qas/op65nNHLNrMA9PtkyTTILXNpe359tgC9m3hCstsXMx41AT4uSPPXCBd3auqlT66d18eylpMVKcijg0aW221dt7Rn9k3JGo+msSo6ZFRhLmV3g0vb6/mwDfBG89YjPOp4mP9O3G/C5M4+3Wlmy5hSkHNtX39oz6pNqAIKU4jqkDO5d2l7NzzYgJe62FWiS/uCTRtP0Z3KoDlizWjcdz7V90Z5RulQFQ0gprkPKdEWXtrfo2ZbrswvoGsFb4Zp2gH0euk1/hvUN6aXcdLwrubYv2jNKl2oAgjVKdUhZDdS17c17tuX87AK6xlBy4ZqmP/ik0bj8DOsb8pRzNa+c21fu7dk1BdonZRrlcq0O2BQpxXVIma4Yq+3l/OwCukbwVrimHWCfhy4P6vLFWPuSquNP+/Ljuu6IdUr9lGIAgjVKdUgV3O///W2unQP6hid34ZqmP/ik0bD2p3yh6TEpU1dybF8lrLFwTYFmnRJiIaW4Hr6p+E3EuI+y0TswH5dB4Zp2gH0eujyoyxcaIKXs+MdsXzE6C6WssXBdd8Q6JcSUstOP8sW6j+Y4uAfkgqugcC7pDz6pDLmv/cHxQtNjUqeuxGhfsdICY6yxaGNtmWu6KempANoSa61a6tROoGQEbxGFdtx8f54AC8cJaR8ldPxjdRZCA9W21pa5rjtinRKAtsQc8KNvA8yWUResbKGpAqWkbCFvsddslZC6EistMHSNRVtry1zTTUl/zlcJaywBF673Ua4BwB1P70hCO24UFUCoFAMAJXT8Yy1sDw1U21xb5rruiHVK+fG9XunsImcu91EGrQE/PMEjCe24UVQAoVINAOTe8Y81OxgaqOZSHY3OfRl8rlc6u8idy32UQWvATz65T4ULXRtUwtoi5K2v++LEXNgessYi1doyl7Ww7OlWDp/rlY2LUYKm99G+PrOAUEFPc2PMe40xXzHG7BljNg997aPGmJeMMV8zxvxo2GHmL3T0v4S1RYcxwp+XXGZ+upDD7GCKFFPXmRZGssvhc72SoYGa9PmZBYQIjQy+LOlvSnp2/4vGmHdIep+k75P0PZI+a4z5L621u4Hvl63Q0f/SyuIywp+frqoKpi6P30b5/VhiV0dznWlhJLscPtcrGRroQqp7MJVwAT9BV5+19qvW2q/N+NJjkj5hrX3dWvsnkl6S9M6Q9ypB6Oh/DrMHTTHCn58uioukXoPT9zU+rjMtuY9kM1t/j8/1WmKGBsqW8h5cQkEsIEep7vhvlfSFff++NXntCGPME5KekKSNjY1Eh9OeHGYJ2jgGRvjz1Pa+OKnX4PR9jY/rTEvOI9l9mq1veg92vV5Ly9BA+XzvwamuAQANZt6MMZ81xnx5xv8eO+7HZrw2s1dvrX3OWrtprd186KGHmh53lmKPUPmMUrc1U5H7CD/akXoNTt/X+LjOtKQcyQ6dNevLbH3qe3BJGRoon889OMa+t8zQA/MtnHmz1r7b4/feknR6379PSfqmx+8pSsxZAt9R6rZmKnIe4Ud7Uq/B6fsaH5+ZlhQj2TFmzfoyW9/32WLUxeceHHIN9GmGHvCV6kr4tKT3GWPuM8a8XdIjkr6Y6L2y4TJCdeX6VW1cflhLTy9p4/LDunL96oGv+45StzVTQa46pPRrcFjjk8dMS4xZs9Jn65vOBvR9thh18bkHh1wDfZmhB0KEbhXwN4wxtyT9V5J+3RjzW5Jkrf2KpE9K+veS/i9JH6q50uRU087JletXdf75J3Vz64asrG5u3dD55588EMD53vza7CDl0KlEt1IH8QwS5CFGQFJyIO6SBlZ6kArs53MPDrkGGPwAFjPW5pOysrm5aa9du9b1YXjbP91/OJVw/41u4/LDurl148jPn17f0I0Pf0OStLU9lJU5kqpgZLW+YJFwk2MAgKZ870eH5VDQyYfL3889GH0Xcg3EutcApTPGvGit3Zz1tfyHPDvk2tFouj7l1tbNmT+//3XfNWVUI0NflBoIpJLy84i1xrXUynIu6/W4B6PvQq4B1tMDixG8zeG7aLZJ5+TU+umZM2+n1u/VeAm5+ZXaQQKaYlH7Qak/jxQBSUnBt2vRBu7B6Dvfa4DBD2AxroY5Ui6avXj2ktZW7j/w2trK/bp49tKB11hTBszGovaDUn0e+4t0bA9HWhusRLkflbb5esnr9YDS0PfpBls0lIMrYo6Ui2bPnXlczzz6rE6vb8jI6PT6hp559FmdO/N40O/lwkPJXNovi9oPSvF5+ARYTc9hacE3hXOAuOiv5KW0AbW+Y9hwjtT7S50783hwsLYfaWQomWv77fv+b4el+Dxc92pyOYc57vm2KI2TVEggDvor+fHZm6+k1Pfa8CnPESNNpunIUowRqNJGsoH9XNsvaWwHpfg8XGfzXM5hG+X0Xe6rjDoD7aG/kh/X+z33zG4RvM0RmibTtGHHugBII0PJXNsvaWwHpfg8XAMsl3OYOvh2va/SmQTaQ38lP673e+6Z3ernMHVDIWkyTaegfaaqZ2krjYxpcqTg0377ksbW9JqL/Xm4lux2OYepK8q53ldjpHFybwSaIe29PU3vS673+xxT3/uESyWRpiNLsUag2kgjY5ocqZAGOVuX15zrbJ7rOUxZUc71vhqaxsm9EWiO+307XO5Lrvf7NlLfMR8fcyJNG3asCyAkbarWCnEoR41pkDWsZXUJsHI6h6731dDOZNfnCShJTveKmrnel1aXl7Q2WNHK0ngAbHs4mvvcIgDvFp9yIk2noF2nqo/jkzZVeoU41KOmNMhY1dRKu+ZyOYeu99XQNM7SzhPQtVzuFTVzvS+5PLfYTL1bBG+JNG3YXV8ALmtDyFMHmiltLWttfO6rIZ1JzhPQDyWtbXW9L7k+twjAu5Nni6tE05Sj/d/3r//wk/rPf+HtWnp6SRuXH9aV61eTHmNOFeKAWviUXZ6VYsk15y/lmrrDOE9A/Upb2+p6X2r63GKD9e4RvGXkyvWrOv/8k7q5dUNWVje3buj8808mDeBc1oaQpw4043JdHdch4Jprn0/HhPME1K+0ta0pipCUFsDWylibT07+5uamvXbtWteH4SXGVPrG5Yd1c+vGkddPr2/oxoe/EetQD9if43x4bQgdD8CPy3W1tT2UlTmS2mJktU5KSqu4HwKY55VXhxrMGIEbjvb04Mny79VN7n88r9pjjHnRWrs562vkdEQQqzjBra2bTq/HCBi7XnMH7FfSeoLjuFxXFLvIR6y1ioCLWu57tat9P90mzy2eV3kgeIsg1gP/1PrpmTNvp9ZPH3ktVsAo1b3olIdi95qeg5htOgdNryuKXeSDjgnaVtt9r2Yxq4PP03Z7mPV8Pm4GjedVHvi4I4i10fbFs5e0tnL/gdfWVu7XxbOXjnxvabnXXSA3u3su56CvbZpiF/lg41m0ra/3vRxcuX5VG5cfblwgbtEastL29vTpI/G8ygOPpAhiPfDPnXlczzz6rE6vb8jI6PT6hp559FmdO/P4ke+NFTDWjIdi91zOQV/bNMUu8hGjY0IlNrjo632va74F4uZVsY01WNxme/DpI/G8ygMFSyLoYpE7i0YXq31xcQlczgFtGjkISbWm4Alccd/rRuwCcbHOY5vtgT5S3o4rWMLTJII2ptIPY+p6MVKguudyDmjTyEHI/nDM9sMV971uuBaIWyTWjFmK9jCvD0ofqVycokhST6XPej+mro/HQ7F7LueANo0chAy2kQIHV9z3ujGrENxxry8SKxCK3R6O64N20UcirTwO7g6JpRyJDRkh7gMeit1zPQe0aXQpdLCNkWz44L7XPpcCcU00CYSaBi4x28NxfdC2+0gUkYuHKYjEKD3drZq3QSgF5wClCN32pY1S4kCpctq6Z1oI7qkXLujW1k2dWj+ti2cvzSwQ18SiPdK62hJiUR+0zecz+2jGQ/CWGHtioE9yejgDrkIH21w2Zwf6JMf97M6deXxusObzLDsuEOoqcMmpD8pkRjw8URJj3RX6gpQIlC5G2uPhlCdJrPFA75VUzCfFs6yr9bA5bX9CWnk8fGSJse4KfVHSwxmYJfZgGwMawFhJxXxSPMu6ClxC+6Ax72FMZsTDJ9YC1vygD0iJQOlipz2yxgMYyyl9b5EUz7Iu18OG9EFj3sNIK4+H4A1AFCU9nIF5Yg62MaABjJVUzCfFs6zUwCX2PYzJjDgI3tAbFNNIq6SHM9AGBjRQk5BnaFvBS4znfKpnWYmBC/ewPPHxoxdYe5Ie6zuBg1jj0T42AU4jxjM09X52sZ7zPMvu4R6WJz599EKKtSfM5B3lM7LI54gcpGiHpaZKlSrHcvSHj6/Ue10J6zdjr8/K5e/qkss9rOT2XRqCN/RC7Lzt3DsJpeBzRA5StkM6ge1pK8Dw6aSWfq8rYf1mCceYC5c23OQeVnr7Lg2fKHohdpleyuLHweeIHNAO69BGOXrf1LzS21gJe3T5HGMf02xTLCMpvX2XJqPLDkgndt52SXvW5IzPETkIaYd97Pzlqo0Aw7eTWvq9roS1T67H2MZa+BzvDykCrdLbd2kI3tALsRcglzAKWQI+R+TAtx1SCCkvbQQYvp3UkHtdygCg6e8uoYiH6zGmni3K9f6QItDiWd6ufIZMgMRirj2hLH4cffkcWcidN992WEIRhz5po0CMb+l03zaWci2R6+8uYf2myzGmXiOX6/0hRfn/vjzLc0HvAfBQwihkCfrwOeY6+op7fNshqUL5SV2O3nd2z7eNpZwd6vs6pdSzRbneH1LMUPfhWZ4TZt4ATyWMQpag9s8x19FXHOTTDtnANh++FSBdf6bt7R9Szg71vTpj6tmiXO8PTduw6/VR+7M8JzxiACChXEdfEa6EIg594DO7HTIj7jO75/t+KWeH+r5OKfVsUcz7Q+x1j4vaMBkjeeMJAxzC+iTElOvoK8KVvgl3rHtd1/dMn9nt2DPiiz4D3/dLOTvEOqW0s0Wx7g9d7KFGxkjeynjCAC1htAmxMTtTt9RrrFKJda/L4Z7pM7sdc0a8yWfg+34pZ4dYpxTXrNmxGPeHLtYmkjGSN3oPKF7MUV+X0aauR5tRhtJnZ1CnWCPrOYzQ+8xux5wRb/IZhLxf6tkhZlLCpZwd62JtIhkjeeM0oGixR32bjjblMNqMcsSanclxw9f9cj8+3BNrZN3n98RuJz6z2zFnxJt8BszA1y3l7FgXaxNpr3kjeEPRYt8wm94k+15iGe3LfcAg9+PDQbE6hK6/J0U78Un/i5ky2OQzcHk/BkHKkzLNsItAipTavBFCo2g+6QTHpTs2XcDd9xLLaF8O6WnHyf34cFCsYhWuvydVO/FJ/4uVMtj0M2jyfjHT70jtb0/KNMOuUu9Jqc0XVzGKFnvUt+loU99LLKN9LiO7XYzcs8C9LLFG1l1/T43tJOYsRaysjpgFaZgFXCz17FiphZGQBjNvKFqKUd8mo02UWEbbmo7s+o7ch47Ss8C9PMfd61zag8sIfYntpMln0XRWbdHviZXVEWOGs4sS9aXqS2EqZnPzwCeOonU16ks+ONrWdGTXZ+Q+xig9C9zrkXL9YmntpO0tFWJldcR41rG2203o7Fjus5ysa84HPU0Uz+WGGTPdkTQGtKnpgIFPpy1GJ40BjXo0bQ8+nc3S2kmsAKbp74kV3MZ41tWY4pqrEgIjgvl85DnUBTTgM31PuiNK1iQ1yyctLVaqFgvc69CkPbim1JWabhXr2mj6e2Kl38V41pWY4lqqEgo+UagtH1yCyE6T0VzfUarSRn0BVz4j9xTgqUOstKsm7cFlFL6EWYV5uthSIUZWR4xnXWkpriUrYZaT50Q++MiRlaYP+ZDpe9IdUTOfThudtPLFDJCatAeXzmbJ6Vaxro2u9uoKedYx2NmeEgIjnhP5yKhZAM0f8iWMUgFdce200UkrX8wAqUl7cOlslny/7mpLhVww2NmOEgKjUttwjfJpFYCa51S75OKHrrUoda0G4qq9HbBerWyx16Msag8ua6pKXzsV69rgGsM8pWw1QBvOA8EbstL0Id+04xC6Tw373ECqtx3UHpD2ScwAqem+Zk07mxSKAo6adZ2tExihAZ7SyErT1IGm0/ehqUQlr9VAPDW2g5KLSOCoWGlXLu2iaUqdS7pVqr2uct9DC/3C/RchCN6QFZeHfJOOQ+hai5LXaiCeGttBjQFpn8Vaj5KqXTS5X6fq0NJRRm64/yJEUNqkMebnJf11SUNJfyTpp6y135l87aOSPihpV9Lfsdb+VuCxomAu6Vkxc6pDU4ma/jzpZ3UrYc2Oaxtkz576xLh3+raLGPfAVHtdlbCHFvqF+y9ChHY9PiPp+621PyDpDyV9VJKMMe+Q9D5J3yfpPZL+mTFmOfC9UKguRz1DU4ma/DyjuvXLvRKYTxtMWZq6xBQ1n2Mu8e9cxKddxLoHpprhrnHmHGUrYWsA5CuomVhrf9taO53j/YKkU5P/fkzSJ6y1r1tr/0TSS5LeGfJeKFeX6QGhqURNfp70h/rlXiLZpw2mCkjbHsyIEUD5HHOtgzY+fGd8QAAAFf9JREFU7SLWPTBVh5aOMnKT+4Ag8hbz1vW3Jf3m5L/fKunmvq/dmryGHup61DPGRqXH/XzXfx/akfN+Rz5tMFVA2uZgRqwAyueYax208WkXse6BqTq0dJSRm9wHBJG3hXcuY8xnJf2nM750wVr7qcn3XJA0kvTx6Y/N+P6ZibzGmCckPSFJGxsbDQ4ZpSlhvVCI2v8+5M+3DabYs6fNtRyx1jL5HHPNa1Zc20Wse2Cqva5K2UML/cKeafC1MHiz1r77uK8bYz4g6VFJP2KtnT61bkk6ve/bTkn65pzf/5yk5yRpc3Oz/Kcejoi5x08XhUEWvSd7GPVXLoVqfNpgqmNvczAjVgDlc8wM2twT8x6YqkNLRxlALYIeM8aY90j6GUk/bq19bd+XPi3pfcaY+4wxb5f0iKQvhrwXyhUrPaCLNSZN3pP0h3RyLgiR05on1zaY8tjbTFGLtZbJ55ib/kzObTgW7oEA0J7Qp+k/lXSfpM8YYyTpC9ba89barxhjPinp32ucTvkha+1u4HuhYDFGPbso99z0PWf9fbnMypRqGmAsLy1psDIezb99ZyebTmFu5cddrrGUx940RS3G9RFrxscnra7Jz+TehmNiZgsA2hEUvFlr/4tjvnZJ0qWQ3w/s18Uak5A9j/rSaUslt+DosJLXPKU+9kUd+VjXR8y1TD7Bx6Kfyb0NAwDKQ6klFKOLNSa+70mnLVxXwVHTGaGS1zx1fewxr4+cZ3xKDvABAHkqoJsBjHVR7tn3Pdk+IJzvZsEh64tc1oKVXH6862Mv/fpo2s7YXwwAEBuPEBTDd1F8aIfeWqvvvDbUK39xR6Pd3UbvSactnGuAEaMIh8veXSUXaej62Eu+PvoS4AMA8sQTBEVxTZEKWVsz/dmV5WX9pTes3O14NcH2AeFc1zPFSMVzTXPLOWVvkS6PveTrw6Wdsb8YACA2gjdULaRDH/KzdNricAkwYqwv6notWF+UfH30KcAHAOSH4A1VC+nQhwYDdNraFSPwKnlGqDSlXh8E+ACALvG4QdVC1taUvC6nj2KsL+p6LRjyxzo2AECXeNqgaiEzKczC5GteOf8YqXilzgihHSWnfAIAykfwhqqFdLTopOVpUREan8Cr6d5ucFfjZ0uADwDoCsEbvJXSKQvpaNFJy0/sDdBDKpI2+d0lXCOppPxsAQDoI56e8BJjTy3AR+wNnl32dnPBNZLuswUAoK8I3uCFThm6EruQTOxgcIprJN1nm9LO7p62tod65dWhtraHvQq2AQD5I3iDlxI6ZSGdMDpw+Ypd7S9VVdESrpHUSqvYymwpACB3mT5CkbvcO2UhnTA6cHmLXc4/NBicF+jnfo20obSy+m3MljIwBAAI0aNuBGLKvVMW0gkj3S1/q8tLWl8b6MGTA62vDYKKX4QEg8cF+rlfI23IZd+8pgFT6tlSBoYAAKEI3uAll07ZPCGdMNLd0NRxgX7u10hbYgbaPlwCptSzpQwMAQBC9WcIGNF1UUa/aen1aSdsfxDWtBMW8rMoT0g5+9GeNFg5GugPR+MAgK0muueytcTaYEW37+xIGn//dLb0gROrUY5lUXsBAGARuqMohssIekjKGuluZfJdSxQyG9LXdW0lrdtymUlPPVva1/YCAIiHRwaK4dLJDumEke5WnpC1RCFpsn0M9Etbt+UaMKVM8+xjewEAxMUTA8VwTTkKSVkj3a0sLqlxh4WkyU4D/e3hSMPR+GdqD/RDPusupE6FdNHH9gIAiIvgDcVgLRrmCVlLFNq571ugX9q6rdwCpr61FwBAXARvKEZbI+hNi6IgH8yetafEQRQCJgBALQjeUIw2OtkhlQdn/S6CwHYwe9Zee8spDREAgL4heEMUbXUcU3eyY63niRkEYrG+z57FaG9Nr+G+f9YAAHSJ4A3BagpUYq3nKa2oQw26nD3repY1tL25XsM1zFQCAFCisnrWyFLIPlm5ibUPU0j5eZQlh9L5oe2tpmsYAICaMfOGYKVVnztOrPU8JRZ1yFXKWa0YvzuHWdbQ9lbTNQwAQM3oSiJYrNmqHMTaoJvNeONIOasV63fnMMsa2t5quoYBAKgZj2Y0srO7p63toV55dait7eGBDm5tgcrq8pLW1wZ68ORA62sDr1meWEFg36VM54v1u3MIfELbW23XMAAAteLJjIUWFTOg+txsFHUIlzKdL9bvzqV0fkh74xoGAKAMBG9YqMmanr4FKl1XF+yLlGsHY/3uWgIfl2uY9g8AQDcI3rAQxQwOCt0agY5vcylntWL+7jYGL3JpNzVtDQIAQGl40mKhHNb05CRkrVQOZeVLknLtYEnrEnNqN2wrAABAd5h5w0K5rOnJRchMZA5l5UuTclarlHTfnNoNM/EAAHQnvyFmZCdkhuK4KpWlCpmJzKGsPMqTU7thJh4AgO7wuEUjPuXzc0r1iimkrDodX/jIqd2wrQAAAN2hy4hkal0bEzITWWPHt8bZVR8pP4ec2k1JawUBAKhNuT1GROdazW7R99e8NsZ3rVQtZeWnqDw4lvpzyK3dlLJWEACA2hC8QZJ757PJ96fco6tkNXV8cyqk0aU2PgffdpPLFgMAACAcwVsmuu5guXY+m3w/VSrb0WXbqXl21UWun4PLoEzX9yAAALAYT+YMhBT2iLXOxrWaXZPvZ21Mel0XhcmpkEaX2vocXK/3putOu25HAACgmZ51sfLkW9gjZofLtfPZ9Pt9qlSiua6LwuRUSKNLbXwOPtd700GZrtsRAABohp50Bnz3cIrZ4XLtfNJpz0PX+38xuzrWxufgc703HWTpuh0BAIBm6GlnwLewR8x1Nq7V7HKrftdXORSFqakAS4hFn0PomjKf673putMc2hEAAFiM4C0DvoU9Yne4Dnc+p+tr5nU293//tGN6m2IHrQopCkOBivbE2ErA53pvOshCcSEAAMpATy0DvilXKVMXXdbXUOygO75th3MWzqV4SIwUZ9/rvcm6U9JfAQAoA0/mTCzqYF25flUblx/W0tNL2rj8sK5cv5q0w+XS2aTYQbd8isJwzsK4Br8x1pSlDrBiFReKVQEXAAAcRdpkAa5cv6rzzz+p7dFrkqSbWzd0/vknJUnnzjyeZL2Ry/qaXPe46huXNEjOWRjXfRFjpTjnvr4wRnooAACYj6dpAZ564cLdwG1qe/SannrhQrL3dNk6gL2+uuc6E8Q5C+M6k9aX6qy5zegyCwgAqA1dtQLc2rrp9HoMLp3NvnRMc+baaeachXENfvuypiynLQdY1wkAqFFdPYdKnVo/7fS6FD7i7NLZ7EvHNGeunWbOWRif4LcPG9a3NaPb5P6W2ywgAAAx1Nd7qNDFs5e0tnL/gdfWVu7XxbOXZn5/rBFnl85mHzqmXVrUWfXpNHPO/BH8ztbGjG7T+1tOs4AAAMTS755GIc6deVzPPPqsTq9vyMjo9PqGnnn0WZ078/jM7+9yxJk1JvE16aySBtm+Pga/i67vNoLapvc31nUCAGrEYyxTh7cGkKQbH/6G9n52Tzc+/I25gZvU3Yiz74zfrG0QcE+TziozQUit6fWdOqhten9jQAMAUCOeYhlatDXAIj5lyV3KzM/jWj5dCv9bYxx37pqW9Q8tI9+HzxL+fK7vFJre36YDGtvDkYaj8dcZ0AAAlI6nWIZCtwZwHXGOtUbOZ8Yv5G+NWU0u53TPNtK/qMzXXM5tJaVc1pC53N/6mNoKAKgbT7IMhW4N4JpCF2uNnE+QEfK3xjru3AOXNtK/qMzXTO5tJaWu1pAdDpYlkSIMAOgtnnYZ8tka4DCXEedYI+o+QUbI3xrruHMPXNpYz5bLrErucm8rKXWxhmxesCyJGTUAQC/xxMuQ69YAoWKNqPsEGSF/a6zjLiFwcQnGfdL6qMzXTAltJZVF13eKdNI+B8sAAMwS1DUzxvxDY8zvG2O+ZIz5bWPM90xeN8aYXzDGvDT5+l+Jc7j94Lo1QKiYI+qzgozjOnUhf2us464pcPFN66MyXzM1tRUf8wYRUqWT9jlYBgBgFmOtXfxd837YmHVr7dbkv/+OpHdYa88bY35M0v8o6cck/bCk/81a+8OLft/m5qa9du2a9/HAX6pKg9NO3XT0fBoULJqRa3o8MY7b9xhztLU9lJU5UonPyGp9QUVAqk0uVlNbiSmk3XXxewEAyJkx5kVr7easrwUNq08Dt4mTkqaR4GOSrthxZPgFY8wbjTFvsdZ+K+T9kE5omfl5fMqL7+8gD1bGHeTbd3ZmdpBjHHdNJcWbbiswS6o2UJOa2kpMi9qd78DA2mBFr7z6unZ27fjpYqTVZaMHT96X4s8AACB7wTlRxphLks5J+nNJ/+3k5bdK2l8u8NbktSPBmzHmCUlPSNLGxkbo4SAzPsFEF/tJ1RK4+OzxBze1tJUmmgZdx7U7l8GYuayVlZGxkwgOAICeWvjkNMZ81hjz5Rn/e0ySrLUXrLWnJX1c0k9Pf2zGr5rZW7fWPmet3bTWbj700EO+fwcy5bNGiHUu/li7hlhc1rEd1+5Cio5sD0c6sbqiN568T286OdAbT96nE6srFCwBAPTWwh6dtfbdDX/Xr0j6dUk/q/FM2/5a76ckfdP56FC8tcHKpLT30TVC8zB75I+0PsTiMgN+XLu7HZDKG5IGDABAjYKG440xj1hrvz75549L+oPJf39a0k8bYz6hccGSP2e9Wz/5BBM+AR/u6VNaH9IVmnENnOa1u5DBGAZyAAA4KDSX6ueMMX9Z0p6kb0g6P3n9NzSuNPmSpNck/VTg+6BgrsEEs0dAM1HWk80RK3AKGYxhIAcAgINCq03+d3Net5I+FPK70W/MHgGLpSzuEytwChmMYSAHAICDqGIAAIVKuSYsZuAUMhjDQA4AAPcQvAFAoVKvCSNwAgAgL+SeAEChYm8NceX6VW1cflhLTy9p4/LDunL9aszDBQAAgZh5A4BCxUxtvHL9qs4//6S2R69Jkm5u3dD555+UJJ0783jU4wYAAH6YeQOAgq0uL2l9baAHTw60vjbwLubx1AsX7gZuU9uj1/TUCxcW/uzO7p62tod65dWhtraHMzfyBgAA4QjeAAC6tXXT6fWp6XYFVkaDlSVZGd2+s0MABwBAAgRvAACdWj/t9PrUrO0KlpfG2xUAAIC4CN4AALp49pLWVu4/8Nrayv26ePbSsT832tOBapfS+N8jJt4AAIiO4A0AoHNnHtczjz6r0+sbMjI6vb6hZx59dmGxkul2BfvF3K4AAADcY6wN38w1ls3NTXvt2rWuDwMA0NB0zds0dXK6XYFv1UsAAPrOGPOitXZz1td4sgIAvE23KzCyGo72ZGQJ3AAASIR93gAATnZ297Q9HGm0N06bXBusaH1t0PVhAQBQPYZGAQCNsTUAAADdIXgDADTG1gAAAHSH4A0A0BhbAwAA0B2CNwBAY2wNAABAd3jcAgAaWxusaHdv724AN90aYG1A/SsAAFIjeAMANMbWAAAAdIehUgCAk9XlJa2yNQAAAK1jqBQAAAAACkDwBgAAAAAFIHgDAAAAgAIQvAEAAABAAQjeAAAAAKAABG8AAAAAUACCNwAAAAAoAMEbAAAAABSA4A0AAAAACkDwBgAAAAAFIHgDAAAAgAIQvAEAAABAAQjeAAAAAKAABG8AAAAAUACCNwAAAAAoAMEbAAAAABSA4A0AAAAACmCstV0fw13GmJclfaPr48BMb5b0H7s+CDTG+SoL56ssnK+ycL7KwvkqC+crjYettQ/N+kJWwRvyZYy5Zq3d7Po40Aznqyycr7JwvsrC+SoL56ssnK/2kTYJAAAAAAUgeAMAAACAAhC8oannuj4AOOF8lYXzVRbOV1k4X2XhfJWF89Uy1rwBAAAAQAGYeQMAAACAAhC84VjGmH9ojPl9Y8yXjDG/bYz5nsnrxhjzC8aYlyZf/ytdH2vfGWN+3hjzB5Pz8X8aY96472sfnZyrrxljfrTL48SYMea9xpivGGP2jDGbh77G+cqQMeY9k3PykjHmI10fDw4yxvyyMebbxpgv73vtQWPMZ4wxX5/8/5u6PEbcY4w5bYz5vDHmq5N74d+dvM45y5Ax5oQx5ovGmOuT8/X05PW3G2N+d3K+/oUxZtD1sdaO4A2L/Ly19gestT8o6XlJ//Pk9b8m6ZHJ/56Q9H90dHy45zOSvt9a+wOS/lDSRyXJGPMOSe+T9H2S3iPpnxljljs7Skx9WdLflPQ7+1/kfOVpcg5+UeN73zskvX9yrpCPf67xNbPfRyR9zlr7iKTPTf6NPIwk/X1r7fdKepekD02uKc5Znl6XdNZae0bSD0p6jzHmXZL+kaTLk/P1Z5I+2OEx9gLBG45lrd3a98+TkqaLJB+TdMWOfUHSG40xb2n9AHGXtfa3rbWjyT+/IOnU5L8fk/QJa+3r1to/kfSSpHd2cYy4x1r7VWvt12Z8ifOVp3dKesla+8fW2qGkT2h8rpAJa+3vSHrl0MuPSfrY5L8/JuknWj0ozGWt/Za19t9O/vu2pK9Keqs4Z1ma9Pf+YvLP1cn/rKSzkv7l5HXOVwsI3rCQMeaSMeampP9e92be3irp5r5vuzV5DXn425J+c/LfnKuycL7yxHkp03dba78ljYMFSd/V8fFgBmPM2yT9kKTfFecsW8aYZWPMlyR9W+Nsnz+S9J19A8fcF1tA8AYZYz5rjPnyjP89JknW2gvW2tOSPi7pp6c/NuNXUbo0sUXnavI9FzROR/n49KUZv4pz1YIm52vWj814jfPVPc4LkIAx5g2Sfk3S3zuU7YPMWGt3J8toTmmcjfC9s76t3aPqn5WuDwDds9a+u+G3/oqkX5f0sxqPrpze97VTkr4Z+dBwyKJzZYz5gKRHJf2IvbcPCOeqIw7X1n6crzxxXsr0p8aYt1hrvzVJ7f921weEe4wxqxoHbh+31v6rycucs8xZa79jjPk3Gq9VfKMxZmUy+8Z9sQXMvOFYxphH9v3zxyX9weS/Py3p3KTq5Lsk/fk0zQHdMMa8R9LPSPpxa+1r+770aUnvM8bcZ4x5u8ZFZr7YxTGiEc5Xnn5P0iOTymoDjYvKfLrjY8Jin5b0gcl/f0DSpzo8FuxjjDGSfknSV621/2TflzhnGTLGPDStYm2MWZP0bo3XKX5e0t+afBvnqwVs0o1jGWN+TdJflrQn6RuSzltr/8PkpvtPNa7s9Zqkn7LWXuvuSGGMeUnSfZL+v8lLX7DWnp987YLG6+BGGqem/Obs34K2GGP+hqT/XdJDkr4j6UvW2h+dfI3zlSFjzI9J+l8lLUv6ZWvtpY4PCfsYY35V0l+V9GZJf6pxlsi/lvRJSRuSbkh6r7X2cFETdMAY899I+r8l/TuN+xiS9A80XvfGOcuMMeYHNC5Isqzx5M8nrbX/izHmP9O4gNODkv5fSY9ba1/v7kjrR/AGAAAAAAUgbRIAAAAACkDwBgAAAAAFIHgDAAAAgAIQvAEAAABAAQjeAAAAAKAABG8AAAAAUACCNwAAAAAoAMEbAAAAABTg/wd3wuEQgCigGwAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 1080x648 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "# This and several other helper functions are defined in a code cell above. Hit the \"code\"\n",
    "# button above if you're curious about how they're implemented.\n",
    "plot_by_title_pattern('Bloomington', figsize=(15, 9), bg_alpha=.05, text=False);"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.7.4"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}