{ "nbformat": 4, "nbformat_minor": 0, "metadata": { "colab": { "name": "7) The Impact of Legalizing Prostitution on Crime.ipynb", "provenance": [] }, "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.7.9" } }, "cells": [ { "cell_type": "markdown", "metadata": { "id": "mfxvKsxwpZjN" }, "source": [ "# 7) The Impact of Legalizing Prostitution on Crime" ] }, { "cell_type": "markdown", "metadata": { "id": "kRSC3oJZQSQ2" }, "source": [ "[Vitor Kamada](https://www.linkedin.com/in/vitor-kamada-1b73a078)\n", "\n", "E-mail: econometrics.methods@gmail.com\n", "\n", "Last updated: 11-2-2020" ] }, { "cell_type": "markdown", "metadata": { "id": "iP3_k0wEcYSY" }, "source": [ "In the Netherlands, there are legal prostitution zones that in Dutch is called tippelzones. Bisschop et al. (2017) reported that the opening of a tippelzone reduces sexual abuse and rape by about 30-40%.\n", "\n", "Let's open the dataset of Bisschop et al. (2017). Each row is a city in the Netherlands. The same city is observed between 1994–2011." ] }, { "cell_type": "code", "metadata": { "id": "3wG1YQm6B55d", "outputId": "fa1285da-5df1-413e-9ca7-41ed341635c1", "colab": { "base_uri": "https://localhost:8080/", "height": 309 } }, "source": [ "import numpy as np\n", "import pandas as pd\n", "pd.set_option('precision', 3)\n", "\n", "# Data from Bisschop et al. (2017)\n", "path = \"https://github.com/causal-methods/Data/raw/master/\" \n", "df = pd.read_stata(path + \"CBSregist2015.dta\")\n", "df.head(5)" ], "execution_count": 1, "outputs": [ { "output_type": "execute_result", "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
cityyearopeningclosingcity1logpopdensopeningRegmayorCDAmayorCUmayorD66mayorVVDlnmaltreatNlnweaponsNlndrugsNlnrapeNlnsexassaultNlnrapesexaNopening12mmayorLibmayorChrmayorSocinkhhinsurWWAOlogpopmale1565nondutchpceduchpcmalfemratioopeningNoReg12openingReg12openingRegP12mopeningRegPeveropenNoRegeveropenNoReg38popeningRegAopeningReg38pAclosespillNoRegregspillNoRegbrothelbanopenspillRegAopenspillNoRegopenspillRegPrapepcNsexassaultpcNrapesexapcNdrugspcNmaltreatpcNweaponspcNpopul_100popmale1565_100pop_dens_100insurWWAO_pcrobbembpcNlnrobbembNrobbpcNlnrobbNsimtheftpcNaggtheftpcNlnsimtheftNlnaggtheftNtheftpcNlntheftNpordoffpcNpordviolpcNlnpordoffNlnpordviolN
0Amsterdam1994-01-010.00.01.08.3810.00.00.00.00.06.7866.5077.4734.8284.1745.2471.00.00.01.029.069650.012.4590.1610.2650.9560.00.00.00.00.00.00.00.00.00.00.00.01.00.01.7260.8982.62424.30612.2229.2537240.962576.0743.640.096132.1659.16612.6366.81959.24657.9348.3648.342117.1819.0467.5965.1106.3105.914
1Amsterdam1995-01-010.00.01.08.3790.00.00.00.00.06.8466.1847.2234.5003.9124.9421.00.00.01.029.070165.012.4630.1530.2810.9590.00.00.00.00.00.00.00.00.00.00.00.01.00.01.2460.6921.93818.96913.0156.7157222.302585.7843.530.097109.4538.97512.9466.84150.81543.8238.2088.06094.6378.8307.0614.3616.2345.753
2Amsterdam1996-01-011.00.01.08.3730.00.00.00.00.06.8516.2447.1114.4433.6894.8280.00.00.01.029.069730.012.4680.1550.3140.9651.00.00.00.01.00.00.00.00.00.00.00.01.00.01.1840.5571.74117.05813.1597.1727181.192599.2143.280.09793.6478.81412.6026.80842.33337.1118.0207.88879.4448.6497.5205.4316.2925.966
3Amsterdam1997-01-011.00.01.08.3690.00.00.00.00.06.8625.9917.1194.0073.4014.4430.00.00.01.029.066790.012.4610.1340.3150.9611.00.00.00.01.00.00.00.00.00.00.00.01.00.00.7690.4191.18917.26913.3545.5937151.482581.2343.100.09393.0588.80311.7466.73346.84332.8608.1177.76279.7048.6486.8524.1956.1945.704
4Amsterdam1998-01-011.00.01.08.3730.00.00.00.00.06.9235.7847.1934.2483.8074.7450.00.00.01.029.063260.012.4680.1400.3200.9620.00.00.00.01.01.00.00.00.00.00.00.01.00.00.9750.6271.60118.52014.1344.5267181.512598.4143.280.08894.0618.81813.0206.84145.25533.9078.0867.79879.1628.6466.1274.5956.0875.799
\n", "
" ], "text/plain": [ " city year opening ... pordviolpcN lnpordoffN lnpordviolN\n", "0 Amsterdam 1994-01-01 0.0 ... 5.110 6.310 5.914\n", "1 Amsterdam 1995-01-01 0.0 ... 4.361 6.234 5.753\n", "2 Amsterdam 1996-01-01 1.0 ... 5.431 6.292 5.966\n", "3 Amsterdam 1997-01-01 1.0 ... 4.195 6.194 5.704\n", "4 Amsterdam 1998-01-01 1.0 ... 4.595 6.087 5.799\n", "\n", "[5 rows x 65 columns]" ] }, "metadata": { "tags": [] }, "execution_count": 1 } ] }, { "cell_type": "markdown", "metadata": { "id": "AWoRmkPF3BMW" }, "source": [ "Let's break the cities into 3 groups. Big cities and medium cities had a tippelzone in at least one year, whereas the other cities in the sample hadn't a tippelzone." ] }, { "cell_type": "code", "metadata": { "id": "Q8FitWN-Di6_" }, "source": [ "big_cities = [\"Amsterdam\", \"Rotterdam\", \"Den Haag\"]\n", "medium_cities = [\"Utrecht\", \"Nijmegen\", \"Groningen\",\n", " \"Heerlen\", \"Eindhoven\", \"Arnhem\"]\n", "\n", "# Classify cities\n", "def classify(var):\n", " if var in big_cities:\n", " return \"Big Cities\"\n", " elif var in medium_cities:\n", " return \"Medium Cities\" \n", " else: \n", " return \"No tippelzone\"\n", "\n", "df['group'] = df[\"city\"].apply(classify)" ], "execution_count": 2, "outputs": [] }, { "cell_type": "markdown", "metadata": { "id": "DAHRi2vD4dW8" }, "source": [ "Below are the annual crime reports per 10,000 inhabitants. Overall, big cities have higher crime rates. The exception is a drug-related crime. " ] }, { "cell_type": "code", "metadata": { "id": "d1hJyks6IO_9", "outputId": "fe717e35-a57c-442a-f9e0-20532405b814", "colab": { "base_uri": "https://localhost:8080/", "height": 202 } }, "source": [ "outcome = ['sexassaultpcN', 'rapepcN', \n", "\t'drugspcN', 'maltreatpcN', 'weaponspcN']\n", "\n", "df.groupby('group')[outcome].mean().T" ], "execution_count": 3, "outputs": [ { "output_type": "execute_result", "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
groupBig CitiesMedium CitiesNo tippelzone
sexassaultpcN0.7750.6260.664
rapepcN1.0320.8460.691
drugspcN14.92115.59912.779
maltreatpcN21.25918.66517.864
weaponspcN5.6354.3854.207
\n", "
" ], "text/plain": [ "group Big Cities Medium Cities No tippelzone\n", "sexassaultpcN 0.775 0.626 0.664\n", "rapepcN 1.032 0.846 0.691\n", "drugspcN 14.921 15.599 12.779\n", "maltreatpcN 21.259 18.665 17.864\n", "weaponspcN 5.635 4.385 4.207" ] }, "metadata": { "tags": [] }, "execution_count": 3 } ] }, { "cell_type": "markdown", "metadata": { "id": "hqQfUR3X6h7Z" }, "source": [ "Tippelzones cities are more populated and densely populated than cities without tippelzones. Average household income (\"inkhh\") measured in 1,000 euros is similar in the 3 groups. Tippelzones cities also have individuals with a higher level of education. The proportion of immigrants is higher in big cities (11.4%). The share of social insurance benefits (\"insurWWAO_pc\") is similar to the 3 groups." ] }, { "cell_type": "code", "metadata": { "id": "BZ9l-Y7yKQx0", "outputId": "6a3ef5cb-3d85-4f11-c20d-259fbcf6e4f1", "colab": { "base_uri": "https://localhost:8080/", "height": 263 } }, "source": [ "demographics = ['popul_100', 'pop_dens_100', 'popmale1565_100',\n", " 'inkhh', 'educhpc', 'nondutchpc', 'insurWWAO_pc']\n", "\n", "df.groupby('group')[demographics].mean().T " ], "execution_count": 4, "outputs": [ { "output_type": "execute_result", "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
groupBig CitiesMedium CitiesNo tippelzone
popul_1005974.8861724.1911131.138
pop_dens_10043.25822.97719.560
popmale1565_1002101.446617.019392.255
inkhh29.05228.98930.502
educhpc0.3000.3170.245
nondutchpc0.1140.0590.052
insurWWAO_pc0.0740.0810.078
\n", "
" ], "text/plain": [ "group Big Cities Medium Cities No tippelzone\n", "popul_100 5974.886 1724.191 1131.138\n", "pop_dens_100 43.258 22.977 19.560\n", "popmale1565_100 2101.446 617.019 392.255\n", "inkhh 29.052 28.989 30.502\n", "educhpc 0.300 0.317 0.245\n", "nondutchpc 0.114 0.059 0.052\n", "insurWWAO_pc 0.074 0.081 0.078" ] }, "metadata": { "tags": [] }, "execution_count": 4 } ] }, { "cell_type": "markdown", "metadata": { "id": "iwmQm89XBL_V" }, "source": [ "The Christian Union party has more mayors (31%) in cities without tippelzone. It is worth mentioning that this party opposes the opening of the tippelzone." ] }, { "cell_type": "code", "metadata": { "id": "osAzJJntM8fE", "outputId": "d1561569-d704-4b60-fcf6-6e6465db53e4", "colab": { "base_uri": "https://localhost:8080/", "height": 141 } }, "source": [ "political_party = ['mayorSoc', 'mayorLib', 'mayorChr']\n", "df.groupby('group')[political_party].mean().T " ], "execution_count": 5, "outputs": [ { "output_type": "execute_result", "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
groupBig CitiesMedium CitiesNo tippelzone
mayorSoc0.4810.5560.410
mayorLib0.2590.3240.278
mayorChr0.2590.1200.312
\n", "
" ], "text/plain": [ "group Big Cities Medium Cities No tippelzone\n", "mayorSoc 0.481 0.556 0.410\n", "mayorLib 0.259 0.324 0.278\n", "mayorChr 0.259 0.120 0.312" ] }, "metadata": { "tags": [] }, "execution_count": 5 } ] }, { "cell_type": "markdown", "metadata": { "id": "0NSIJJNb12a2" }, "source": [ "The dataset is a balanced Panel Data. It is necessary to declare the indices in order: the unit of analysis and time unity." ] }, { "cell_type": "code", "metadata": { "id": "fw285SLhc9tN" }, "source": [ "df['year'] = pd.DatetimeIndex(df['year']).year\n", "df['Dyear'] = pd.Categorical(df.year)\n", "\n", "# Set Panel Data\n", "# Set city as the unit of analysis\n", "df25 = df.set_index(['city1', 'year'])" ], "execution_count": 6, "outputs": [] }, { "cell_type": "markdown", "metadata": { "id": "HLjOUaocpP2i" }, "source": [ "Let $Y_{ct}$ be the crime in city $c$\n", "in year $t$. Let $D_{ct}$ = 1 if city $c$ had an open tippelzone in year $t$; and 0 otherwise. Let's estimate the following model:\n", "\n", "$$ln(Y_{ct})=\\alpha_c+\\rho D_{ct}+\\beta X_{ct}+\\gamma_t + \\epsilon_{ct}$$\n", "\n", "where $\\alpha_c$ is the city fixed effects, $X_{ct}$ is the vector of control variables, $\\gamma_t$ is the year fixed effects, and $\\epsilon_{ct}$ the usual error term. " ] }, { "cell_type": "code", "metadata": { "id": "r0CEG--xh1no", "outputId": "cdde9d3a-c74a-41a7-9ad3-63d359d47788", "colab": { "base_uri": "https://localhost:8080/" } }, "source": [ "import statsmodels.formula.api as smf\n", "\n", "Ys = [\"lnsexassaultN\", \"lnrapeN\", \"lndrugsN\",\n", " \"lnweaponsN\", \"lnmaltreatN\"]\n", "\n", "base = \"~ 1 + opening\"\n", "fe = \"+ C(city) + C(Dyear)\"\n", "\n", "controls = ['logpopmale1565', 'logpopdens', 'inkhh', \n", "\t'educhpc', 'nondutchpc', 'insurWWAO', 'mayorCDA',\n", " 'mayorCU', 'mayorD66', 'mayorVVD']\n", "\n", "Xs = \"\"\n", "for var in controls:\n", " Xs = Xs + '+' + var\n", "\n", "columns = []\n", "for Y in Ys:\n", " result = smf.ols(Y + base + fe + Xs, df25).fit(cov_type='cluster',\n", " cov_kwds={'groups': df25['city']})\n", " columns.append(result)" ], "execution_count": 7, "outputs": [ { "output_type": "stream", "text": [ "/usr/local/lib/python3.6/dist-packages/statsmodels/tools/_testing.py:19: FutureWarning: pandas.util.testing is deprecated. Use the functions in the public API at pandas.testing instead.\n", " import pandas.util.testing as tm\n" ], "name": "stderr" } ] }, { "cell_type": "code", "metadata": { "id": "PUbxgAwQ-lsb", "outputId": "d4116195-4f14-4297-8c86-e98986dd3824", "colab": { "base_uri": "https://localhost:8080/" } }, "source": [ "# Library to print professional publication\n", "# tables in Latex, HTML, etc.\n", "!pip install stargazer" ], "execution_count": 8, "outputs": [ { "output_type": "stream", "text": [ "Requirement already satisfied: stargazer in /usr/local/lib/python3.6/dist-packages (0.0.5)\n" ], "name": "stdout" } ] }, { "cell_type": "markdown", "metadata": { "id": "RgjGIUuYGv37" }, "source": [ "Column 1 indicates that an open tippelzone decreases sexual abuse by 26% ($e^{-0.302}-1$). In the other columns, the coefficients of tippelzone are not statistically significant. It looks that legalizing prostitution decreases sexual abuse, but not other crimes as rape, assault, illegal weapon, and drug-related. " ] }, { "cell_type": "code", "metadata": { "id": "wq2Qmjyn-sU5", "outputId": "b18dbc49-93c7-4442-a100-0ca7ea9a2b00", "colab": { "base_uri": "https://localhost:8080/", "height": 535 } }, "source": [ "# Settings for a nice table\n", "from stargazer.stargazer import Stargazer\n", "stargazer = Stargazer(columns)\n", "\n", "stargazer.title('The Impact of Tippelzone on Crime')\n", "\n", "names = ['Sex Abuse', 'Rape', 'Drugs', 'Weapons', 'Assault']\n", "stargazer.custom_columns(names, [1, 1, 1, 1, 1])\n", "\n", "stargazer.covariate_order(['opening'])\n", "\n", "stargazer.add_line('Covariates', ['Yes', 'Yes', 'Yes', 'Yes', 'Yes'])\n", "\n", "stargazer.add_line('City Fixed Effects', ['Yes', 'Yes', 'Yes', 'Yes', 'Yes'])\n", "stargazer.add_line('Year Fixed Effects', ['Yes', 'Yes', 'Yes', 'Yes', 'Yes'])\n", "\n", "stargazer" ], "execution_count": 9, "outputs": [ { "output_type": "stream", "text": [ "/usr/local/lib/python3.6/dist-packages/statsmodels/base/model.py:1752: ValueWarning: covariance of constraints does not have full rank. The number of constraints is 52, but rank is 24\n", " 'rank is %d' % (J, J_), ValueWarning)\n", "/usr/local/lib/python3.6/dist-packages/statsmodels/base/model.py:1752: ValueWarning: covariance of constraints does not have full rank. The number of constraints is 52, but rank is 24\n", " 'rank is %d' % (J, J_), ValueWarning)\n", "/usr/local/lib/python3.6/dist-packages/statsmodels/base/model.py:1752: ValueWarning: covariance of constraints does not have full rank. The number of constraints is 52, but rank is 24\n", " 'rank is %d' % (J, J_), ValueWarning)\n", "/usr/local/lib/python3.6/dist-packages/statsmodels/base/model.py:1752: ValueWarning: covariance of constraints does not have full rank. The number of constraints is 52, but rank is 24\n", " 'rank is %d' % (J, J_), ValueWarning)\n", "/usr/local/lib/python3.6/dist-packages/statsmodels/base/model.py:1752: ValueWarning: covariance of constraints does not have full rank. The number of constraints is 52, but rank is 24\n", " 'rank is %d' % (J, J_), ValueWarning)\n" ], "name": "stderr" }, { "output_type": "execute_result", "data": { "text/html": [ "The Impact of Tippelzone on Crime
\n", "
Sex AbuseRapeDrugsWeaponsAssault
(1)(2)(3)(4)(5)
opening-0.302***-0.042-0.057-0.074-0.017
(0.091)(0.060)(0.076)(0.093)(0.060)
CovariatesYesYesYesYesYes
City Fixed EffectsYesYesYesYesYes
Year Fixed EffectsYesYesYesYesYes
Observations450450450450450
R20.7270.7940.8990.9070.957
Adjusted R20.6920.7670.8860.8950.952
Residual Std. Error0.485 (df=397)0.449 (df=397)0.330 (df=397)0.289 (df=397)0.165 (df=397)
F Statistic26.328*** (df=52; 397)65.173*** (df=52; 397)45.273*** (df=52; 397)10.007*** (df=52; 397)252.545*** (df=52; 397)
Note:\n", " *p<0.1;\n", " **p<0.05;\n", " ***p<0.01\n", "
" ], "text/plain": [ "" ] }, "metadata": { "tags": [] }, "execution_count": 9 } ] }, { "cell_type": "code", "metadata": { "id": "tiIWOkWXHz4o", "outputId": "e0b469a1-e110-4654-984b-10473c2dbc2b", "colab": { "base_uri": "https://localhost:8080/" } }, "source": [ "import math\n", "math.exp(-0.302) - 1 " ], "execution_count": 10, "outputs": [ { "output_type": "execute_result", "data": { "text/plain": [ "-0.2606619351104681" ] }, "metadata": { "tags": [] }, "execution_count": 10 } ] }, { "cell_type": "markdown", "metadata": { "id": "_LhtUSOZE5s3" }, "source": [ "## Exercises" ] }, { "cell_type": "markdown", "metadata": { "id": "tx6eyoPl3yWU" }, "source": [ "1| In the introduction of the paper, Bisschop et al. (2017: 29) state: \"Our study is one of the first to provide causal evidence for the connection between the regulation of prostitution and crime.\" In the discussion section, Bisschop et al. (2017:44) state: \"The opening of a tippelzone with or without a licensing system is correlated with a short-run decrease of 30-40 percent in sexual abuse and rape, and the results are robust to different specifications.\" Why do Bisschop et al. (2017) use the word \"causal\" in the introduction and \"correlated\" in the discussion section? Do you consider the main results of Bisschop et al. (2017) as \"causal\" or \"correlation\"? Justify." ] }, { "cell_type": "markdown", "metadata": { "id": "MCiUud-ijFP9" }, "source": [ "2| Bisschop et al. (2017: 29) state: \"We perform several empirical tests to assess endogenous crime trends around the moment of opening a tippelzone.\" Why do they do that? What is the logic of that? Are there or not endogenous crime trends? Justify and be specific in your answer. " ] }, { "cell_type": "markdown", "metadata": { "id": "BkDjoB43t15x" }, "source": [ "3| Bisschop et al. (2017: 36) state: \"... the time trend $\\mu_t$ is modeled using year fixed effects\". What are the other ways to model the time trend? Write different pieces of code to create time trends under different assumptions. Hint: Remember that this is a Panel Data. A proper code that works in Cross-Sectional Data will create the wrong variable in a Panel Data Structure." ] }, { "cell_type": "markdown", "metadata": { "id": "jU9MHjsxy6K2" }, "source": [ "4| Bisschop et al. (2017: 36) state: \"We use a difference-in-difference specification to study the effect of the presence of a tippelzone on various types of crime.\" What is the key assumption to deploy a difference-in-difference estimator?" ] }, { "cell_type": "markdown", "metadata": { "id": "77lK90TzNuVg" }, "source": [ "5| Replicate the table \"The Impact of Tippelzone on Crime\" without Amsterdam, Rotterdam, and Den Haag. Additionally, replace the variable \"opening\" with the following four variables:\n", "\n", "i) \"everopenNoReg\": 1 if city $c$ ever has open tippelzone in or before year $t$ which opened without licensing, and 0 otherwise.\n", "\n", "ii) \"openingRegP\": 1 if city $c$ has open tippelzone and introduced licensing ex-post in or before year $t$, and 0 otherwise.\n", "\n", "iii) \"openingRegA\": 1 if city $c$ ever has open tippelzone in or before year $t$ which opened with licensing, and 0 otherwise.\n", "\n", "iv) \"closing\": 1 if city $c$ closes tippelzone in or before year $t$, and 0 otherwise.\n", "\n", "Interpret the results." ] }, { "cell_type": "markdown", "metadata": { "id": "i-MPIny2O8Om" }, "source": [ "## Reference" ] }, { "cell_type": "markdown", "metadata": { "id": "ZOwX1_TNpbMt" }, "source": [ "Bisschop, Paul, Stephen Kastoryano, and Bas van der Klaauw. (2017). [Street Prostitution Zones and Crime](https://github.com/causal-methods/Papers/raw/master/Street%20Prostitution%20Zones%20and%20Crime.pdf). American Economic Journal: Economic Policy, 9 (4): 28-63.\n", "\n" ] } ] }