topic-keywords.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Topic keywords\n",
    "\n",
    "Check which words are typical for topic tweets in different time frames. Uses the Ttest."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {},
   "outputs": [],
   "source": [
    "import csv\n",
    "import datetime\n",
    "import math\n",
    "import matplotlib.pyplot as plt\n",
    "import os\n",
    "import pandas as pd\n",
    "import random\n",
    "import re\n",
    "import sys\n",
    "import time\n",
    "from IPython.display import clear_output\n",
    "from nltk.tokenize import TweetTokenizer\n",
    "from collections import Counter\n",
    "sys.path.append(\"/home/erikt/projects/newsgac/fasttext-runs\")\n",
    "import tscore"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "DATADIR = \"/home/erikt/projects/puregome/data/text/\"\n",
    "ID = \"id_str\"\n",
    "REPLYID = \"in_reply_to_status_id_str\"\n",
    "TEXT = \"text\"\n",
    "TOKEN = \"token\"\n",
    "USER = \"user\""
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "def squeal(text=None):\n",
    "    clear_output(wait=True)\n",
    "    if not text is None: print(text)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Count tweets with topic words"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [],
   "source": [
    "def countTweets(datePattern, query):\n",
    "    count_query = 0\n",
    "    count_all = 0\n",
    "    fileList = sorted(os.listdir(DATADIR))\n",
    "    for inFileName in fileList:\n",
    "        if re.search(datePattern, inFileName):\n",
    "            squeal(inFileName)\n",
    "            df = pd.read_csv(DATADIR+inFileName,compression=\"gzip\",index_col=ID)\n",
    "            count_query += len(df[df[TEXT].str.contains(query, flags=re.IGNORECASE)])\n",
    "            count_all += len(df)\n",
    "    return(count_query, count_all)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [],
   "source": [
    "#FILEPATTERN = \"20200522\"\n",
    "\n",
    "#for query in \"corona|covid\".split():\n",
    "#    count_query, count_all = countTweets(FILEPATTERN,query)\n",
    "#    print(count_query, count_all, query)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Find other relevant words in topic tweets"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [],
   "source": [
    "import warnings\n",
    "warnings.simplefilter(\"ignore\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [],
   "source": [
    "AT = r\"@\"\n",
    "HASH = r\"#\"\n",
    "\n",
    "def get_tokens_of_matched_tweets(file_pattern, query, query2, spy=True):\n",
    "    file_list = sorted(os.listdir(DATADIR))\n",
    "    match_tokens = {}\n",
    "    non_match_tokens = {}\n",
    "    for file_name in file_list:\n",
    "        if re.search('^' + file_pattern, file_name):\n",
    "            if spy: squeal(file_name)\n",
    "            df = pd.read_csv(DATADIR+file_name, index_col=ID)\n",
    "            if query2 == \"\":\n",
    "                matching_rows = df[TEXT].str.contains(query, flags=re.IGNORECASE)\n",
    "            else:\n",
    "                matching_rows = df[TEXT].str.contains(query, flags=re.IGNORECASE) & df[TEXT].str.contains(query2, flags=re.IGNORECASE)\n",
    "            df_match = df[matching_rows]\n",
    "            df_non_match = df[matching_rows == False]\n",
    "            text = \"\"\n",
    "            for i in range(0, len(df_match)):\n",
    "                text += \" \" + re.sub(\"\\\\\\\\n\", \" \", str(df_match.iloc[i][TEXT]))\n",
    "            for token, count in Counter(TweetTokenizer().tokenize(text.lower())).items():\n",
    "                if not token in match_tokens: \n",
    "                    match_tokens[token] = count\n",
    "                else: \n",
    "                    match_tokens[token] += count\n",
    "            text = \"\"\n",
    "            for i in range(0, len(df_non_match)):\n",
    "                text += \" \" + re.sub(\"\\\\\\\\n\", \" \", str(df_non_match.iloc[i][TEXT]))\n",
    "            for token, count in Counter(TweetTokenizer().tokenize(text.lower())).items(): \n",
    "                if count > 1:\n",
    "                    if not token in non_match_tokens: \n",
    "                        non_match_tokens[token] = count\n",
    "                    else: \n",
    "                        non_match_tokens[token] += count\n",
    "    match_tokens_clean = { t:match_tokens[t] for t in match_tokens if not re.search(\"^http\", t, flags=re.IGNORECASE) and \n",
    "                                                                      not re.search(\"^@[a-z]\", t, flags=re.IGNORECASE) }\n",
    "    non_match_tokens_clean = { t:non_match_tokens[t] for t in non_match_tokens if not re.search(\"^http\", t, flags=re.IGNORECASE) and \n",
    "                                                                                  not re.search(\"^@[a-z]\", t, flags=re.IGNORECASE) }\n",
    "    return(match_tokens_clean, non_match_tokens_clean)\n",
    "\n",
    "\n",
    "def readData(fileName):\n",
    "    return(pd.read_csv(fileName,index_col=TOKEN).to_dict())\n",
    "\n",
    "def writeData(data,fileName):\n",
    "    pd.DataFrame(data).to_csv(fileName,index_label=TOKEN)\n",
    "    \n",
    "def findKeysStartingWithChar(data,char):\n",
    "    keysStartingWithChar = []\n",
    "    for key in data:\n",
    "        try:\n",
    "            if re.search(r\"^\"+char+r\"\\w\",key):\n",
    "                shortKey = key[1:]\n",
    "                if shortKey in data: keysStartingWithChar.append(shortKey)\n",
    "        except: pass\n",
    "    return(keysStartingWithChar)\n",
    "\n",
    "def combineKeysStartingWithChar(data,char,keysStartingWithChar):\n",
    "    for key in keysStartingWithChar:\n",
    "        data[key] += data[char+key]\n",
    "        del(data[char+key])\n",
    "    return(data)\n",
    "\n",
    "def combineInitialHashAt(data):\n",
    "    keysStartingWithHash = findKeysStartingWithChar(data,HASH)\n",
    "    data = combineKeysStartingWithChar(data,HASH,keysStartingWithHash)\n",
    "    keysStartingWithAt = findKeysStartingWithChar(data,AT)\n",
    "    data = combineKeysStartingWithChar(data,AT,keysStartingWithAt)\n",
    "    return(data)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [],
   "source": [
    "NBROFEXAMPLES = 20\n",
    "\n",
    "def dictTopN(dictionary,n=NBROFEXAMPLES):\n",
    "    return([(x[1],x[0]) for x in dictionary.items()][0:n])\n",
    "\n",
    "def dictBottomN(dictionary,n=NBROFEXAMPLES):\n",
    "    return([(x[1],x[0]) for x in dictionary.items()][-n:])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [],
   "source": [
    "NBROFTOKENS = \"totalFreq\"\n",
    "NBROFTYPES = \"nbrOfWords\"\n",
    "WORDFREQS = \"wordFreqs\"\n",
    "\n",
    "def makeTscoreData(tokenList):\n",
    "    data = { NBROFTOKENS:0, NBROFTYPES:0, WORDFREQS:{} }\n",
    "    for token in tokenList:\n",
    "        if not math.isnan(tokenList[token]):\n",
    "            data[WORDFREQS][token] = tokenList[token]\n",
    "            data[NBROFTYPES] += 1\n",
    "            data[NBROFTOKENS] += tokenList[token]\n",
    "    return(data)\n",
    "\n",
    "def sortTscores(tscores):\n",
    "    return({token:tscores[token] for token in sorted(tscores.keys(),key=lambda t:tscores[t],reverse=True)})"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [],
   "source": [
    "QUERY_DISTANCE= r\"1[.,]5[ -]*m|afstand.*hou|hou.*afstand|anderhalve[ -]*meter\"\n",
    "QUERY_TOPIC = \"corona|covid|mondkapje|rivm|blijfthuis|houvol|huisarts|flattenthecurve\"\n",
    "QUERY_TOPIC_SMALL = \"corona|covid\"\n",
    "QUERY_TRAVEL = \"reis|reizen\"\n",
    "QUERY_WASH_HANDS = \"hand.*was|was.*hand|pfizer|astra\"\n",
    "HAPPY_QUERY = r'\\b(geluk|gelukkig|gelukkige|blij|happy)\\b'\n",
    "LONELY_QUERY = r'eenza|alleen.*voel|voel.*alleen|lonely|loneli'\n",
    "IK_QUERY = r'\\b(ik|mij|mijn|me|mn|m\\'n|zelf|mezelf|mijzelf|i)\\b'\n",
    "# 20220127 added: herstelbewij|^cst$\n",
    "# 20220127 removed: apartheid\n",
    "CTB_QUERY = r'(ctb|qr|toegangsbewij|testbewij|coronapas|vaccinatiepas|vaccinpas|\\bcodes\\b|2g|3g|1g|herstelbewij|^cst$)'\n",
    "# 20211202 added: booster|vax|mrna|inject\n",
    "QUERY_VACCIN = \"vaccin|ingeënt|ingeent|inent|prik|spuit|bijwerking|-->|💉|pfizer|moderna|astrazeneca|astra|zeneca|novavax|biontech|booster|vax|mrna|inject\"\n",
    "TOPICQUERY = \"corona|covid|huisarts|mondkapje|rivm|blijfthuis|flattenthecurve|houvol\"\n",
    "TESTQUERY = r'\\btest|getest|sneltest|pcr'\n",
    "LOCKDOWNQUERY = \"lock.down|lockdown\"\n",
    "DISTANCEQUERY = \"1[.,]5[ -]*m|afstand.*hou|hou.*afstand|anderhalve[ -]*meter\"\n",
    "# 20220127 added: delta|omicron|omikron|antistof|serologi|immuun|immuni|^ema$|^lareb$|^ifr$|aerosol\n",
    "PANDEMICQUERY = \"|\".join([TOPICQUERY, TESTQUERY, QUERY_VACCIN, LOCKDOWNQUERY, DISTANCEQUERY, CTB_QUERY,\n",
    "                                      r'virus|besmet|ziekenhui|\\bic\\b|intensive.care|^zorg|vaccin|[^ad]arts|uitbraak|uitbrak|pandemie|ggd|'+\n",
    "                                      r'mondkapje|quarantaine|\\bwho\\b|avondklok|variant|verple|sympto|e.golf|mutant|^omt$|umc|hcq|'+\n",
    "                                      r'hydroxychloroquine|virolo|zkh|oversterfte|patiënt|patient|intensivist|🦠|ivermectin|'+\n",
    "                                      r'delta|omicron|omikron|antistof|serologi|immuun|immuni|^ema$|^lareb$|^ifr$|aerosol'])\n",
    "\n",
    "def show_top_tokens(file_pattern, query, query2=\"\", n=200, spy=False):\n",
    "    queryTokens = {}\n",
    "    nonQueryTokens = {}\n",
    "    tscoresDataQuery = {}\n",
    "    tscoresDataNonQuery = {}\n",
    "    start_time = time.time()\n",
    "    query_tokens, non_query_tokens = get_tokens_of_matched_tweets(file_pattern, query, query2=query2, spy=spy)\n",
    "    tscores_data_query = makeTscoreData(query_tokens)\n",
    "    tscores_data_non_query = makeTscoreData(non_query_tokens)\n",
    "    print(file_pattern, end=\" \")\n",
    "    top_tokens = dictTopN(sortTscores(tscore.computeTscore(tscores_data_query, tscores_data_non_query)), n=n)\n",
    "    for x in top_tokens: \n",
    "        if re.search(query, x[1]):\n",
    "            print(x[1], end=\" \")\n",
    "        else:\n",
    "            print(f'\\x1b[31m{x[1]}\\x1b[30m', end=\" \")\n",
    "    print(f\"\\ncomputation took: {int(time.time() - start_time)} seconds\")\n",
    "    return top_tokens"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 22,
   "metadata": {},
   "outputs": [],
   "source": [
    "def save_top_tokens(top_tokens, month):\n",
    "    data_file = open(f\"csv/ttest_{month}.txt\", \"w\")\n",
    "    for data in top_tokens:\n",
    "        print(int(data[0]), data[1], file=data_file)\n",
    "    data_file.close()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 24,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "20210131-23.out.gz\n",
      "202101 corona covid \u001b[31m19\u001b[30m #corona #coronamaatregelen #coronavirus #covid19 #coronadebat coronavirus \u001b[31m-\u001b[30m coronamaatregelen coronacrisis \u001b[31mvaccin\u001b[30m \u001b[31mmaatregelen\u001b[30m \u001b[31m#avondklok\u001b[30m \u001b[31m#lockdown\u001b[30m \u001b[31mvirus\u001b[30m #covid19nl \u001b[31mtegen\u001b[30m \u001b[31maantal\u001b[30m #coronavaccin \u001b[31m%\u001b[30m \u001b[31mmensen\u001b[30m \u001b[31mbritse\u001b[30m coronavaccin \u001b[31mpatiënten\u001b[30m \u001b[31mde\u001b[30m \u001b[31mic\u001b[30m coronabeleid \u001b[31mgriep\u001b[30m \u001b[31moverleden\u001b[30m \u001b[31mvariant\u001b[30m \u001b[31mpositief\u001b[30m \u001b[31mdoden\u001b[30m \u001b[31mgetest\u001b[30m \u001b[31m#fvd\u001b[30m \u001b[31mdoor\u001b[30m \u001b[31mvaccinatie\u001b[30m \u001b[31mbesmettingen\u001b[30m \u001b[31mtesten\u001b[30m \u001b[31mziekenhuis\u001b[30m \u001b[31m#vaccinatie\u001b[30m coronatijd \u001b[31mna\u001b[30m #coronaprotest #covid \u001b[31mzorg\u001b[30m \u001b[31min\u001b[30m \u001b[31mziekenhuizen\u001b[30m #coronacrisis \u001b[31mcrisis\u001b[30m \u001b[31m#vaccin\u001b[30m coronadoden \u001b[31mbesmet\u001b[30m \u001b[31mopgenomen\u001b[30m \u001b[31mvanwege\u001b[30m \u001b[31mziek\u001b[30m \u001b[31mvaccineren\u001b[30m \u001b[31mcijfers\u001b[30m coronapatiënten \u001b[31mtest\u001b[30m coronaregels \u001b[31mover\u001b[30m \u001b[31mvaccins\u001b[30m \u001b[31mrivm\u001b[30m coronabesmettingen \u001b[31mlockdown\u001b[30m \u001b[31m“\u001b[30m \u001b[31mhuisarts\u001b[30m coronavariant \u001b[31mggd\u001b[30m #coronabeleid \u001b[31mtotaal\u001b[30m \u001b[31mverspreiding\u001b[30m \u001b[31m”\u001b[30m coronavaccins \u001b[31msterven\u001b[30m \u001b[31mpositieve\u001b[30m \u001b[31moverlijden\u001b[30m \u001b[31mpandemie\u001b[30m \u001b[31martsen\u001b[30m \u001b[31mhuidig\u001b[30m \u001b[31mkabinet\u001b[30m \u001b[31mouderen\u001b[30m \u001b[31mivermectine\u001b[30m \u001b[31mpfizer\u001b[30m \u001b[31mlijden\u001b[30m \u001b[31mbesmetting\u001b[30m \u001b[31mwe\u001b[30m coronatest \u001b[31mduitsland\u001b[30m #covid_19 \u001b[31mivm\u001b[30m \u001b[31mhet\u001b[30m \u001b[31mgevaccineerd\u001b[30m \u001b[31mkrijgen\u001b[30m \u001b[31mtijdens\u001b[30m \u001b[31maanpak\u001b[30m \u001b[31m\"\u001b[30m \u001b[31mmondkapjes\u001b[30m \u001b[31m2020\u001b[30m \u001b[31moversterfte\u001b[30m \u001b[31mtot\u001b[30m \u001b[31m#alleensamen\u001b[30m \u001b[31monderliggend\u001b[30m \u001b[31m#rutte\u001b[30m \u001b[31mminder\u001b[30m \u001b[31m⤵\u001b[30m \u001b[31mbedden\u001b[30m \u001b[31m#stopdelockdown\u001b[30m \u001b[31m#avondklokdebat\u001b[30m \u001b[31marts\u001b[30m \u001b[31mangst\u001b[30m coronaproof \u001b[31mgevolgen\u001b[30m coronasystemen \u001b[31mklachten\u001b[30m \u001b[31mdood\u001b[30m \u001b[31muitbraak\u001b[30m coronacijfers \u001b[31mjonge\u001b[30m \u001b[31mziekte\u001b[30m \u001b[31mgolf\u001b[30m \u001b[31mmaanden\u001b[30m \u001b[31mafdeling\u001b[30m \u001b[31mbijwerkingen\u001b[30m \u001b[31mliggen\u001b[30m coronavaccinatie \u001b[31mnederland\u001b[30m \u001b[31mzuid-afrikaanse\u001b[30m \u001b[31mlockdowns\u001b[30m \u001b[31mvarianten\u001b[30m \u001b[31mprivégegevens\u001b[30m \u001b[31minfectie\u001b[30m \u001b[31m#persconferentie\u001b[30m \u001b[31mbeleid\u001b[30m \u001b[31msterfte\u001b[30m \u001b[31mbuild\u001b[30m \u001b[31mpatiënt\u001b[30m \u001b[31m#rivm\u001b[30m \u001b[31mreset\u001b[30m \u001b[31m✔\u001b[30m \u001b[31mzorgpersoneel\u001b[30m \u001b[31mgezonde\u001b[30m corona-aanpak \u001b[31mcontrole\u001b[30m \u001b[31m'\u001b[30m \u001b[31mitalië\u001b[30m \u001b[31mondanks\u001b[30m corona-maatregelen \u001b[31msterfgevallen\u001b[30m \u001b[31mzich\u001b[30m \u001b[31mhoax\u001b[30m covid-patiënten \u001b[31m#vaccinatiestrategie\u001b[30m corona-uitbraak \u001b[31mmedische\u001b[30m \u001b[31mvaccinaties\u001b[30m \u001b[31mvd\u001b[30m \u001b[31mernstige\u001b[30m \u001b[31m💉\u001b[30m \u001b[31m📽\u001b[30m coronapandemie \u001b[31mhugo\u001b[30m \u001b[31msterft\u001b[30m \u001b[31mintensive\u001b[30m \u001b[31mgriepje\u001b[30m \u001b[31mavondklok\u001b[30m \u001b[31mprik\u001b[30m \u001b[31mpatienten\u001b[30m \u001b[31mbeschermen\u001b[30m \u001b[31mbetter\u001b[30m \u001b[31mlanden\u001b[30m \u001b[31mpcr\u001b[30m \u001b[31mscholen\u001b[30m \u001b[31mquarantaine\u001b[30m \u001b[31mcare\u001b[30m \u001b[31mmiljoen\u001b[30m \u001b[31m70\u001b[30m \u001b[31mwho\u001b[30m \u001b[31m23\u001b[30m \u001b[31monderliggende\u001b[30m \u001b[31moverlijdt\u001b[30m \u001b[31meconomie\u001b[30m \u001b[31m#vaccineren\u001b[30m coronawet \u001b[31mbestrijden\u001b[30m \u001b[31mimmuunsysteem\u001b[30m \u001b[31mnevenschade\u001b[30m \u001b[31m+\u001b[30m \u001b[31mdissel\u001b[30m #coronarellen \u001b[31mmeer\u001b[30m \u001b[31mblijkt\u001b[30m \u001b[31mmedewerkers\u001b[30m \u001b[31mstatistieken\u001b[30m \u001b[31mgeen\u001b[30m \u001b[31mzorgen\u001b[30m \u001b[31mmaatregel\u001b[30m \u001b[31m‘\u001b[30m \n",
      "computation time: 5977 seconds\n"
     ]
    }
   ],
   "source": [
    "for month in \"202003 202004 202005 202006 202007 202008 202009 202010 202011 202012 202101\".split():\n",
    "    top_tokens = show_top_tokens(month, QUERY_TOPIC_SMALL, spy=True)\n",
    "    save_top_tokens(top_tokens, month)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 26,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "20210601-23.out.gz\n",
      "20210601 vaccin corona vaccinatie vaccineren mondkapjes prik covid vaccins gevaccineerd pfizer virus moderna pandemie \u001b[31mmensen\u001b[30m ic \u001b[31m%\u001b[30m rivm \u001b[31mmiljoen\u001b[30m vaccinaties testen prikken \u001b[31mkrijgen\u001b[30m who \u001b[31m19\u001b[30m besmettingen bijwerkingen \u001b[31mna\u001b[30m test \u001b[31mspike\u001b[30m mondkapje coronavirus huisarts #vaccinatie \u001b[31mkans\u001b[30m \u001b[31mlienden\u001b[30m mrna besmet ziekenhuis \u001b[31mgeen\u001b[30m besmetting ggd astrazeneca \u001b[31masteroïde\u001b[30m gevaccineerden \u001b[31mdoden\u001b[30m lockdown \u001b[31martsen\u001b[30m mondkapjesdeal variant \u001b[31mlaten\u001b[30m 💉 \u001b[31mafspraak\u001b[30m \u001b[31mcijfers\u001b[30m \u001b[31marts\u001b[30m patiënten \u001b[31mziek\u001b[30m \u001b[31mdr\u001b[30m #mondkapjes immuunsysteem \u001b[31mstaatssteun\u001b[30m gevaccineerde \u001b[31musa\u001b[30m \u001b[31mzonder\u001b[30m \u001b[31mzich\u001b[30m \u001b[31meiwit\u001b[30m \u001b[31m1981\u001b[30m \u001b[31mwetenschap\u001b[30m verpleegkundigen \u001b[31mbrengt\u001b[30m \u001b[31m2020\u001b[30m getest \u001b[31mbreaking\u001b[30m \u001b[31mnemen\u001b[30m \u001b[31mkreeg\u001b[30m #coronavirus \u001b[31mjanssen\u001b[30m vaccinatiegraad \u001b[31m28\u001b[30m #corona \u001b[31mpropaganda\u001b[30m \u001b[31mgehad\u001b[30m \u001b[31mdissel\u001b[30m \u001b[31morder\u001b[30m virologen \u001b[31mgolf\u001b[30m #mondkapjesplicht \u001b[31mtv\u001b[30m \u001b[31manderen\u001b[30m \u001b[31mlicht\u001b[30m \u001b[31mseptember\u001b[30m \u001b[31mduitse\u001b[30m \u001b[31mwauw\u001b[30m niet-gevaccineerden virussen \u001b[31mvanwege\u001b[30m \u001b[31m2e\u001b[30m \u001b[31mnetwerk\u001b[30m \u001b[31mhoren\u001b[30m \u001b[31mgaven\u001b[30m \u001b[31msars\u001b[30m \u001b[31mwaarheid\u001b[30m \u001b[31moverleden\u001b[30m \u001b[31mmedische\u001b[30m \u001b[31msolidariteit\u001b[30m \u001b[31mmiljoenen\u001b[30m \u001b[31mwegzetten\u001b[30m \u001b[31maerosolen\u001b[30m \u001b[31mvanaf\u001b[30m uitbraak \u001b[31monbruikbare\u001b[30m \u001b[31mmaatregelen\u001b[30m \u001b[31mname\u001b[30m \u001b[31msteeds\u001b[30m besmettelijk \u001b[31mwuhan\u001b[30m \u001b[31mgriep\u001b[30m \u001b[31menkel\u001b[30m \u001b[31mtoegang\u001b[30m \u001b[31moverdragen\u001b[30m \u001b[31mrisico\u001b[30m \u001b[31m.\u001b[30m \u001b[31mtenpenney\u001b[30m viroloog #pfizer \u001b[31mafstand\u001b[30m \u001b[31mgeboren\u001b[30m #covid19 \u001b[31mcashen\u001b[30m \u001b[31mwetenschappers\u001b[30m injecteren \u001b[31m0,002\u001b[30m \u001b[31m0,0002\u001b[30m \u001b[31msymbolisme\u001b[30m \u001b[31mlogischer\u001b[30m \u001b[31m☄\u001b[30m #coronamaatregelen \u001b[31maz\u001b[30m quarantaine spuit \u001b[31m1,5\u001b[30m bijwerking #vaccinaties \u001b[31m11:30\u001b[30m \u001b[31mopzetten\u001b[30m \u001b[31moverwegingen\u001b[30m \u001b[31mdood\u001b[30m vaccinatieafspraak \u001b[31mbeurt\u001b[30m \u001b[31mlareb\u001b[30m \u001b[31mvoorkomen\u001b[30m \u001b[31m#beau\u001b[30m varianten \u001b[31mhelm\u001b[30m ziekenhuizen biontech \u001b[31m@0987654321id\u001b[30m \u001b[31m@010wijk11\u001b[30m spuiten astra \u001b[31m23.000\u001b[30m \u001b[31mgommers\u001b[30m \u001b[31mweken\u001b[30m coronacrisis #astrazeneca \u001b[31mopgerichte\u001b[30m \u001b[31mkinderen\u001b[30m \u001b[31mone\u001b[30m \u001b[31mifr\u001b[30m \u001b[31mnobelprijswinnaar\u001b[30m \u001b[31mdrukte\u001b[30m \u001b[31mopgenomen\u001b[30m \u001b[31mvertrouwt\u001b[30m #vaccineren \u001b[31moverheidsgeld\u001b[30m \u001b[31mzak\u001b[30m #vaccin \u001b[31mbijstandfraude\u001b[30m \u001b[31mondermijnend\u001b[30m \u001b[31maanbiedingen\u001b[30m \u001b[31mmails\u001b[30m \u001b[31monethisch\u001b[30m \u001b[31mscheppen\u001b[30m \u001b[31monrecht\u001b[30m \u001b[31mtegen\u001b[30m #coronavaccinatie coronamaatregelen geprikt #rivm \u001b[31mgroter\u001b[30m \u001b[31mtermijn\u001b[30m \u001b[31mvideo\u001b[30m \u001b[31mduitsland\u001b[30m \u001b[31minfectie\u001b[30m \u001b[31melens\u001b[30m #coronavaccin \u001b[31mzegt\u001b[30m \u001b[31mleveren\u001b[30m #covid19nl \u001b[31m\"\u001b[30m \u001b[31mstudie\u001b[30m \n"
     ]
    }
   ],
   "source": [
    "show_top_tokens(\"20210601\", PANDEMICQUERY, spy=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 25,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "20211201-23.out.gz\n",
      "20211201 corona #coronadebat gevaccineerd covid ongevaccineerden vaccins ic vaccin vaccineren ziekenhuis vaccinatie #vaccinatieplicht \u001b[31m%\u001b[30m gevaccineerden qr ziekenhuizen vaccinatieplicht \u001b[31mde\u001b[30m virus \u001b[31m19\u001b[30m 2g testen prikken besmettingen patiënten prik #coronatoegangsbewijs pandemie test oversterfte vaccinaties booster bijwerkingen #qrsamenleving \u001b[31mpositieve\u001b[30m \u001b[31mmaatregelen\u001b[30m \u001b[31maantal\u001b[30m boosters pfizer rivm variant \u001b[31mgesloten\u001b[30m gevaccineerde \u001b[31martsen\u001b[30m \u001b[31m20\u001b[30m \u001b[31mwerken\u001b[30m #corona vaccinatiegraad \u001b[31mnovember\u001b[30m mondkapjes mondkapje lockdown #2g ongevaccineerde \u001b[31mlaten\u001b[30m \u001b[31mmensen\u001b[30m #coronamaatregelen ggd \u001b[31mkomen\u001b[30m \u001b[31mwonen\u001b[30m \u001b[31m⚠\u001b[30m \u001b[31mdoden\u001b[30m \u001b[31mgolf\u001b[30m besmet \u001b[31mwe\u001b[30m zorg coronapas \u001b[31mvon\u001b[30m \u001b[31mleyen\u001b[30m \u001b[31m“\u001b[30m \u001b[31m#codezwart\u001b[30m \u001b[31m”\u001b[30m \u001b[31mvolledig\u001b[30m who \u001b[31mkuipers\u001b[30m \u001b[31mdata\u001b[30m \u001b[31mjonge\u001b[30m #coronavirus #coronapas \u001b[31maandeel\u001b[30m \u001b[31mopgenomen\u001b[30m \u001b[31mdissel\u001b[30m #qrcode \u001b[31monwettig\u001b[30m getest \u001b[31marts\u001b[30m \u001b[31mrechtbank\u001b[30m \u001b[31mstraks\u001b[30m pcr immuunsysteem \u001b[31m70\u001b[30m boosterprik \u001b[31msportschool\u001b[30m \u001b[31meerste\u001b[30m \u001b[31mcode\u001b[30m verpleegkundigen immuniteit \u001b[31mkinderen\u001b[30m coronapatiënten huisarts \u001b[31muitsluiting\u001b[30m patiënt \u001b[31mtegen\u001b[30m \u001b[31mmaanden\u001b[30m \u001b[31mcijfers\u001b[30m omicron \u001b[31mtoename\u001b[30m \u001b[31moostenrijk\u001b[30m \u001b[31meu\u001b[30m \u001b[31mwooow\u001b[30m boosteren \u001b[31mverplichte\u001b[30m 💉 #covid19 besmettelijk ongevaccineerd \u001b[31mmeer\u001b[30m omt \u001b[31mziek\u001b[30m \u001b[31mbezette\u001b[30m \u001b[31messentiële\u001b[30m \u001b[31mlanden\u001b[30m quarantaine \u001b[31m21\u001b[30m coronavirus \u001b[31mdringt\u001b[30m qr-code \u001b[31mwallonië\u001b[30m \u001b[31m80\u001b[30m \u001b[31mboetes\u001b[30m #vaccinatiedwang #2gnee \u001b[31m1,5\u001b[30m \u001b[31m-\u001b[30m \u001b[31millusies\u001b[30m #ongevaccineerden mrna \u001b[31mseh\u001b[30m \u001b[31malle\u001b[30m \u001b[31mgriep\u001b[30m \u001b[31mmachtsgreep\u001b[30m zorgsalarissen \u001b[31mrenden\u001b[30m \u001b[31mmotie-wilders\u001b[30m \u001b[31muitgezochte\u001b[30m omikron \u001b[31mcwz\u001b[30m \u001b[31mvoorstelt\u001b[30m \u001b[31m30\u001b[30m varianten coronadoden \u001b[31mziekte\u001b[30m \u001b[31mnu\u001b[30m besmetten \u001b[31mwinkels\u001b[30m \u001b[31msnelle\u001b[30m \u001b[31mernstige\u001b[30m \u001b[31mkans\u001b[30m coronabeleid \u001b[31mgezonde\u001b[30m \u001b[31mschuld\u001b[30m \u001b[31mhugo\u001b[30m \u001b[31mpolitieagenten\u001b[30m coronamaatregelen \u001b[31m🇦\u001b[30m \u001b[31mverplicht\u001b[30m \u001b[31mhoreca\u001b[30m \u001b[31msterfte\u001b[30m \u001b[31mrutte\u001b[30m \u001b[31mmijnwerkers\u001b[30m lockdowns \u001b[31mverzet\u001b[30m geprikt #boosterprik qr-codes \u001b[31mgeen\u001b[30m \u001b[31mperth\u001b[30m besmetting ic-bedden spuit \u001b[31mdwingen\u001b[30m \u001b[31monverklaarbare\u001b[30m \u001b[31mzonder\u001b[30m \u001b[31mwerkuniformen\u001b[30m \u001b[31mzwemles\u001b[30m \u001b[31mstigmatisering\u001b[30m cst \u001b[31mstudies\u001b[30m \u001b[31mniet\u001b[30m \u001b[31montnomen\u001b[30m #omicron \u001b[31mpositief\u001b[30m \u001b[31maustralische\u001b[30m \u001b[31mpostcodegebieden\u001b[30m \u001b[31mherhaalt\u001b[30m #oversterfte \u001b[31mstijging\u001b[30m \u001b[31mexperimenteel\u001b[30m verpleegsters \u001b[31mbedden\u001b[30m \n"
     ]
    }
   ],
   "source": [
    "show_top_tokens(\"20211201\", PANDEMICQUERY, spy=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 23,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "20220114-23.out.gz\n",
      "20220114- 2g qr corona lockdown covid besmettingen virus booster omicron \u001b[31mcode\u001b[30m vaccin quarantaine vaccinatie ic ziekenhuis mondkapjes testen ctb rivm \u001b[31m#nederlandgaatopen\u001b[30m vaccins #qrcode coronapas gevaccineerd prikken gevaccineerden \u001b[31mmaatregelen\u001b[30m pandemie ziekenhuizen \u001b[31momikron\u001b[30m \u001b[31m#persconferentie\u001b[30m ongevaccineerden mondkapje vaccineren \u001b[31mmeer\u001b[30m \u001b[31mmodellen\u001b[30m test lockdowns #coronamaatregelen \u001b[31mgeen\u001b[30m bijwerkingen #2g \u001b[31mopen\u001b[30m \u001b[31mzonder\u001b[30m qr-code vaccinaties besmet \u001b[31mgriep\u001b[30m \u001b[31mimmuunsysteem\u001b[30m prik 3g boosters \u001b[31m19\u001b[30m \u001b[31msamenleving\u001b[30m #stopdelockdown delta \u001b[31m\"\u001b[30m \u001b[31mziek\u001b[30m #lockdown #geenqr variant \u001b[31mtegen\u001b[30m oversterfte #coronapas \u001b[31maantal\u001b[30m immuniteit besmettelijk omt coronabeleid \u001b[31mpositieve\u001b[30m coronavirus boosterprik ziekenhuisopnames \u001b[31muitsluiting\u001b[30m coronamaatregelen \u001b[31mminder\u001b[30m \u001b[31mema\u001b[30m \u001b[31mnu\u001b[30m \u001b[31m#horeca\u001b[30m \u001b[31mmemo\u001b[30m \u001b[31mverwijderde\u001b[30m getest \u001b[31mverplichte\u001b[30m mondkapjesplicht gevaccineerde \u001b[31mduitsland\u001b[30m \u001b[31mnatuurlijke\u001b[30m #2gregel \u001b[31m‼\u001b[30m #corona \u001b[31mhoreca\u001b[30m vaccinatieplicht \u001b[31mdoor\u001b[30m besmetting pfizer \u001b[31mde\u001b[30m \u001b[31mchinese\u001b[30m \u001b[31mondernemers\u001b[30m \u001b[31mpolio\u001b[30m \u001b[31mcijfers\u001b[30m \u001b[31mmilitaire\u001b[30m \u001b[31mpositief\u001b[30m \u001b[31mengeland\u001b[30m \u001b[31m#horecaopen\u001b[30m \u001b[31m%\u001b[30m zorglast patiënten \u001b[31mstop\u001b[30m \u001b[31mper\u001b[30m \u001b[31mmensen\u001b[30m \u001b[31mwetenschappers\u001b[30m \u001b[31mamerikaanse\u001b[30m \u001b[31m#allesopen\u001b[30m \u001b[31marts\u001b[30m \u001b[31mdalen\u001b[30m boosteren \u001b[31montwikkeld\u001b[30m bijwerking \u001b[31m#kuipers\u001b[30m \u001b[31mschokkende\u001b[30m \u001b[31mhoger\u001b[30m #oversterfte \u001b[31montwrichting\u001b[30m \u001b[31m🔹\u001b[30m \u001b[31mna\u001b[30m \u001b[31mwerkt\u001b[30m \u001b[31m2018\u001b[30m qrcode besmetten \u001b[31m1,5\u001b[30m \u001b[31mmeldingen\u001b[30m pcr \u001b[31mjongeren\u001b[30m ongevaccineerde #stoplockdownnu \u001b[31mcafé\u001b[30m \u001b[31mhet\u001b[30m \u001b[31martsen\u001b[30m \u001b[31monderbouwing\u001b[30m #eindelockdownnu \u001b[31montwricht\u001b[30m \u001b[31mweg\u001b[30m \u001b[31mopnames\u001b[30m \u001b[31mlaten\u001b[30m \u001b[31mgehad\u001b[30m \u001b[31m#stopdewaanzin\u001b[30m #booster \u001b[31mwe\u001b[30m #omicron 1g \u001b[31mrechtmatigheid\u001b[30m \u001b[31mhanteert\u001b[30m \u001b[31mrestaurant\u001b[30m ggd \u001b[31mmet\u001b[30m \u001b[31mrecord\u001b[30m \u001b[31mgrafiek\u001b[30m \u001b[31mlager\u001b[30m \u001b[31msteeds\u001b[30m \u001b[31mgezondheid\u001b[30m \u001b[31minbreuk\u001b[30m \u001b[31mlareb\u001b[30m \u001b[31mverwerpt\u001b[30m #vaccinatieplicht \u001b[31mbuitensporige\u001b[30m \u001b[31mberoep\u001b[30m \u001b[31mfilmpje\u001b[30m \u001b[31mpats\u001b[30m \u001b[31mstoppen\u001b[30m #qrankzinnigheid \u001b[31mafstand\u001b[30m #coronatoegangsbewijs geboosterd \u001b[31mdragen\u001b[30m \u001b[31mdicht\u001b[30m \u001b[31mbiden\u001b[30m niet-gevaccineerden \u001b[31mdata\u001b[30m \u001b[31mwerkvloer\u001b[30m \u001b[31mondanks\u001b[30m #covid19 \u001b[31mbezoeken\u001b[30m \u001b[31mziekmakend\u001b[30m \u001b[31mdecember\u001b[30m \u001b[31mlanden\u001b[30m \u001b[31mernstige\u001b[30m \u001b[31mziekte\u001b[30m \u001b[31mdus\u001b[30m \u001b[31m#steundehoreca\u001b[30m \u001b[31mconclusie\u001b[30m 💉 \u001b[31mmaatregel\u001b[30m \u001b[31mrechtvaardigen\u001b[30m \u001b[31mlanger\u001b[30m ongevaccineerd \u001b[31minvoeren\u001b[30m besmette \u001b[31mkuipers\u001b[30m \u001b[31mrapport\u001b[30m \u001b[31m#omikron\u001b[30m \n"
     ]
    }
   ],
   "source": [
    "show_top_tokens(\"20220114\", PANDEMICQUERY, spy=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "20211029-23.out.gz\n",
      "20211029 qr coronapas de data code #coronamaatregelen #coronapas qr-code qrcode werk echte toegang tonen worden te ontzeggen vaccinatiegraad #apartheidspas tegen gevaccineerd maatregel qr-codes vaccinatiedwang #qrcode staatsgeheim informeren facto buitenlandse afpakken qr-pas directe keuzes families inkomen aanval hugo normale omt gevaccineerden actuele totale hun miljoenen gebruikt crisis grondrechten ollongren restricties nepnieuws-verwijt qr-discriminatie bevragen gratuit behelst geschonden manipuleert nederlanders punt blijkbaar verantwoordelijke tevoren uk alsof #fvd doel bereiken overlijdens maken pandemie % tijdens hand 75 ziekenhuisopnames ziekenhuisdata gedaan patiënten beste 80 reactie lancet nut kappen kans buitensluiten onderscheid prullenbak hieruit om #nooitmeerlockdown dus gezondheidsraad gevolgtrekkingen besmettingen discriminatie eindelijk slaat e wisten nergens horeca ziekenhuizen staan wanneer krijgt geven maart nooit allemaal zien “ #qrterreur bezwoer #hugodeqrankzinnige sorteren via enkel rechters 2020 aparte herhaling niet spanje haar test 17 meer plekken minister verboden basismaatregelen december invoering 🤥 #qrsamenleving woorden omdat apartheid en strijd #coronatoegangsbewijs positieve af weer #hugodejonge reden laat cijfers grondwet hardwerkende besmetten #geenqr terecht codes #vaccinatieplicht ” testen - regels uitbreiding advies weken konden #vaccinatiedwang sinds 1,5 overal pas horecaondernemer verlinken wordt ongevaccineerden vaccinatiepas komen groen verplichte ingevoerd bulan 🤦 🇹 geen tot is coronapaspoort geopenbaard sportscholen 🇮 vinkje opgelopen verbergt europees wil #corona een liegt uitbreiden mét #medischeapartheid parlement feest democratisch \n"
     ]
    }
   ],
   "source": [
    "show_top_tokens(\"20211029\", CTB_QUERY, spy=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "20211130-23.out.gz\n",
      "20211130 qr code 2g #coronatoegangsbewijs #qrsamenleving winkels eerste qr-code tegen #qrcode #coronapas cu achterban coronapas #2g zwart de essentiële ctb leden kamer motie uw niet-essentiële qr-codes #coronatoegangsbewijzen inzet gestemd voor bewijst stemmen #geenqr uitzien wíj dermate uitbreiden #coronamaatregelen vrijheidsberoving florida onwettig jl christenunie gehersenspoeld gevaccineerd partijcongres kamerfractie helaas leiders onbegrijpelijk 😮 heldere 🧵 zulke heeft onderstaande 3g nodig dienstverlening velen denken partij stemde schande ongekend knip mogelijk uitspraak wetsvoorstel europa aangegeven segers verlaten onbetrouwbaar waard rechtbank ziel ongevaccineerden steunt vrijheid #fvd stel controleren neus vaccins wallonië uitsluiting kapper om #qr straks gemaakt coronatoegangsbewijs zonder gevaccineerden € leven anders hadden weken opzeggingen qrcode aangenomen probleempje kamerfracties kamp overtreding binnen vandaag 1ste #qrankzinnigheid wetsontwerp werkstraf 3.600 weigeraars 7.200 omgezet weigering ontstaan lidmaatschap ' #2gregel vaccinatieplicht coronatoegangsbewijzen #coronavirus inperking voorwendselen gaat oostenrijk god discriminatie 20 beleid bestrijden kleren zoveelste weg uitbreiding opzeggen fysiek ingang verklaard apartheid credit ‼ systeem #christenunie maatschappij controle #qrindekliko #ctb vaccinatiepaspoorten vrijheidsbeweging gedaan het bijwonen events corona massaal protesteren uitsluiten boete besmettingen - bijzonder #hugodeqrankzinnige valse vaccinatiestatus te kost ernstige #qrterreur verbergen discriminerende fase #knettergek hakken klaar zwijgen codes november social inefficiënte introduceren bij invoeren dwang geldbulkwagens prik-automaat milieubewuste autobezitters laadpaal betekent zand daarover herstelbewijs #eerstekamer rooskleurig onzorgvuldig stop 3/2 \n"
     ]
    }
   ],
   "source": [
    "show_top_tokens(\"20211130\", CTB_QUERY, spy=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "202105 corona covid 19 - #coronamaatregelen #covid19 #corona #coronavirus coronavirus vaccin coronamaatregelen tegen coronacrisis vaccinatie griep #covid19nl maatregelen vaccins virus mensen #coronavaccin % overleden patiënten ic 2020 doden aantal na gevaccineerd positief #coronapaspoort coronapaspoort testen vaccineren coronavaccin ziekenhuizen ziek #vaccinatie ivermectine getest #covid ' ziekenhuis besmettingen artsen gehad door coronabeleid coronavaccins coronapatiënten rivm coronaregels bijwerkingen cijfers pandemie india via opgenomen klachten test tijdelijke vaccinaties zorg coronatijd besmet overlijden onderzoek besmetting #testsamenleving #coronadebat prik risico long miljoen wuhan wet pfizer worden coronavaccinatie #longcovid #coronavaccinatie jonge krijgen #coronacrisis #ivermectine sterven totaal #covid_19 paspoort angst huidig coronabesmettingen \" mondkapjes sterfgevallen de gevaccineerden lab #vaccinatiepaspoort gevaccineerde behandeling ( coronacijfers kinderen ifr vanwege coronaprik verspreiding #ikweiger wetenschappers ) crisis positieve miljard kans maanden geen ziekte uitbraak ernstige februari variant medicijn coronadoden coronapas #lockdown lockdown #vaccin #vaccineren ernstig amerikaanse pcr aerosolen sterfte jaar dissel immuunsysteem gevolgen gevallen gezonde coronatest covid-vaccins medische medicijnen dr overheid miljarden coronabesmetting weken alle geneesmiddel volledig groter texas coronapandemie opnames #testmaatschappij indirecte coronajaar #coronahoax infectie geregistreerd arts #fvd slechts versoepelingen huisarts oversterfte ivm beleid mogelijk coronasteun maart bewijs #vaccinatieplicht reizen behandelen besmettelijk ruim griepje landen #rivm wetenschappelijke influenza over minder geneesmiddelen immuniteit levensjaren aanpak who effect ziekenhuisopnames symptomen 65 #vaccinaties bereik 28 coronaproof \n",
      "202106 corona covid 19 #covid19 - #coronavirus #corona coronavirus vaccin #coronamaatregelen % kinderen vaccineren coronacrisis vaccinatie virus #coronavaccinatie ic #coronadebat gevaccineerd coronamaatregelen vaccins #covid19nl tegen na maatregelen overleden mensen coronaregels aantal doden griep positief gehad patiënten bijwerkingen ziek #vaccinatie prik 2020 getest testen crisis coronavaccin besmettingen ziekenhuis vaccinaties coronabeleid coronavaccins de ivermectine #voorweinigcovid rivm risico kans #coronacheckapp besmet artsen cijfers gezonde krijgen opgenomen \" who long overlijden ifr klachten pandemie #covid door fauci vanwege test variant sterven ziekenhuizen coronaprik besmetting jonge totaal jongeren maanden coronatijd coronabesmettingen jaar prikken coronavaccinatie positieve infectie antistoffen wereldwijd #vaccin #coronapaspoort huidig #coronacrisis #longcovid gevaccineerde coronadoden amerikaanse ( vrijheid dissel #mondkapjesplicht #fvd ernstige ziekte #coronavaccin volledig verspreiding tijdens coronapandemie angst waren wuhan coronapaspoort steeds ggd toegang app normaal coronapatiënten #ivermectine zorg jun vaccinatiegraad alle miljoen middel lockdown #vaccineren voorkomen arts delta ) het #deltavariant geen #wuhanlableak 0,15 coronacijfers koppelt groter onder effect mrna chinese per hugo behandeling coronasteun gekoppeld inenten miljoenen versoepelingen ernstig lab a-status gestorven vorig ic-capaciteit ' longcovid sterfte dood overheid nu 1,5 sterfgevallen gevaccineerden bewijs deltavariant india blijkt september worden #versoepelingen onderzoek paspoort afgelopen immuniteit ziekenhuisopnames afspraak #rivm via ouderen senaat + msm najaar ambtsmisdrijf experimenteel gezondheid #lableaktheory gevallen voorbij mondkapje studie coronacheck digitale \n",
      "202107 corona covid 19 #coronamaatregelen #corona #covid19 - #coronadebat #coronavirus coronavirus #coronavaccinatie tegen coronamaatregelen #covid19nl ic % mensen gevaccineerd maatregelen virus testen vaccin griep long vaccinatie positief ziekenhuis #vaccinatie besmettingen test aantal coronabeleid vaccins coronacrisis positieve getest ziek a-status #vaccinatieplicht doden de overleden patiënten vaccineren jongeren opgenomen coronabesmettingen kabinet #fvd ' ziekte klachten 2020 krijgen ifr kamer ziekenhuizen gevaccineerden na #longcovid status besmet over #covid longcovid rivm hugo jonge alle #persconferentie medische coronapatiënten #voorweinigcovid cijfers coronavaccin #nooitmeerlockdown lockdowns ebola gehad huidig ✅ debat pandemie vaccinaties #rutte coronatest ziekenhuisopnames coronacheck rutte 0,15 jul deltavariant door overlijden beleid + a bijwerkingen #astatus coronaregels sterven besmetting liggen a-lijst coronapas coronavaccins risico coronadoden #coronavaccin coronavaccinatie coronapaspoort pcr-test prik vws \" coronagevallen totaal volledig kans #coronahoax antistoffen #bvnl ministerie coronaprik #hugodejongekanniks geen vanwege app apartheid volksgezondheid pcr #hugodejonge gevaarlijk 1,5 #deltavariant infectie frankrijk varianten ernstig crisis #coronacrisis ( coronatijd #covid_19 ) long-covid 17.474 reset ventilatie who coronakaart demonstraties sinds het coronabesmetting #coronabesmettingen nieuwe ernstige ggd opnames verplicht #rivm experimentele ivm nu gevaccineerde leugen gevolgen artsen paniek stemming variant kamerleden cdc #frankrijk gevallen corr perc kinderen symptomen lijst coronadebat coronacijfers uk langdurige @2ekamertweets #testenvoortoegang stijging #vaccin #ikprikniet vrijheid mondkapjes tweede toegang speech 0 negatieve overlijdens influenza klaus \n",
      "202108 corona covid 19 - #coronavirus #coronamaatregelen #corona #vaccinatieplicht #covid19 coronavirus #coronabeleid gevaccineerd tegen % vaccinatie vaccin coronamaatregelen ziekenhuis ic virus mensen #covid19nl maatregelen 2020 patiënten overleden #coronadebat gevaccineerden griep vaccins #vaccinatiedwang testen ziek coronabeleid door #coronavaccinatie ) positief coronacrisis ( vaccineren ongevaccineerden coronapas #vaccinatie doden opgenomen ' ziekenhuizen #covid coronavaccins bijwerkingen gehad cijfers na artsen geen getest \" aantal coronapatiënten – besmettingen volledig test prik #vaccinaties ivermectine coronavaccin gevaccineerde rivm besmet #covid_19 long pandemie ernstige coronaregels ziekte infectie media worden #nooitmeerlockdown overlijden injectie gemeld ernstig beleid kans vaccinaties coronadoden krijgen sterven lockdown antistoffen alle risico “ huidig vanwege volksgezondheid #voorweinigcovid zorg klachten nasser #vaccinatieschade dood mohammad jonge #corona_origin_untold_secret positieve pfizer israël #longcovid ” kinderen meer imam gezonde #coronavaccin urk 1,5 voorkomen #bvnl coronapaspoort immuunsysteem arts pcr beschermd per denemarken jongeren overheid ziekenhuisopnames besmetting who australië #coronacrisis coronavaccinatie medicijn longcovid reguliere coronabesmetting sinds sterfgevallen het aug onderzoek influenza crisis coronabesmettingen miljoen immuun deltavariant #belangvannederland géén #medischeapartheid covid-patiënten medische hoax doodsoorzaak variant bpoc straf coronaprik gezond coronatijd totaal duitse #vaccin coronatest ivm lareb moderna prikken beschermt spoedwet #coronapaspoort bewijs 20 gezeik immuniteit bescherming opnames massaal cdc #vrijheid al-yemeni lockdowns verplicht sars-cov wereldwijd bild vóór duizenden overleed natuurlijke — 12 ongevaccineerd uitbraak lab \n",
      "202109 #coronapas corona covid coronapas 19 #coronamaatregelen tegen #coronadebat #coronapaspoort #vaccinatiedwang #coronatoegangsbewijs horeca #corona de coronamaatregelen #coronavirus coronatoegangsbewijs - maatregelen #medischeapartheid gevaccineerd vaccinatie invoering #qrcode virus ic mensen #persconferentie #covid19 griep ziekenhuis coronabeleid overleden coronavirus keijzer vaccin handhaven gevaccineerden vaccineren jonge % testen na #covid19nl kabinet coronacrisis qr-code medische prik ongevaccineerden #coronabeleid restaurant patiënten zorg alle geen #horeca #vaccinatieplicht géén 1,5 mona ongevaccineerd vaccins test volksgezondheid ziekenhuizen #qrsamenleving coronapaspoort coronaregels coronacheck verplichte #ikweiger verplicht samenleving ontslagen #hugodeqrankzinnige ziek vaccinatiegraad #fvd coronabewijs pas invoeren #apartheidspas staatssecretaris stemde vrijheid besmetten tweedeling #bvnl denemarken ingevoerd ) rivm doden #covid app coronatoegangsbewijzen positief noorwegen coronapatiënten controle zich apartheidsmaatschappij besmettingen zonder coronavaccin #nooitmeerlockdown getest beleid ( #monakeizer ' #vaccinatie culturele #covid_19 gehad long klachten vanaf apartheid terras dwang vaccinaties opgenomen antistoffen 25 zweden meer gevaccineerde vaccinatiedwang immuniteit risico sterven bijwerkingen kritiek controleren besmet vrijheden kinderen toegang worden #hugodejonge discriminatie vanwege geding 2020 coronacheck-app sector binnen toegangsbewijs ziekte artsen maes voeren inzake bewijs dus drang demissionair middel ivermectine opnames meerderheid systeem overheid overlijden maatregel hugo restaurants #rutte vinkje weigeren #vrijheid verworpen cijfers laten qr #longcovid ziekenhuisopnames paspoort patienten ernstig pandemie wordt grondrechten ‘ door #testenvoortoegang omt coronasteun #coronacheckapp meter stemming ondernemers september #coronavaccinatie infectie gezondheid maatschappij crisis \n",
      "202110 corona covid #coronapas #coronamaatregelen coronapas 19 ic gevaccineerd griep % ziekenhuis #corona #coronavirus coronapatiënten patiënten tegen maatregelen - de coronamaatregelen #covid19 #vaccinatiedwang vaccin ziekenhuizen testen gevaccineerden mensen coronabeleid zorg vaccinatie ongevaccineerden aantal opgenomen #qrcode overleden test virus liggen positief coronacrisis 2020 coronavirus #covid19nl ⁩ cijfers hugo #coronatoegangsbewijs gevaccineerde geen ⁦ coronatoegangsbewijs #persconferentie volledig na besmettingen ongevaccineerd vaccins #vaccinatieplicht ziek #coronapaspoort \" vaccineren bedden ” doden getest artsen gehad “ jonge ' besmet #covid volksgezondheid #qrsamenleving horeca bijwerkingen isala maanden druk ongevaccineerde vaccinatiegraad opnames positieve vaccinaties het waanzin huidig door crisis coronavaccins krijgen patienten rivm coronavaccin operaties probleem coronabesmettingen ziekte prik 1,5 umc totaal omt #ongevaccineerden #bvnl coronadoden sinds kans #hugodejonge klachten invoering coronaprik bewijs ernstige ggd oversterfte beleid steeds coronabewijs coronapatienten onderliggend ivermectine #covid_19 sterfgevallen vorig blijkt #vaccinatie schijnveiligheid overlijden long zweden pandemie angst laten covid-patiënten immuniteit vergelijken safe besmetting vrijheid maastricht infectie #zorginfarct grondwet #vaccinatieschade lijden ziekenhuisopnames okt lockdowns coronacijfers zonder risico influenza afdeling ticket patiënt werkvloer coronagevallen bron golf coronacheck coronapaspoort nu verplicht stop discriminerende ‘ grondrechten bescherming gevallen bruls lockdown tekort personeel ic's gezondheidszorg levi afgelopen coronaregels dus italië overheid besmetten #coronabeleid terwijl #hugodeqrankzinnige negatief bevolking vinkje uitbraak bezet strengere media extra discriminatie qrcode #longcovid toegang coronastatus \n"
     ]
    }
   ],
   "source": [
    "for month in \"202105 202106 202107 202108 202109 202110\".split():\n",
    "    show_top_tokens(month, QUERY_TOPIC_SMALL)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "20211031-23.out.gz\n",
      "20211031 #coronapas qr code coronapas bruls werkvloer qr-code werk #bruls #coronamaatregelen italië voor #qrcode #corona mijn toegang #coronavirus apartheid 👏 #geenqr verplicht #vaccinatiedwang geen tegen leerlingen wiskundedocent massaal leerlingenenquete bied skills gegevens #coronapaspoort gezonde men restaurant waarderen #vaccinatieplicht volksgezondheid pleit please maandag uitleg humor vaccinatiedwang beschikbaar om ingevoerd momenteel dus bleek retweet willen online dm meerderheid qr-pas overheid ingezet zonder welnee qrcode krijgen verbazen sprake de thuis 🇨 retweeten werkend coronamaatregelen besmettelijk mensen leven demonstreren miljoenen vaccinatiepas immuunsysteem verplichte werkgever ondertussen vanwege besmetten discriminatie protesten #qrsamenleving bevorderen restaurants krabbé #welevennietinchina inzetten #persconferentie weken blijf weigert #buitenhof zal wanneer zorginfarct accepteert controle breed durven jeroen leidt beschermen qrpas #jinek bereik tot gevaccineerden gezond we teugels gisteravond 2min strakker medische segregatie aantrekt problemen 30kg ranst afgevallen testmaatschappij vitaliteitsdeskundige econometrist toeten bankrekeningen ach steunpakketten algehele geblokkeerd daarbij steeds inzien maatschappij blazen reguliere onvoorstelbaar merkt #qrterreur werkgevers gebruik 82 keus . australië opnames werkelijk alsnog maak china wordt zomaar europa mogen ellende 1,5 #coronapassen certificaat beelden #coronatoegangsbewijs #verplichtevaccinaties covid-gezondheidspakket 200 maart zermatt helpt vragen pas 🇭 pijnstillers #ivermectine salvador invoering groter gestopt straat absolute veiligheidsberaad aard komen weigeren zink coronatoegangsbewijs heer #freechoice bewijzen zorg invoeren overlijdt onrechttv polarisatie walliserkanne graag dit plaatsen genomen vitamine horeca \n"
     ]
    }
   ],
   "source": [
    "show_top_tokens(\"20211031\", CTB_QUERY, spy=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "20210531-23.out.gz\n",
      "202105 ik alleen eenzaam ! voel me eenzaamheid gevoel mij voelen ben voelt mijn zo maar eenzame heb eigen n , #maatjegezocht wereld voelde als mag #eenzamejongeren #walgelijk moraalstoel anne zittende vergelijkingen maanden zelf recept pleindansjes #keeponmoving vijfjarig stilzitten pyjama dansje #vrijheid enz niveau ouderen meedoen frank onmenselijk dat kracht chinese kan mensen meisje mis leidt hoog slingeren piano lonely pest omdat doodeng bezoek #durftevragen hobby speel erg het corona i . beste joden duurt soms ziet 22 en muziek bejaarden leegte depressie steeds woon quarantaine m zit achtergelaten ookal vanavond niet vertrouwen lieve opgesloten erge denken.zat vermagering genas je schoonmoeder 5x overleed politiek bezig onderstaand 92 pure 😪 iedereen gestorven gevoeld rondom sommige doet 7 wensen drama vrienden bepaalde pijn herdenk triest vele vooral limburg-brabant boswandelingen doodsstrijd ongeventileerde kamers guur nederweert aangemaakt gevoelens gamen metal #eenzaamheid verstrekken chillen gepast gevoelig waarbij klein heel leven overleden festivals vrij contact duo meer deed sluipenderwijs indrukwekkender begrijp jarige bestaande jezus vaak medicijnen speciaal dodenherdenking verdriet prima heen invoelbaar ofzo verhaal gefrustreerd mezelf lessen #hemels boenwas 992 werklui 4s zwembadje sotheby crispy ayn lijdensweg picasso herlezen wereldleed speakers bach draadje limburg vergeten lege opsluiten dagen daarna zoek blauwe hoe terecht wzc janker liefde weekje \n"
     ]
    }
   ],
   "source": [
    "show_top_tokens(\"202105\", LONELY_QUERY, query2=IK_QUERY, spy=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "20210531-23.out.gz\n",
      "202105 ik blij gelukkig ben mijn heb dat geluk zo mij me mee happy maar heel word was zelf had gehad wel erg daar maakt niet echt ook ❤ vind wens moeder hoop gelukkige prik veel nog mn weer m'n kon nu ontzettend er 🥰 denk hoef met mezelf toen oprecht jou woon eindelijk leven kreeg voel . super krijg soms trots 🥳 mag dochter want beetje kan oudste ken positiefs maken vandaag wees profvoetballer lieve papieren depressie doorstaan 🙏 moderna weet even hem vagina 😊 injecties okay voetballen kwam hier zie chun dankbaar aardig excited verjaardag charcuterie kaasland altijd schotels kei srry montagnier bestellingen nobelprijswinnaar dingen support school vond fijn lief al combineren hardop voelt bestelling vrienden zó uitleg mondkap n werk 😅 twitter luc moest injectie zoon allernieuwste last i verblijden snap kaas zou kind pais _ vree goed creaties zag experimentele :) like persbedeling follow gekregen 875 voortzet meemaak verhult insta verlegenheid mocht dief ging verdedigen appie ✨ birthday muilkorf ouders 7,50 knijpen voetbal ze gerekend betekenen ga favoriet gestart zijn mrna afgelopen daaromtrent kracht gevoel enige gestemd 💖 merendeel medelijden sektarische action beschilderd geweest enorm keer flats serum anders dacht vader sta kinderen hierna pijn #happy zinvol \n"
     ]
    }
   ],
   "source": [
    "show_top_tokens(\"202105\", HAPPY_QUERY, query2=IK_QUERY, spy=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "20210131-23.out.gz\n",
      "20210131 thuis thuiswerken dicht werken , scholen winkels avondklok reizen 🤔 meer 👇 maar ... gewoon 80 gaat werkt werkgevers ouders procent bewijzen niemand wel rutte grote wat vrijheid #thuiswerken kantoor zegt werk minder #nieuwsuur er bedrijven kinderen ? thuiswerk kan door kantoren na niet corona hun kantoorleven kantoorgebouwen noodgedwongen vlucht coronacrisis ná toekomst thuiswerkers structureel thuiswerkende genomen school kantoorplannen files hetzelfde ✔ gevolgen 25 zo'n denken kids sluiting grootste schrappen werkgever informatie ✅ kunnen kantoorruimte verkeer % mogelijk thuisonderwijs online thuiswerkt nodig ontmoeten wifi besmettingen zoom #coronamaatregelen #alleensamen kpn nederlandse ⤵ ziggo #scholendicht tips werkplek verantwoordelijk snelweg gaan verplichte werknemers contacten belangrijke gewerkt zwaar #corona verplicht verspreiding zal onderdrukking blijf afstand thuiswerker kantoorpersoneel effectief fijn knopen infrastructuur voorkomen lockdown thuiswerkplek #werkgevers #bnrmobility gevolg investeren tijdens meeting sluiten #coronavirus mogen morgen combineren nieuwsuur schoolsluiting file asfalt evt huishouden winkel ontmoet oneens janken ipv houden extra leerlingen volhouden bezorgen prioriteit blijven veel energie bso #a27 belangrijk #kantoorruimte #amelisweerd min klagende werkomgeving naderende ⬆ stelling coronavirus achterhoofd bevorderen weer vergaderen teams conclusie werknemer coronapandemie leerkracht minister meet drukker moesten voordeel :/ hoef wegen iedereen etc zolang langs oproepen beperken helder klas houding woningen nergens #longziekte #ergotherapeut stoplichtmodel lepeltheorie thuiswerktips nicolaashoeve thuiswerkvakantie #hotel #benb \n"
     ]
    }
   ],
   "source": [
    "show_top_tokens(\"20210131\", \n",
    "                r'(thuiswerk|thuis.werk|werk.thuis)', \n",
    "                spy=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "20210331-23.out.gz\n",
      "20210331 avondklok uur 22.00 #avondklok later gaat langer vanaf buiten informatie licht 4.30 31 maart #alleensamen blijft om woensdag ⤵ per tussen steeds omdat de maatregelen ga vanavond dodelijker én 22:00 meer uurtje ️ naar tijdens vandaag lockdown 21.00 besmettelijker strengere gemodelleerd #britsevariant vk kuipers bewijs gek #coronamaatregelen < eerder gehouden straat geldt minder 20 #lockdown zoals we paar ; zei alcohol schenkt terug dagen tot 22 allemaal 04:30 verboden 21:00 woensdagavond #corona boete dus zonder 21 beroep dicht frankrijk #pvda nachts 10 npo ingang effect variant mondkapjes groepjes net avonds andere winkels politie • patrouilleren klok ronde s drank overtreders handhaven hapjes rechter rusten 04.30 ramadan zomertijd 🚨 verlenging marcheert opeenvolgende boekarest roemenië britse #museumplein strenge tij protesten t politiestaat afschaffen verschoven wet ochtends menselijk ingrijpen lanceerde liveblog overal jongeren partijvoorzitter avond ingaan volgt april álle misleiding uitspraak mogen pro hoger kleine vonnis 19:00 regels demonstratie burgemeester % verzetten #mondkapjes masker overdag handhaving #zomertijd thuis trivialiseert daarom 🕙 #hogerberoep holocaust terrassen staatskas vergeet belgische graaien vergelijkt #op1npo werkgeversverklaring directe beboeten geldige #politie vooruit verklaring maatregel verzet wo #sp uitlaten 5.000 19u derde tijdstip spoedwet breaking aanvechten haarlemmermeer tien ii jullie #denhaag gebeurt nergens eraf basis betaalt pasen 22u negeren ergens \n"
     ]
    }
   ],
   "source": [
    "show_top_tokens(\"20210331\", \n",
    "                r'avondklok', \n",
    "                spy=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "202104 corona covid 19 #coronamaatregelen #corona #coronavirus #covid19 #coronadebat - ic coronacrisis coronamaatregelen mensen vaccin coronavirus maatregelen coronabeleid patiënten griep #covid19nl virus aantal % doden tegen ziekenhuis door overleden coronapaspoort ziekenhuizen besmettingen #coronavaccin coronapatiënten risico positief vaccinatie de #lockdown testen kans #vaccinatie na india #covid vaccineren getest opgenomen vaccins coronaregels gevaccineerd overlijden 2020 cijfers rivm #zerocovid #stopdelockdown #versoepelingen coronatijd test bijwerkingen #fvd jonge #coronacrisis sterven liggen coronadoden coronavaccin lockdown totaal #coronapaspoort tijdens jaar coronabesmettingen beleid miljoen zorg ziek #rivm #astrazeneca klachten besmet kabinet meer ifr positieve crisis maanden sterfgevallen avondklok pandemie huidig coronawet #longcovid sterfte vanwege coronavaccins #avondklok ivermectine besmetting ernstige ziekenhuisopnames wappies variant coronatest apr bevolking artsen alle open coronacijfers 60 #zerocoviddayofaction in influenza coronavaccinatie #vaccin mogelijk stap coronapandemie bedden ic's brazilië gehad #coronavaccinatie angst ziekte coronaproof \" medische krijgen terrassen we opnames landen covid-patiënten dood long + wereldwijd verspreiding versoepelingen #testsamenleving nauwelijks duitsland effect golf hcq ivm inwoners coronaprik britse 🚜 regels dissel griepje ) lockdowns patiënt #ivermectine 📽 per infectie #horeca immuniteit #terrassen medicijn overlijdens aan sinds ondanks voorstander gevolgen ernstig demonstratie versoepelen buitenlucht aanpak omt vrijheid weken ( onderzoek vorig trombose #covid_19 verspreiden vaccinaties tijdelijk hydroxychloroquine overheid 70 intensive zweden behandeling strategie gestorven minder ⤵ gevaccineerden oversterfte \n",
      "202103 corona covid #coronamaatregelen 19 #corona #coronavirus coronamaatregelen griep #covid19 coronavirus maatregelen virus coronacrisis tegen - coronabeleid #coronadebat mensen vaccin #lockdown vaccineren #covid19nl doden #fvd #coronapaspoort aantal volstrekt testen #coronavaccin coronapaspoort #avondklok ic door dissel overleden na 2020 risico positief % 70 besmettingen bijwerkingen stierven getest patiënten vaccinatie vaccins de 151 lockdown immuunsysteem coronavaccin coronaregels < besmet ziekenhuis #stopdelockdown rivm ziekenhuizen coronatijd test ifr ziek verwaarloosbaar onbekend vergelijkbaar eraan coronabesmettingen heftig ) #vaccinatie duidelijkheid kabinet #covid vaccinaties cijfers zorg jaar termijn-effect lockdowns beleid het coronacijfers ( coronavaccins opgenomen crisis coronadoden \" angst lange #coronacrisis absoluut overlijden artsen klachten #persconferentie #zerocovid sterven gevaccineerd jonge #stemnederlandterug mrt #longcovid vanwege jaap huidig ggd regels sterftecijfer pandemie ivermectine coronapatiënten coronaproof effect golf coronatest ' positieve hugo quarantaine motie #rivm vrijheid corona-uitbraak who sterfgevallen alle sterfte gehad variant “ in besmetting tijdens ivm corona-maatregelen dood meer #vaccin ” verspreiding ondanks laten overheid vrijheden 📽 totaal groep grootte duitsland via medicijnen haga ernstig weken volksgezondheid fragment modellen ernstige reguliere coronaprik avondklok vitamine houden waanzin aanpak + oversterfte kans coronapandemie heftige we uitbraak coronasteun long #ikspreekmeuit #mondkapjesplicht huisartsen #vaccinatiepaspoort ondernemers griepje impact d belgië versoepelingen massaal waarin lijden onderzoek strenge gezonde zware wybren medewerkers bevolking zich griepprik tijdelijke \n",
      "202102 corona covid #coronamaatregelen 19 #corona #coronavirus #covid19 coronamaatregelen coronacrisis coronavirus - maatregelen #lockdown #coronadebat coronabeleid virus vaccin overleden #covid19nl doden tegen #avondklok de besmettingen griep na patiënten vaccinatie door #fvd ic aantal mensen testen getest #coronavaccin coronatijd ziekenhuis positief #vaccinatie coronaregels crisis jongeren #coronaregels variant % cijfers kabinet \" lockdown #stopdelockdown ziekenhuizen #persconferentie zorg coronavaccin besmet #stemnederlandterug onderliggende ' overlijden #coronacrisis vaccineren artsen wetenschappers angst vaccins kinderen #zerocovid ouderen rivm britse ziek coronavaccins 2020 gedwongen duitse coronapatiënten #covid pandemie medicijn sterven we coronabesmettingen opgenomen test klachten dood vanwege #longcovid worden omt #vaccin positieve golf beleid lijden aanpak coronadoden regering coronaproof totaal #rivm coronaprik oversterfte het huidig maanden tijdens uitbraak #horeca bijwerkingen isolatie coronapandemie onderliggend varianten corona-uitbraak coronacijfers coronawet via jaar coronatesten jonge gevaccineerd kluytmans zweden ⤵ #coronabeleid krijgen ) israël mondkapjes #vanangstnaarvertrouwen #alleensamen rutte ivm onder verpleeghuis who corona-maatregelen schade gevolgen coronatest feb verspreiding huisarts vaccinaties quarantaine prik coronavaccinatie reguliere ggd influenza avondklok verlengd gezonde overheid ivermectine long ernstige pfizer #omt open alle ziekte minder #ivermectine #rutte over #persco maart aan groter scholen 352 hugo weken griepje bewoners #covid_19 geregistreerd #lockdownnl besmetting winkels overbelast umcg #kutcorona medisch arts omt-lid sterfgevallen coronabesmetting corona-aanpak #scholenopen overlijdens intensive bewijs gemeentes mutatie wereldwijd 70 \n",
      "202101 corona covid 19 #corona #coronamaatregelen #coronavirus #covid19 #coronadebat coronavirus - coronamaatregelen coronacrisis vaccin maatregelen #avondklok #lockdown virus #covid19nl tegen aantal #coronavaccin % mensen britse coronavaccin patiënten de ic coronabeleid griep overleden variant positief doden getest #fvd door vaccinatie besmettingen testen ziekenhuis #vaccinatie coronatijd na #coronaprotest #covid zorg in ziekenhuizen #coronacrisis crisis #vaccin coronadoden besmet opgenomen vanwege ziek vaccineren cijfers coronapatiënten test coronaregels over vaccins rivm coronabesmettingen lockdown “ huisarts coronavariant ggd #coronabeleid totaal verspreiding ” coronavaccins sterven positieve overlijden pandemie artsen huidig kabinet ouderen ivermectine pfizer lijden besmetting we coronatest duitsland #covid_19 ivm het gevaccineerd krijgen tijdens aanpak \" mondkapjes 2020 oversterfte tot #alleensamen onderliggend #rutte minder ⤵ bedden #stopdelockdown #avondklokdebat arts angst coronaproof gevolgen coronasystemen klachten dood uitbraak coronacijfers jonge ziekte golf maanden afdeling bijwerkingen liggen coronavaccinatie nederland zuid-afrikaanse lockdowns varianten privégegevens infectie #persconferentie beleid sterfte build patiënt #rivm reset ✔ zorgpersoneel gezonde corona-aanpak controle ' italië ondanks corona-maatregelen sterfgevallen zich hoax covid-patiënten #vaccinatiestrategie corona-uitbraak medische vaccinaties vd ernstige 💉 📽 coronapandemie hugo sterft intensive griepje avondklok prik patienten beschermen better landen pcr scholen quarantaine care miljoen 70 who 23 onderliggende overlijdt economie #vaccineren coronawet bestrijden immuunsysteem nevenschade + dissel #coronarellen meer blijkt medewerkers statistieken geen zorgen maatregel ‘ \n"
     ]
    }
   ],
   "source": [
    "for month in \"202104 202103 202102 202101\".split():\n",
    "    show_top_tokens(month, r'corona|covid', spy=False)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 21,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "20201130-23.out.gz\n",
      "20201130 ik je heb ben mij me jij jullie dat jou ga vind als bent hebt jouw mijn weet kan denk niet doe wil kun dan moet maar krijg kunt snap ook had wel ken wat zou hoop dacht zie voel hou kom hoe begrijp bedoel vond zeggen zeg of echt stel iets geloof mee hoor iemand wilt neem lijkt benieuwd vraag doorgeleerd mag al hooggeleerd bijltjes-dag daar achten smetten wens blij sta hoef zo jezelf hem zag soms mezelf word dank denken ? mn hoogleraren nogal maak aangifte dm zelf m'n woon belgie-japan dus gewoon doen even nou wist probeer verwelkomen verteld gevraagd vermoed meenemen ooit sorry uitgemaakt zit meld kon gestemd altijd vreugde graag gwn eens volgens stuur gehad wilde mis bang haat durf zal zei houd toen gelezen racistische fantastische noem ergens geef las dingen vrees volg loop geleerd gezien kinderspel volgers kreeg verwacht app want oogst gevoelige voorstellen helemaal vrij vader ogen omdat ommetje gedaan gevoel #detijdloze hoorde vaak haal examens persoonlijk bed types missen hoi bewijs keer oprecht gronde kies draag lachwekkend officieren nooit merk tweet opiniemakers meen gemist zíj da fractie spreek messenger schrijf twitter belasting-benaming idee tweets niks ́ kmoet uitblazen snel 19u45 betaal \n"
     ]
    }
   ],
   "source": [
    "show_top_tokens(\"20201130\", \n",
    "                r'\\b(ik|mij|me|jij|je|jou|jouw|jullie|heb|ben)\\b', \n",
    "                spy=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "20210331-23.out.gz\n",
      "20210331 gelukkig blij ben ik dat zo word mijn ⛔ mee had goed weer mag heel met 💐 gehad 🥊 buurman beleving js werkend gekomen immuunsysteem maakt zes ️ arrogantie 🥰 sprak kus geven nog er gebeurd heb pieter schuld erg 🧚 eruit ♥ 🤭 💞 ander geleden 🧡 bundel boerenverstand #eenvandaag #nosjournaal macht 🥳 publiek gehouden maar ✨ wel arm slechts gezond mus m'n #pvv-team ! zijn 🍀 wees vanaf 4e ziet prive klant kiesraad 🌞 waanzinnig iemand digitale onderzoeksjournalistiek geïnformeerd paar eer ons superblij gestemd eindelijk oom opvallende niet eigen uitstaat flitsen reageren #ziekvoorgezond #groenlinks meuk speciaal noodgedwongen #vergeetonsniethugo langskomen situatie strot dode 😀 bijdruk geldsysteem pro-eu vergis hefboom super aangesloten week opblazen munt rammen ♀ aanvallen waardoor maatregelen 🏻 versneld moedwillig oh #ollongren bank geleid #vivaldi #pandemiewet blijft slag grafiek oprecht ontzettend ze x tevreden trots echt centrale bepalen daar stond #moslimbroederschap 👊 50.000 hij mensen haar hebben net schikt gestelden documenten gestrekt eisenjij vonnis luistert lot huishoudelij europese was :) gezicht #zorg tv autosector bent 1000x 20.35 passion eerlijk wilders 1e opgeslagen allen onzin medisch patiëntendossier z'n bende even gister verandermanagers wacht aangeven bestuurskundigen #klimaat grenzen 😁 geruststelling yess buren wens #onzeboerderij zorgt been 😃 we steun \n"
     ]
    }
   ],
   "source": [
    "show_top_tokens(\"20210331\", \n",
    "                r'\\bgelukkig\\b|blij\\b', \n",
    "                spy=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "20210331-23.out.gz\n",
      "20210331 alleen maar niet kijken uur 22.00 breder lockdowngekkies regeringsnarratief water praten uitpraten #alleensamen nog nat ruimte heen nam besmettelijkheid kunt hem . licht actieve onmogelijk bewijst pcr besmetting alsnog mensen eens 4.30 avondklok anders blijven langer haga klaar ⤵ krijgt #fvd testen stellen negatief om waarom staat door paar eenzaam doctrine steunt dbrink informatie gaat als ollongren kulldorff omdat test #vvd want eenzame dat arib blijft later vast vaccinaties chaos héle ellende stuurt rede scoringsdrang halfwassen rondvliegen stukslaat stijf pilletje adrenaline aangericht nederland #gaswinning onacceptabel miljard daling verwijten virus plussers mafkezen #landverrader 24u rondliep klungelig schokkende totale blaadje verzoek #groningen zwaar rest kiezers smsjes opvallend mutatie viel buitenspel eenzaamheid mag rustig gebruiken kosten gespreksverslagen ratten vuurwerk eventjes ambtenaren notities aanpakken 1,5 werkt vuile incl rutte kreeg meer havardwetenschapper risico-ouderen doortweeten akeligs kijkt elites instemming opvatting “ geen tweette 60 zette appjes buiten betalen ook bekendheid fvd woensdag niemand acht fysiek kabinet wij martin waarheid lockdown kant kunnen #verkenners gesnoerd 31 baas dierenkwakzalver 18-55 > verhaal gehuild afgewezen openbaar klachten gaande kwetsbaren moeten ab alleenstaande opsluiting leeftijdsgroep brief voel gevaar gebruikt min steeds kan andere besluiten omtzigt spel vertegenwoordigers hoekstra afzeikt #ahmadrezadjalali #executie #saveahmadreza #freedjalali ogen mevrouw denken fysieke alstu betuttelracisme \n"
     ]
    }
   ],
   "source": [
    "show_top_tokens(\"20210331\", \n",
    "                r'eenzaam|eenzam|alleen', \n",
    "                spy=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "20210331-23.out.gz\n",
      "20210331 bang angst paniek moe stress ben brazilië corona te doden osterhaus ab jullie angstporno = pandemie lege tijdens paniekzaaier 2009 man podium miquel kostelijk ekkelenkamp buitenproportioneel #coronavirus pompen schoot jan-en-alleman verbergen veranderd samenleving buitenlucht lichaamstaal borstademhaling duidend psychoanalyse tong bijten 318500 2110000 0,0015 17.500000 0,00094 inwoneraantal 16.500 leg 4000 lippen 00 snapt eerdmans boetes blik heden houding oké #functieelders gedeeld zaaien mensen delen overleden hou depressie plek regels gesloten tent lezen gespreksverslag maanden al iedereen deel hallen ontvangst niets gelijk angstig 100.000 openbaring word #op1 hiervan belangstelling maken jennifer publicatie #meidenvanholland #gekkepetertje return enig vóór inzien aanpraten nederland cijfers meiden holland eerlijkheid depressief #depressief geven #functieeldersdebat bange #positieomtzigt deze bent angstzaaier wybren zogenaamd graaiende en medew binnenkort super gesmeekt schoppen maatregels aankaartte campagne ollongren onnodig diep opsluit zorgen #ploumen #klaver tij vaccinatielijn gewaardeerd spuiten flag false voor nep verliezen paniekerig golf #kaag herhaal duimen mexicaanse hebt viroloog geheugenverlies angstoornissen 👿 virologen #arib banger leven jagen keren armen peeuw houden liegen paniekerige paniekzaaien bevolking kennelijk opheldering geschorst langst draaien b12 lopen bizarre angsthazen vaccinatie vegen om curve leggen moves gemotiveerde extreem langste ofzo heeft kostbare ranst blubber imaginaire oppoetsen oprakelen griep prikkers raken relaxed tb slangenkuil imploderen controle 900.000 ratje intro \n"
     ]
    }
   ],
   "source": [
    "show_top_tokens(\"20210331\", \n",
    "                r'moe\\b|depress|depri|bang|angst|paniek|stress', \n",
    "                spy=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "20201130-23.out.gz\n",
      "20201130 mijn ouders familie vader moeder vrienden collega zoon vriend haar dochter vriendin collega's hij broer zoontje > kinderen oma vond opa damien school zus vriendelijk neef tante hun ’ (8 verandert vriende boze kind buurman compensatie n zíj belasting-benaming ́ ogen vriendjes onzin aaahhh messenger cringe altij zooooo vibe neefje ic denkt waarbij #trotsedochter taakstraffer #rommel mn nice erin ruimt bedankje gij families vriendinnen gedeactiveerd erom hahahaha bernhard vrijwillig vriendje lijn mooier ontvangen schouders moeders anonieme gevraagd ze lacht gooien uitgemaakt ma da das absurd altijd vriendelijke idioten ge diegene haat irritante broertje leo oom belasting genoemd twitter soms ⃣ ;) ni 🤬 wezenlijke melden getreiterd stereotypen lucassen bericht uitgekeerd hiermee anti-joodse nicht naar najaarsnota uitvoeringskosten prinsjesdag zusje sturen dm 235 boreale . . . . . hem ik dezelfde stukje kerst ♥ m'n micha straat naam unfocussed emine middelvinger en samen wist account kinderopvangtoeslagschandaal kijken dochters rest focussen paniek kleinkinderen \" ideologie kleinzoon 🤍 chef vaders herenigd yazidi-meisje isis-terroristen gevangenschap christina vastgehouden kutkarakter inhalige werd ziet vieren emotioneel functie machteloosheid vriendschap miljoen mag beschrijven neefs dofheid 🇱 darth die momenten waarin man lichtje . . . . van staren slavernij uw vaderland broers dikke én bevriend #yazidigenocide s 🇮 leuk lichtjes statushouders meermalen zei jack wapens \n"
     ]
    }
   ],
   "source": [
    "show_top_tokens(\"20201130\", \n",
    "                r'\\bopa|\\boma|vader|moeder|ouders|broer|\\bzus|zoon|dochter|kleinkind|\\boom|tante|neef|nicht|zwager|schoonzus|familie|'\n",
    "                'vriend|\\bmaat\\b|\\bmaatje|\\bmaten\\b|collega|\\bbaas|chef|buurman|buurvouw|\\bkennis', \n",
    "                spy=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "202002 coronavirus #coronavirus corona virus het china #corona in besmet #covid19 italië nederland mensen quarantaine @rivm - griep #coronavragen besmetting chinese besmettingen over covid doden patiënt uitbraak 19 rivm paniek ggd ziekenhuis #covid2019 eerste #covidー19 huisarts vragen vanwege via wuhan verspreiding #coronavirusnederland maatregelen gevallen : #rivm vastgesteld corona-virus #coronavirusoutbreak informatie coronapatiënt tilburg getest #virus contactonderzoek patiënten chinezen handen #china weken aantal dagen besmette symptomen iran cruiseschip grens update ziek @coronanederland mondkapjes met contact voorkomen viroloog wereldwijd wassen landen #covid #wuhan blok iraanse #coronavirius @telegraaf nieuwe europa nieuws @nos bruins epidemie positief gestart overleden breaking 😷 mondkapje pandemie testen who grenscontroles incubatietijd besmettelijk gevolgen bang nederlanders coronabesmetting gezondheid minister chinees japan #covid_19 bevestigd passagiers geval besmetten lombardije nu grenzen aziaten ziekte arts limburg dood 2 zorgen bekend volksgezondheid diemen belgië mogelijk advies italiaanse bereikt #iran hoest hoesten test #tilburg griepje nieuwsblik verspreid kruis voorbereid vaccin kwestie @nunl virussen syriërs vanavond tenerife geweest angst #ncov uitgebroken reisadvies erdogan verspreiden zich beurzen labonderzoek genezen diagnose loon zand hotel #patiënt dokter nos https://t.co/80e0vcvakc https://t.co/ugqpidj7tl erger duitse getroffen frankrijk nl uitzending #schengen dodental overheid #sarscov2 italie autoriteiten uitzonderingsclausules aanknopingspunten antwoorden coronagebied tegen doorgang noord-italië risico https://t.co/dqwxcypmon westerdam heeft meldt 🦠 #coronavirusus uit #italie dode verspreidt zieke \n",
      "202004 corona #corona #coronavirus coronavirus coronacrisis #covid19 covid #coronacrisis 19 de - door in #coronadebat tijdens #coronavirusnederland #covid19nl maatregelen #coronanederland #coronavirusnl over tegen coronamaatregelen patiënten coronapatiënten overleden aantal virus #coronamaatregelen doden #covidー19 #covid_19 #samentegencorona getest vanwege crisis van via coronatijd ziekenhuis gevolgen besmet #covid2019 strijd ic china corona-app testen verspreiding #rivm tijden \" kabinet aan #blijfthuis het aanpak #alleensamen #coronaapp nederland #covid19be we #coronanl rivm samen coronadoden #rutte ouderen getroffen ' zorg corona-crisis positief mensen #covid2019nl impact cijfers controle onderzoek zorgmedewerkers beschermingsmiddelen #pvv app huisarts #wilders griep opgenomen corona-apps onze verpleeghuizen rutte ziekenhuizen coronatijden #mondkapjes uitbraak besmettingen pandemie #nieuwsuur #covid #lockdown patiënt artsen bestrijding na informatie vaccin krijgen wetenschappers economie test bestrijden gemeld intensive miljoen lees overlijden genezen onder #persconferentie mondkapjes mogelijk sterven wereldwijd besmetting verband minder motie ziek 👇 ivm zich ondanks hulp minister nl overledenen sterfgevallen personeel oproep zorgverleners #samensterk medicijn care jongeren corona-virus totaal blijkt v februari update corona-maatregelen weken apps verpleeghuis risico klachten ziekenhuisopnames coronabonds coronapatiënt arts & worden zorgen alle symptomen afgelopen economische legt tekort #china ondernemers voorkomen #op1 sinds bedrijven 👉 medische inwoners sterfte coronaregels corona-patiënten corona-doden april nieuwe onderzoekers @rivm wuhan kwetsbare experts thuis landen 1,5 #zorg privacy ggd coronapandemie + kamer zorgpersoneel #nederland \n",
      "202005 corona #corona covid coronacrisis coronavirus 19 #coronavirus #covid19 #coronacrisis - de door #coronamaatregelen in maatregelen virus tijdens coronamaatregelen crisis doden coronatijd #coronadebat #covid19nl aantal patiënten tegen #covid_19 #coronavirusnl testen vanwege besmet kabinet via rivm #coronavirusnederland over #rivm #covidー19 besmettingen #coronanederland onderzoek coronadoden getest coronapatiënten we vaccin gevolgen #samentegencorona overleden tijden verspreiding aanpak #lockdown griep ziekenhuis 👉 na impact versoepeling ic rutte @rivm huisarts lees stap corona-crisis lockdown zorg minister china #persconferentie uitbraak #alleensamen meter test artsen #covid coronabeleid positief pandemie chinese #mondkapjes inwoners nieuwe getroffen hydroxychloroquine miljoen elens ouderen #rutte strijd coronabesmettingen \" besmetting coronatijden controle #covid__19 van klachten coronaregels mei afstand weken mensen 1,5 mogelijk coronaproof cijfers ' @hugodejonge #covid19be het ondernemers verpleeghuizen nl economie zink #stopdelockdown extra ziekenhuisopnames anderhalve corona-uitbraak nederland bedrijven vanaf ziekenhuizen ⤵ ivm overheid ondanks blijkt corona-maatregelen update medicijn who genezen ziekte onze onafhankelijk corona-app #hydroxychloroquine ziek oversterfte verband patiënt informatie corona-aanpak personeel arts wuhan persconferentie versoepelen duitsland #anderhalvemetersamenleving coronagolf nauwelijks gevolg per golf dissel onder symptomen medewerkers meer @minpres krijgen juni ggd 🇪 landen @kellyovercorona samen wereldwijd risico voorwaarden gezondheid slachthuizen maart versoepeld gates premier #nieuwsuur chloroquine maanden piek #wilders economische horeca 27 met verspreidt veehouderij overlijden sinds zich worden #anderhalvemeter middel sterfgevallen #blijfthuis \n",
      "202006 corona covid de #corona 19 coronacrisis coronavirus #coronavirus #covid19 in coronamaatregelen maatregelen #coronacrisis door virus tijdens #coronamaatregelen coronatijd - #coronadebat coronawet doden vanwege #coronawet tegen crisis coronaregels demonstratie dam aantal vaccin regels rivm het testen #halsema @minpres halsema #covid19nl kabinet besmettingen nieuwe tweede van besmet verspreiding voor over @hugodejonge maanden meter ouderen patiënten golf zorg onderzoek 1,5 #covid_19 weken #coronavirusnl beroep afstand per lockdown houden #amsterdam aan alle 👉 jonge overleden aanpak @rivm op coronadoden uitbraak na china ic economie overheid coronabesmettingen ziekenhuis getest #rivm gevolgen klachten meer ondernemers ivm nu coronapatiënten test via tijden #anderhalvemeter we wet spoedwet nederland pandemie ⤵ hugo lees impact corona-crisis cijfers demonstreren femke juni corona-app coronabeleid ' minister #rutte positief coronaproof risico rutte ggd corona-maatregelen tijdelijke wereldwijd ondanks griep mensen verslaan miljoen getroffen #spoedwet klemmend verboden #viruswaanzin coronagolf samen corona-uitbraak zorgmedewerkers burgermeester controle handhaven juli #wilders burgers collectief corona-aanpak #alleensamen corona-regels gericht maakten artsen #lockdown coronatest keuzes #pvv en vliegtuigen horeca #demonstratie blijkt angst economische medicijn verpleeghuizen #blijfthuis burgemeester hoax zich volksgezondheid extra zware mondkapjes #covidー19 om volle boete ruim demonstraties gezondheid besmetting coronavaccin coronatijden duizenden mogelijk #coronanederland moeilijke #covid euro gebleken onder arts ziekenhuizen coronapandemie grondrechten personen mag laars zorgen beleid #samentegencorona umc \n",
      "202007 corona covid 19 #corona coronavirus #coronavirus coronacrisis #covid19 de - aantal coronamaatregelen besmettingen maatregelen #coronamaatregelen in virus rivm testen #coronacrisis door doden coronabesmettingen #covid19nl tijdens vaccin het golf we tweede 👉 patiënten @rivm coronatijd mensen getest besmet artsen positief griep test maanden crisis vanwege #covid_19 #rivm overleden \" over klachten afstand cijfers 1,5 verspreiding gevallen coronaregels tegen nieuwe coronagolf ic houden #covidー19 ivm lockdown meter coronawet coronadoden coronapatiënten #viruswaanzin stijgt mondkapjes @hugodejonge ziek meer controle via #covid__19 coronavaccin uitbraak #mondkapjes #coronawet onderzoek na pandemie stijging @viruswaanzin gevolgen regels @coronaquatsch coronatest #coronavirusnl ondanks angst medicijn landen corona-app aanpak #covid quarantaine onder aan #lockdown #hugodejonge besmetting belgië overheid ziekenhuis ggd ziekenhuizen economie impact who totaal ' wereldwijd #rutte dagen coronaproof @coronaconnectie kabinet zorg coronabeleid getroffen app opgenomen #covid19be tot @demasquecorona #spoedwet rutte voorkomen jonge gestorven toename arts #samentegencorona basisregels #alleensamen #covid-19 nu corona-maatregelen symptomen ouderenzorg ouderen hoax nederland #nieuwsuur huisarts lees per coronacijfers coronarichtlijnen maart brandbrief blijkt hydroxychloroquine europese hugo @burgerfront1 regering gates leidden #anderhalvemeter verpleeghuizen corona-uitbraak milde ✅ vermijd @coronavaccinnee besmette volksgezondheid neemt @nieuwsuur stierf tests tijden ‘ griepje 👇 wéér medische china meldde mogelijk onveiligheid #antwerpen zich economische horeca coronafonds 2e @minvws stijgen ettema bezoek corr behandeling afgelopen verpleegkundigen weken \n",
      "202008 corona covid #corona 19 coronavirus #coronadebat #covid19 #coronavirus coronacrisis #coronamaatregelen maatregelen - coronamaatregelen virus aantal de doden besmettingen #covid19nl testen griep test door mensen #coronacrisis coronabesmettingen vanwege coronaregels besmet @hugodejonge #covid_19 getest quarantaine kabinet ivm positief rivm tegen vaccin kamer tweede regels over coronatest coronatijd verspreiding klachten houden coronamelder ic ggd coronabeleid tijdens golf jonge we #persconferentie crisis totaal aan @rivm maanden #wilders cijfers #coronamelder patiënten artsen corona-app #rivm app motie hugo #grapperhaus coalitie overleden rutte besmetting #coronavirusnl grapperhaus ziekenhuis miljoen ventilatie \" minister in ' zich opgenomen worden pandemie #coronawet ziekenhuizen coronadoden via economie afstand coronapatiënten pcr coronaproof positieve ziek verplicht controle #spoedwet blijkt het #covid mondkapjes ⁦ coronawet overheid verplichte corr meer mrd scholen zorgpersoneel + overtreden ⁩ sijbesma hoax coronavaccin gevolgen #viruswaanzin dissel alle #coronaregels milde sluiten huisarts 1,5 coronasteun geen opnames uitbraak coronadebat onder #rutte zorg #alleensamen salaris aanpak hcq corona-maatregelen tests persconferentie berlijn #hugodejonge vaccins 👇 lockdown #lockdown coronagevallen zorgen vóór beleid #mondkapjes boetes ‼ rent coronabonus coronagolf #coronahoax coronagezant ziekte gevaar volksgezondheid oversterfte financiële #anderhalvemeter zorgmedewerkers per landelijke spoedwet coronabesmetting krijgen astrazeneca ➔ twentse horeca arts hydroxychloroquine besmettelijk ⤵ pcr-test ziekenhuisopnames sterven #covid__19 onderzoek bestrijding besmette stijging salarisverhoging 2e toename @minpres ouderen getroffen wet \n",
      "202009 corona covid #corona 19 coronavirus #covid19 #coronadebat #coronamaatregelen maatregelen #coronavirus coronacrisis coronamaatregelen aantal virus de besmettingen - coronaregels testen #covid19nl doden ic griep test getest patiënten #grapperhaus @hugodejonge coronabesmettingen mensen positief kabinet coronabeleid totaal coronatest #coronacrisis besmet #persconferentie rivm vaccin #fvd ziekenhuis grapperhaus rutte vanwege we jonge #covid_19 het maanden opgenomen klachten regels coronatijd ivm cijfers minister golf zorg in overleden coronaproof crisis verspreiding horeca coronavaccin @rivm coronawet positieve ziekenhuizen lockdown pcr-test @minpres scholen door ggd coronapatiënten hugo beleid #spoedwet coronadoden #mondkapjes #rivm maart aanpak #coronawet #coronaregels tweede besmetting tijdens overheid over #covid pcr ' opnames quarantaine corr #ikdoenietmeermee statistieken artsen mondkapjes #pcrgate ziekenhuisopnames nu per angst gevolgen #coronatest ziek tests krijgen coronatests #ikdoewelmee economie coronatesten hoax tegen houden zorgmedewerkers #alleensamen meer corona-maatregelen gevallen overtreden weken #covidー19 pandemie corona-app zich #rutte #coronavirusnl dissel strafblad onderzoek 1,5 huidig reguliere medewerkers @wybrenvanhaga bewezen uitslag alle coronagolf worden strengere coronacijfers lab sterfgevallen coronabesmetting symptomen coronaboetes coronadebat grappersmaus laten ouderen moeten arts overlijden coronaboete controle uitbraak @thierrybaudet terwijl persconferentie who spoedwet wegens voorkomen bruiloft broekers verdacht economische getroffen #tweedegolf april afstand snol nieuwe onder dhabi abu #covid__19 regio's huisarts % hcq kamer #coronavaccin risico stijging informatie testcapaciteit boete impact wereldwijd via https://t.co/csg3rreysl afgelopen #hugodejonge \n",
      "202010 corona covid #corona 19 #coronamaatregelen #covid19 coronavirus #coronavirus maatregelen coronamaatregelen #covid19nl - coronacrisis #coronadebat virus patiënten griep de ziekenhuizen zorg ziekenhuis positief aantal ic besmettingen + kabinet coronapatiënten mensen getest coronabesmettingen door besmet coronaregels #lockdown totaal #coronacrisis #spoedwet test doden testen coronatijd ivm #fvd positieve golf coronabeleid opgenomen #covid_19 overleden app #coronawet vanwege in % maanden rivm krijgen we #covid coronamelder tegen coronawet verspreiding lockdown cijfers #coronamelder coronatest #persconferentie crisis #pcrgate @hugodejonge controle tweede artsen vaccin reguliere klachten #covid__19 coronaproof rutte regels tijdens #coronaregels #alleensamen quarantaine coronagolf #mondkapjes extra nieuwe horeca corona-app ziek minder #rutte gisteren corr huidig #wilders alle #samentegencorona #rivm #horeca aanpak covid-patiënten coronadoden @minpres who duitsland 7d het coronacijfers okt opnames niet-westerse pcr-test \" onder coronavaccin zorgmedewerkers verspreiden #coronavirusnl meer @rivm pcr overlijden over besmetting statistieken strengere beleid worden angst ) nederland #coronapandemie informatie aan houden ggd gevolgen pandemie hugo trump ⤵ ( weken kamer onderzoek huisarts jonge hoax getroffen #covid19be nu uitgesteld helplijn afdeling gezin ouderen dagen personeel afgelopen griepje economie ifr coronabesmetting uitnodigen bedden wetenschappers sinds scholen #pvv coronagevallen gevallen corona-maatregelen wet coronatests sterfte verdacht oversterfte inwoners patiënt medicijn mortaliteit opgelopen patienten operaties hydroxychloroquine arts 1,5 hcq samen gehad uitslag snol ondanks ziekte https://t.co/csg3rreysl wegens \n",
      "202011 corona covid 19 #corona #coronamaatregelen #covid19 coronavirus - coronacrisis #covid19nl vaccin #coronavirus maatregelen coronamaatregelen aantal + virus besmettingen griep door #coronadebat coronavaccin coronatijd de testen mensen coronabesmettingen doden getest test positief patiënten vanwege kabinet in overleden @hugodejonge ic zorg positieve tegen coronaregels cijfers ziekenhuizen #coronacrisis vaccinatie we opgenomen besmet ziekenhuis tijdens coronabeleid #pcrgate sinds lockdown #vaccin 7d 󠁧 huidig coronapatiënten #coronavaccin crisis weken veilig pandemie coronaproof normaal rivm #covid_19 coronavaccins #covid % golf ivm gisteren besmetting vaccins who * artsen @rivm nieuwe ' nov totaal ziek #vaccinatie pcr-test #lockdown jonge minder #persconferentie okt pcr #mondkapjes #pcrtest gevolgen coronatest #alleensamen @minpres \" daling vaccineren #vaccinatieplicht maanden inwoners quarantaine #samentegencorona nertsen dagen onder geen per #mondkapjesplicht denemarken ) economie verspreiding werkzaam hugo regels 󠁢 󠁿 coronagolf overlijdens horeca 🏴 extra bemerk 100.000 coronacijfers ⤵ ( pfizer coronadoden ziekte statistieken ondanks mondkapjes aanpak mogelijk oversterfte getroffen ggd coronawet na lcps hoax opleveren genou opnames bocht #ditismijnzorg vaccinatieplicht 󠁥 󠁮 risico @edwinveldhuizen impact klachten gelden #rivm coronatesten over coronagevallen #coronavirusnl regering duurt coronapandemie arts ervoor corr cbd rutte stijging tweede #fvd via patiënt #spoedwet omgeving coronatijden bijwerkingen veiligheid wereldwijd overlijden schande #blijfthuis informatie meer duitse 🟣 feesten tijden sterven beleid personeel het aangenomen angst \n",
      "202012 corona covid 19 #coronamaatregelen #corona #covid19 - coronavirus #coronavirus coronacrisis #covid19nl maatregelen griep coronamaatregelen vaccin #coronadebat #lockdown kabinet virus #fvd de aantal #coronavaccin coronavaccin zorg patiënten mensen besmettingen cijfers ic doden overleden positief % door ziekenhuizen getest coronabeleid + test tegen ziekenhuis coronatijd rivm #coronacrisis lockdown ' vanwege artsen opgenomen tijdens vaccinatie #covid_19 coronabesmettingen het verspreiding coronaregels in testen #covid horeca #vaccinatie #vaccin besmet positieve variant pcr crisis vaccineren huidig pfizer who ziek #coronabeleid we vrijheid #persconferentie aanpak oversterfte bedden jonge coronatest coronaproof golf #rutte @hugodejonge ivm \" nederland totaal coronavaccins rutte coronacijfers zweden feiten mutatie sterven coronapatiënten #alleensamen ziekte vaccins #lockdownnl corona-maatregelen gevaarlijk over #mondkapjesplicht pandemie hoax personeel capaciteit klachten lockdowns coronadoden onacceptabel weken 17 scholen maanden @wybrenvanhaga geen ⤵ mrd reguliere * beleid inwoners patienten verspreidt update britse mondkapjes #codezwart #horeca coronatesten bijwerkingen extra @rivm #vaccinatieplicht ziekenhuisschip gevolgen zich overlijden angst wantrouwen ouderen regering nederlanders ivermectine na #pcrgate intensive beschermen maart informatie #avondklok opgelopen corona-beleid totale besmetting ernstig voorkomen #rivm nieuwe #mondkapjes volksgezondheid coronapandemie alle #blijfthuis januari terwijl influenza huisarts media contacten economie stageplek coronagevallen bron statistieken zieke motie afgelopen ernstige corr gestorven #hugodejonge wetenschappers covid-patiënten ondanks miljoen @minpres tests coronavariant hugo via #longcovid dood #coronahoax cbd dicht \n"
     ]
    }
   ],
   "source": [
    "for month in \"202002 202004 202005 202006 202007 202008 202009 202010 202011 202012\".split():\n",
    "    show_top_tokens(month, r'corona|covid', spy=False)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "20210131-23.out.gz\n",
      "^202101 corona covid 19 #corona #coronamaatregelen #coronavirus #covid19 #coronadebat coronavirus - coronamaatregelen coronacrisis vaccin maatregelen #avondklok #lockdown virus #covid19nl tegen de aantal % mensen #coronavaccin britse coronavaccin patiënten ic coronabeleid griep overleden variant positief doden door getest #fvd vaccinatie besmettingen testen ziekenhuis in #vaccinatie na coronatijd #coronaprotest zorg #covid ziekenhuizen #coronacrisis crisis #vaccin coronadoden besmet vanwege opgenomen over vaccineren ziek cijfers test coronapatiënten coronaregels vaccins rivm lockdown coronabesmettingen “ huisarts ggd coronavariant ” totaal #coronabeleid verspreiding coronavaccins sterven positieve het pandemie overlijden artsen we kabinet huidig ouderen ivermectine pfizer lijden besmetting \" coronatest @hugodejonge duitsland krijgen #covid_19 ivm tijdens gevaccineerd 2020 tot aanpak mondkapjes oversterfte minder #alleensamen #rutte onderliggend ⤵ bedden #stopdelockdown #avondklokdebat arts angst gevolgen coronaproof coronasystemen klachten dood jonge uitbraak ziekte coronacijfers maanden nederland golf liggen ' afdeling bijwerkingen coronavaccinatie lockdowns zuid-afrikaanse varianten beleid #persconferentie privégegevens infectie sterfte build @rivm patiënt reset zich gezonde #rivm ✔ zorgpersoneel controle corona-aanpak ondanks italië corona-maatregelen sterfgevallen @who hoax #vaccinatiestrategie vd covid-patiënten medische vaccinaties corona-uitbraak ernstige 💉 avondklok hugo 📽 coronapandemie sterft intensive meer griepje prik beschermen landen geen patienten scholen better pcr quarantaine miljoen care 70 who ‘ 23 economie onderliggende + overlijdt #vaccineren bestrijden coronawet kinderen immuunsysteem blijkt dissel nevenschade zorgen \n"
     ]
    }
   ],
   "source": [
    "show_top_tokens(\"^202101\", r'corona|covid', spy=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "20200331-23.out.gz\n",
      "^202003 corona #coronavirus coronavirus #corona #covid19 #coronavirusnederland #coronacrisis coronacrisis covid #coronavirusnl #coronanederland maatregelen #covid19nl 19 #coronadebat #covid_19 het - vanwege in virus nederland verspreiding patiënten #covidー19 #covid2019nl door tegen over overleden #rivm #covid2019 besmet aantal ziekenhuis coronapatiënten getest #coronahulp crisis aanpak besmettingen kabinet rivm italië #covid19nederland de #samentegencorona #covid19be doden #lockdown gevolgen mensen we strijd #coronanl testen scholen corona-virus om #rutte tijdens ivm zorg informatie coronamaatregelen ziekenhuizen uitbraak positief alle update ic #blijfthuis corona-crisis bestrijden van via #wilders rutte nieuwe #pvv #coronalied besmetting vragen griep totaal mogelijk klachten tijden ouderen #covid-19 bestrijding verband intensive voorkomen rondom afstand #coronavirusnetherlands #coronalul onze ziek corona-patiënten gevallen #hamsteren getroffen china te #socialdistancing ggd test patiënt gesloten huisarts beademing #coronavid19 houd #lockdownnl impact beschermen thuis #scholendicht sterven @rivm nederlandse corona-maatregelen collega wegens miljoen #nieuwsuur overleven #fvd huisartsen debat care minister coronadoden hulp roche #coronavirusupdate eu lees coronabesmettingen opgenomen ondernemers #covid extra jongeren nemen 1,5 koorts milde anesthesist #zorg verdere beperken zorgen bedden offer controle coronapatiënt miljard landen #coronapandemie evenementen verliezen symptomen onwaarschijnlijk farmaceut kroeg #groepsimmuniteit #samensterk zorgverleners #houdafstand sinds weken cijfers risico kwetsbare ⤵ oproep krijgen intensieve bevolking #stayathome blijf helpen inwoners dicht ➡ ontwikkelingen #eu persconferentie vastgesteld \" #socialdistance ( 450 italiaanse medewerkers afgelast #nederland \n"
     ]
    }
   ],
   "source": [
    "show_top_tokens(\"^202003\", r'corona|covid', spy=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "20210131-23.out.gz\n",
      "^20210131 ic - p op mensen ( % bedden totaal ) 100k treft / 🇳 🇱 + zkh \" ziekenhuis ic's covid 60 personeel artsen lockdown opgenomen liggen aantal 700 evenementen opnames etc jan 2400 c19 sportscholen p= 1500k 1000k leerl 600k 25k https://t.co/lvnm2bq4ya 10k 70k 4m 250k wv kermis theater student 90 kapper clubs ca ziekenhuizen https://t.co/sgrrmonmf9 horeca ): huidig winkel 27 zorg accepteren verpleegkundigen positief https://t.co/fbckeomuog patiënten doodgaan ernstig kind #stopdelockdown overlijden patiënt bed leer scholen ic-verpleegkundige belangrijke oproep besmettingen overlijdt daalt verplaatst repatriëring https://t.co/3uwyolioe8 42 ondersteunt pubers ondersteuning onvoldoende #covid19 6 vader 🇩 646 maatregelen fig 75 #rutte #covid19nl 72 🇪 142 vanaf update getest ic-opnames 1.573 omt risico ic-bedden 30 capaciteit 14 +3.700 13.998 978.475 vanuit óók #goedzo https://t.co/7rfppio2ca liegende yess #nosjournaal #npo1 overleden 20:00 vaccin #msm zender 3.714 @annstrikje corr perc griep 226.298 11,2 https://t.co/mvoqen2aoh reëel toegediend ic-artsen positieve ivm #coronaspreekuur lagen ● r 25 ic-verpleegkundigen ' verhalen corona haalt #buildbackbetterdebat momenteel 10 40 #ic leeg man inwoners vpk 19 1766 #vvd overgewicht nou strategie ligt er 2 😎 verdwenen beeld #scholenopen dec zh waarvan verpleegafdeling 4w pd risicofactoren 29 @spnl covid-unit signaalwaarde publicaties vh bijwerking werking heropening jaar mannen veiliger gezondheidswinst 56 ic-capaciteit patienten portugal \n"
     ]
    }
   ],
   "source": [
    "show_top_tokens(\"^20210131\", r'\\bic\\b', spy=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "20210131-23.out.gz\n",
      "^20210131 % long ziekte corona ( griep ) patiënten ; beschermen kanker diabetes 92 hartziekte obesitas nier gehospitaliseerde https://t.co/vdez5lc8k0 sciensano focus zwakke 32 bestaande 23 75 rapport leer 28 virus 11 groep 9 feiten rest normaal longcovid + leven shot beter aalter logistiek 58 @survivalitis ultieme papa job zoekt brugge vaste regio tip walgelijk iemand longen china had / omt voorkomen covid longontsteking al longread #longcovid menselijkheid ollongren faciliteren https://t.co/wu6eavrxus coronaschade misdaden medeplichtig herhalen doden risico toenemen hualong 1161 kernreactor https://t.co/q6gj83iy86 gebouwde commercieel line ontworpen mw klimaat @2ekamertweets ruimte vele one on tot waarin bouwen onderweg generatie , – vijf serie krijgt variant betekent ziekenhuizen zal kabinet ِ along to bedrijf britse zouden later eerste ruim genomen @longnetherlands advies neemt jaar overheid َ blijkt druk @kajsaollongren mijn volgende @long_goldsilver @bosveldklong was ♡ laten ganga aartis bijnor ballia 1,100 5-10 river 5x ziek soon sites long-keel #delaydiscounting https://t.co/dra9cq1z4b blaaskanker vruchtbaar beroerte @ntisec from hartaanval sigaret hart @gracelongollo longkanker besteed daarvan @elongated_musk @toilucknow at belangen zeker @longtallerne #dtv #omt @rroy9944 chronisch @openrotterdam blind omiya jongeren @cdavandaag komen mensen 50 miljard ass 3 kwetsbaar 4x إ ّ ل my @christenunie say isoleren vh kwetsbare verraadster collaborateur bilderberger ek nederlanden druppels gezondheidsrisico zieke risicogroepen #prinsesbeatrix \n"
     ]
    }
   ],
   "source": [
    "show_top_tokens(\"^20210131\", r'long', spy=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 22,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "20210131-23.out.gz\n",
      "^20210131 corona virus vaccin covid avondklok lockdown vaccins gevaccineerd vaccineren vaccinatie 19 besmettingen scholen ic variant testen maatregelen aantal de prikken ziekenhuis #avondklok vaccinaties pandemie mensen #coronamaatregelen getest positief na who % test #corona cijfers #lockdown rivm golf ziekenhuizen @hugodejonge coronamaatregelen pcr britse kinderen griep #coronaspreekuur #covid19 quarantaine het zorg doden ggd #coronavirus beu besmet afstand ziek @rivm prik besmetting week pfizer hcq mondkapjes mondkapje bewoners 1,5 artsen basisscholen patiënten 💉 @margabult overleden + virologen test-schandaal nu #covid19nl coronavirus #covid we hydroxychloroquine waren economie vanwege #scholendicht astrazeneca lockdowns blijkt #stopdelockdown 📽 arts kloppen overlijdens kregen dood opruiers worden tegen verspreiding dalen positieve ouderen klas coronabeleid leraren engelse beleid @who corona-uitbraak coronacrisis 30-1- varianten men omt #vaccinatie \" besmettelijk coronadoden minder reguliere 👊 besmetten er derde p #scholenopen crisis 2020 milaan alle winkels horeca @nos sneller @arieslob https://t.co/1qzacge4ga inenting kunnen #rutte medische overlijden houden groep dicht avondklokrellen spuit bestrijden openen personeel zwitserland tegenkomen maanden mogen medewerkers bedden besloten hugo virussen bevolking ziekte effect gaan viroloog meer ziekenhuisopnames krijgen wordt ziektes 🙄 anderhalf risico klachten eu umcg lansingerland #pcrtest totaal ontdooide . binnenkomst protesten bewaard jonge oversterfte uitbraak verplicht #stemzeweg #vaccinatiestrategie huisarts 🇩 adviezen leerkrachten aarhus weken ministerie kwetsbaren nexit #vaccin ontkomen in patiënt toe onderzoeken https://t.co/u9rerqs0vz 100k gemeld rond media sneltesten hebben nederlandse #wnlopzondag zult verpleegkundigen veilig ongevaccineerd dan binnen intensivisten @viruswaarheid 1 per logica etc #rivm geen terwijl gevaar s-avonds opnames effectief opgenomen https://t.co/imwujbreil afdeling treft nauwelijks vk sterven bootjes inentingen zelfquarantaine hun zkh intensive opgehokt uk onvoldoende bijwerkingen 🦠 blijven floreren betogen verpleeghuizen maatschappij berekening gevolgen dagen care #omt gezondheidswinst vervolging onomwonden huisartsen 8) laten staakt vanaf #avondklokrellen lak ( pleit bonden @rikkleinentink r geprikt ouders doses ingezet pleiten israel #ggd duizenden staken bescherming verpleegkundige gezette farmaceuten gezondheid symptomen ingeënt steeds omgekeerde 400 reizen juf gestraft voorkomen niet aangezien onderzoek begrip hoeveel @hommel_b begonnen @minpres #vaccineren @borgerpieter toestand enschede https://t.co/xxl2cj8pzl @mgdbertelli coronavaccins gates veiligheid gezet zich besmettelijker uitbraken @2ekamertweets liggen bureaucratische kabinet covid-patiënten 90 gooien pcr-test tvl-steun rijschoolhouder tennistrainer https://t.co/sfyrmmwf0o sars-cov verpleeghuis 14.000 aanpak tests mrna slechts reëel #coronaprotest toegediend heesch miljoen gevaccineerden c19 protest heelwijk ivermectine 75 @blckbxnews voorjaar 🇮 besmettingscijfers mutanten verlenging medicijnen lijden #vvd 60 waarvan onbewust wrang bron zweden pas vervolgd immuunsysteem kwetsbare gezonde medicijn 28 wél schade wuhan 3e indicatie coronatijd eumaffiabende gevolg #opendebat massaal kind 200k moeten gooit medicatie zitten onderliggende hoogste hoe hartstikke bezwaar héél waarin kantoor coronawandelingen 🇰 videochats chatgesprekken werken europa \n"
     ]
    }
   ],
   "source": [
    "show_top_tokens(\"^20210131\"   , r'corona|covid|huisarts|mondkapje|rivm|blijfthuis|flattenthecurve|houvol|virus|besmet|ziekenhui|\\bic\\b|'+\n",
    "                                r'intensive.care|^zorg|vaccin|[^ad]arts|uitbraak|pandemie|ggd|mondkapje|quarantaine|\\bwho\\b|avondklok|'+\n",
    "                                r'variant|verple|sympto|e.golf|pfizer|moderna|astrazeneca|astra|zeneca|novavax|biontech|\\btest|getest|'+\n",
    "                                r'sneltest|pcr|vaccin|ingeënt|ingeent|inent|prik|spuit|bijwerking|-->|💉|lock.down|lockdown|1[.,]5[ -]*m|'+\n",
    "                                r'afstand.*hou|hou.*afstand|anderhalve[ -]*meter|mutant|^omt$|hcq|hydroxychloroquine|umc|virolo|zkh|'+\n",
    "                                r'oversterfte|patiënt|patient|intensivist|🦠|uitbrak|ivermectin', \n",
    "                                spy=True, n=400)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "20210131-23.out.gz\n",
      "^20210131 vaccin pfizer astrazeneca bewaard bewoners koud #coronaspreekuur binnen uur prikken ontdooide binnenkomst toestand overlijdens anderhalf kregen gevaccineerd begonnen zorg 3 na werk worden wij met moeten miljoen ze .. zijn hoe vaccins moderna dit nu * eu kan levert @astrazeneca mijn doses overleden in @ema_news @gezondheidsraad @hugodejonge effectief viroloog / % 9 55 👇 toegediend astra en goedgekeurd 48 kwartier #astrazeneca 40 65 italië kwamen leden bevat effectiviteit zeneca the heen onder boechtvandenaldi vaccine coronavaccin mrna vertrouwen astrazeneca-vaccin getest resultaten bijwerking articles betrefd twijfels looking below plussers vaccinoorlog not it's see #pfizer gegevens zaterdag kwartaal artikelen blunder rij zuid-afrika good janssen https://t.co/bjdoqm6gza astrazenica 60 @demorgen dosis virus vaccinatie noemt all hadden pfizer-vaccin @pfizernl leert farmaceuten 66 bijsluiter mocht beschermt werkelijkheid injectie ziekenhuisopname for europese kortom minder ergste toe az lezen bedrijf europa tonen meldingen definition estimates https://t.co/garbdry6vm astra-zeneca ruimer sars-cov contract experimentele lauw covid ve beloofde adenovirus doorgaans ouderen duwt vrees #gentech #chimpansee #gmo #menselijke foetuscellen #hek293 gentechnologie #denkzelfna https://t.co/e6k6fcj3gm enthousiast gemodificeerd gekweekt #experimentele plaatsen afspraken géén extra 100 #ikprikhetnooit case genetisch bel omlaag @geertbourgeois die https://t.co/dfnbzujgex levering draadje italiaanse coronastrategie productie oa vk #coronavaccin keurt nogal #coronavirus duitsers groep openheid ceo betwisten coronavaccins bescherming biontech bierevenement génen nutshell pestgedrag \n"
     ]
    }
   ],
   "source": [
    "show_top_tokens(\"^20210131\", r'pfizer|moderna|astrazeneca|astra|zeneca|novavax|biontech', spy=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "20201231-16.out.gz\n",
      "^20201231-16 vaccin vaccinatie vaccineren vaccins indirecte vaccinatieplicht nl probleem liberale ziekenhuizen d66 gommers zover @ddstandaard partijen kuipers vvd zouden maandag \" verwacht gevaccineerd pfizer week directe liggen 200.000 vrijwillige https://t.co/zl1fej1o3u tegen #fvd zit forum democratie jonge qua koelkast eventbranche wereldtop #coronavaccin organisatie vrijwel enkel ziekenhuispersoneel vaccinaties terwijl @joostniemoller prik moeite @thierrybaudet gaan slechter @hugodejonge #vaccinatiestrategie @wybrenvanhaga registratie leg trek lijkt beginnen niet ervaring sector bijwerkingen zich events crowd deadlines #takeaction oproep ticket werk communicatie management uitnodigen verbinding nederland hugo covid slecht #vaccinatie sterft zonder aanstaande ingeënt vrijwilli neer niets virus zijn it begrijp personeel 26 zullen bevolking coronavaccin het tunnel getest is doden ‼ te @shossontwits kabinet steeds luister artsen ingewikkeld % vrijwillig spuiten @irenriquecortes eve #covid19 symbolisch #vaccin gevaccineerde --> @benwoldring landen 💉 https://t.co/evfn2ys4do inenten én met @isayael2020 - corona sterven israël verpleeghuispersoneel levensreddend eerst gezin #coronavirus zichzelf @rosenthal_uri zorgpersoneel 24 voorrang duitsland #vaccincovid #vaccineren prikken zorginfarct termijn willen onderliggend veerkracht collega aantal 🇩 weigeren @nos @vogelvrijeharts ggd naast massaal https://t.co/obgyp1atyl 99,935 https://t.co/2wodxacliz verpleeghuizen licht lijden toegang weiger dinsdag zwaar @nieuwsuur dreigt #corona mouwen opstropen gevaccineerden pfizer-vaccin 37,273 78.109 🇪 advies bekend ondertussen #covid hen kennis bevestigen plan basis tegemoet wie genoeg leggen koortsachtige dadeloosheid spuit spuitje vaccineer vaccine vaccinatieprogramma \n"
     ]
    }
   ],
   "source": [
    "show_top_tokens(\"^20201231-16\", r'vaccin|pfizer|ingeënt|ingeent|inent|prik|spuit|bijwerking', spy=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 23,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "20210124-23.out.gz\n",
      "^20210124 #eindhoven eindhoven in #amsterdam station amsterdam politie geplunderd #rellen rellen winkels ondernemers geplaagde riem https://t.co/bko9hkxpzr naam demonstranten #urk vrijheid hart beelden geweld uit steken waterkanon tuig dit jumbo urk plunderen geslagen alle demonstraties worden zojuist relschoppers onder en #tuig #coronaprotest plunderingen #museumplein #schilderswijk vrouw #demonstratie slopen protest #politie werkelijk @gem_eindhoven herkent vernielingen #avondklokprotest tientallen vreedzame demonstratie traangas vandaag binnenstad #tilburg ontgelden treinverkeer wie nu door gekste rt paard vanwege politiepaarden @brabostudent #venlo onderuit haag optreden hard https://t.co/wkp71cnjjs protesteren respect politiepaard afgesloten confrontatie https://t.co/u9bih0l5vf probeerden kappen https://t.co/eokknu0xw6 auto gestaakt #musemplein demonstreren tegen dame vuurwerk bebloede naïef misbruik #doeermeewatjewil https://t.co/dfnuke4fxp stad https://t.co/lhndn5tzxh treinstation leest den collega buitenlandse geschrokken nadat https://t.co/avi8a7fk3r agent vierluik paarden https://t.co/inq27yda1l zwakzinnige emoties noemen 💙 eén #denhaag schandalige eindhovense @hartvnl hoofd schreef golfballen gezicht met zie #enschede gesloten @samhvnl tilburg onrust #coronamaatregelen jezus actie raad centrum livestream oké ongehoord @cwbbakker #stein https://t.co/gtkhdgayui venlo #helmond jongen messen @jjorritsma040 vreedzaam maken @omroepbrabant wapens mathildelaan heet zooitje #video grijpt muur stein ingezet jongeren 18 gif'jes https://t.co/yaa6lwnlfd losgebroken https://t.co/bwc0jcohr9 vreselijk trieste vanuit bedoel centraal #avondklok https://t.co/nmidfvdxcb beëindig onze brand liefdevolle @quintenvdmolen septemberplein #pleurop gewond erachter 💔 het lof stenen rellend wist s #iksteundepolitie moest sterkte @gelehesjesnl begonnen horecaondernemer terreur loopt daarom chaos hulp \n"
     ]
    }
   ],
   "source": [
    "show_top_tokens(\"^20210124\", r'eindhoven', spy=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 21,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "20210124-23.out.gz\n",
      "^20210124 #eindhoven eindhoven in #amsterdam station amsterdam politie geplunderd #rellen rellen winkels ondernemers geplaagde riem https://t.co/bko9hkxpzr naam demonstranten #urk vrijheid hart beelden geweld uit steken waterkanon tuig dit jumbo urk plunderen geslagen alle demonstraties worden zojuist relschoppers onder en #tuig #coronaprotest plunderingen #museumplein #schilderswijk vrouw #demonstratie slopen protest #politie werkelijk @gem_eindhoven herkent vernielingen #avondklokprotest tientallen vreedzame demonstratie traangas vandaag binnenstad #tilburg ontgelden treinverkeer wie nu door gekste rt paard vanwege politiepaarden @brabostudent #venlo onderuit haag optreden hard https://t.co/wkp71cnjjs protesteren respect politiepaard afgesloten confrontatie https://t.co/u9bih0l5vf probeerden kappen https://t.co/eokknu0xw6 auto gestaakt #musemplein demonstreren tegen dame vuurwerk bebloede naïef misbruik #doeermeewatjewil https://t.co/dfnuke4fxp stad https://t.co/lhndn5tzxh treinstation leest den collega buitenlandse geschrokken nadat https://t.co/avi8a7fk3r agent vierluik paarden https://t.co/inq27yda1l zwakzinnige emoties noemen 💙 eén #denhaag schandalige eindhovense @hartvnl hoofd schreef golfballen gezicht met zie #enschede gesloten @samhvnl tilburg onrust #coronamaatregelen jezus actie raad centrum livestream oké ongehoord @cwbbakker #stein https://t.co/gtkhdgayui venlo #helmond jongen messen @jjorritsma040 vreedzaam maken @omroepbrabant wapens mathildelaan heet zooitje #video grijpt muur stein ingezet jongeren 18 gif'jes https://t.co/yaa6lwnlfd losgebroken https://t.co/bwc0jcohr9 vreselijk trieste vanuit bedoel centraal #avondklok https://t.co/nmidfvdxcb beëindig onze brand liefdevolle @quintenvdmolen septemberplein #pleurop gewond erachter 💔 het lof stenen rellend wist s #iksteundepolitie moest sterkte @gelehesjesnl begonnen horecaondernemer terreur loopt daarom chaos hulp \n"
     ]
    }
   ],
   "source": [
    "show_top_tokens(\"^20210124\", r'eindhoven', spy=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "show_top_tokens(\"^2021012[34]\", r'avondklok', spy=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "20201231-23.out.gz\n",
      "^(20200[7-9]|20201) corona covid #corona 19 #coronamaatregelen #covid19 coronavirus de #coronavirus - coronacrisis maatregelen coronamaatregelen #coronadebat virus aantal #covid19nl besmettingen griep in mensen door doden het vaccin testen patiënten kabinet positief getest test ic coronabesmettingen we #coronacrisis zorg coronaregels rivm besmet ziekenhuis vanwege ziekenhuizen coronatijd tegen @hugodejonge overleden cijfers coronabeleid #covid_19 tijdens + golf coronavaccin #lockdown totaal over ivm lockdown opgenomen maanden crisis coronapatiënten ' coronatest artsen verspreiding positieve \" #fvd klachten tweede % #persconferentie rutte regels @rivm . aan jonge #covid meer quarantaine coronaproof ziek horeca #spoedwet #rivm nieuwe houden coronawet coronadoden besmetting pandemie #coronawet aanpak ggd #mondkapjes pcr #coronavaccin geen krijgen hugo who nu #alleensamen @minpres gevolgen weken #rutte alle beleid coronagolf mondkapjes onder met controle economie corona-app via worden #coronavirusnl coronamelder overheid angst vaccinatie corr app zich scholen en hoax huidig #pcrgate pcr-test #grapperhaus afstand minister per coronacijfers 1,5 #vaccin na kamer extra ) corona-maatregelen gevallen opnames #coronamelder onderzoek statistieken ( nederland #samentegencorona huisarts maart ziekte #coronaregels #covid__19 ouderen uitbraak vaccins minder reguliere oversterfte tests sinds arts miljoen ondanks overlijden ⤵ zorgmedewerkers getroffen stijging #covidー19 sterven grapperhaus afgelopen van gisteren coronatesten personeel impact ziekenhuisopnames #mondkapjesplicht bedden voorkomen dagen informatie coronagevallen wereldwijd coronatests blijkt zweden vaccineren risico #hugodejonge #horeca symptomen \n"
     ]
    }
   ],
   "source": [
    "show_top_tokens(\"^(20200[7-9]|20201)\", r'corona|covid', spy=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 25,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "20201231-23.out.gz\n",
      "^(2020) corona #corona covid coronavirus #coronavirus 19 #covid19 coronacrisis maatregelen #coronacrisis #coronamaatregelen virus - #coronadebat #covid19nl de coronamaatregelen aantal in door patiënten doden #coronavirusnederland besmettingen #coronavirusnl testen besmet crisis vanwege #covid_19 rivm mensen getest griep #coronanederland het kabinet ziekenhuis verspreiding tijdens ic overleden we zorg positief tegen vaccin coronapatiënten coronatijd #rivm ziekenhuizen #covidー19 test over @rivm coronabesmettingen #lockdown nederland coronaregels aanpak rutte cijfers ivm china gevolgen #samentegencorona lockdown weken italië klachten afstand besmetting coronabeleid tijden opgenomen totaal artsen uitbraak ziek via coronadoden #rutte scholen ouderen @hugodejonge coronavaccin #persconferentie #covid alle pandemie golf #covid2019 #covid19be #blijfthuis #covid2019nl economie ggd mondkapjes aan #alleensamen coronatest thuis getroffen quarantaine huisarts controle nieuwe 1,5 houden #mondkapjes nu mogelijk corona-crisis informatie + @minpres impact corona-app #fvd positieve maanden regels overheid coronaproof maart coronawet minister #coronanl tweede who miljoen strijd sterven gevallen krijgen patiënt #coronawet corona-maatregelen extra zorgmedewerkers voorkomen horeca risico april landen \" ondernemers corona-virus symptomen meter jonge bestrijding arts intensive onder #spoedwet #wilders bestrijden bedden personeel zich #coronahulp inwoners coronagolf lees ondanks wereldwijd ziekte overlijden economische gezondheid ⤵ onderzoek coronatijden angst medische medicijn verband update #covid-19 zorgen verpleeghuizen sterfgevallen % worden onze volksgezondheid hoax medewerkers persconferentie #coronavaccin gesloten tests samen per coronacijfers #covid19nederland na #nieuwsuur vragen \n"
     ]
    }
   ],
   "source": [
    "show_top_tokens(\"^(2020)\", r'corona|covid', spy=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "show_top_tokens(\"^202001\", r'corona|covid', spy=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "show_top_tokens(\"^202002\", r'corona|covid', spy=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "show_top_tokens(\"^202003\", r'corona|covid', spy=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "show_top_tokens(\"^202004\", r'corona|covid', spy=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "show_top_tokens(\"^202005\", r'corona|covid', spy=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "show_top_tokens(\"^202006\", r'corona|covid', spy=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "show_top_tokens(\"^202007\", r'corona|covid', spy=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "show_top_tokens(\"^202008\", r'corona|covid', spy=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "show_top_tokens(\"^202009\", r'corona|covid', spy=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "show_top_tokens(\"^202010\", r'corona|covid', spy=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "show_top_tokens(\"^202011\", r'corona|covid', spy=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "show_top_tokens(\"^202012\", r'corona|covid', spy=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "show_top_tokens(\"20200831\", r'\\b(je|jij|jou|jouw|jezelf|u|uw|ge|gij|aub)\\b', spy=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "show_top_tokens(\"20200630\", r'\\btest|getest|sneltest|pcr', spy=True,\n",
    "                query2=(\"verkoud|snot|ziek|sympto|griep|koorts|verhoging|pijn|hoest|nie[sz]|klacht|loopneus|benauwd|vermoeid\"+\n",
    "                        \"|(verlies|verlo).*(reuk|smaak)|(reuk|smaak).*(verlies|verlo)\"))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "show_top_tokens(\"20210114\", r'avondklok', spy=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "show_top_tokens(\"20200831\", r'\\b(geluk|gelukkig|gelukkige|blij|happy)\\b', spy=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "show_top_tokens(\"20201031\", r'eenza|alleen|isolatie|isolee?r', spy=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "show_top_tokens(\"20200831\", r'\\b(ik|mij|mijn|me|mn|m\\'n|zelf|mezelf|mijzelf|i)\\b', spy=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "show_top_tokens(\"20200831\", r'\\btest|getest|sneltest|pcr', spy=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "show_top_tokens(\"20201130\",\n",
    "    (r'shit|fuck|klote|kanker|kut|lul|dom|reet|gvd|strot|waanzin|nep|rotzooi|idio|schaamtelo|fake|strot|'+\n",
    "     r'zeik|stom|onbeschoft|oprot|bekrompen|ongelo|bah|kwalijk|onbekwa|achterlijk|jat|schand|triest|puinho|'+\n",
    "     r'immore|kleptocratie|nepotism|absurd|misselijkma|\\bbek\\b|\\brot\\b|\\bkk\\b'),\n",
    "    spy=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "MONTH = \"202012\"\n",
    "\n",
    "for day in range(1,32):\n",
    "    show_top_tokens(MONTH+str(day).zfill(2), QUERY_TOPIC_SMALL)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "top_tokens = {}\n",
    "infile = open(\"202012-tokens.txt\",\"r\")\n",
    "for line in infile:\n",
    "    tokens = line.strip().split()\n",
    "    date = tokens.pop(0)\n",
    "    top_tokens[date] = tokens\n",
    "infile.close()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "top_tokens_freq = {}\n",
    "for date in top_tokens:\n",
    "    top_tokens_freq[int(date)] = {}\n",
    "    for i in range(0,len(top_tokens[date])):\n",
    "        token = top_tokens[date][i]\n",
    "        if i < 10: top_tokens_freq[int(date)][token] = 3\n",
    "        elif i < 100: top_tokens_freq[int(date)][token] = 2\n",
    "        else: top_tokens_freq[int(date)][token] = 1"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "top_tokens_freq7 = {}\n",
    "for date in top_tokens_freq:\n",
    "     for d in range(date-3,date+4):\n",
    "            if d in top_tokens_freq:\n",
    "                for token in top_tokens_freq[date]:\n",
    "                    if not token in top_tokens_freq7: top_tokens_freq7[token] = {d:top_tokens_freq[date][token]}\n",
    "                    elif d not in top_tokens_freq7[token]: top_tokens_freq7[token][d] = top_tokens_freq[date][token]\n",
    "                    else: top_tokens_freq7[token][d] += top_tokens_freq[date][token]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "interesting_tokens = [\"virus\",\"griep\",\"vaccin\",\"overleden\",\"besmettingen\",\"doden\",\"kabinet\",\"ic\",\"test\",\"pfizer\",\"lockdown\",\"variant\"]\n",
    "\n",
    "plt.subplots(figsize=(24,6))\n",
    "for token in top_tokens_freq7:\n",
    "    if sum(top_tokens_freq7[token].values()) > 100 and token in interesting_tokens:\n",
    "        plt.plot_date([datetime.datetime.strptime(str(d),\"%Y%m%d\") for d in top_tokens_freq7[token].keys()], \n",
    "                      [x-0.4*random.random()+0.1 for x in top_tokens_freq7[token].values()], fmt=\"-\", label=token)\n",
    "plt.legend(loc=\"upper left\")\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Month vocabulary analysis"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 26,
   "metadata": {},
   "outputs": [],
   "source": [
    "month_data = {}\n",
    "data_file = open(\"month-topics.txt\", \"r\")\n",
    "for line in data_file:\n",
    "    fields = line.strip().split()\n",
    "    month = fields.pop(0)\n",
    "    month_data[month] = fields\n",
    "data_file.close()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 27,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "^202001 virus china corona wuhan chinese mensen uitbraak doden nederland besmet besmettingen frankrijk rovirus rvirus gevallen europa maatregelen patiënten duitsland eerste \n",
      "^202002 virus corona china besmet italië nederland mensen  quarantaine rivm griep vragen besmetting chinese besmettingen covid doden patiënt uitbraak paniek \n",
      "^202003 corona virus nederland crisis covid maatregelen debat  verspreiding patiënten overleden aantal rivm besmet ziekenhuis getest mensen hulp besmettingen kabinet \n",
      "^202004 corona virus crisis covid  debat nederland maatregelen patiënten aantal overleden doden samentegencorona getest ziekenhuis tijd gevolgen besmet ic china \n",
      "^202005 corona covid crisis virus  maatregelen doden tijd debat aantal patiënten testen kabinet nederland besmet rivm onderzoek besmettingen getest vaccin \n",
      "^202006 corona covid crisis virus maatregelen  tijd doden debat wet demonstratie dam aantal regels vaccin rivm halsema minpres testen nieuwe \n",
      "^202007 corona covid virus crisis  aantal besmettingen maatregelen rivm testen doden mensen vaccin tweede golf 👉 patiënten tijd getest besmet \n",
      "^202008 corona covid virus debat crisis maatregelen  aantal doden besmettingen mensen testen griep test hugodejonge besmet regels getest kabinet quarantaine \n",
      "^202009 corona covid virus debat maatregelen crisis aantal  besmettingen testen regels doden mensen ic griep test getest hugodejonge grapperhaus patiënten \n",
      "^202010 corona covid maatregelen virus  crisis debat patiënten griep zorg ziekenhuizen ziekenhuis aantal positief mensen ic besmettingen kabinet getest besmet \n",
      "^202011 corona covid maatregelen  virus crisis vaccin aantal besmettingen griep mensen debat tijd testen doden getest test positief patiënten kabinet \n",
      "^202012 corona covid maatregelen  virus crisis griep vaccin kabinet debat lockdown fvd mensen aantal zorg patiënten besmettingen cijfers ic doden \n",
      "^(20200[7-9]|20201) corona covid maatregelen virus  crisis debat aantal besmettingen griep mensen doden vaccin testen patiënten kabinet positief getest test ic \n",
      "^202101 corona covid maatregelen virus debat  crisis vaccin avondklok lockdown aantal britse mensen variant patiënten ic beleid griep overleden positief \n",
      "^(2020) corona covid virus crisis maatregelen  debat aantal patiënten doden nederland besmettingen testen besmet rivm mensen getest griep kabinet ziekenhuis \n"
     ]
    }
   ],
   "source": [
    "MAX_WORDS = 20\n",
    "STOP_WORDS = [ \"de\", \"het\", \"in\", \"over\", \"via\", \"tijdens\", \"tegen\", \"vanwege\", \"door\", \"voor\", \"op\", \"van\", \"aan\", \"om\", \"en\", \"meer\", \"per\", \"na\", \"sinds\", \"tot\",\n",
    "               \"we\", \"alle\", \"zich\", \"huidig\", \"extra\", \"meer\", \"minder\", \"ivm\", \"geen\", \"virusnl\", \"virusnederland\", \n",
    "               \":\", '-', '+', '%', '\"', \"'\", \".\", \"(\", \")\",'19', '19nl', '_19', 'ー19', '2019nl', '2019', '__19', '19be', 'sen' ]\n",
    "\n",
    "for month in month_data:\n",
    "    seen = {}\n",
    "    counter = 0\n",
    "    print(month, end=\" \")\n",
    "    for i in range(0,len(month_data[month])):\n",
    "        word = month_data[month][i]\n",
    "        word = re.sub(\"^#\", \"\", word)\n",
    "        word = re.sub(\"^@\", \"\", word)\n",
    "        word = re.sub(r'^corona(.)', r'\\1', word)\n",
    "        word = re.sub(r'^covid(.)', r'\\1', word)\n",
    "        word = re.sub(r'^virus(.)', r'\\1', word)\n",
    "        word = re.sub(r'^nl$', 'nederland', word)\n",
    "        word = re.sub(\"^-\", \"\", word)\n",
    "        if word not in seen and word not in STOP_WORDS:\n",
    "            print(word, end=\" \")\n",
    "            seen[word] = True\n",
    "            counter += 1\n",
    "            if counter >= MAX_WORDS:\n",
    "                break\n",
    "    print()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Analysis frequencies \"de\" and \"ik\" in corona|covid tweets"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 28,
   "metadata": {},
   "outputs": [],
   "source": [
    "def count_words(file_pattern, query1, query2, spy=True):\n",
    "    file_list = sorted(os.listdir(DATADIR))\n",
    "    nbr_of_matches1 = 0\n",
    "    nbr_of_matches2 = 0\n",
    "    for file_name in file_list:\n",
    "        if re.search(file_pattern,file_name):\n",
    "            if spy: squeal(file_name+\" \"+file_pattern+\" \"+query2)\n",
    "            df = pd.read_csv(DATADIR+file_name, index_col=ID)\n",
    "            matching_rows1 = df[TEXT].str.contains(query1, flags=re.IGNORECASE)\n",
    "            matching_rows2 = df[TEXT].str.contains(query2, flags=re.IGNORECASE)\n",
    "            try:\n",
    "                nbr_of_matches1 += matching_rows1.value_counts()[True]\n",
    "            except:\n",
    "                pass\n",
    "            try:\n",
    "                nbr_of_matches2 += (matching_rows1 & matching_rows2).value_counts()[True]\n",
    "            except:\n",
    "                pass\n",
    "    return(nbr_of_matches1, nbr_of_matches2)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 29,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "20201231-23.out.gz 20201231 \\bik\\b\n"
     ]
    }
   ],
   "source": [
    "base_pattern = \"corona|covid\"\n",
    "\n",
    "results = {\"de\":{}, \"ik\":{}}\n",
    "for month in range(1, 13):\n",
    "    for day in range(1,32):\n",
    "        file_pattern = \"2020\"+str(month).zfill(2)+str(day).zfill(2)\n",
    "        results[\"de\"][file_pattern] = count_words(file_pattern, base_pattern, r'\\bde\\b')\n",
    "        results[\"ik\"][file_pattern] = count_words(file_pattern, base_pattern, r'\\bik\\b')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 30,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXoAAAEICAYAAABRSj9aAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/Il7ecAAAACXBIWXMAAAsTAAALEwEAmpwYAABplElEQVR4nO2dd3gc1dWH36PeiyVZsiXLveBuYxswzYRmemihJBAghBAgJF8aNQmBEEgIhARIHAIJkNB7Nw6YYjCuuOBuucuyLVm9S6u93x93Znd2tSutZLWV7/s8elY7c2fmzs7MmXN/99xzRSmFwWAwGPovEb1dAYPBYDB0L8bQGwwGQz/HGHqDwWDo5xhDbzAYDP0cY+gNBoOhn2MMvcFgMPRzjKE/zBGRH4rIARGpEZGMDm77lIj8rrvq1hlE5CoR+TzEsutFZI71/10i8t9OHO/bIrKgo9uFCyIyR0Q+6eZjvC8i3w2ybpiIKBGJ6s469HeMoe8B+qJBBBCRaOAh4DSlVJJSqtRv/TAR2dlFx9opIsO6Yl+HUAef81FKTVBKfRLCdneJyF2B1imlnlVKndZllTwMUUqdoZR6ujPbWvfVKV1dp3aO2SmnoDcxht7iMPUYsoE4YH1vV8TQOURjnmNDm4T9DWK90W8TkQ0iUi4i/xaROMf6s0VktYhUiMhiEZnst+0tIrIWqBWRKBE5zipXISJ7ROQqq2ysiPxJRHZbUsc8EYm31s0RkUIR+ZmIFIvIPhG52lp3HfBt4JeWPPK2tfxWEdkmItVW3c931CtSRB4UkYMiskNEbnI2X0UkVUSetI6zV0R+JyKRQX6fWBF5WESKrL+HrWVjgM1WsQoRWRjCbz1NRL6y6vwi+iXhXB/0t25nv1eLyEZrv9tF5AeOdUF/W2t9hoi8JSJVIrIMGBnKMa1tA3qDIhItIs+LyKsiEtPOPnykIus6XS8iW6378TERkSDbRorI7Y77YKWIDLHWzRaR5SJSaX3Odmz3iYjcKyJfAHXAiBDK3yMiX1jHWSAimY71L4vIfmvbz0RkQhvnG2/dm7us8p87noNzRcthFdYxj7CW3yoir/jt5y8i8ldH/a51/CZ/su797cBZbdTlP0A+8Lb1bP1SRJ4WkZ9Z63Ot63GD9X2UiJTZ16Ot+1VEBlvXv0T0M3iztXwucDtwiXXMNdbyq6x7t9oq/+1g9e4VlFJh/QfsBNYBQ4ABwBfA76x104Fi4CggEviuVT7Wse1qa9t49E1TDVwGRAMZwFSr7MPAW9YxkoG3gfusdXMAF3C3td2Z6Acw3Vr/lF0nR70vBgajX7aXALXAIGvd9cAGIA9IBz4EFBBlrX8D+AeQCAwElgE/CPL73A0sscplAYuBe6x1w5z7bed3jgF2Af9nneNFQHOov3U7+z4LbaAFONH67aaH+Nu+ALxk/RYTgb3A5x24d06x/r8L+K91H7xrXbPIEPZxlfN41u/5DpBm3U8lwNwg2/4C+BoYa537FPQ9NwAoB64AotD3YzmQYW33CbAbmGCtzw6h/DZgjHV+nwD3O+pxDfqejkXf56sd6+YAnzi+P2Ztn2td59nWdmPQ9/Cp1nX6JVBg3TdDrWuWYu0jEtgHHO2o37WOe38T3uf5Y9q4R53X0HEub1v/X26d94uOdW+2d7+in8mVwK+t+o8AtgOnO+8VxzETgSpgrPV9EDCht22jz+/U2xU45BPQF+d6x/czgW3W/3/HMmqO9ZuBEx3bXuNYdxvweoBjiHUTj3QsOwbY4XgY6p03o3UT2TfyU/gZ+gDHWA2cZ/2/EIfhBk6xb3b0Q90IxDvWXwZ8HGS/24AzHd9PB3Za/w9r6yHy288JQBEgjmWL8Rr6Nn/rDl7TN4Aft/fbWg9oMzDOse73dN7QvwV8CvzVeZ7t7OMqWhv64xzfXwJuDbLtZvua+y2/Aljmt+xL4Crr/0+AuztY/k7HuhuA+UHqlGadQ6rj9//E+j/CuhZTAmz3K+Alx/cI9Et3jvX9c+BK6/9TsZ5RR/1sQ78Q3+f5tLbuUVob+pFAhXX8ecAPgEJr3dPAT9u7X9HGf7ffutuAfzvuFX9DXwFciOO57Et/YS/dWOxx/L8L7SmD9iR+ZjXNKkSkAu0pDA6y7RC0YfQnC0gAVjr2M99ablOqlHI5vtcBScEqLCJXOpqNFWhv1G5OD/arl/P/oWiPaZ9j23+gPfZADEb/JjbO36cjDAb2KuvOduzLWa/2fuuAiMgZIrLEalZXoF/WmY4iwX7bLPTLz//6d5ajgclob/dQsv3td/zf1n0Q7H7zv2ZY33Md353nHEr5gHWypJL7LfmoCm04wff3x7EsLpQ6K6XcVh3tOjyHdkhAe9rPBdiHvZ9OX0+l1DagBpgKHI9uXRWJyFi0Ef/UKtrW/ToUGOy37na0kxXomLXoVvn16OfyXREZ15F6dzf9xdAPcfyfj/Y8Qd8w9yql0hx/CUqp5x3lnQ/0HgJrvAfRnswEx35SlVJBDbkfPkZDRIYC/wRuQjev09Dyk63l7kPLNoHObw/ao8901CVFKRVMVy1C37g2zt+nI+wDcv305ny/erX3W7dCRGKBV4E/AdnWb/Ee3t+iLUrQso7/9e8sC4D7gI9EJOBD3cUEu9/8rxno89rr+O68p0IpH4zLgfPQrcZUdCsPAv/+B4GGUOps3SdDHHV4GZgjInnA+QQ39Pvo2PUM9EL+FC0txiil9lrfr0TLoKutMm3dr3vQrXXnumSl1JnBjqmU+kApdSpattmEfr77DP3F0N8oInkiMgD95n3RWv5P4HoROUo0iSJylogkB9nPs8ApIvIt0R2zGSIy1fJO/gn8WUQGgqej5/QQ63cArfPZJKJvlhJrX1ejPXqbl4AfW8dIA26xVyil9qEN0oMikiIiESIyUkRODHLs54E7RSTL6oD7NVqL7ihfoo3qzdZvcwEwy7G+o7+1TQxaFy0BXCJyBrq53i5KqRbgNeAuEUkQkfForbXTKKX+iDZCH4mjw7KbeAK4R0RGW7/ZZNFjGd4DxojI5dZvfQkwHu2dBqKj5Z0kox2HUnSr9ffBClrPwb+Ah6zOykgROcZ6Wb8EnCUiJ4sO2/2Ztd/F1rYlaInm32gjujHIYV5C32N5IpIO3NpO/f2fLdCG/SbgM+v7J8CP0BJbi7Wsrft1GVAlOlAj3jrPiSIy03HMYWJFO4lItuiO6ETrnGsA+zh9gv5i6J9DG7/t1t/vAJRSK4DvA4+iO6cK0JpqQJRSu9Gywc+AMvTbf4q1+hZr+yVWE/dDdCdaKDwJjLeagW8opTYAD6KN5wFgEroT2eaf1vmsBVahH2QX3pvnSrSB3GCd1ytoTyIQvwNWWPv6GvjKWtYhlFJNwAXo368c3VR9zbG+Q7+1Y7tq4Gb0A16O9jDf6kDVbkLLEPvRfSH/7sC2wep0D7qf4EPLeeguHkKf9wJ0Z96TaI23FDgbfR+Wojs2z1ZKHQxS3w6V9+MZtDyyF30/LWmn/M/R99Fy9DPyByBCKbUZ+A7wCNrzPwc4x7pvbJ5DtxyCefOg7/0PgDXoe/W1NsqCboHdaT1bP7eWfYp+gdmG/nP0S8z+3ub9ar0MzkHLPzus83kC3eIB3ToBKBWRr9B29GfoVk0ZWiK6oZ169yhyaFJk7yN6AMy1SqkPe7su3YXl5c5TSvk3zw2GbkX0yOG7lFJzercmhkOhv3j0/QqruXim1QzPBX4DvN7b9TIYDOGJMfR9EwF+i25SrgI2orV1g6Gn2YmWxAxhTNhLNwaDwWBoG+PRGwwGQz+n1xJ5ZWZmqmHDhvXW4Q0GgyEsWbly5UGlVFb7Jb30mqEfNmwYK1as6K3DGwwGQ1giIh0e/W2kG4PBYOjnGENvMBgM/Rxj6A0Gg6GfczjOqmQwGA4TmpubKSwspKGhober0mHi4uLIy8sjOjr6kPdlDL3BYOi3FBYWkpyczLBhw5DAE331SZRSlJaWUlhYyPDhww95fyFJNyIyV0Q2i0iBiLTKJid6aru3RWSN6KnErg60H4PBYOhJGhoayMjICCsjDyAiZGRkdFlLpF1DL3ou0seAM9CpTy+z0sE6uRHYoJSagp6R5kFpZ65Ng8Fg6AnCzcjbdGW9Q/HoZwEFSqntVsrRF9ATFThRQLI12UASOlWnix5AKcWrKwtpaO5T6Z8NBoOhzxCKoc/Fd2qvQnynKAOd0/kIdD7mr9Hzfbr9dyQi14nIChFZUVJS0skq+/LZ1oP87OU13P/+pi7Zn8FgMHQXd911F3/60596/LihGPpA7Qf/TGinoyfpGIxO1v+oiKS02kipx5VSM5RSM7KyOjSCNyhV9c0AlFQ3dsn+DAaDob8RiqEvxHcOxzxazzl6NfCa0hSgZ2Xp2clxw1OGMxgM/Zx7772XsWPHcsopp7B582YAtm3bxty5cznyyCM5/vjj2bSpexWJUMIrlwOjRWQ4erqxS9HTvTnZDZwMLLImVR6LntKv27GbFsbOGwyGtvjt2+vZUFTVpfscPziF35wzIej6lStX8sILL7Bq1SpcLhfTp0/nyCOP5LrrrmPevHmMHj2apUuXcsMNN7Bw4cIurZuTdg29UsolIjeh53GMBP6llFovItdb6+cB9wBPicjXaJt7S4jzVR4ydj79cO1ZNxgM/ZdFixZx/vnnk5CQAMC5555LQ0MDixcv5uKLL/aUa2zsXuk5pAFTSqn30BNUO5fNc/xfBJzWtVULDXveFGPmDQZDW7TleXcn/k6o2+0mLS2N1atX91gdwj7XjcL26Hu5IgaDweDHCSecwOuvv059fT3V1dW8/fbbJCQkMHz4cF5++WVAqxJr1qzp1nqEv6E3Hr3BYOijTJ8+nUsuuYSpU6dy4YUXcvzxxwPw7LPP8uSTTzJlyhQmTJjAm2++2a31CPtcNx5Db1x6g8HQB7njjju44447Wi2fP39+j9Uh/D1669OYeYPBYAhM+Bt6o90YDAZDm4S/obc+xbL0ja4WmltaZV8wGAyGw5awN/R4NHr9+b2nVnDvuxt7rz4Gg8HQxwj/zlg7vNL6vq+ynsTYyN6rkMFgMPQxwt6jV34evVvpP4PBYDBowt/QW5+2Ru9WCrex9AaDoQ8we/ZsAD755BPOPvvsXqtH+Bt6P4++xa1wK2PoDQZD77N48eLergLQHwy9XwoEt1vRYuy8wWDoAyQlJbVatnz5cqZNm8b27T2S4BfoD52xHqNuSzeO2HqDwWCwef9W2P911+4zZxKccX/IxRcvXsyPfvQj3nzzTfLz87u2Lm0Q/obe+vRIN0rRYjR6g8HQx9i4cSPXXXcdCxYsYPDgwT167LA39LZLL56vRqM3GAwB6IDn3R0MGjSIhoYGVq1aZQx9R2nl0bsVbjMw1mAw9DHS0tJ48sknOe2000hMTGTOnDk9duyQOmNFZK6IbBaRAhG5NcD6X4jIautvnYi0iMiArq9ua7ypbrSlN1E3BoOhr5Kdnc3bb7/NjTfeyNKlS3vsuO169CISCTwGnIqeKHy5iLyllNpgl1FKPQA8YJU/B/g/pVRZ91TZF7fyjbpRSuv0BoPB0NvU1NQAMGfOHI8Hn5+fz/r163u0HqF49LOAAqXUdqVUE/ACcF4b5S8Dnu+KyoWC3fEaYVn6FqXMyFiDwWBwEIqhzwX2OL4XWstaISIJwFzg1SDrrxORFSKyoqSkpKN1DYi/TGNGxhoMBoMvoRj6QJneg1nSc4Avgsk2SqnHlVIzlFIzsrKyQq1jm/jbdLe7tfE3GAyHL+E6rqYr6x2KoS8Ehji+5wFFQcpeSg/KNuCVbkwcvcFg8CcuLo7S0tKwM/ZKKUpLS4mLi+uS/YUSXrkcGC0iw4G9aGN+uX8hEUkFTgS+0yU1CxFbpnEmNQuza2owGLqJvLw8CgsL6SqpuCeJi4sjLy+vS/bVrqFXSrlE5CbgAyAS+JdSar2IXG+tn2cVPR9YoJSq7ZKahUiLI+pGWUbeRN0YDAaA6Ohohg8f3tvV6HVCGjCllHoPeM9v2Ty/708BT3VVxULF2fFq/2s0eoPBYPAS9tkrbe/d7dDmTdSNwWAweAk7Q7++qJLfvr2e4qoGAOx5wJXyevLGzhsMBoOXsDP0Ow/W8e8vdlJe1wx4jbszmZmRbgwGg8FL2Bn6CDvVAb4yjVthpBuDwWAIQNgZerEC5u0MlU6N3l5m7LzBYDB4CTtDb3v0HpnG4dHby0x4pcFgMHgJQ0OvLb1ty1scGr3zf4MhEPPX7Wfz/moWbzvInAc+ZsuB6t6uksHQ7YSfobdq7PHe3d7v3mXG0PdX/vHpNh77uKDT29/62lqeWLSdL7eVsrO0ju89vbwLa2cw9E3CztA7Ux2AV7pRCqPR91P2Vzbgst7o972/iQc+2NyqzOurClm3t7LN/bS4FZX1zZTXNVPT6AJgT1k9JdWNXV/pLmT+un1s3FfV29UwhDHhZ+g9Gr3+bHHEzrcoE3XT39hVWsvR933EPxft8FnulOfqmlz834trOPuRz5n78GcUltcF3FdVfTNKQUVdEwdrmjzL1xe1/YLoDLWNri65DzcUVXH9f7/ixme/6oJaGbqa8tqm9gv1AcLO0NsaPfh79MrRMWsMfX/hsy06GdWm/VU+hvOfi7ZTZj1kK3eVe5Zv2l/Nur2Bvd/yOl2+or6Zg9WNjB6YBMD6Il2+obmFlbsOfWK0mkYXE37zAQ/9b8sh7+uxT7RMFRkRKFu4oTdZur2UGfd+yK7SHk3v1SnC1tD757Xx0eiNoe83LN2hDe+g1HiPoQb4/Xub+N07ejbLpdvtMjqla32zK+C+Kur1IDvt0TcyIiuRvPR4NliG/s3Ve7lo3pccsEZdB2PLgWpufn4Vja6WgOsrreO8sHx3SOfYFjsPaiNS1xT4WP6s21vJtpKaQz5uX6C5xc3Hm4u7bf+uFjcPLtjscRg6ypbiGlrcipufX8W1T6/o4tp1LWFo6PWn7d21eAy+0/j3QsXCnIq6rmmCfralhD1lgaWTzmAb+tpGFyU1vlp6o0vr9st2ljFlSBpv/+g4AKobAhv6yjrb0DdzsKaRzKRYJgxO8Ug3+ysbUQoKy+vbrNPv3t3IW2uK+HJbacD1Ndbx7fodCvZL50BVQ7tBBvVNLZz9yOec/OCn3Pz8qrCXMN9aXcTV/15OQXH3REZ9uPEAjyws4IEPNrVbduWuck8/kY2dhmVNYSWfbS3p00EgYWfoxd+jd8g1ZmRs51iyvZSpd/+PF5Ydugd65b+WcfKDn9LiVuwuPTSD3+hq8XSUVjc0c7Da72Uk+lqv31vJ1LxUkuN0MtYqy6P2x24RuNyK8rpmy9CnsrO0juqGZs/69jz6FOs4Ow4GbrJXN+jjN7ccmqFvbnFzsKaJzKQYXG5FaU0jZbVNrNpdHrB8mVX/mKgI3lpTxJZuMpA9hS2p7a/sus7yguIa/v2F7u9x9tPYlNY0trpuu0prufDvi3n4w6088MEmT/+Q8z5pcrnZ246D0JuEnaH3pEDwC6X0TWpmDL0/P31pNVc8uZSDNa0fmvnr9gNw62tf88Si7QG3f+h/W3jqix1BX6Ird5Xz7NJdADS1uLnzjXWc8MDHHi+6MxRXeeta0+hqVfcDlQ3sLK2ltqmFCbmpxEZFEhsVQVUQj77Cry6ZydqjB7jxuVV8vdf27Ns29DYrdpVTH0BSqbYieppbDu0+tF9yU4ekAVBU2cCFf1/M+X9bHHCsiN0xeNXsYUDo59FT3PbaWqb8dgFvrQk2QZ0vm/ZrQ19a23WG/pSHPuW3b2+grLbJI7GlxEUDWso5+aFPefwz32dgb4U24I9+XMBjH2/zSD37q3zr5ZTM3l5TxH3vb+yyeh8qIeWj70v4e/SeQVIon+yVSilP2cOd0ppGXvtqLwCPfVzAb86ZwNLtpbyxuohzpgxi8baDTMpNZUBiDPe9v4nzp+WyZHsZJ47NIik2itKaRv760VYA6ppbiImM4HvHDff5fS/8+2KfYz5vtQ72VtSTmqAfpKqGZlbsLOMb47JDqvd+h8dU1eDyGL5lt5/M3e9sYPWeCtZZXt/EwakApMRHezxqf/zlqaykGCZY29mdvtC+R19s1ePdtfvYeqCa7JQ4fvSN0cwaPgDwSkeH2pS36zElL40PNxazv7Le04qorG8mLSHGp7xtuMZmJ4d0Hp1h5a5yBiTGMDwzEYBlO8qoaWwO6Zq+s3Yf1Q0uFqzfz2njs1EK4mMiA5ZVSnlCSgN53p3BKQEVFNd4XozRkdrf3VVWR0Vdc6sWk3/47b7KBjKSYj3Sjc22khpOGjcQgJ+/vIZGl5tjRmQwZ+zALqn/oRCSRy8ic0Vks4gUiMitQcrMEZHVIrJeRD7t2mp6CZrUzO37YBmn3sv/NhwAICYygk379M0+79NtPL9sN7e8upYtB2o4e/IgbjtzHC1uxZ8WbObG577iv0u0h758p/fG/+P8zfzu3Y0+N7/teQViX6W3Ofvc0t1c89SKVg9Ii1vR0NzaM7Y90vwBCdQ0aI8+JjKCrORYctPjKa5q5OvCCmIiIxidrSNoUuKiqKrXhvbHL6zi8c+28dKKPdzzzgY+Lzjos/9hmYlkp8S2imjZ346BLKluRAQSYyLZcqCGRVsPsmD9fs/6YC+ajnLA8hinWB793gpvvQLF/tvS09gcbei7UvIAbXwv/PtiTn1IP95ut+Jb//iSa54K3hG5anc5NY0uKuubPS/Ad9buY/yv53PEr+fzi5fXBNyupLrRk6E2UCu0M6zZ4w2j3VpczT7r/qpt0vXaekB75Jv2e18IG/dVsWp3hc9+bA/f/0W63SHlDc1IAPQ97+SD9ftZvM33PuwJ2jX0IhIJPAacAYwHLhOR8X5l0oC/AecqpSYAF3d9VT3HAgLF0fvOFWsib7wsKjjI4NQ4zps6mK3FNSilWFuob/o9ZfqmnZyXxricFI4YlMLzy/YA3rDF5TvLiI2KYGJuimefFQ4dfMH6A0GPvc8hH2yxHqBdfp21Dy7YzJl/WcQXBQc58YGP+cIyyPaDNGpgEjWNLlbtqWBEViIiQk5KHE0tbj7aVMyE3BSPV5YcF01VQzNKKd5cXcTv39vEL19Zy1OLd/KV3wM7NjsZEWHp7Sczygq1hPYlj+KqBq6ePZz1d8/lwYunAL4dwM7/mw6hQ7a4Wtdj3KBkMhJjWFvorX9JdSN1TS5ueWUtpZYhtKWprORYMhJjgr6wmlvcLN520DNoLBgvrdjDYsfLcafV5+KyHr5lO72hqFV+L7eiinqW7Sjj/L8t5pjff8Sdb6wDvDKUW0FsVISns92fjQ5jW9pFht4pARUU11BkOSG11u9gSy+F5fWe3+aMvyziqcU7ffZTVFFPo6uF8rpmUuOjiYmKYFxOsidCCvDIO+V+rcgf/Gcll/9zaZecT0cIxaOfBRQopbYrpZqAF4Dz/MpcDrymlNoNoJTqtpioVknNHJE2To/+cNTp95TV8cSi7a302837qxk/OJUx2ckcrGlkfVEVpbVNPsZtnOUFnj9tsGfZV7vKUUrx6ZYSpuWnMSk31bOurLaJpdtLKaluDNo5CL5Gs8B6kPyjcr7YVsr2g7V8+4ml7Cqt492v9wHa0MdGRZCbFs/usjpW7Czj1PFaIshJ0aGU20tqmTE03bOvlPhoqhpcPjr9PedNYNEvT/I55pFD0z1OQ2ZSLMMsDwy0R9/c4uaO179u1VqprGumtqmFgSmxAFx4ZB4Tc1M4UO09T6dH7/+g2yzfWcbnW9v27A5UNRAZIWQmxjIxN5V31+7zrCuubmT1ngpeXLGHL6zoH1uaSo2PJjsljgNVDbjdiuYWN59sLubRhVtpaG7xGJtnvtwJaC/zpeV7PPt2Wy2se97ewCMLvekm/K/zm6v3ev7fst+34/fKfy3jWiu9RN6ABN62dPnTJ+QAcNGReXxrxpCgL5tNlmyTmxZP6SFIN3e8/jVn/XURAKW1TcRERjApN1Ubesszr7X6WQqKvRr7hxsOBJXeiirqPff1T08dw7s/Oo4x2cmeaC27Ex28cpq9vLcIRaPPBfY4vhcCR/mVGQNEi8gnQDLwF6XUM/47EpHrgOsA8vPzO1NfR1Kz1gOmnF68uxd+01++soazJg/mxDFZPXrc6oZmvvHgp5TXNuFyK2aPzGREViI/+M9KfnrqGHYcrGXuhBxGWfKGrddfMD2XP87fTHZKLOmJWu89b2ou972/idioCEprm3j1q70UFNdw7XGTfOSVu9/ewIZ9VVx0ZB6r91SQmxbvadJ+5+h8bjppNBf87QuP16SUYpv1IO12GPrmFner4f1fbitl/rp97K9qJCc1zhNN41Z4DP34wd7WxZFDB3j+T46LorCsztPc//MlUzh/Wh4AT109k8gIobnFzazhGT7HzErWhltEv5w+33qQZ5fuZldpHf+99ih+9Pwqdhys8QzGGmiVB/3SccoqNY6XTGlNE9nWS8nJxfO+BGDn/We1WmdTXtdMekI0ERHC5LxUPnX0I/zkxdWe6B9bCquoayY+OpK46EhyUuPYX9nAuY99zvqiKk9rd9XuCpZs1y8G2wP9wX9WAvCtmUMAuO/9jZ6RyGsLK2hxKyIjxCN9JcRE0uJWLFh/gGn5aazaXcHmA9XMGKavQ0NzC9tKajzH/NVZR3D5E9qLvXhGHplJMZw9eTB/XbiVaqv1JSIs2lrCC8v38JdLprJpfzWDUuMYkZUYULopqqinttHFaKs/wub6/6xkeFYit8wdB8CzlnSilKK0pomMpBgm5aXy6spCmizDa3v0m/ZXMy4nmT1ldfzkxdX8aUHrVBv62A2ejvtp+WmMzk4mLz2e977eR4tbefpwAI+MaNe5twjFow/Uo+n/qosCjgTOAk4HfiUiY1ptpNTjSqkZSqkZWVmdM4YR/vno3U7ppvc8eqUUL68s9MgOPcmm/dWUVDd6mtRX/msZf/pgM59uKeHJz3fQ4laMyUlm/CBtHJ9ftpvk2CjOmay997E5XqOZnRLH3789nQcvngrAH+ZvIjEmknOmDOaoEV7juMEyzp9vPUh5XTMnjfNez5tPHk1OapzH2ICWcGzPaU9ZPQXFNazbW8noO973kTeuPnYYOw7Wcv1/v+KrXeXkpMSRZBm0CPF2ug7NSOS7xwwFYOYwh0dvSTcHrYctM8lrkOeMHcjxo7P4xrhskmJ9fZyEGP09KymWRpfb05ls6/dvrynyGXE7xmFgBqbE+fQ7OKUb/9j/jlBV3+yJCJnoaE151lvHsfX6inr9YgB9HTfsq2Ld3iqGpCfw6OXT+N5xw/loU7Fn8NXeIIbnpRWFnv9rm7TRfvLzHR4Hoa6phReX76G0tomrjx1OUmwUmx0e/a7SOo+RH5AY4+mkBshIjOHiGUOIj4kkKTaK5hblGW/w8IdbeXftPt5YXcTGfVWMy0kmMynWpzO2xa1Ysr2U2fcv5NQ/f8ZDCzb75Dj6YttBn5HSNqv2VLCvsp4BiTF877jhNLe4SYqJYlxOMrWNLvZW1LNxXxXnTc1lye0nc/2JI4OOpyisqGflrnLioiM4wnqmhgxIwOVW7K9q8NzzowYm+Uhauw4x3PhQCMXQFwJDHN/zAP/4qEJgvlKqVil1EPgMmNI1VfRF/KQbn1w3Di++pw29y637CA5Fk+0sNX7hhAdrGnnic+2RfbhR6+djs5PJTonjzEk51De3cMnMIeSlx5M/IIGjHA8iwNyJgzj5iIFEiDYi04emkxgbxRGDUth0z1yfsrYOfJIjssA2ToNS49lf2YBSij/O14NSkmOj2FlayykPfcrZj3zu2eZ7xw3nnm9OZHKe16DtrahnRFYSyZZRzk2PJ8LRcXrXuRNY/etTyXAY85T4KKoaXB7j4DT0bREXrR+F1Hhd9wVWB3agwV+Lb/2Gj+HNTo6jtLbJM1K2qsHFYGuUbnuJ1tqisr6ZFKs+x43K5LvHDG0lQYE3CqiirolUKxLnG+O81+Otm47l7MmDfX5baD0wTCnFwZpGxmQn+Sx/dsku7n13A3Mn5HDfBZMAuP31r5mYm8IpRwxkUGqcj9e93RFmOGRAAlGREZ7WhzNSy15myzeZSbruDy7YTEFxDUcMSiEzKYaSmkaPLPXEou1c+vgSzz7+urCAHz2/Sv9edbrD17+zH+CCvy3mi4JSMpJiGZmVxJ1njefeCyaRlx5PTWOLpzP99AnZJMdFc/bkQa32AVgjqSv5YN1+puSlefqG8tLj9W9aVuc5/uiBSdQ1tXgkG2dLtqcHV4Vi6JcDo0VkuIjEAJcCb/mVeRM4XkSiRCQBLe10SxCpf1Izj3SDr3HvaenGNvDBhsUfKiXVjT4aopPi6uCdh3VNOhzSDof7+WljmT0yg+8dr8MjP/zpifzwxJGttouLjmSYtY3TqMVFe8PhBlhyT0JMpMfbj42K8JQZlBpHUWU9eyvqeWN1EVfNHsapE7JbeVy/Ons8vzp7PFccPZQzJg7i2uOGe9aNGpjk8eiHpCf4bCcirUIMU+KiaXK5uevt9UDohv4oS8o5bYI3TDAmKoJdZXU+L++BybEMTov32Tbb0uudg7vy0hMYkZUYsP/COcKyrcF9VQ0uj6FPjI3it+dNZMiAhFbl7OtfUddMmlX+1PHZ/HLuWH5yymjPb+TcdsbQdE+nos2yHWXM+N2HPlFWkRHCG6uLcCv4zbnjPecK8LfLjyQhJor4mEif8QTO6JN865iLbvkGy+84xafe9nU9WNPInz7YzOb91cRFR7C/qoH46Ei+ffRQrekr+P4zK2hobuGf1jiPb83I8+zH7vi0DWlxdSNKqYD6f4Z1z15z3HDOnTKYxNgoCoqreXRhAZNyUxmRpV9yzv4rJzedNIrmFkVRZQPHjPS2cO17c095vafT2paVquqb2VNW5xNt01WRWaHSrkavlHKJyE3AB0Ak8C+l1HoRud5aP08ptVFE5gNrATfwhFJqXXdU2D+pmXOyEXcvdsZ6DX33vGHO+usiiqsbfTTdL7eV8vhn29i4z9tsvmHOSFLjo7nvfe+w7iMGJRMTpd/pI7KSeO77R3vW2csDMTY7me0ltT6dsE6OHZXJ22uKmJSbSlJsFPHRkR49HSAnNY6GZjfbSmo95Z0x/aA96WuOHeb4HskdZx3Biyv2UN3gYtTAJE/SKH9DHwjbs7KNrv0yao8TxmSx9PaTaXK5eezjbYA2lu+u3cf2g94XbCAJxdbgC4pryEtPoLrBxaDUOPIzEli4qbjVmI4yRwdtTZPL0wIC3Umbnay3ra5vZki670sF4PvHD/fJ5llc1Yjb0oadkVE3zBnls53z95sxbAArdpV7wm2BVlFJvz57PH//dJvnt8xMiiUrydvfMGSArltcdCT1jv6bbcU1pMRFUd3o8nRy2y0lJ0mxetmjCwt4x+povmB6LieOyWJgchy5afHkpsXz2/MmcNtrX/PIwq0crGli3nemM3fiIMrrmvlkczGVliHdY2UtrWtqoabR5aOV22T43Q+JlnxUUd/MS9cf41nudGhW3HkKNQ0ufvfuRs6aPIiv91bS4lZc73CQBqfFI6LHYyzYsJ8pQ9IYnqnP/X8bDnDra1/7HLeirvU4iO4kpDh6pdR7SqkxSqmRSql7rWXzlFLzHGUeUEqNV0pNVEo93E31bZ3UzJkCwWHcezq80u7Y6S7pxr5pnb34j368lY83l/iE0eWmx3tieG2cHZcdwY7HDmbojxulPZpp+VojH5AY4/NAD0rVhsDOJTMgMYaTj/AdWDMoNb7VwDYR8fGsbI330llDaI+zJw/mkcumeb53JOtjdkqcT8epLUfZncU5KXFcc+zwVttNGZJGVnIsP3lxNRV1TdQ0ukiOi2J6fjpltU2ttFlnFIl/uoaL533JCQ98DOhrHchA3nHWeB64aLLn+9biGm5+YRW7y+o4blTwvq/MpBjiLQM2PT8NgNV7KjzrnZ3tl80awjXHDfcYxrSEaKIjI8hM9hon+7ppQ++moLiGq/69jE+3lDBz2AD+fdVMzyjdQNhOwTuOaKKclDjOm5rr4y1ffGQeg1Pj+Odn+uVmtzT/eeUM3rjxWADWFFb4yGwHqhpDGjCWaA3YGjoggZFZvl68HYmWmRTLsMxEnvjuDJLjorn3/Encf+Fkn5dBTFQEwzISeXttEQ3Nbh69bJqPDCgC7958HE9cOQPwfY57grBNgdBKo3f7JjPraY++sbl7Db3N79/dyMZ9VTS3uPlqV0Wr9blp8eQPSPRZlhWifOHPFUcP5S+XTg0oF4A2hOdNHcz503IBbQx8DH2aNpp2zpKMxBiykmN95BSnFOBk9MAkkmKjGJwax7icFHbef5bnhdIWkRHCGRNzQjvBADhbONMsY7jDapH8/PSxHDc6s9U2AxJjmPed6VTUNesRrFUNDEyJY/pQvf1Xu8tpcSuueWo5H28q9jH0zgfeKTUUltdR1eDV6P3JSfWN5Hln7T5+cspoLmvjZSgi5KXHk5kU42mZLNjgHejlHNxme5sZlm5ut4zsa3fF0UM9ZeOjI2hsbuHZpbv4ZHMJpbVNHDViAHPGDvTpP/HHv0M80HkBREVGMCkv1eNM5TqkM/verqhr9nj0oCOR7BQab954LL/75kTAG0ppk2jVIT+j9T3+xo3H8vVdpwWtvz9js5NRSt/neenxnmfhi4KDjMhMZMLgVNIT9bKKHjb0YZ8CwTurlJ9009MafUuL9dk9B46NiqDR5ebFFXt4ccUeXrjuaOqbWzhzUg7vfa0f1rMmDWL60HRPq2dMdpIe9TplcFu7DkpGUiznTc0Nuj4rOZa/XOr1nq87YYSncwq8aYPtNMADLKPx+S0n8faaIn7xylqP1+/PT04ZzcVH5nUqjUVUZARXHD3U8xB3Ftu732ZpzoG8a5tpQ9LJSo7l6cU7aXK5GZWVxOiBySTFRvHV7nJOHJPFwk3FLNxU7OnQBN/wO2de87fWFNHconxkHSf2OILkuCiqG1z88aLJfGtG+y2eIwalcKCqgcFp8UwYnMIXBd4MnM7O2QG2oU+MtT719+jICNbedRqJMd7fNt6SbpwvyaNH+IavBiLQuQ1Mbm3owevFp8RFkezYzn4RVjU0s6esnuRYLRltPlDtkVFHDkxibE4yG/ZVcfM3RvvsN8Hy6LMDHDfOClUNlXGDkpm/fj/jB6cgIp7za3S5Pa3i1Hj9O3ZVtthQCTtD75/UzO3R6P06Y3vao3d1r0efHBdNoyOywR5p+N1jhnkM/WPfnu5Zf/GReZw2IccTd96VPHPNLNYVVbYywv4vhaykWCJEZ3mMjhRP9ExcdCQjsvSDGyjGHCAvPYG8EDT5YNxjeXCd4c6zjqChuYXEmEhioiI8Hn1aQnBDHxEhnDxuIC9YA49GDkwkMkKYMiSVr3ZV+OSTf/Jzr77+/WdW8N9rj2LqkDSfbJ92uGJKfOBHNCc1jqgI4YY5o5gzNssT5tcev79gkifi47TxOZ7WFvgaevtcbY/eNvjQ2kDbnbF2WGF6QrQnlLctkhz9OQOTYymubgya+2aEZehz/e6JuGh9jWydfmp+Gou2HuS3b2/wtF7slsPvz5/Uar92KOzAIC3LjmBLPfa1cLbG7BaU/bv2tHQThobed85YZxy9M2Spp8OXursztqbRe2NkJsVQUFzD0IwEZlqDVM6b6uu1P3Bxt0S3Arrj8oQQBoVFRUaQnRLHvsoG0hNifF4MQwYkECG06k/oC1x7/AjP/xmJMZ5EYmltePQApxyR7TX0lt47PT+dxz4u8OkYdEZP1TS6+M4TS7nuhBGe/oSUuCjPcPxgrYjkuGhe/eFsRg1M6lDLxSmXXHHMUBZuOkBlfTM7S+t8pmC0pRrbk7dbY4GwO2P3VTYwa/gAXvrBMUHLBqvL09fM4r9LdnH0iAEByw7P1L9nblrrFmBKXDSVdc0Ultdz6vhshmcm8syXuygsr2f2yLZbFrahD9ay7AiT89KIiYxglvVMOq+dPSLYXuafSbW7CTtD7wmvtAdMebJX+nrxPZ0BoambPHo9AnE/Dc1uvnvMUH5yyhhueXUtCzYcYOawAURECKt/fapnwE9fY3BaPPsqG1pFvwxMjuONG4/1dPj2VdITYjz5elLb8OhBRxXFRkWQHBfl0bin56fjVnhGo84Yms6KXeWcN3Uwb67Ww1FqGl089L8txEVHkJ4QzdCMRK9HH0S6AW+ys84yIDGGN2/Sk7WMvfN9HyfFlixsjd0/WsVJfHQkDc3ao7f7NULBKfWMy0nm3gAet80wK4IlL0AUUkp8FAXFNTS1uMkbkMAVRw9l1e4Kvt5b6TOwLRA3njSK5hY3F0wPLlGGyuC0eJbfcYqnFRZrnd/MYemefq7oyAgGp8YFHIXanfRN69AGnhQI1ndn1I3Tie8vUTcfbjzAD62JoUcOTCI9MYYZw9JZsOGAx3PoyTCtjjJz2ABW7ir33PROJuel9XyFOkiGw5NtS6MHLWGcN3Wwzz1gG75FW3X6gp+eOgYFHDV8gMfQ2zQ0u5mSl0ZyXDSrmysAgnbGdjVpCdGebJng7ST39+wDERcdSXOLorC8jjMmda4jvL2+mKykWC6ZMYS5ATraU+OjPRKUHY46MTc1JEOflRzb5gumozidARFhxZ2ntGoJLr7t5C47XqiEXdRN2yNjez+Ovqs7Y52TGdiG5rTxOUzPT2POuPblk97GHqG5wS+fTbhgG7r46Ehio9rvmPvjRVN42NFBnZYQw4isRJZY89qmJkRz7KhMoiJbP3pHDk3ngYumePLugHf0aHdjH3NsdjILf3YiowbaoYW2dBNcw7ZDNt0KBneBBBIIEeEPF00O2MmbEhftieO3PWd7FPDYnMADn3qKzKTYgNe6pwlfj94RVml/9x0Z20uGvos9+l0Hvbqp7d0Ny0zktRuO7dLjdBd2vPbciYGHlPd1bEM/yS99QEcYmZXEdqtD1xmtkpMS5xkDcf60XP58yVRAT4gCOlQ0N4BU0R3kpsWzbm8VKfFRnjEMoCc9+empYzhpbHCnIs7RgRqscz0Y7//4eJ9Bdp3Bfi5EvBr+N6fmIuhoKEMYG3pPPnof6cbp0fdsvWxPvqs7Y3c6Qu7a0mv7KlGREaz61akkxIYeptaXsDMbzmnD0LWHs+nu/B2W3H4yLyzbza2vfc2wDO/YB9u7HpQaF1IroiuwI5z8O3ajIiO4+eTRgTbxEO8IQcxso9M2EKFGC7VFqqWJ5w9I8PQtxMdEcumszmXI7Y+EoaHXn/7zw+oZprzlejrqxhte2bW5bpyjKtvTiPsq6SGmIeiLjLYkjFOO6HyYqjMsM9Gv09z22O3ORl3eyk1zCOGlHcX2hF2dmOc2PkD+o57EdoDsa2VoTdgZ+lYDplQwjz78wyvrm1p80hsEi6k2dB/XHDecuRNzgo4ODgVnZ3m83wCcWcMHcPPJo31SQ9jx5e2FBnYl9gunor7jA3niY7wadG8YetuLzz+Ea9Tf6f1egg4ifgOmfDpgHf/vLK1tNdNSd+LsjO2q49r5wu3IjbT48PWMw5XICDkkIw++LbEIv9w7sVGR/PTUMT4x5XPGZPHkd2dww0m+Scm6EztssTMDeeIc8lJvyIv2VIM5qYc+6Km/EnaG3tsZq7+3BElqdtNzq/isnanauhJbo1fKO6fmoWJnDfzF6WPZef9ZbWaaNPRd2hpRGwgR4eQjsjuUkO1QsaWbQKkA2sPZGev/IusJ7DDdUNIuHK6EnRbQ1pyx/va1qyYVDgVntE2Ty+2T86Wz2LMTdTYpmaFvEA4tsbQEnZht+tCOR6n4y1E9zQXTczl+dCYDOxjxczgRdobeX6P3ibrxs/Sd6VjqLP6GPrELbLM9HZ4zrtoQfnTUo+8tOhsC29uGXkSMkW+HkNxOEZkrIptFpEBEbg2wfo6IVIrIauvv111fVY1/UjPvxCOtO2C7K5NkIJzH6qrjltQ0Eh0pYRttY9D09+sXLBGZoe/QrkcvIpHAY8Cp6Llhl4vIW0qpDX5FFymlzu6GOvrgn9TMHSSpGeCZq7En8Pfou4KS6kYyk2I7larX0HdoL0dOuNORVL6G3iEUj34WUKCU2q6UagJeAM7r3moFxzln7IaiKk/Hp394JXSvod+8v9pnhndnWOWhzBurlPLMKVpS3Whkm35A8iHmxe/r2NLNlEMYPWzoXkIx9LnAHsf3QmuZP8eIyBoReV9EJgTakYhcJyIrRGRFSUlJJ6rr69H/7OU1nuUqQGdsczdp9K4WN6c//BnXPbMCgJ++uJrnl+32rD+UWPrHP9vOqDve59ZX1/LplhLTEdsP6O8tspioCJ77/lE8fc2s3q6KIQihuBqB7lJ/C/oVMFQpVSMiZwJvAK3GTSulHgceB5gxY0anrLA3jl7Ptzl7ZAaDUuNZsr20x6QbexDTip3lALy2aq/P+kORbv7xmZ7l3s5rXlrbszPRGLqH+y+YxKiBvZtgqzuZPbL1FIuGvkMoHn0h4JyjLA/wya+qlKpSStVY/78HRItIt1x5Z1KzphY3QzMSiYzQHr7/QKXuMvT2TDzBUsi2Z+h/9cY6PtsSuEWT4Nex9e2jTL6O/sCls/KZMSzwpBoGQ3cTiqFfDowWkeEiEgNcCrzlLCAiOWK1T0VklrXf0lZ76gKcSc0amluIi44gQsTqjPUt213SjW3og0VTtBV1U9fk4j9LdnHlv5a1Wud2Kw46Yv+X3HYyF4cwD6jBYDC0RbvSjVLKJSI3AR8AkcC/lFLrReR6a/084CLghyLiAuqBS1U35R9wDphqbHYTFx1JQ7Nb56P3D6/swrwz6/ZWcsWTS5n/kxM8U64FmsU+2HGbXG6+2l1OvWPuULdb+Ywk3FNeR0Ozm7joCIZlJJKTamKDDQbDoRNSOIAlx7znt2ye4/9HgUe7tmqBsTu2WtxauomNiiBCtJTTndLNxn1VlNc1s2l/tcejb7AmO0iOi/LMPQm6M7aoop6KumbGD9ZpWP+7ZBd3v+Mbkfrl9lIO1jRyzuTBREQIWw7oSUaeueYoz6zxBoPBcKiEZdyXiNfIxkVHItJ6hinoWkNvJ3s6UNng8eht4+7fdmlyuZl9/0IAdt5/FgBFVoIyJ1c8uRS3gq92lfPb8yayy8o9P3pgkhmEYjAYuoywzJIVIUKdJYFoj16sGaZ8y3VlCoQqy9Dvq2zwJBurbmjG7VbUNmmDb8cR1zW1jqMvq/ONnrlgeq6nvm+v3Qfo3PN6Yun+PcDGYDD0LGFq6PHMERkXHWl1xnZvCoQKy9Dvr2qgsl4b9tqmFqobXSgFt50xjv9cexSgO1z9Ka9tYmiGN93tbWccwQ9OGMEPThxBWW0Tja4WdpbWMiwjsd/HXRsMhp4lLA29iHg6NeOiIyzppntTINjSzf7KeqoamomxslMesGLqk+KiSLBGCNY2ej16e5RrWV0zQ/2mi7vtzCMYkamXlVQ3srusjvwMM3mCwWDoWsJSo48QHNJNpCXdtPbouzK80jb0u8rqaHK5GZGZyPaDteyrtAx9bBRRkRHERUd4pBzQXn9qfATltU0Mz0jgtRtm+8TK21n39pbXs7e8nrMnh+ck2gaDoe8Slh59hIhDutFRN4HSFHelR19Rpw39zoO6w9Seeu27Vjy8PRdoYkwUNY1eQ//Gqr0UltdRXttEemIM0/PTGZfjnRDZnuhhyfYyXG7lM0m0wWAwdAXha+ht6SYqErEGTPl3xnYmjv7/XlzNJ5uLWy23O2PtY9hTr9nY83wmxEZSVuPteP3NW+t58vMdVDe6GJDQegKK7BSdy+alFTrlwTE9OE+owWA4PAhLQy94O2NjLY2+odnNf5bs8inX0Sn9mlxuXl+1l483+Rr6zfurKaqs95lgISfFz9DHej36fY4JvcHbCkgPMHFyekIM0ZHC3op6xmYnk5duNHqDwdC1hKehFzweva3RB6Kj0o0tuZQ40hDYmSobmt2MyfYmpZowOMVn20Tb0MdGsc8vZn5nqY67HxDA0EdEiKcv4fQJ2R2qr8FgMIRCWBr6iIjWGn0gOirdVFv55UuqG3lnbRHT7l7A13srPetHZyd7/h+WmciO+870fE+M1d5+QkwkxdW+c9XutAZCpQeQbgDOmTKYo0cM4IaTRnWovgaDwRAKYRp1I55Y9UPx6L/cVsqY7CQyrJzv9kjXkupGPt5UQnldM997eoWnfGyU972YEh+FiJCXHk9heT3JsXqQU6D8N3YwUCCPHuCvl041sfMGg6HbCE+P3tLkwU6BEMzQB9foV+ws47J/LuGBDzZ7llU5PHo7BLLMygc/a5j2uBOt5Slx2rC/cv1s/nLpVE/KgoSY4O/OYLNFGSNvMBi6k7A09E7DGNuGdONqw6P/w/xNgG9+HNujr21q8QyEsvn31TPJTYsnJzWOmKgIzzyZOalxnDfVO+GWLeH4ExkhpPXzSaINBkPfJCwNvdOwx0VFIgEnwYImP49+64FqTn7wE0prGj3auzMwx5mB0tbVbezO1pzUOI83Hwjbo/d30gckxvikJDYYDIaeIiwNvW3YIwSiI8XH8DvlEX+N/m+fbGNbSS3/XbLbI/3UNHon+K52TPa9s7SOKUPSWh37mBEZzByWHrRuSZZHP3SAb5ikMzTTYDAYepKwNPS2YY+1BkvZnvLE3BQW/fIkTzl/Qx8dqctt3FflWebMS1Pj8OibXG4yrc7TTMcE3Td9YzR//86RQetme/RD/Ay9syPXYDAYepKQrI+IzBWRzSJSICK3tlFupoi0iMhFXVfFgMcBdGil/q6XR0VEEOlw7/0NfZSViGzTfm3oJ+Wm+qQrqG70zTqZFBfFsjtO5qOfntjhOvoPfIoxht5gMPQS7VofEYkEHgPOAMYDl4nI+CDl/oCecrBbibBqHRul5RA7vDImMsIn1LK5xXfWqWjrJbCztI4BiTEMGRBPTaPLM4lJdUMzGYkxnhZDclwUA5PjSO1Afvj9Vieuf4qE6Ehj6A0GQ+8QivWZBRQopbYrpZqAF4DzApT7EfAq0DpRTBcT4efR24Y5OkpaReA40yA4je3A5FgSY6IoKK5h3K/ms6GoiqoGF6kJ0QxI1FJNUmzHo2S+O3sYx4/O5PJZ+Sz65Uk8fMlUwEg3BoOh9whlwFQusMfxvRA4yllARHKB84FvADOD7UhErgOuA8jPz+9oXT3Yht423M7vIsKZk3Ioq21iyfYymlvcnnJRDkOfGBvliaQB+HRLCdUNLpLjoomNiuRgTSPJcR0fT5abFs9/vqd/nvTEGI+Bv2r2sI6fqMFgMHQBobiZgWIC/UciPQzcopRqPYeecyOlHldKzVBKzcjKygqxisErFOUnh0RZms7fvn0kp43PAaDZ5fTovaeSGBvlY8gLy+s4WN1ISlyUJ3In0CjXjjIwJY6d95/FGZNMnnmDwdA7hGLJCoEhju95QJFfmRnAC1YnaSZwpoi4lFJvdEUl/bFleNtwezT6KK8hj7Y8aed0gk7pJjEm0sejf3bpbgB+OGckxVU6V01XGHqDwWDobULx6JcDo0VkuIjEAJcCbzkLKKWGK6WGKaWGAa8AN3SXkQevYbcjbGxdPsZhyO2OV2fkjVO/T4jxlW4AfnbqGH5x2lgyk3VYZVInpBuDwWDoa7RryZRSLhG5CR1NEwn8Sym1XkSut9bP6+Y6tsI29FG2oY/w1eyd/7sco2OduW8SYyM9g5sAbpk7jh+cMIKICCHLipvvjEZvMBgMfY2QLJlS6j3gPb9lAQ28UuqqQ69W2zjj5iGwZh9IunHmtUmMjfJE1YzLSeaHc0Z61tnzuLaV6sBgMBjChbCM+fN49JZGb4dQxjg6W+3/nTnpm93e/xNjIj2ZKONjfNMTnDY+m999cyLjB/lOLmIwGAzhSFhqE/aAqUg/Hd4p3diDqZwevVPGSYiJotl6QST4Gfq46Ei+c/TQrq+4wWAw9AJh6dHbSc1s6cbW3qMdg5LslAONzd6IT1/pJpKRWYkAXDKz8zH9BoPB0NcJT4/eo9H7yjO+Hn1rjd4ZgZMYG0VeegLbf3+mSR9sMBj6NeHp0dvhlZG+0o2PRu/x6ANLN4lWlklj5A0GQ38nLA29J7dNRzV6t1OjN/nhDQbD4UGYGnp7wJSfRh8ZQKN3eTV6l9tXujEYDIbDgbA29FH+Hn1Ua40+qHRjDL3BYDhMCEtD7xkw5afRR0e01uh9pRvfOHqDwWA4HAhLt1b8om4CSTdOj766oZnoyAjfOHrj0RsMhsOEsPTo/TX6CYP1CNbhVlw8+Hr0k+5awHmPfuEZIDUiM5EEM1m3wWA4TAhLt9Y70Yj+vObY4Rw/OouxOcmeMnYmy1prHtjNB6rJTI5hxtB0Xvnh7B6uscFgMPQeYenR29JNpCN7pdPI6zJCbFQEW4trPMuaW5TP5OEGg8FwOBCWht4/6iYYMVERbCiq8nx3OaYVNBgMhsOFsLR6dpeq/1SC/sRGRbK3oh6A7JRYWtzKE6ljMBgMhwvhaeiVNvXtyTCxjrj6CBGaW1S7rQCDwWDob4Rk6EVkrohsFpECEbk1wPrzRGStiKwWkRUiclzXV9WLnYWyPaPtNPR1TS243G5PxkuDwWA4XGg36kZEIoHHgFPRE4UvF5G3lFIbHMU+At5SSikRmQy8BIzrjgqDw9C3I93EOAx9fVMLLiPdGAyGw5BQ3NtZQIFSartSqgl4ATjPWUApVaNsPQUS8cro3YJ9pI549E0tbhqb3Ua6MRgMhx2hGPpcYI/je6G1zAcROV9ENgHvAtcE2pGIXGdJOytKSko6U18A3CFq9E6PHqCqobndVoDBYDD0N0KxeoGsaSuPXSn1ulJqHPBN4J5AO1JKPa6UmqGUmpGVldWhijppUXbKg9AMfUqcVqiqG1ztbmMwGAz9jVAMfSEwxPE9DygKVlgp9RkwUkQyD7FuQbHTyke207FqSzyDUuM9y8yAKYPBcLgRiqFfDowWkeEiEgNcCrzlLCAio8Sa9klEpgMxQGlXV9bG7g5oT29vsOaLzU6N8ywzUTcGg+Fwo92oG6WUS0RuAj4AIoF/KaXWi8j11vp5wIXAlSLSDNQDlzg6Z7scb9RN24a+3spFn5MS61lmpBuDwXC4EVJSM6XUe8B7fsvmOf7/A/CHrq1acEKNo2+0PPocH+nGePQGg+HwIiytngpRo7elm0EO6cZ49AaD4XAjLA29HV7ZnnTT4LKlG6PRGwyGw5ewtHotIXbG1jfZ0o3D0BuP3mAwHGaEpaH3SjdtG+1zpgwCINvHozeG3mAwHF6E5QxTdmdse7nl7z1/EredcQSJjvlhzchYg8FwuBHWhr49jz46MoL0xBifZcajNxgMhxth6d6GOmDKSV66DrE0Gr3BYDjcCEtD7/Zkrwy9+qMHJgEQKcbQGwyGw4uwNPQtIYZXOhmTrScPL65u7JY6GQwGQ18lLA19qFMJOpmUlwpAkxVbbzAYDIcLYdkZa0s30R2Qbs6aNIiaC1ycM2VwN9XKYDAY+iZhaeg9UTcdkG5EhEtn5XdXlQwGg6HPEpbSjTvEpGYGg6GXaaiCJfO8oxwNvUJ4GvpOhFcaDGGPUlBX1tu16Bgb34L5t0Dxht6uyWFNmBp6/WkSlBkOK7Z8AA+Og9qDvV2T0Kk5oD/rum0eIkMIhGQpRWSuiGwWkQIRuTXA+m+LyFrrb7GITOn6qnqxwys7otEbDGFPaQG0NELF7t6uSejUlOhPY+h7lXYNvYhEAo8BZwDjgctEZLxfsR3AiUqpyeiJwR/v6oo66czIWIMh7Km3ZJtwMpoejz7MJKd+Rige/SygQCm1XSnVBLwAnOcsoJRarJQqt74uQU8g3m14pZsghr6qCFxmYJShn2Eby9qS3q1HR7DrWm8MfW8SiqHPBfY4vhday4LxPeD9QCtE5DoRWSEiK0pKOn+znj4hGwgyYEopeOgIePnq0HfYUNnpuhgMPUa95UuFlUZfrD+NR9+rhGLoA7nNAWOlROQktKG/JdB6pdTjSqkZSqkZWVlZodfSjz9fMpWlt5+MBMpbo6yRr5vfDW1nu76EP46Aij3tlzUYehOPdBNOht5IN32BUAx9ITDE8T0PKPIvJCKTgSeA85RS3SoixkZF+kwm4oPb1bGdVRbqbSp2HXrFDIbupM726MNEo3c1QUOF/j+c+hVCRSk4EB5ho6EY+uXAaBEZLiIxwKXAW84CIpIPvAZcoZTa0vXV7AAtzR0sb2n59RVdXhWDoUupDzON3lnP/qjRL/4r/P0Y2PtVb9ekXdo19EopF3AT8AGwEXhJKbVeRK4XkeutYr8GMoC/ichqEVnRbTVuj4569Hanra1/9lV2fg6r/tvbtehdvvgr/OPE0Mr+79fwwR3dW5+exr5Hw0W6qbX0+ejE/indrHxKf5bv6NVqhEJIuW6UUu8B7/ktm+f4/1rg2q6tWidxt3SsfEuT/uzrhv6ps/TntO94l31yP+xYBFeH2B8R7vzvV/qzxQWR7dy62z8F1cF7oS/T3ADNdfr/cOmMteWazNFQ1veNYYeoPQhl2/X/5Tt7tSqh0P+Glrod0s2uxfDaD9rOs2F79LaW2BdxNXn/dzvSLBethr0rD788ItX72i9TX+7VtPsDtvQRm6oHTH36QO/WJxTs33/geGishKba3q1PV+Fq0o5XlNVPaAx9L+CUbv5zAax9AZpqgpfvSx69UrD/69bL9632/u+sZ20JuOoPv/DQyhAipBoqelcX3vhO5+QKVyM8c56W6pzY+zr5V5B/DCz/56HXsbuxPfq8I/VnGBjEkFj7IpRsggufgLyZHW+tvPdL2NSzrfD+behjEvVnWx2trj7UGbv9E5h3HBxY713matQSjY2zg8vWau0Qts5SfQCev6znddTN86GxuuPbVRa2vd7dol9+zXVa8nAeb+M7HT9eR6krgxe/DS9+p/2y/uz8XN8HH/7Wd7md9mDwNBh3lr7mvSHhuFvgvxfqvDvtUVcKEqHrDFC6rXvrFgotrs5FLTlbzSufguyJMO5sSB8O5R2I2GuohGX/gOKNHa/DIdD/DH2Lw9DH6nli2/TWW/pQZ6zt8Tg9n03vwraP9E0FvobevmGr9x/acQuXw+b3YPeSQ9tPR6gqgucvgTUv6O/ulvb7VyKi9Wd7Hr2zheO8rs9fog1wd2Mfv7ATMQnbFurPgeN8l5ds0p+ZYyDbykDidAi6kxYXlGzW/1fugYIP4blvtb9dXSnEp0PGKP3d1rRB37s7Puv6ugajoRKeOhueuxgeGNExB6NwJdw3xFvfqr0waCqIQPowqCoMHO339Lnw5k2+y/at0Z/2y6+H6H+G3sej1/PEtqm/u/qQdGOPInQa7v1rISIKTviF/l5bor2+xmposm7WQ/Xo7Zu+J6MHqiyd3dbbX/0evH598PItLm//S3sevfNaBpJvOtph31Eaq/Sn7UR8dLf+a4/mBp3WF3wdFtCGNiUX4lJg4AS9LNTUv1VF8MDoznuRG9+Cx2bBZw+07ZVX7dOyhP1M1ZVCQgbEpUJCpq+hX3AHPPNNfa1qS+G5S33XV+3z7Zs6VHYvgZ2LvC/SnV+Etp3bDS9fpZ+1Te9qz772ICRm6PXJOXqQpn/rqnAl7PgU1r6kc/LbFK3Sn8bQHyIBpZsQPPqOdMY2N8D9+V0f7mgbbKfh3r8Ossbphxxg95e6I+iLv3jLhNI52Ra2oe/JyAg79M4+1/3r2jZEzmvY3ihmpwwXSI7q7uyPzge7shAWPaj/2uPzh7x18+93KdmkvXmApIHagIZquDe8qX/v5U8EXr/xbfjqmeDb28Z94e9gz1Lvcn8JZOPbWpbYv1Z/tw09wIARXkNeXwHr39BRUbu+hA9/A1ve97bummrh0Zk6Tv1QaG7w3tO2J22z/WPf71/+Dd76kf7/0wfgkSO143VgHVRa1+TAev2suJv1iwsg0RrhX+uX6uGrp0AitX1x6vFFqyAtHxIGHNq5dZB+aOgdTSiPdFMRvHxnPPq1L+oHcfEjHa5em9hGz/boD2zQzeTsidaNIV5tdMt873bVh+rRW0alJz16z0vNkqJqi9vuPHWuK2tH623Po+9urdhppJc/2X75+gr46j/a+I04CYYe5+t4uN1wcKt+4YOWDJIHe1uA7dFcrz+jgowm//g+WPAr34guJ7WO43zmiPbZ5ecV+4cb1pV5DX3ukfolUVcG617RQQQSoWXJ9W/oMrbDsvcr7UFv/6Tt81r9PBR8FHz9h3fBX6dC5d7Wht4pG7W44PM/65fd6ufh49/plND713pbAOPO1i2CTVYfT6Jl6JMG6s/aEn1f/WkMrHtV1+uIs7VR//olXQelYM8y/Vv0MP3Q0Dua5R3y6CtDb9KvflZ/Zo1ru1wwlj8RWF+tcXi5jTV61B0KciZCRKR+aOxUDc7onJr9umnvrP9zl8CqZ0OrT3d69OU7A3vQ9rnWFusO54bKtjuD7XV5M/U+nZ2s/jiNZKB9lm5tp9KHSKPDo1/miI4JJkUs/Qe8dRMc3AxDjtJSh/NlUb0Pmmshc5R3WWJG6AOn7PqUbW99T9SXawmooQJK/FoIix/RnddVRZB1hHY4lFt7sYlZ2oCBvu+ePheW/l1/t6+3rdGDHv/R0qSdpK+egZxJMPRY/SzYEqTdQtlj9RUVrvDVvusrfF9GH92tvf7Sbb6dpUpZ6QnW6e9L58G+td71STn6Xt88X79ktn/sfZm970jTVVmoHa2B4719ZG/8UH/6e/Tzb4enz9GO5kd3ax1/xEkw6WL9svjzeFj5b718+An0NP3Q0Dukm8gY/dmWoXemMw41TLHKSvXjn7/jle/pm6ct6srg3Z/B5w+3XueUbmyPI382TLlM/2/fVE6SB+uyj83STWvQXvKW+fDmDSGdjsfQV+zuev36zZsCa++el1qJ939Xvdf79Mf2zPOP1samLa/ex6O3/ndq3t3u0TsMfVM15EzW/9+Xq+WNxmpdh+Z6bci2OiJYBk+F+DTfVqjtIacP9y5LyNS68KKH4NmLoanOtw47FnlHBtt9Gpvf0/eEsyNyz3I8OQqdurVSsOBO3XldvQ9SBsEZf9Drhh4Lky/R93pdme6g3fGpd9uKXda0hw7pJmeiNu6LH9He9bQrvcYTtEEs3qi3223JQ656r4GuK4OHJ8GTp2ppx9Wk67X9Ey2zfP2KtU0TPDAKPv2jbvkArPiXll9GnarllAnf1Pt+/hJ4+bs64iwhEyZ9S7duU/N1uV1f6v6wcWfDxAvhiHO89bU1evuZPLhZG3Hn9RoxR/9ONu/8n/4cHuLo7i6k/xl6pwdgG/229PcWh5dVXw5PnKI1zbbwDEV3eIvNDbpJWvBh29sWWXkxvn4J/jLV+xAq5eiMPaD3E50IV77hbSZOv0J/2no9wLDjvJE4ay0Pq9jRWrDrqJROH7BkHq2wH3x3s/dmBV+D5cT2mEKhaq8e2OX/ArFfarXFvhJEMK/eXj7kaP15sI2USvb1iYj2viCaHMbtUDuv28P2oLOO0AEBM6yU2S1N2nj853x4ZDrcm6PvN2eulEFTW3v0tqQ2wGnoM/Tyj34LWxd4o3LcbljzIjx9Nnz5qH4Z+Hde71sLz1+uHYS9K7SEkjjQ60mDb79PZaF2KIYdBz9eC2f/WRs/d7PuMzro10Iq36VbD+5mHZViM2KO9/464hyYcL533dDZerzLXqsT034JHPha32tL/q5/170r4MvHrP3Y96CCjW/q1spTZ+qWzie/h5ItOpDBHkdz5h/h51t0PWyGHa9bHd96BmZbGv3wEyBlsB6DA7o1EhUDJzlSatgvsNhk33M/9scw6hQ44wF9vbLGwi+2wxl/1OvTh+v+ih4mpBQIYYXTo7eNfpsevUMCKFqlQw1fuw7GnxekfJP3xnE2nUOd/adwpff/8h2w/nU45iZtbF31EBUP1UWw4kk44lyIivWWP+ZGyJ6gm9BFq7SBb673NqFrS7RBdcpC2xbCpIu07r9vtd7f0ZaHvf51iEny9fDKdmhdsXAFPHEyXPAETL5Yd7zFJunt3/k/3dy//CX4xwlw7iMwIoiXUlemZYeyHb7Sg/1yamnylVLqyyA1wHQHdr/F0NmA6Ie4tlQ3/dOHwRSH51Rfrs8rLs0bDeEclel/jT77EwyaAqNPDXwOHaWhUl/HG5fo6+E05BW79D1mYw+Gm/UDbVBTBmlD31TtTfVQtkN7mKmOJLL2y9+mcg/kTofCZfD6dd7lRatbG/qvX9ZpvPcshTFztZQxaLLuE1r/Oow9yxtOCfrFmDJI/58+VH8OnqZbzC9crl8UTip2eUN184/xLh9+ovbocyZ793f8z7Xhyxipv7/7U/3cnnIXbP2fbvm893N9nYccraOOlv6jddRKwULdWnJSsx9mXqtfAGn5XgOblu8tc9ZDOkWD7f2f+4huRZdt07/p0GO95+w00LZ040yVfvlLMOZ0WpGYAUf9AEafpuXkQOnVu5l+aOgdnqNt9NvrjI2M1Vq93fHivBH8sV8asanaYCilL5xtPNrTTfeu0NvaHaAL7tQyjn3jzvye9sTA14Owsb0R2yg5m9stjVqXPLAB4gfoB2bn59rQ217+3pVa/wcdNgb6Zk7L19JN+Q7gRG8Y2GvXwvhzdX/BpIu1F7by33rdhjf1Q/3hXXDVu7D+Nd1vkTfDqo/L25rav9bX0Ncc0N6W26U9NJtgHn35TkgepDuls8bqjrGoWO25AUz+lvcBOrBeG46IaK9nap8z+I5FaG6Ahffo/+86hBHGVft0qJ2I9jzjUvTyiEhvhx14PW8nMUlw+u+9+Xvi0vTn32fDFa/pc0/Ng8ho7za2R2ljG3P/fpbCZV6p0caOFotN1sYwaaC+blvm63vi3EdaS2jJg3y/R8fpe6xmv3cOiLFn6euz6j9aIoqI8kYKgZbdYpK0A2NzspW/yG7B7FujO6MzR2uPuHSbfmbShsKlz+r78tmLYPVzvvVpdrzIY5K9LbjMsXDWn3SLxcb5whwwwtfwTr9Sf9p2ZNxZ3nVOp8vu/3OSNbb1MifOFlkP0/+kG2fUjcejb6OTr6URkvWMVR5D77wp/LENfcZIbaTsZrptoNobXVq8Sb/1v/8xTLjA2uYgFPxP34Cn3AXXfQrfebX1oJlADDxCf2ZaN1nJFm3UcybBkFm6aQ3evNlul/bk1r3i3UdjlZYZIqK9HVtOeWPj2/r7+je8Rh5gmTU1cFQcfHA7vHkjvPZ973qnZObsPFZKe+i51gvBDscDra0GiiIq3+mVAY44R0d8OKMy7OvidmtjkDtDe412S8BuhaXl+xp6p8zV2eilotXw0Dh44du6ldFQBbEp3vVOQ19a0Hr7IUf5JmmLS9WfBzfrmdLKd7Q2Ek5DH53oDTn17/j+/GEdxhjjkBjsZyQiSl/X5Bzd4WhT8KGWxmKStdEF3Yr0Z/y5vt8ve05LFHNu19+HnwgRDhMTmww/WgnH/aT1vuJSvec9dLb+HDBSh13WHIDjf6ZbMUOP1S2J9a/5bh9tGd5jfwI/+FRvC7p1OO07MOY0x7Ec1yZYcryYBP05dm7g9YG88tQ2HMReph8aeod0Y9/Q/h1VTlxNuukKXu+vrU5Z26BkjtaftjRge/RtDUt3NelRdAOG62b2affAOX+BX5XCaffCZS9or23wVK3zhULCAN1RdMLPAdEPaPFG/WAOna09yKp9ukkeP0A/3BvehIX3WjsQ/XKKT9NGcPFf4ZVrfGPV7XSsVYXaGxxtPTR2VENpAWy2Zo8s2w5/OwbuztBD5W32f63r9cEd+sFtroOJF8D1fiF6uz73tmiclO/UXh3A+G9qL9IZC217tAe36BdX7pHaC/V49JaHlz5cX8OGKtj2sW/YXcH/gv/Oe5Z55Qh3i69BtaWYze/Ck6f4evQA0fHB9xsR5asZg9fggfbI9631vshtnNJN2hDvaGH/CXTsqLIhM/XnlMu96yp26xdh0kBfp2LTe7D2Zd0yu/pduK0Q8o9qXfdT74GfbtROwuyb9bKYBJhzC9y4HM4P0B+UnOPbMnFiy2v5Vj9Mxghva8GOVIlJ0JFXoO9n0NEtI0/S/8+8VjthEy0nyr8lYpN7pO58DcZ5f4MLn2ytp3//Y/2cOjn7YTjyat+XWh+j79ass/ho9Nb/zg5Xf1oafT0uaFtn93j0lgzh8eRLvZ+BOiordmstXbm9nmlqHhx5lfYqZt+kIxM6w0X/0tJF6hBtrFwNepi8Le88OlN3AucfA4Onw1dP63oefSOgdMdWbLLXWKx/TXuT2ZP0952LfCN+Jl8CKdb87ym5Vofqfu/DXrxBd3LZ+nPyIG3o/3a0NuLbrQiN9GH6nG9cpjvDbOyXhk1Dla6j/btlT/D+/rbHWVmoZRg7lXHeTG1UGir1i9726O19vPMT+M83tT4fm6o1dedsQa4m+PeZupXXWK2jPf5l6a/Ln9QRIPZ5HFinO/Tm3q9fdPvX+Xr0oPth7A7G9GHacOVMhpuWw9E/9C0bEak/h1jG1d0Mo072LWNrxIi+j+wXnfMFZF+PCed75cjpV8BRP9T3Qkujbt0k5XhfJIOn6eM1VsKcW/Uy/w5Hm+g43Wl54xLttDjJGtP6uWqPkd/Qn7Yht6/V4Om+nbp29MuxP4YffgmX/Fd78sfcpH8L0K2KaxZohyoQ318IF7aRGC41V0ue/uROh7Fn+C6bcTWc83DwffUB+o9G31yvvZOWAB69M4TSH1dT6xvZqb37Y8tAHkNvGXiPdNBseXQOr+zABism3iK9m7S6jJFeL3fgeN3BOP27Om757If1jbvwd9pLHHWKbpYusfTx2GStE79/i27qF63SHk/FLn0+Y8/UaRi2zNcd1RkjtRGNSdKdtgCzvq8N3bDjIG+Wty7DT/RGMIBXcrEf3qyxvvpm6VYdybH7Sx2lsuY53/Ii2qtf9CdtAHd9rl8EW+brCJQ5t2lDs3ewLl+9z6vR2/tY96r+rNqrz6dsh28kT/kOLQ89txxOut273NXo3fa9n2tDvX+dbkHZhrm2GOJm+16b0+/VoaR1pboVZ/dPBIrAGHmybuEdeRX8dZr+nYcd71vGHlkZn6Zf8IXLdd2cCbayxsItO7W05nbp+uUfo1t62xbq6B/QBjk6Dm7fp1sfe5bp62571j3FhU/qlpvdGhr/TR2RNftHvs/iUdfDjO/pSBibITO9rRbQ3nWgVshhSkgevYjMFZHNIlIgIrcGWD9ORL4UkUYR+XnXVzMElj8J84739d5tjb6lDUPf0uiNtwftTbQ0wqd/0GmO/UcL2gY928o3YjeVna0A/xaBf8pZp3fSldgvH/AO5jrrIbj5K+11xCZ7Q9qO+z/fTqnYZG2o7zzgfUklDfT2QUz9tpYIZn3fkpemackhbwbctAKunq+9xkuf1R6qLW1Ba2nCNvT+nd4X/Vv/Abx4hR6Svuld3d8AvpEW076tPeKJF+rrV7lHdzRHxsBxP9Vlkm1Jbr/Xo3dq3ZO+BSfdCd+cpzsNnYbe9oxbmmDl097lhSu8xzm4BQ4WeKUy+55wHttJUhZcM18b4IyRwTvv7BZebBLMuk7LEbZmbJOQqY3dFW/oEZgNlfC/3/iGxyYO1C2N6Hh9fade7jWYg6a2rmtMgl6ff5RuIfY0cSk6+scmMVO/ZP0dMRFfI29ol3YNvYhEAo8BZwDjgctEZLxfsTLgZuBPXV7DUKkt0T3t9iw84A2dbCs5kqvRtzfdbjZ+cp8enl1a4BuxUV+uQ90yx0BStjd0zmnc/XOAlG3Tnq9NICPQFQw7Vn8OP9FrGCKjfL3G/KPhzmIYeoxvPL4tNURG65cD6FbBKb/Vn0NmBT9u5mi9PydOjXrkN/TL7fx/aANZXaTlHH/teuIFOiIjMlaP0sybCbfugus/1/0YTh15wAi4fpEOfUsZrIeY7/1KG3/bCNj67Pu/9F4f50v2jD/Aib/Qv1XmGG3c7WgTZwbR8h3el8fKp3Sr7URrBOWjR+qw2DGn+95HR4c4WK09TvwFzP196+UREXD2Q97+nCmX65GpqkVHv4B+sQTDmWulreADQ78gFI9+FlCglNqulGoCXgB8gsyVUsVKqeVAB2fm7kJsecY5yYht9Nv06Jv8PPoZvuvnHQtv/9j7vb5ce0kiWjsschh6O42uv0dfWqA9uBNvhYkXdV8c7YTz4dfl8N232i5nG6SoGK1zjjzZd7TepIvg5tVaiz/uJ3DDl4dW56SB8OM1MOVSbysimHwVGeU16PnHeI/b1tSBaUN1p3PRKl9NNi1fx3gfWKfD8STCG40x9Tu+xi5zNKC8g39sjz46Ucsmx/2f1vFt2cY54nHu/d7OwO++o+Op7djrnuIbd+j7eNBUXdfMsd5zDcbZf9afvRj2Z+gZQtHocwFnusBCoFPil4hcB1wHkJ/fxaFItjF3Doyxo22UO/g8o/4evX8vfUuT1ixB67gb3oQMS5bIPVKHfzVUav01c7TuiFz3io4lvvhpfczSAl32pNu65lzboqM9/6cGSZ/bFQ//jO/pJFDOl4RtQCdfHHy77Ek6GiZvZvAyToYf703/MNShjcckwB0H4IGRWmKLTdUtjZ9vbZ1OYsgs/aJ+5/+0l7z8CS2F3bTCW/+cSbp/IzbV2+kHWkt31qU3SM2Db7+s79+ssXDTsva3mXGNfuEZGaTfE4pVCOTKhTj+3W8jpR5XSs1QSs3IymqjWdkZXAEMvVPGCeTVu1t0Uzcy1hs2ZRtxW74AnSdj/zrdqVlfAd/8m15ua8a7l2gDP9rKpfH1y9rA7fxMy0YVu33188OFsx/SQ86dzLJGbU5tYwKQvCO10W1LLnIyyjGidYxfRERUjFfjt73upIGtWyipebr/Ye8Kbzrf+nLfcvaLZ8AwvfymFTrqw+ko9CYj5rQ/aMcfY+QPC0Lx6AsBR68deUBRkLK9h63HO4fz++j1ja1Hs9kvh6gY3yb/XZU66ubdn3qXzbP07/Th3uHadufbV8/oF8awE3S8uu21rn9De1jKfXga+kCcfi+c8pu2jeO0K7SUFGpfRs5kLauMP09Hj/iTN0NH0BxzY9v7Ofk3OopnwEid3naSX6vDM+LXUiidHc4GQx8mFEO/HBgtIsOBvcClwOVtb9ILBPLoAd0gUYFj6W0vPzKA0QmmSTsljeQcPWhj0zv6OENm6tFxtqFf+yJEW52igUYWHo6ItO8BR0Z7X6ahEBGhwwiDDcSZc5vuoGyvhRAd5x2odntR6/ztdh5xZ3SNwRAGtCvdKKVcwE3AB8BG4CWl1HoRuV5ErgcQkRwRKQR+CtwpIoUikhJ8r91AMENvR7v4pCOu0pruM1afcnvNV6cc4OxEFPGGMebN1GGJdsjglMt1CNzSv+sXifH+upfoOO9Ao1br4jseUx2T2Hp/6UPhyre8nZgGQ5gQ0oAppdR7wHt+y+Y5/t+PlnR6D1u6aWXoE61MgA6P/n+/tob1W10NEUE8wZnX6tzYl7+gh/MXfNi6k9L2Io+1RiHahj5vhvZKF96jOwCDeZuG8CJYlk6DoQ/Tf1IgBAqvBK8u7/ToD27V8eSnWZEawVIknPUg/NAa7GQ34/2jNU7/vU64ZMcup1ndGQOGe/Nt+A+HNxgMhh6k/6RAaMujB9+om5oDWmc9+kY96Gnsme3vf9oVWov3n+8xZ6JvjppRp8DkS3UKgNgk/bLI9xsObzAYDD1I/zH0tlfeyqO3NHrnzFM1xXq0ZkRE6EO9x87VozPbGrgDuoP2gn94v8+8NrT9GwwGQzfRj6SbYB69FfViSzfNDTozX1vDw4PRnpE3GAyGPkg/MvRBom7s8Ebb47dne0/K7pl6GQwGQy/Tjwy95dErv0mo/cMra4yhNxgMhxf9yNAHiZzx74y1p8jr6KQIBoPBEKb0H9HZ9uj9sQ19Y7XOb26nZDUevcFgOEzoH4a+xeUr2SRkeFMF29JN4XI932liFiCOqdgMBoOhf9M/pBt/bz4hw/u/7dGXbteftSV6wg2Ttc9gMBwm9A9D7z+y1ZmkzDb0Zdu9y8xECwaD4TCifxh6f4/emYzKNvQ1+73LumvOVoPBYOiD9FNDH6WnjQNvHL0T49EbDIbDiH5i6P2lm2iIsWaOj4zxnRMWfCfLNhgMhn5OPzH0ATz6WCvaJjLKa+jtkMrMDk63ZjAYDGFM/wivdPnNBxsRCbGWRx8R7U10dtT1eoq67PE9Wz+DwWDoRULy6EVkrohsFpECEbk1wHoRkb9a69eKyPRA++k2bI8+wnpvRUR74+fFcYr5R+uJpw0Gg+Ewol1DLyKRwGPAGcB44DIR8XeJzwBGW3/XAX/v4nq2jR1eaU/w4ZRumh1JzgZP69FqGQwGQ18gFOlmFlCglNoOICIvAOcBGxxlzgOeUUopYImIpInIIKXUvi6vccGH8MEdvssaq/VnXCrUl2npxp4JytlRGx3f5dUxGAyGvk4ohj4X2OP4Xgj4z7QcqEwu4GPoReQ6tMdPfn5+R+uqiU2BrACdqfGnwujTYe0LMPU7es7WlMEw+jQ4/3FIGNC54xkMBkOYE4qhlwDLVCfKoJR6HHgcYMaMGa3Wh8SQWTDkmeDrxzmmBTz1bv055ZJOHcpgMBj6A6F0xhYCQxzf84CiTpQxGAwGQy8QiqFfDowWkeEiEgNcCrzlV+Yt4Eor+uZooLJb9HmDwWAwdJh2pRullEtEbgI+ACKBfyml1ovI9db6ecB7wJlAAVAHXN19VTYYDAZDRwhpwJRS6j20MXcum+f4XwE3dm3VDAaDwdAV9I8UCAaDwWAIijH0BoPB0M8xht5gMBj6OcbQGwwGQz9HdD9qLxxYpATY1cnNM4GDXVid3qK/nAeYc+nL9KfzMecCQ5VSWR3ZoNcM/aEgIiuUUjN6ux6HSn85DzDn0pfpT+djzqVzGOnGYDAY+jnG0BsMBkM/J1wN/eO9XYEuor+cB5hz6cv0p/Mx59IJwlKjNxgMBkPohKtHbzAYDIYQMYbeYDAY+jtKqW7/Q+eq/xjYCKwHfmwtHwD8D9hqfaZby08FVgJfW5/fcOzrSGt5AfBXLPkpwDEDlgNOAL4CXMBFYXwe11vLVwOfA+PD+FyuAkqsc1kNXBvG5/Jnx3lsASrC/HkZCnwErAU+AfLC4FzuRc94V+O3vK88+wHr14Fr0uHz6FKD3kaFBwHTrf+TrQdgPPBH4FZr+a3AH6z/pwGDrf8nAnsd+1oGHIOe1ep94IwgxwxYDhgGTAae6cTF7kvnkeIocy4wP4zP5Srg0f5wf/mV+RE6rXfYng/wMvBd6/9vAP8Jg3M52jquv6EfRt949gPWrwPXpMPn0akH61D/gDfRb7zNwCDHD7k5QFkBSoFYq8wmx7rLgH8EuShtlgOe6ujF7ovn4Vj+frieC4do6PvSufiVWwycGs7ng/Ze8xz7rurL5+K3fUBDSi8++6HUL9R7rCPn0eMavYgMQ7/tlgLZypqJyvocGGCTC4FVSqlG9ITjhY519iTk/oRartP0hfMQkRtFZBvaq7g5nM8FuFBE1orIKyIyhE7SR84FERkKDAcWdupEvPsZRu+ezxprnwDnA8kiktGHz6XbOcTzCJUuPd8eNfQikgS8CvxEKVUVQvkJwB+AH9iLAhRTgTYNsVyn6CvnoZR6TCk1ErgFuLO9egSpW184l7eBYUqpycCHwNPt1SNI3frCudhcCryilGpprx5t1K8vnM/PgRNFZBVwIrAXrQ13iB48l26lC84j5EMFWNbp8+0xQy8i0egf6Fml1GvW4gMiMshaPwgodpTPA14HrlRKbbMWF6InHrfJA4pEJFJEVlt/dwcr14/P4wXgm+F6LkqpUoe38090J1RYnouDS4HnO3oefe18lFJFSqkLlFLTgDusZZV9+Fy6jS46j2D77l4bdig6VQf0LEF3HDzst/wBfDsy/mj9n4bVZAywr+Xozgy7g+LMIMdssxyd0On60nkAox1lzgFWhPG5DHKUOR9YEq7nYq0bC+wkSFRIOJ0POsNihPX/vcDdff1cHOW7TKPvyvNor36hnm9HzqPDN2Enb9zj0M2OtXhDz84EMtChW1utzwFW+TuBWkfZ1cBAa90MYB2wDXg02MMUrBwwE/22rEV3kKwP0/P4C7qjbDU67GtCGF+T+6xzWWOdy7hwPRdr3V3A/f3kebnIOt4W4An8OhT76Ln8Ef2Mu63Pu/rYsx+wft1pw0wKBIPBYOjnmJGxBoPB0M8xht5gMBj6OcbQGwwGQz/HGHqDwWDo5xhDbzAYDP0cY+gNBoOhn2MMvcFgMPRz/h/drnLNu+HSsAAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "DATE_PATTERN = \"%Y%m%d\"\n",
    "\n",
    "for word in [\"de\", \"ik\"]:\n",
    "    x = []\n",
    "    y = []\n",
    "    for date in results[word]:\n",
    "        if results[word][date][0] != 0:\n",
    "            x.append(datetime.datetime.strptime(date, DATE_PATTERN))\n",
    "            y.append(results[word][date][1]/results[word][date][0])\n",
    "    plt.plot_date(x, y, fmt=\"-\", label=word)\n",
    "\n",
    "plt.title(f\"percentage of 'de' and 'ik' in corona|covid tweets \")\n",
    "plt.legend()\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Old code"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "FILENAME202002 = \"query-tokens-202002-202003.csv\"\n",
    "FILENAME202004 = \"query-tokens-202004-202005.csv\"\n",
    "NON = \"non-\"\n",
    "WEEKS = {\"2020020[2-8]\":\"20200202\",\"20200209|2020021[0-5]\":\"20200209\",\"2020021[6-9]|2020022[0-2]\":\"20200216\",\"2020022[3-9]\":\"20200223\",\\\n",
    "         \"2020030[1-7]\":\"20200301\",\"2020030[89]|2020031[0-4]\":\"20200308\",\"2020031[5-9]|2020032[01]\":\"20200315\",\"2020032[2-8]\":\"20200322\",\\\n",
    "                                   \"20200329|20200330|2020040[1-4]\":\"20200329\",\\\n",
    "         \"2020040[5-9]|2020041[01]\":\"20200405\",\"2020041[2-8]\":\"20200412\",\"20200419|2020042[0-5]\":\"20200419\",\"2020042[6-9]|2020043|2020050[1-2]\":\"20200426\",\\\n",
    "         \"2020050[3-9]\":\"20200503\",\"2020051[0-6]\":\"20200510\",\"2020051[7-9]|2020052[0-3]\":\"20200517\"}"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "queryTokens202002 = readData(FILENAME202004)\n",
    "nonQueryTokens202002 = readData(NON+FILENAME202004)\n",
    "tscoresDataQuery = {}\n",
    "tscoresDataNonQuery = {}\n",
    "tscores = {}\n",
    "for query in queryTokens202002:\n",
    "    tscoresDataQuery[query] = makeTscoreData(combineInitialHashAt(queryTokens202002[query]))\n",
    "    tscoresDataNonQuery[query] = makeTscoreData(combineInitialHashAt(nonQueryTokens202002[query]))\n",
    "    print(query)\n",
    "    tscores[query] = sortTscores(tscore.computeTscore(tscoresDataQuery[query],tscoresDataNonQuery[query]))\n",
    "    for x in dictTopN(tscores[query]): print(round(x[0],1),x[1])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "# writeData(tscores,\"tscores-202004-202005.csv\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "tscores202002 = readData(\"tscores-202002-202003.csv\")\n",
    "tscores202004 = readData(\"tscores-202004-202005.csv\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "TOPICS = \"maatregelen mondkapje anderhalve besmet rivm \".split()\n",
    "\n",
    "data = {}\n",
    "for topic in TOPICS:\n",
    "    for date in tscores202002.keys():\n",
    "        if not topic in data: data[topic] = {}\n",
    "        data[topic][date] = tscores202002[date][topic]\n",
    "    for date in tscores202004.keys():\n",
    "        data[topic][date] = tscores202004[date][topic]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "from datetime import datetime\n",
    "import matplotlib.pyplot as plt\n",
    "from matplotlib.dates import DateFormatter\n",
    "\n",
    "DATEFORMAT = \"%Y%m%d\"\n",
    "\n",
    "fig,ax = plt.subplots(figsize=(8,4))\n",
    "ax.xaxis.set_major_formatter(DateFormatter(\"%d %b\"))\n",
    "for topic in TOPICS:\n",
    "    plt.plot_date([datetime.strptime(WEEKS[date],DATEFORMAT) for date in data[topic]],list(data[topic].values()),label=topic,fmt=\"-\")\n",
    "ax.set(xlabel=\"date (weeks)\",ylabel=\"t-scores\")\n",
    "plt.title(\"t-scores for relevant words, comparing selected topic tweets with unselected tweets\")\n",
    "plt.legend()\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "**Note**: We expected only positive scores for topic words but the graph shows that *anderhalve* had negative scores  in two weeks. This means that most of the tweets containing this word are outside our topic tweets. We checked a sample of these missing tweets and most of them turned out to be on-topic, but they lacked the words *corona* and *covid*."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "data[\"anderhalve\"]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "for date in tscores202002: \n",
    "    print(\"#####\",date)\n",
    "    for x in dictTopN(sortTscores(tscores202002[date]),n=50): print(round(x[0],1),x[1])\n",
    "for date in tscores202004:\n",
    "    print(\"#####\",date)\n",
    "    for x in dictTopN(sortTscores(tscores202004[date]),n=50): print(round(x[0],1),x[1])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "queryTokens202004 = readData(FILENAME202004)\n",
    "nonQueryTokens202004 = readData(NON+FILENAME202004)\n",
    "for query in queryTokens202004:\n",
    "    tscoresDataQuery[query] = makeTscoreData(combineInitialHashAt(queryTokens202004[query]))\n",
    "    tscoresDataNonQuery[query] = makeTscoreData(combineInitialHashAt(nonQueryTokens202004[query]))\n",
    "    print(query)\n",
    "    for x in dictTopN(sortTscores(tscore.computeTscore(tscoresDataQuery[query],tscoresDataNonQuery[query]))): print(round(x[0],1),x[1])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "#writeData(queryTokens,\"query-tokens-202004-202005.csv\")\n",
    "#writeData(nonQueryTokens,\"non-query-tokens-202004-202005.csv\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "query = \"2020051[7-9]|2020052[0-3]\"\n",
    "for x in dictTopN(sortTscores(tscore.computeTscore(tscoresDataQuery[query],tscoresDataNonQuery[query])),n=40): print(round(x[0],1),x[1])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "FILEPATTERN = \"20200601\" \n",
    "QUERY= r\"1.5m|1,5m|afstand.*hou|hou.*afstand|anderhalve\"\n",
    "\n",
    "queryTokens = {}\n",
    "nonQueryTokens = {}\n",
    "tscoresDataQuery = {}\n",
    "tscoresDataNonQuery = {}\n",
    "queryTokens[FILEPATTERN],nonQueryTokens[FILEPATTERN] = getTokensOfMatchedTweets(FILEPATTERN,QUERY)\n",
    "tscoresDataQuery[FILEPATTERN] = makeTscoreData(queryTokens[FILEPATTERN])\n",
    "tscoresDataNonQuery[FILEPATTERN] = makeTscoreData(nonQueryTokens[FILEPATTERN])\n",
    "for x in dictTopN(sortTscores(tscore.computeTscore(tscoresDataQuery[FILEPATTERN],tscoresDataNonQuery[FILEPATTERN])),n=40): \n",
    "    print(round(x[0],1),x[1])"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "**Note**: the query \"zorg\" produces many false positives. There are irrelevant types (\"bezorgd\" and \"bezorgen\"), irrelevant syntactical forms (\"ik zorg dat\" and \"Zorg dat je\") and even the correct sense is not always related to the pandemic topic (\"zorg voor ouderen/gehandicapten\")."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Estimate coverage of crawler by looking back for messages with replies"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "DATEPATTERN = \"20200522\"\n",
    "\n",
    "fileList = sorted(os.listdir(DATADIR))\n",
    "seenIds = {}\n",
    "for inFileName in fileList:\n",
    "    if re.search(DATEPATTERN,inFileName):\n",
    "        df = pd.read_csv(DATADIR+inFileName,compression=\"gzip\",dtype=str)\n",
    "        known = 0\n",
    "        unknown = 0\n",
    "        for i in range(0,len(df)):\n",
    "            idstr = df.iloc[i][ID]\n",
    "            seenIds[idstr] = True\n",
    "            replyParent = df.iloc[i][REPLYID]\n",
    "            if str(replyParent) != \"nan\":\n",
    "                if replyParent in seenIds: known += 1\n",
    "                else: unknown += 1\n",
    "        print(\"{0} {1:.3f}\".format(inFileName,round(known/(known+unknown),3)))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "python37",
   "language": "python",
   "name": "python37"
  },
  "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.3"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}