Gill-Wiehl_et_al_sensitivity_Adopt_KPT.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>protocol_type</th>\n",
       "      <th>protocol_version</th>\n",
       "      <th>protocol_id</th>\n",
       "      <th>monitoring_period</th>\n",
       "      <th>stove_year_cohort</th>\n",
       "      <th>total_verified_credits</th>\n",
       "      <th>type</th>\n",
       "      <th>verified_credits</th>\n",
       "      <th>verified_credits_per_stove</th>\n",
       "      <th>verified_credits_per_stove_isolating</th>\n",
       "      <th>...</th>\n",
       "      <th>0.14</th>\n",
       "      <th>0.15</th>\n",
       "      <th>0.16</th>\n",
       "      <th>0.17</th>\n",
       "      <th>0.18</th>\n",
       "      <th>0.19</th>\n",
       "      <th>0.20</th>\n",
       "      <th>0.21</th>\n",
       "      <th>0.22</th>\n",
       "      <th>0.23</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>Units</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>tCO2e</td>\n",
       "      <td>0</td>\n",
       "      <td>tCO2e</td>\n",
       "      <td>tCO2e/stove/day</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>...</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>CDM AMS II G</td>\n",
       "      <td>AMS-II.G. ver. 3</td>\n",
       "      <td>GS10974</td>\n",
       "      <td>2</td>\n",
       "      <td>domestic</td>\n",
       "      <td>69556</td>\n",
       "      <td>firewood</td>\n",
       "      <td>16258.34792</td>\n",
       "      <td>0.001856</td>\n",
       "      <td>0.002730</td>\n",
       "      <td>...</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>CDM AMS II G</td>\n",
       "      <td>AMS-II.G. ver. 3</td>\n",
       "      <td>GS10974</td>\n",
       "      <td>2</td>\n",
       "      <td>commercial</td>\n",
       "      <td>69556</td>\n",
       "      <td>firewood</td>\n",
       "      <td>8463.65208</td>\n",
       "      <td>0.044424</td>\n",
       "      <td>0.065329</td>\n",
       "      <td>...</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>CDM AMS II G</td>\n",
       "      <td>AMS-II.G. ver. 3</td>\n",
       "      <td>GS10974</td>\n",
       "      <td>3</td>\n",
       "      <td>domestic</td>\n",
       "      <td>69556</td>\n",
       "      <td>firewood</td>\n",
       "      <td>29387.015158</td>\n",
       "      <td>0.001801</td>\n",
       "      <td>0.002648</td>\n",
       "      <td>...</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>CDM AMS II G</td>\n",
       "      <td>AMS-II.G. ver. 3</td>\n",
       "      <td>GS10974</td>\n",
       "      <td>3</td>\n",
       "      <td>commercial</td>\n",
       "      <td>69556</td>\n",
       "      <td>firewood</td>\n",
       "      <td>8294.984842</td>\n",
       "      <td>0.04365</td>\n",
       "      <td>0.064191</td>\n",
       "      <td>...</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>147</th>\n",
       "      <td>CDM AMS II G</td>\n",
       "      <td>AMS-II.G. Version 11.1</td>\n",
       "      <td>GS 2077</td>\n",
       "      <td>MR 2018-2020</td>\n",
       "      <td>2018 vintage</td>\n",
       "      <td>975</td>\n",
       "      <td>firewood</td>\n",
       "      <td>175</td>\n",
       "      <td>0.008561</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>...</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>148</th>\n",
       "      <td>CDM AMS II G</td>\n",
       "      <td>AMS-II.G. Version 1</td>\n",
       "      <td>GS 411</td>\n",
       "      <td>MR5</td>\n",
       "      <td>0</td>\n",
       "      <td>97297</td>\n",
       "      <td>firewood</td>\n",
       "      <td>3016</td>\n",
       "      <td>0.006747</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>...</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>149</th>\n",
       "      <td>CDM AMS II G</td>\n",
       "      <td>AMS-II.G. Version 12</td>\n",
       "      <td>GS 10914</td>\n",
       "      <td>MP2</td>\n",
       "      <td>0</td>\n",
       "      <td>61013</td>\n",
       "      <td>charcoal</td>\n",
       "      <td>18752</td>\n",
       "      <td>0.006303</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>...</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>150</th>\n",
       "      <td>GS Simplified</td>\n",
       "      <td>The Gold Standard Simplified Methodology for E...</td>\n",
       "      <td>GS 10781</td>\n",
       "      <td>0</td>\n",
       "      <td>MP6</td>\n",
       "      <td>2344</td>\n",
       "      <td>firewood</td>\n",
       "      <td>2345</td>\n",
       "      <td>0.0117</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>...</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>151</th>\n",
       "      <td>GS Simplified</td>\n",
       "      <td>The Gold Standard Simplified Methodology for E...</td>\n",
       "      <td>GS 5660</td>\n",
       "      <td>0</td>\n",
       "      <td>MR1</td>\n",
       "      <td>3397</td>\n",
       "      <td>firewood</td>\n",
       "      <td>2771</td>\n",
       "      <td>0.011774</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>...</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>152 rows × 355 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "     protocol_type                                   protocol_version  \\\n",
       "0            Units                                                  0   \n",
       "1     CDM AMS II G                                   AMS-II.G. ver. 3   \n",
       "2     CDM AMS II G                                   AMS-II.G. ver. 3   \n",
       "3     CDM AMS II G                                   AMS-II.G. ver. 3   \n",
       "4     CDM AMS II G                                   AMS-II.G. ver. 3   \n",
       "..             ...                                                ...   \n",
       "147   CDM AMS II G                             AMS-II.G. Version 11.1   \n",
       "148   CDM AMS II G                                AMS-II.G. Version 1   \n",
       "149   CDM AMS II G                               AMS-II.G. Version 12   \n",
       "150  GS Simplified  The Gold Standard Simplified Methodology for E...   \n",
       "151  GS Simplified  The Gold Standard Simplified Methodology for E...   \n",
       "\n",
       "    protocol_id monitoring_period stove_year_cohort total_verified_credits  \\\n",
       "0             0                 0                 0                  tCO2e   \n",
       "1       GS10974                 2          domestic                  69556   \n",
       "2       GS10974                 2        commercial                  69556   \n",
       "3       GS10974                 3          domestic                  69556   \n",
       "4       GS10974                 3        commercial                  69556   \n",
       "..          ...               ...               ...                    ...   \n",
       "147     GS 2077     MR 2018-2020       2018 vintage                    975   \n",
       "148      GS 411               MR5                 0                  97297   \n",
       "149    GS 10914               MP2                 0                  61013   \n",
       "150    GS 10781                 0               MP6                   2344   \n",
       "151     GS 5660                 0               MR1                   3397   \n",
       "\n",
       "         type verified_credits verified_credits_per_stove  \\\n",
       "0           0            tCO2e            tCO2e/stove/day   \n",
       "1    firewood      16258.34792                   0.001856   \n",
       "2    firewood       8463.65208                   0.044424   \n",
       "3    firewood     29387.015158                   0.001801   \n",
       "4    firewood      8294.984842                    0.04365   \n",
       "..        ...              ...                        ...   \n",
       "147  firewood              175                   0.008561   \n",
       "148  firewood             3016                   0.006747   \n",
       "149  charcoal            18752                   0.006303   \n",
       "150  firewood             2345                     0.0117   \n",
       "151  firewood             2771                   0.011774   \n",
       "\n",
       "     verified_credits_per_stove_isolating  ... 0.14 0.15 0.16 0.17 0.18 0.19  \\\n",
       "0                                0.000000  ...    0    0    0    0    0    0   \n",
       "1                                0.002730  ...    0    0    0    0    0    0   \n",
       "2                                0.065329  ...    0    0    0    0    0    0   \n",
       "3                                0.002648  ...    0    0    0    0    0    0   \n",
       "4                                0.064191  ...    0    0    0    0    0    0   \n",
       "..                                    ...  ...  ...  ...  ...  ...  ...  ...   \n",
       "147                              0.000000  ...    0    0    0    0    0    0   \n",
       "148                              0.000000  ...    0    0    0    0    0    0   \n",
       "149                              0.000000  ...    0    0    0    0    0    0   \n",
       "150                              0.000000  ...    0    0    0    0    0    0   \n",
       "151                              0.000000  ...    0    0    0    0    0    0   \n",
       "\n",
       "    0.20 0.21 0.22 0.23  \n",
       "0      0    0    0    0  \n",
       "1      0    0    0    0  \n",
       "2      0    0    0    0  \n",
       "3      0    0    0    0  \n",
       "4      0    0    0    0  \n",
       "..   ...  ...  ...  ...  \n",
       "147    0    0    0    0  \n",
       "148    0    0    0    0  \n",
       "149    0    0    0    0  \n",
       "150    0    0    0    0  \n",
       "151    0    0    0    0  \n",
       "\n",
       "[152 rows x 355 columns]"
      ]
     },
     "execution_count": 1,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "import pandas as pd\n",
    "df = pd.read_excel('9_27_23_python_inputs_charc_firewood_baseline_adjusted_gillwiehl_et_al.xlsx')\n",
    "df"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "\n",
    "VERs_all=pd.read_excel('VERS.xlsx')\n",
    "CDM_I_E_VERs=pd.read_excel('CDM_I_E_VERs.xlsx')\n",
    "CDM_II_G_VERs=pd.read_excel('CDM_II_G_VERs.xlsx')\n",
    "GS_VERs=pd.read_excel('GS_VERs.xlsx')\n",
    "GS_Meter_VERs=pd.read_excel('GS_Meter_VERs.xlsx')\n",
    "\n",
    "\n",
    "bangladesh_VERs=df.loc[df['country']=='Bangladesh','verified_credits_per_stove'].sum()\n",
    "burkinafaso_VERs=df.loc[df['country']=='Burkina Faso','verified_credits_per_stove'].sum()\n",
    "china_VERs=df.loc[df['country']=='China','verified_credits_per_stove'].sum()\n",
    "cameroon_VERs=df.loc[df['country']=='Cameroon','verified_credits_per_stove'].sum()\n",
    "ethiopia_VERs=df.loc[df['country']=='Ethiopia','verified_credits_per_stove'].sum()\n",
    "ghana_VERs=df.loc[df['country']=='Ghana','verified_credits_per_stove'].sum()\n",
    "guatemala_VERs=df.loc[df['country']=='Guatemala','verified_credits_per_stove'].sum()\n",
    "haiti_VERs=df.loc[df['country']=='Haiti','verified_credits_per_stove'].sum()\n",
    "honduras_VERs=df.loc[df['country']=='Honduras','verified_credits_per_stove'].sum()\n",
    "india_VERs=df.loc[df['country']=='India','verified_credits_per_stove'].sum()\n",
    "kenya_VERs=df.loc[df['country']=='Kenya','verified_credits_per_stove'].sum()\n",
    "lesotho_VERs=df.loc[df['country']=='Lesotho','verified_credits_per_stove'].sum()\n",
    "malawi_VERs=df.loc[df['country']=='Malawi','verified_credits_per_stove'].sum()\n",
    "mali_VERs=df.loc[df['country']=='Mali','verified_credits_per_stove'].sum()\n",
    "mexico_VERs=df.loc[df['country']=='Mexico','verified_credits_per_stove'].sum()\n",
    "mozambique_VERs=df.loc[df['country']=='Mozambique','verified_credits_per_stove'].sum()\n",
    "myanmar_VERs=df.loc[df['country']=='Myanmar','verified_credits_per_stove'].sum()\n",
    "nepal_VERs=df.loc[df['country']=='Nepal','verified_credits_per_stove'].sum()\n",
    "nigeria_VERs=df.loc[df['country']=='Nigeria','verified_credits_per_stove'].sum()\n",
    "peru_VERs=df.loc[df['country']=='Peru','verified_credits_per_stove'].sum()\n",
    "rwanda_VERs=df.loc[df['country']=='Rwanda','verified_credits_per_stove'].sum()\n",
    "sudan_VERs=df.loc[df['country']=='Sudan','verified_credits_per_stove'].sum()\n",
    "uganda_VERs=df.loc[df['country']=='Uganda','verified_credits_per_stove'].sum()\n",
    "zambia_VERs=df.loc[df['country']=='Zambia','verified_credits_per_stove'].sum()\n",
    "\n",
    "firewood_VERs=df.loc[df['type']=='firewood','verified_credits_per_stove'].sum()\n",
    "charcoal_VERs=df.loc[df['type']=='charcoal','verified_credits_per_stove'].sum()\n",
    "biomass_VERs=df.loc[df['type']=='biomass','verified_credits_per_stove'].sum()\n",
    "switch_VERs=df.loc[df['type']=='switch','verified_credits_per_stove'].sum()\n",
    "\n",
    "gs_simple_firewood_VERs=df.loc[df['stove']=='gs_simple_firewood','verified_credits_per_stove'].sum()\n",
    "gs_lpg_VERs=df.loc[df['stove']=='gs_lpg','verified_credits_per_stove'].sum()\n",
    "gs_imp_firewood_VERs=df.loc[df['stove']=='gs_imp_firewood','verified_credits_per_stove'].sum()\n",
    "gs_imp_charcoal_VERs=df.loc[df['stove']=='gs_imp_charcoal','verified_credits_per_stove'].sum()\n",
    "cdm_imp_firewood_VERs=df.loc[df['stove']=='cdm_imp_firewood','verified_credits_per_stove'].sum()\n",
    "cdm_imp_charcoal_VERs=df.loc[df['stove']=='cdm_imp_charcoal','verified_credits_per_stove'].sum()\n",
    "\n",
    "\n",
    "CDM_II_G_verified_total=df.loc[df['protocol_type'] == 'CDM AMS II G', 'verified_credits_per_stove'].sum()\n",
    "CDM_I_E_verified_total=df.loc[df['protocol_type'] == 'CDM AMS I E', 'verified_credits_per_stove'].sum()\n",
    "GS_verified_total=df.loc[df['protocol_type'] == 'GS TPDDTEC', 'verified_credits_per_stove'].sum()\n",
    "GS_metered_verified_total=df.loc[df['protocol_type'] == 'GS Metered', 'verified_credits_per_stove'].sum()\n",
    "\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "bangladesh_total=df.loc[df['country'] == 'Bangladesh', 'verified_credits'].sum()\n",
    "burkinafaso_total=df.loc[df['country'] == 'Burkina Faso', 'verified_credits'].sum()\n",
    "cameroon_total=df.loc[df['country'] == 'Cameroon', 'verified_credits'].sum()\n",
    "china_total=df.loc[df['country'] == 'China', 'verified_credits'].sum()\n",
    "ethiopia_total=df.loc[df['country'] == 'Ethiopia', 'verified_credits'].sum()\n",
    "ghana_total=df.loc[df['country'] == 'Ghana', 'verified_credits'].sum()\n",
    "guatemala_total=df.loc[df['country'] == 'Guatemala', 'verified_credits'].sum()\n",
    "haiti_total=df.loc[df['country'] == 'Haiti', 'verified_credits'].sum()\n",
    "honduras_total=df.loc[df['country'] == 'Honduras', 'verified_credits'].sum()\n",
    "india_total=df.loc[df['country'] == 'India', 'verified_credits'].sum()\n",
    "kenya_total=df.loc[df['country'] == 'Kenya', 'verified_credits'].sum()\n",
    "lesotho_total=df.loc[df['country'] == 'Lesotho', 'verified_credits'].sum()\n",
    "malawi_total=df.loc[df['country'] == 'Malawi', 'verified_credits'].sum()\n",
    "mali_total=df.loc[df['country'] == 'Mali', 'verified_credits'].sum()\n",
    "mexico_total=df.loc[df['country'] == 'Mexico', 'verified_credits'].sum()\n",
    "mozambique_total=df.loc[df['country'] == 'Mozambique', 'verified_credits'].sum()\n",
    "myanmar_total=df.loc[df['country'] == 'Myanmar', 'verified_credits'].sum()\n",
    "nepal_total=df.loc[df['country'] == 'Nepal', 'verified_credits'].sum()\n",
    "nigeria_total=df.loc[df['country'] == 'Nigeria', 'verified_credits'].sum()\n",
    "peru_total=df.loc[df['country'] == 'Peru', 'verified_credits'].sum()\n",
    "rwanda_total=df.loc[df['country'] == 'Rwanda', 'verified_credits'].sum()\n",
    "uganda_total=df.loc[df['country'] == 'Uganda', 'verified_credits'].sum()\n",
    "sudan_total=df.loc[df['country'] == 'Sudan', 'verified_credits'].sum()\n",
    "zambia_total=df.loc[df['country'] == 'Zambia', 'verified_credits'].sum()\n",
    "\n",
    "firewood_total=df.loc[df['type'] == 'firewood', 'verified_credits'].sum()\n",
    "charcoal_total=df.loc[df['type'] == 'charcoal', 'verified_credits'].sum()\n",
    "biomass_total=df.loc[df['type'] == 'biomass', 'verified_credits'].sum()\n",
    "switch_total=df.loc[df['type'] == 'switch', 'verified_credits'].sum()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [],
   "source": [
    "ratio_analysis_project_mean_compilation=[]\n",
    "ratio_project_analysis_mean_compilation=[]\n",
    "ratio_analysis_project_lower_compilation=[]\n",
    "ratio_analysis_project_upper_compilation=[]\n",
    "ratio_project_analysis_lower_compilation=[]\n",
    "ratio_project_analysis_upper_compilation=[]\n",
    "\n",
    "# For protocol wide over-estimation\n",
    "GS_Meter_N=0\n",
    "GS_N=0\n",
    "CDM_II_G_N=0\n",
    "CDM_I_E_N=0"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [],
   "source": [
    "bangladesh_N=0\n",
    "burkinafaso_N=0\n",
    "cameroon_N=0\n",
    "china_N=0\n",
    "ethiopia_N=0\n",
    "ghana_N=0\n",
    "guatemala_N=0\n",
    "haiti_N=0\n",
    "honduras_N=0\n",
    "india_N=0\n",
    "kenya_N=0\n",
    "lesotho_N=0\n",
    "malawi_N=0\n",
    "mali_N=0\n",
    "mexico_N=0\n",
    "mozambique_N=0\n",
    "myanmar_N=0\n",
    "nepal_N=0\n",
    "nigeria_N=0\n",
    "peru_N=0\n",
    "rwanda_N=0\n",
    "sudan_N=0\n",
    "uganda_N=0\n",
    "zambia_N=0\n",
    "\n",
    "firewood_N=0\n",
    "charcoal_N=0\n",
    "biomass_N=0\n",
    "switch_N=0\n",
    "\n",
    "ALL_N=0\n",
    "\n",
    "gs_simple_firewood_N=0\n",
    "gs_lpg_N=0\n",
    "gs_imp_firewood_N=0\n",
    "gs_imp_charcoal_N=0\n",
    "cdm_imp_firewood_N=0\n",
    "cdm_imp_charcoal_N=0\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [],
   "source": [
    "bangladesh_VERs_agw=[]\n",
    "burkinafaso_VERs_agw=[]\n",
    "cameroon_VERs_agw=[]\n",
    "china_VERs_agw=[]\n",
    "ethiopia_VERs_agw=[]\n",
    "ghana_VERs_agw=[]\n",
    "guatemala_VERs_agw=[]\n",
    "haiti_VERs_agw=[]\n",
    "honduras_VERs_agw=[]\n",
    "india_VERs_agw=[]\n",
    "kenya_VERs_agw=[]\n",
    "lesotho_VERs_agw=[]\n",
    "malawi_VERs_agw=[]\n",
    "mali_VERs_agw=[]\n",
    "mexico_VERs_agw=[]\n",
    "mozambique_VERs_agw=[]\n",
    "myanmar_VERs_agw=[]\n",
    "nepal_VERs_agw=[]\n",
    "nigeria_VERs_agw=[]\n",
    "peru_VERs_agw=[]\n",
    "rwanda_VERs_agw=[]\n",
    "sudan_VERs_agw=[]\n",
    "uganda_VERs_agw=[]\n",
    "zambia_VERs_agw=[]\n",
    "\n",
    "firewood_VERs_agw=[]\n",
    "charcoal_VERs_agw=[]\n",
    "biomass_VERs_agw=[]\n",
    "switch_VERs_agw=[]\n",
    "\n",
    "gs_simple_firewood_VERs_agw=[]\n",
    "gs_lpg_VERs_agw=[]\n",
    "gs_imp_firewood_VERs_agw=[]\n",
    "gs_imp_charcoal_VERs_agw=[]\n",
    "cdm_imp_firewood_VERs_agw=[]\n",
    "cdm_imp_charcoal_VERs_agw=[]\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [],
   "source": [
    "#compiling our results on a protocol basis\n",
    "all_VERS_agw=[]\n",
    "GS_Meter_VERs_agw=[]\n",
    "GS_VERs_agw=[]\n",
    "CDM_II_G_VERs_agw=[]\n",
    "CDM_I_E_VERs_agw=[]\n",
    "\n",
    "\n",
    "# GS \n",
    "GS_Meter_ratio_analysis_project_mean_compilation=[]\n",
    "GS_Meter_ratio_project_analysis_mean_compilation=[]\n",
    "GS_Meter_ratio_analysis_project_lower_compilation=[]\n",
    "GS_Meter_ratio_analysis_project_upper_compilation=[]\n",
    "GS_Meter_ratio_project_analysis_lower_compilation=[]\n",
    "GS_Meter_ratio_project_analysis_upper_compilation=[]\n",
    "\n",
    "# GS \n",
    "GS_ratio_analysis_project_mean_compilation=[]\n",
    "GS_ratio_project_analysis_mean_compilation=[]\n",
    "GS_ratio_analysis_project_lower_compilation=[]\n",
    "GS_ratio_analysis_project_upper_compilation=[]\n",
    "GS_ratio_project_analysis_lower_compilation=[]\n",
    "GS_ratio_project_analysis_upper_compilation=[]\n",
    "\n",
    "#CDM II G\n",
    "CDM_II_G_ratio_analysis_project_mean_compilation=[]\n",
    "CDM_II_G_ratio_project_analysis_mean_compilation=[]\n",
    "CDM_II_G_ratio_analysis_project_lower_compilation=[]\n",
    "CDM_II_G_ratio_analysis_project_upper_compilation=[]\n",
    "CDM_II_G_ratio_project_analysis_lower_compilation=[]\n",
    "CDM_II_G_ratio_project_analysis_upper_compilation=[]\n",
    "\n",
    "#CDM I E\n",
    "CDM_I_E_ratio_analysis_project_mean_compilation=[]\n",
    "CDM_I_E_ratio_project_analysis_mean_compilation=[]\n",
    "CDM_I_E_ratio_analysis_project_lower_compilation=[]\n",
    "CDM_I_E_ratio_analysis_project_upper_compilation=[]\n",
    "CDM_I_E_ratio_project_analysis_lower_compilation=[]\n",
    "CDM_I_E_ratio_project_analysis_upper_compilation=[]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [],
   "source": [
    "# bangladesh \n",
    "bangladesh_ratio_analysis_project_mean_compilation=[]\n",
    "bangladesh_ratio_project_analysis_mean_compilation=[]\n",
    "bangladesh_ratio_analysis_project_lower_compilation=[]\n",
    "bangladesh_ratio_analysis_project_upper_compilation=[]\n",
    "bangladesh_ratio_project_analysis_lower_compilation=[]\n",
    "bangladesh_ratio_project_analysis_upper_compilation=[]\n",
    "\n",
    "# burkinafaso \n",
    "burkinafaso_ratio_analysis_project_mean_compilation=[]\n",
    "burkinafaso_ratio_project_analysis_mean_compilation=[]\n",
    "burkinafaso_ratio_analysis_project_lower_compilation=[]\n",
    "burkinafaso_ratio_analysis_project_upper_compilation=[]\n",
    "burkinafaso_ratio_project_analysis_lower_compilation=[]\n",
    "burkinafaso_ratio_project_analysis_upper_compilation=[]\n",
    "\n",
    "# cameroon \n",
    "cameroon_ratio_analysis_project_mean_compilation=[]\n",
    "cameroon_ratio_project_analysis_mean_compilation=[]\n",
    "cameroon_ratio_analysis_project_lower_compilation=[]\n",
    "cameroon_ratio_analysis_project_upper_compilation=[]\n",
    "cameroon_ratio_project_analysis_lower_compilation=[]\n",
    "cameroon_ratio_project_analysis_upper_compilation=[]\n",
    "\n",
    "\n",
    "# china \n",
    "china_ratio_analysis_project_mean_compilation=[]\n",
    "china_ratio_project_analysis_mean_compilation=[]\n",
    "china_ratio_analysis_project_lower_compilation=[]\n",
    "china_ratio_analysis_project_upper_compilation=[]\n",
    "china_ratio_project_analysis_lower_compilation=[]\n",
    "china_ratio_project_analysis_upper_compilation=[]\n",
    "\n",
    "# ethiopia \n",
    "ethiopia_ratio_analysis_project_mean_compilation=[]\n",
    "ethiopia_ratio_project_analysis_mean_compilation=[]\n",
    "ethiopia_ratio_analysis_project_lower_compilation=[]\n",
    "ethiopia_ratio_analysis_project_upper_compilation=[]\n",
    "ethiopia_ratio_project_analysis_lower_compilation=[]\n",
    "ethiopia_ratio_project_analysis_upper_compilation=[]\n",
    "\n",
    "# ghana \n",
    "ghana_ratio_analysis_project_mean_compilation=[]\n",
    "ghana_ratio_project_analysis_mean_compilation=[]\n",
    "ghana_ratio_analysis_project_lower_compilation=[]\n",
    "ghana_ratio_analysis_project_upper_compilation=[]\n",
    "ghana_ratio_project_analysis_lower_compilation=[]\n",
    "ghana_ratio_project_analysis_upper_compilation=[]\n",
    "\n",
    "# guatemala \n",
    "guatemala_ratio_analysis_project_mean_compilation=[]\n",
    "guatemala_ratio_project_analysis_mean_compilation=[]\n",
    "guatemala_ratio_analysis_project_lower_compilation=[]\n",
    "guatemala_ratio_analysis_project_upper_compilation=[]\n",
    "guatemala_ratio_project_analysis_lower_compilation=[]\n",
    "guatemala_ratio_project_analysis_upper_compilation=[]\n",
    "\n",
    "# haiti \n",
    "haiti_ratio_analysis_project_mean_compilation=[]\n",
    "haiti_ratio_project_analysis_mean_compilation=[]\n",
    "haiti_ratio_analysis_project_lower_compilation=[]\n",
    "haiti_ratio_analysis_project_upper_compilation=[]\n",
    "haiti_ratio_project_analysis_lower_compilation=[]\n",
    "haiti_ratio_project_analysis_upper_compilation=[]\n",
    "\n",
    "# honduras \n",
    "honduras_ratio_analysis_project_mean_compilation=[]\n",
    "honduras_ratio_project_analysis_mean_compilation=[]\n",
    "honduras_ratio_analysis_project_lower_compilation=[]\n",
    "honduras_ratio_analysis_project_upper_compilation=[]\n",
    "honduras_ratio_project_analysis_lower_compilation=[]\n",
    "honduras_ratio_project_analysis_upper_compilation=[]\n",
    "\n",
    "# india \n",
    "india_ratio_analysis_project_mean_compilation=[]\n",
    "india_ratio_project_analysis_mean_compilation=[]\n",
    "india_ratio_analysis_project_lower_compilation=[]\n",
    "india_ratio_analysis_project_upper_compilation=[]\n",
    "india_ratio_project_analysis_lower_compilation=[]\n",
    "india_ratio_project_analysis_upper_compilation=[]\n",
    "\n",
    "# kenya \n",
    "kenya_ratio_analysis_project_mean_compilation=[]\n",
    "kenya_ratio_project_analysis_mean_compilation=[]\n",
    "kenya_ratio_analysis_project_lower_compilation=[]\n",
    "kenya_ratio_analysis_project_upper_compilation=[]\n",
    "kenya_ratio_project_analysis_lower_compilation=[]\n",
    "kenya_ratio_project_analysis_upper_compilation=[]\n",
    "\n",
    "# lesotho \n",
    "lesotho_ratio_analysis_project_mean_compilation=[]\n",
    "lesotho_ratio_project_analysis_mean_compilation=[]\n",
    "lesotho_ratio_analysis_project_lower_compilation=[]\n",
    "lesotho_ratio_analysis_project_upper_compilation=[]\n",
    "lesotho_ratio_project_analysis_lower_compilation=[]\n",
    "lesotho_ratio_project_analysis_upper_compilation=[]\n",
    "\n",
    "# malawi \n",
    "malawi_ratio_analysis_project_mean_compilation=[]\n",
    "malawi_ratio_project_analysis_mean_compilation=[]\n",
    "malawi_ratio_analysis_project_lower_compilation=[]\n",
    "malawi_ratio_analysis_project_upper_compilation=[]\n",
    "malawi_ratio_project_analysis_lower_compilation=[]\n",
    "malawi_ratio_project_analysis_upper_compilation=[]\n",
    "\n",
    "# mali \n",
    "mali_ratio_analysis_project_mean_compilation=[]\n",
    "mali_ratio_project_analysis_mean_compilation=[]\n",
    "mali_ratio_analysis_project_lower_compilation=[]\n",
    "mali_ratio_analysis_project_upper_compilation=[]\n",
    "mali_ratio_project_analysis_lower_compilation=[]\n",
    "mali_ratio_project_analysis_upper_compilation=[]\n",
    "\n",
    "# mexico \n",
    "mexico_ratio_analysis_project_mean_compilation=[]\n",
    "mexico_ratio_project_analysis_mean_compilation=[]\n",
    "mexico_ratio_analysis_project_lower_compilation=[]\n",
    "mexico_ratio_analysis_project_upper_compilation=[]\n",
    "mexico_ratio_project_analysis_lower_compilation=[]\n",
    "mexico_ratio_project_analysis_upper_compilation=[]\n",
    "\n",
    "# mozambique \n",
    "mozambique_ratio_analysis_project_mean_compilation=[]\n",
    "mozambique_ratio_project_analysis_mean_compilation=[]\n",
    "mozambique_ratio_analysis_project_lower_compilation=[]\n",
    "mozambique_ratio_analysis_project_upper_compilation=[]\n",
    "mozambique_ratio_project_analysis_lower_compilation=[]\n",
    "mozambique_ratio_project_analysis_upper_compilation=[]\n",
    "\n",
    "# myanmar \n",
    "myanmar_ratio_analysis_project_mean_compilation=[]\n",
    "myanmar_ratio_project_analysis_mean_compilation=[]\n",
    "myanmar_ratio_analysis_project_lower_compilation=[]\n",
    "myanmar_ratio_analysis_project_upper_compilation=[]\n",
    "myanmar_ratio_project_analysis_lower_compilation=[]\n",
    "myanmar_ratio_project_analysis_upper_compilation=[]\n",
    "\n",
    "# nepal \n",
    "nepal_ratio_analysis_project_mean_compilation=[]\n",
    "nepal_ratio_project_analysis_mean_compilation=[]\n",
    "nepal_ratio_analysis_project_lower_compilation=[]\n",
    "nepal_ratio_analysis_project_upper_compilation=[]\n",
    "nepal_ratio_project_analysis_lower_compilation=[]\n",
    "nepal_ratio_project_analysis_upper_compilation=[]\n",
    "\n",
    "# nigeria \n",
    "nigeria_ratio_analysis_project_mean_compilation=[]\n",
    "nigeria_ratio_project_analysis_mean_compilation=[]\n",
    "nigeria_ratio_analysis_project_lower_compilation=[]\n",
    "nigeria_ratio_analysis_project_upper_compilation=[]\n",
    "nigeria_ratio_project_analysis_lower_compilation=[]\n",
    "nigeria_ratio_project_analysis_upper_compilation=[]\n",
    "\n",
    "# peru \n",
    "peru_ratio_analysis_project_mean_compilation=[]\n",
    "peru_ratio_project_analysis_mean_compilation=[]\n",
    "peru_ratio_analysis_project_lower_compilation=[]\n",
    "peru_ratio_analysis_project_upper_compilation=[]\n",
    "peru_ratio_project_analysis_lower_compilation=[]\n",
    "peru_ratio_project_analysis_upper_compilation=[]\n",
    "\n",
    "\n",
    "# rwanda \n",
    "rwanda_ratio_analysis_project_mean_compilation=[]\n",
    "rwanda_ratio_project_analysis_mean_compilation=[]\n",
    "rwanda_ratio_analysis_project_lower_compilation=[]\n",
    "rwanda_ratio_analysis_project_upper_compilation=[]\n",
    "rwanda_ratio_project_analysis_lower_compilation=[]\n",
    "rwanda_ratio_project_analysis_upper_compilation=[]\n",
    "\n",
    "# sudan \n",
    "sudan_ratio_analysis_project_mean_compilation=[]\n",
    "sudan_ratio_project_analysis_mean_compilation=[]\n",
    "sudan_ratio_analysis_project_lower_compilation=[]\n",
    "sudan_ratio_analysis_project_upper_compilation=[]\n",
    "sudan_ratio_project_analysis_lower_compilation=[]\n",
    "sudan_ratio_project_analysis_upper_compilation=[]\n",
    "\n",
    "# zambia \n",
    "zambia_ratio_analysis_project_mean_compilation=[]\n",
    "zambia_ratio_project_analysis_mean_compilation=[]\n",
    "zambia_ratio_analysis_project_lower_compilation=[]\n",
    "zambia_ratio_analysis_project_upper_compilation=[]\n",
    "zambia_ratio_project_analysis_lower_compilation=[]\n",
    "zambia_ratio_project_analysis_upper_compilation=[]\n",
    "\n",
    "# uganda \n",
    "uganda_ratio_analysis_project_mean_compilation=[]\n",
    "uganda_ratio_project_analysis_mean_compilation=[]\n",
    "uganda_ratio_analysis_project_lower_compilation=[]\n",
    "uganda_ratio_analysis_project_upper_compilation=[]\n",
    "uganda_ratio_project_analysis_lower_compilation=[]\n",
    "uganda_ratio_project_analysis_upper_compilation=[]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [],
   "source": [
    "\n",
    "\n",
    "# gs_lpg\n",
    "gs_lpg_ratio_analysis_project_mean_compilation=[]\n",
    "gs_lpg_ratio_project_analysis_mean_compilation=[]\n",
    "gs_lpg_ratio_analysis_project_lower_compilation=[]\n",
    "gs_lpg_ratio_analysis_project_upper_compilation=[]\n",
    "gs_lpg_ratio_project_analysis_lower_compilation=[]\n",
    "gs_lpg_ratio_project_analysis_upper_compilation=[]\n",
    "\n",
    "# gs_imp_firewood\n",
    "gs_imp_firewood_ratio_analysis_project_mean_compilation=[]\n",
    "gs_imp_firewood_ratio_project_analysis_mean_compilation=[]\n",
    "gs_imp_firewood_ratio_analysis_project_lower_compilation=[]\n",
    "gs_imp_firewood_ratio_analysis_project_upper_compilation=[]\n",
    "gs_imp_firewood_ratio_project_analysis_lower_compilation=[]\n",
    "gs_imp_firewood_ratio_project_analysis_upper_compilation=[]\n",
    "\n",
    "# gs_imp_charcoal\n",
    "gs_imp_charcoal_ratio_analysis_project_mean_compilation=[]\n",
    "gs_imp_charcoal_ratio_project_analysis_mean_compilation=[]\n",
    "gs_imp_charcoal_ratio_analysis_project_lower_compilation=[]\n",
    "gs_imp_charcoal_ratio_analysis_project_upper_compilation=[]\n",
    "gs_imp_charcoal_ratio_project_analysis_lower_compilation=[]\n",
    "gs_imp_charcoal_ratio_project_analysis_upper_compilation=[]\n",
    "\n",
    "# cdm_imp_firewood\n",
    "cdm_imp_firewood_ratio_analysis_project_mean_compilation=[]\n",
    "cdm_imp_firewood_ratio_project_analysis_mean_compilation=[]\n",
    "cdm_imp_firewood_ratio_analysis_project_lower_compilation=[]\n",
    "cdm_imp_firewood_ratio_analysis_project_upper_compilation=[]\n",
    "cdm_imp_firewood_ratio_project_analysis_lower_compilation=[]\n",
    "cdm_imp_firewood_ratio_project_analysis_upper_compilation=[]\n",
    "\n",
    "# cdm_imp_charcoal\n",
    "cdm_imp_charcoal_ratio_analysis_project_mean_compilation=[]\n",
    "cdm_imp_charcoal_ratio_project_analysis_mean_compilation=[]\n",
    "cdm_imp_charcoal_ratio_analysis_project_lower_compilation=[]\n",
    "cdm_imp_charcoal_ratio_analysis_project_upper_compilation=[]\n",
    "cdm_imp_charcoal_ratio_project_analysis_lower_compilation=[]\n",
    "cdm_imp_charcoal_ratio_project_analysis_upper_compilation=[]\n",
    "\n",
    "# gs_simple_firewood\n",
    "gs_simple_firewood_ratio_analysis_project_mean_compilation=[]\n",
    "gs_simple_firewood_ratio_project_analysis_mean_compilation=[]\n",
    "gs_simple_firewood_ratio_analysis_project_lower_compilation=[]\n",
    "gs_simple_firewood_ratio_analysis_project_upper_compilation=[]\n",
    "gs_simple_firewood_ratio_project_analysis_lower_compilation=[]\n",
    "gs_simple_firewood_ratio_project_analysis_upper_compilation=[]\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [],
   "source": [
    "#import relevant packages\n",
    "%matplotlib inline\n",
    "import pandas as pd\n",
    "import numpy as np\n",
    "import matplotlib.pyplot as plt\n",
    "import time\n",
    "\n",
    "from scipy import stats as stats\n",
    "from scipy.stats import triang\n",
    "\n",
    "import seaborn as sns\n",
    "\n",
    "#triangular distributions\n",
    "\n",
    "np.random.seed(4)\n",
    "N = 10000                # random numbers to generate\n",
    "\n",
    "    \n",
    "#Adoption \n",
    "min_adoption=.51\n",
    "mean_adoption=.59\n",
    "max_adoption= .65\n",
    "n_adoption=np.random.triangular(left=min_adoption, mode=mean_adoption, right=max_adoption, size=(N))\n",
    "    \n",
    "req_adoption=np.random.uniform(0,mean_adoption,N)      \n",
    "#Usage\n",
    "min_usage=0.16\n",
    "mean_usage=0.52\n",
    "max_usage=0.85\n",
    "n_usage=np.random.triangular(left=min_usage, mode=mean_usage, right=max_usage, size=(N))\n",
    "req_usage=np.random.uniform(0,mean_usage,N)\n",
    "                               \n",
    "#Stacking\n",
    "min_stacking=.193\n",
    "mean_stacking=.673\n",
    "max_stacking=.99\n",
    "n_stacking=np.random.triangular(left=min_stacking, mode=mean_stacking, right=max_stacking, size=(N))\n",
    "req_stacking=np.random.uniform(0,mean_stacking,N)\n",
    "\n",
    "# Renewable GWP from Whitman & Lehman\n",
    "\n",
    "gwp_ch4_renewable=28 #tCO2e/tonne CH4 i.e., instead of typically 25\n",
    "\n",
    "# Whitman & lehmann 2011 also cite a value for CO \n",
    "# but we do not include CO in our analysis\n",
    "    \n",
    "gwp_ch4_agw=30\n",
    "gwp_n2O_agw=298\n",
    "    \n",
    "eff_trad_firewood=0.14\n",
    "eff_improv_firewood=.21\n",
    "eff_trad_charcoal=0.2\n",
    "eff_improv_charcoal=.22\n",
    "eff_kerosene=.47\n",
    "eff_LPG=.5\n",
    "eff_ethanol=.47\n",
    "eff_pellets=.5 #only two pellet projects is the mimi motto\n",
    "\n",
    "ncv_firewood_agw=.016\n",
    "ncv_charcoal_agw=.027\n",
    "ncv_kerosene_agw=.043\n",
    "ncv_LPG_agw=.045\n",
    "ncv_ethanol_agw=.027\n",
    "ncv_pellets_agw=.017\n",
    "\n",
    "kg_per_liter_ethanol_agw=0.783 #kg/liter\n",
    "    \n",
    "ef_combustion_co2_firewood=95.1 #tCO2/TJ updated good \n",
    "ef_combustion_ch4_firewood=0.58 #tCH4/TJ updated good\n",
    "ef_combustion_n2O_firewood=0 #tN2O/TJ\n",
    "\n",
    "ef_combustion_co2_charcoal=114.75 #tCO2/TJ updated good\n",
    "ef_combustion_ch4_charcoal=0.255 #tCH4/TJ updated good\n",
    "ef_combustion_n2O_charcoal=0 #tN2O/TJ good\n",
    "\n",
    "\n",
    "ef_total_not_combustion_co2_charcoal=87.725 #tCO2/TJ updated good\n",
    "ef_total_not_combustion_ch4_charcoal=2.06 #tCH4/TJ updated good\n",
    "ef_total_not_combustion_n2O_charcoal=0.0075#tN2O/TJ good\n",
    "    \n",
    "ef_total_co2_kerosene=79.9 #tCO2/TJ good\n",
    "ef_total_ch4_kerosene=.107 #tCH4/TJ good\n",
    "ef_total_n2O_kerosene=0.000141 #tN2O/TJ good\n",
    "    \n",
    "ef_total_co2_LPG=82.8#tCO2/TJ good\n",
    "ef_total_ch4_LPG=0.117 #tCH4/TJ good\n",
    "ef_total_n2O_LPG=0.000188 #tN2O/TJ good\n",
    "    \n",
    "ef_total_co2_ethanol=60.45#tCO2/TJ good\n",
    "ef_total_ch4_ethanol=0.61 #tCH4/TJ good\n",
    "ef_total_n2O_ethanol=0.05 #tN2O/TJ good \n",
    "    \n",
    "ef_total_co2_pellets=90.05#tCO2/TJ good\n",
    "ef_total_ch4_pellets=0.105#tCH4/TJ good\n",
    "ef_total_n2O_pellets=0 #tN2O/TJ good\n",
    "    \n",
    "ef_total_zeroed_co2_ethanol=-9.78#tCO2/TJ good\n",
    "ef_total_zeroed_ch4_ethanol=0.61 #tCH4/TJ good\n",
    "ef_total_zeroed_n2O_ethanol=0.05 #tN2O/TJ good\n",
    "    \n",
    "ef_total_zeroed_co2_pellets=0 #tCO2/TJ good\n",
    "ef_total_zeroed_ch4_pellets=0.105 #tCH4/TJ good\n",
    "ef_total_zeroed_n2O_pellets=0 #tN2O/TJ good\n",
    "    \n",
    "sum_VERS_not_per_stove_per_day_agw=0  "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "7.522792138375865e-05\n",
      "We find that GS10974 project reports emissions which are 24.67 [ -132.84  ,  182.19 ] times overcredited.\n",
      "Our analysis finds emissions which are 4.05 [ 4.05  ,  4.05 ] % of the verified credits in the PDD.\n",
      "0.030081576300695405\n",
      "We find that GS10974 project reports emissions which are 1.48 [ 1.48  ,  1.48 ] times overcredited.\n",
      "Our analysis finds emissions which are 67.72 [ 67.71  ,  67.72 ] % of the verified credits in the PDD.\n",
      "7.35921782099292e-05\n",
      "We find that GS10974 project reports emissions which are 24.47 [ -135.61  ,  184.55 ] times overcredited.\n",
      "Our analysis finds emissions which are 4.09 [ 4.09  ,  4.09 ] % of the verified credits in the PDD.\n",
      "0.029570674125931264\n",
      "We find that GS10974 project reports emissions which are 1.48 [ 1.48  ,  1.48 ] times overcredited.\n",
      "Our analysis finds emissions which are 67.75 [ 67.74  ,  67.75 ] % of the verified credits in the PDD.\n",
      "0.00017404169818050857\n",
      "We find that VCS1216 project reports emissions which are 44.9 [ 43.48  ,  46.31 ] times overcredited.\n",
      "Our analysis finds emissions which are 2.23 [ 2.23  ,  2.23 ] % of the verified credits in the PDD.\n",
      "0.0001735346887718854\n",
      "We find that VCS1216 project reports emissions which are 45.03 [ 43.56  ,  46.5 ] times overcredited.\n",
      "Our analysis finds emissions which are 2.22 [ 2.22  ,  2.22 ] % of the verified credits in the PDD.\n",
      "0.00017548861800918564\n",
      "We find that VCS1216 project reports emissions which are 44.72 [ 43.27  ,  46.17 ] times overcredited.\n",
      "Our analysis finds emissions which are 2.24 [ 2.24  ,  2.24 ] % of the verified credits in the PDD.\n",
      "0.00017542809819886268\n",
      "We find that VCS1216 project reports emissions which are 44.74 [ 43.32  ,  46.15 ] times overcredited.\n",
      "Our analysis finds emissions which are 2.24 [ 2.24  ,  2.24 ] % of the verified credits in the PDD.\n",
      "0.00017286341043794602\n",
      "We find that VCS1216 project reports emissions which are 45.06 [ 43.62  ,  46.51 ] times overcredited.\n",
      "Our analysis finds emissions which are 2.22 [ 2.22  ,  2.22 ] % of the verified credits in the PDD.\n",
      "0.00017278437290405376\n",
      "We find that VCS1216 project reports emissions which are 45.08 [ 43.66  ,  46.51 ] times overcredited.\n",
      "Our analysis finds emissions which are 2.22 [ 2.22  ,  2.22 ] % of the verified credits in the PDD.\n",
      "0.00017729119972823625\n",
      "We find that VCS1216 project reports emissions which are 44.47 [ 43.01  ,  45.92 ] times overcredited.\n",
      "Our analysis finds emissions which are 2.25 [ 2.25  ,  2.25 ] % of the verified credits in the PDD.\n",
      "0.0001775149000570274\n",
      "We find that VCS1216 project reports emissions which are 44.41 [ 43.0  ,  45.82 ] times overcredited.\n",
      "Our analysis finds emissions which are 2.25 [ 2.25  ,  2.25 ] % of the verified credits in the PDD.\n",
      "0.00017504701510788023\n",
      "We find that VCS1216 project reports emissions which are 44.74 [ 43.31  ,  46.16 ] times overcredited.\n",
      "Our analysis finds emissions which are 2.24 [ 2.24  ,  2.24 ] % of the verified credits in the PDD.\n",
      "0.00017468097581905918\n",
      "We find that VCS1216 project reports emissions which are 44.83 [ 43.4  ,  46.26 ] times overcredited.\n",
      "Our analysis finds emissions which are 2.23 [ 2.23  ,  2.23 ] % of the verified credits in the PDD.\n",
      "0.000171225087806943\n",
      "We find that VCS1216 project reports emissions which are 45.34 [ 43.93  ,  46.75 ] times overcredited.\n",
      "Our analysis finds emissions which are 2.21 [ 2.21  ,  2.21 ] % of the verified credits in the PDD.\n",
      "0.0001711796327391189\n",
      "We find that VCS1216 project reports emissions which are 45.35 [ 43.92  ,  46.78 ] times overcredited.\n",
      "Our analysis finds emissions which are 2.21 [ 2.2  ,  2.21 ] % of the verified credits in the PDD.\n",
      "0.00014381488386633643\n",
      "We find that VCS1216 project reports emissions which are 27.84 [ 26.93  ,  28.75 ] times overcredited.\n",
      "Our analysis finds emissions which are 3.59 [ 3.59  ,  3.59 ] % of the verified credits in the PDD.\n",
      "0.00014376178785705386\n",
      "We find that VCS1216 project reports emissions which are 27.85 [ 26.97  ,  28.74 ] times overcredited.\n",
      "Our analysis finds emissions which are 3.59 [ 3.59  ,  3.59 ] % of the verified credits in the PDD.\n",
      "0.00013911990828437783\n",
      "We find that VCS1216 project reports emissions which are 28.44 [ 27.53  ,  29.35 ] times overcredited.\n",
      "Our analysis finds emissions which are 3.52 [ 3.52  ,  3.52 ] % of the verified credits in the PDD.\n",
      "0.00013902488034517538\n",
      "We find that VCS1216 project reports emissions which are 28.46 [ 27.55  ,  29.36 ] times overcredited.\n",
      "Our analysis finds emissions which are 3.51 [ 3.51  ,  3.51 ] % of the verified credits in the PDD.\n",
      "0.00014568332432113017\n",
      "We find that VCS1216 project reports emissions which are 27.64 [ 26.72  ,  28.56 ] times overcredited.\n",
      "Our analysis finds emissions which are 3.62 [ 3.62  ,  3.62 ] % of the verified credits in the PDD.\n",
      "0.000145688498712371\n",
      "We find that VCS1216 project reports emissions which are 27.64 [ 26.76  ,  28.52 ] times overcredited.\n",
      "Our analysis finds emissions which are 3.62 [ 3.62  ,  3.62 ] % of the verified credits in the PDD.\n",
      "0.00014632057250420888\n",
      "We find that VCS1216 project reports emissions which are 27.57 [ 26.71  ,  28.44 ] times overcredited.\n",
      "Our analysis finds emissions which are 3.63 [ 3.63  ,  3.63 ] % of the verified credits in the PDD.\n",
      "0.0001464492447135454\n",
      "We find that VCS1216 project reports emissions which are 27.55 [ 26.67  ,  28.44 ] times overcredited.\n",
      "Our analysis finds emissions which are 3.63 [ 3.63  ,  3.63 ] % of the verified credits in the PDD.\n",
      "0.00014625108557972852\n",
      "We find that VCS1216 project reports emissions which are 27.59 [ 26.69  ,  28.49 ] times overcredited.\n",
      "Our analysis finds emissions which are 3.62 [ 3.62  ,  3.63 ] % of the verified credits in the PDD.\n",
      "0.00014629168206307277\n",
      "We find that VCS1216 project reports emissions which are 27.58 [ 26.69  ,  28.47 ] times overcredited.\n",
      "Our analysis finds emissions which are 3.63 [ 3.63  ,  3.63 ] % of the verified credits in the PDD.\n",
      "0.00015840075837407635\n",
      "We find that VCS1216 project reports emissions which are 26.26 [ 25.43  ,  27.09 ] times overcredited.\n",
      "Our analysis finds emissions which are 3.81 [ 3.81  ,  3.81 ] % of the verified credits in the PDD.\n",
      "0.00015831638362769253\n",
      "We find that VCS1216 project reports emissions which are 26.27 [ 25.45  ,  27.09 ] times overcredited.\n",
      "Our analysis finds emissions which are 3.81 [ 3.81  ,  3.81 ] % of the verified credits in the PDD.\n",
      "We find that GS2094 project reports emissions which are 19.36 [ 18.77  ,  19.94 ] times overcredited.\n",
      "Our analysis finds emissions which are 5.17 [ 5.17  ,  5.17 ] % of the verified credits in the PDD.\n",
      "We find that GS2094 project reports emissions which are 19.35 [ 18.77  ,  19.93 ] times overcredited.\n",
      "Our analysis finds emissions which are 5.17 [ 5.17  ,  5.17 ] % of the verified credits in the PDD.\n",
      "We find that GS2094 project reports emissions which are 19.94 [ 19.35  ,  20.54 ] times overcredited.\n",
      "Our analysis finds emissions which are 5.01 [ 5.01  ,  5.01 ] % of the verified credits in the PDD.\n",
      "We find that GS2094 project reports emissions which are 19.95 [ 19.35  ,  20.55 ] times overcredited.\n",
      "Our analysis finds emissions which are 5.01 [ 5.01  ,  5.01 ] % of the verified credits in the PDD.\n",
      "We find that GS2094 project reports emissions which are 20.39 [ 19.77  ,  21.0 ] times overcredited.\n",
      "Our analysis finds emissions which are 4.91 [ 4.9  ,  4.91 ] % of the verified credits in the PDD.\n",
      "We find that GS2094 project reports emissions which are 20.38 [ 19.77  ,  20.99 ] times overcredited.\n",
      "Our analysis finds emissions which are 4.91 [ 4.91  ,  4.91 ] % of the verified credits in the PDD.\n",
      "We find that GS2094 project reports emissions which are 22.0 [ 21.34  ,  22.66 ] times overcredited.\n",
      "Our analysis finds emissions which are 4.55 [ 4.54  ,  4.55 ] % of the verified credits in the PDD.\n",
      "We find that GS2094 project reports emissions which are 22.0 [ 21.34  ,  22.66 ] times overcredited.\n",
      "Our analysis finds emissions which are 4.55 [ 4.55  ,  4.55 ] % of the verified credits in the PDD.\n",
      "We find that GS2094 project reports emissions which are 24.05 [ 23.33  ,  24.78 ] times overcredited.\n",
      "Our analysis finds emissions which are 4.16 [ 4.16  ,  4.16 ] % of the verified credits in the PDD.\n",
      "We find that GS2094 project reports emissions which are 24.05 [ 23.32  ,  24.77 ] times overcredited.\n",
      "Our analysis finds emissions which are 4.16 [ 4.16  ,  4.16 ] % of the verified credits in the PDD.\n",
      "We find that GS2094 project reports emissions which are 29.84 [ 28.95  ,  30.74 ] times overcredited.\n",
      "Our analysis finds emissions which are 3.35 [ 3.35  ,  3.35 ] % of the verified credits in the PDD.\n",
      "We find that GS2094 project reports emissions which are 29.83 [ 28.93  ,  30.72 ] times overcredited.\n",
      "Our analysis finds emissions which are 3.35 [ 3.35  ,  3.35 ] % of the verified credits in the PDD.\n",
      "0.00022915039436547532\n",
      "We find that VCS1719 project reports emissions which are 25.27 [ 24.51  ,  26.04 ] times overcredited.\n",
      "Our analysis finds emissions which are 3.96 [ 3.96  ,  3.96 ] % of the verified credits in the PDD.\n",
      "0.0003609146862268365\n",
      "We find that VCS1721 project reports emissions which are 23.87 [ 23.15  ,  24.59 ] times overcredited.\n",
      "Our analysis finds emissions which are 4.19 [ 4.19  ,  4.19 ] % of the verified credits in the PDD.\n",
      "We find that GS7312 project reports emissions which are 8.48 [ 8.44  ,  8.52 ] times overcredited.\n",
      "Our analysis finds emissions which are 11.79 [ 11.79  ,  11.79 ] % of the verified credits in the PDD.\n",
      "We find that GS7312 project reports emissions which are 8.48 [ 8.44  ,  8.52 ] times overcredited.\n",
      "Our analysis finds emissions which are 11.8 [ 11.8  ,  11.8 ] % of the verified credits in the PDD.\n",
      "We find that GS500 project reports emissions which are 38.03 [ 37.75  ,  38.31 ] times overcredited.\n",
      "Our analysis finds emissions which are 2.63 [ 2.63  ,  2.63 ] % of the verified credits in the PDD.\n",
      "We find that GS10884 project reports emissions which are 3.57 [ 3.54  ,  3.6 ] times overcredited.\n",
      "Our analysis finds emissions which are 28.02 [ 28.02  ,  28.02 ] % of the verified credits in the PDD.\n",
      "We find that GS10884 project reports emissions which are 0.97 [ 0.96  ,  0.97 ] times overcredited.\n",
      "Our analysis finds emissions which are 103.56 [ 103.56  ,  103.56 ] % of the verified credits in the PDD.\n",
      "We find that GS447 project reports emissions which are 8.68 [ 8.64  ,  8.72 ] times overcredited.\n",
      "Our analysis finds emissions which are 11.52 [ 11.52  ,  11.52 ] % of the verified credits in the PDD.\n",
      "We find that GS447 project reports emissions which are 8.82 [ 8.78  ,  8.87 ] times overcredited.\n",
      "Our analysis finds emissions which are 11.33 [ 11.33  ,  11.33 ] % of the verified credits in the PDD.\n",
      "We find that GS2564 project reports emissions which are 12.07 [ 12.02  ,  12.13 ] times overcredited.\n",
      "Our analysis finds emissions which are 8.28 [ 8.28  ,  8.28 ] % of the verified credits in the PDD.\n",
      "We find that GS7438 project reports emissions which are 5.11 [ 5.09  ,  5.14 ] times overcredited.\n",
      "Our analysis finds emissions which are 19.57 [ 19.56  ,  19.57 ] % of the verified credits in the PDD.\n",
      "We find that GS7438 project reports emissions which are 46.09 [ 45.87  ,  46.3 ] times overcredited.\n",
      "Our analysis finds emissions which are 2.17 [ 2.17  ,  2.17 ] % of the verified credits in the PDD.\n",
      "0.00035540638564643217\n",
      "We find that GS913 project reports emissions which are 21.18 [ 20.53  ,  21.83 ] times overcredited.\n",
      "Our analysis finds emissions which are 4.72 [ 4.72  ,  4.72 ] % of the verified credits in the PDD.\n",
      "0.0004459253932837502\n",
      "We find that GS6212 project reports emissions which are 13.08 [ 12.69  ,  13.48 ] times overcredited.\n",
      "Our analysis finds emissions which are 7.64 [ 7.64  ,  7.64 ] % of the verified credits in the PDD.\n",
      "0.0004460847325234239\n",
      "We find that GS6212 project reports emissions which are 14.55 [ 14.11  ,  14.99 ] times overcredited.\n",
      "Our analysis finds emissions which are 6.87 [ 6.87  ,  6.87 ] % of the verified credits in the PDD.\n",
      "We find that GS2758 project reports emissions which are 7.51 [ -0.64  ,  15.66 ] times overcredited.\n",
      "Our analysis finds emissions which are 13.32 [ 13.32  ,  13.32 ] % of the verified credits in the PDD.\n",
      "We find that GS2758 project reports emissions which are 7.54 [ -2.27  ,  17.35 ] times overcredited.\n",
      "Our analysis finds emissions which are 13.26 [ 13.26  ,  13.26 ] % of the verified credits in the PDD.\n",
      "We find that GS2758 project reports emissions which are 7.66 [ -1.18  ,  16.51 ] times overcredited.\n",
      "Our analysis finds emissions which are 13.05 [ 13.04  ,  13.05 ] % of the verified credits in the PDD.\n",
      "We find that GS2758 project reports emissions which are 8.22 [ 0.15  ,  16.29 ] times overcredited.\n",
      "Our analysis finds emissions which are 12.17 [ 12.17  ,  12.17 ] % of the verified credits in the PDD.\n",
      "We find that GS2758 project reports emissions which are 8.31 [ -1.61  ,  18.23 ] times overcredited.\n",
      "Our analysis finds emissions which are 12.04 [ 12.03  ,  12.04 ] % of the verified credits in the PDD.\n",
      "We find that GS2758 project reports emissions which are 8.34 [ -1.95  ,  18.64 ] times overcredited.\n",
      "Our analysis finds emissions which are 11.98 [ 11.98  ,  11.99 ] % of the verified credits in the PDD.\n",
      "We find that GS3112 project reports emissions which are 20.2 [ 19.59  ,  20.81 ] times overcredited.\n",
      "Our analysis finds emissions which are 4.95 [ 4.95  ,  4.95 ] % of the verified credits in the PDD.\n",
      "We find that GS3112 project reports emissions which are 19.72 [ 19.12  ,  20.31 ] times overcredited.\n",
      "Our analysis finds emissions which are 5.07 [ 5.07  ,  5.07 ] % of the verified credits in the PDD.\n",
      "We find that GS3112 project reports emissions which are 19.71 [ 19.12  ,  20.3 ] times overcredited.\n",
      "Our analysis finds emissions which are 5.07 [ 5.07  ,  5.07 ] % of the verified credits in the PDD.\n",
      "We find that GS3112 project reports emissions which are 19.39 [ 18.81  ,  19.97 ] times overcredited.\n",
      "Our analysis finds emissions which are 5.16 [ 5.16  ,  5.16 ] % of the verified credits in the PDD.\n",
      "We find that GS5642 project reports emissions which are 3.51 [ 3.5  ,  3.53 ] times overcredited.\n",
      "Our analysis finds emissions which are 28.47 [ 28.47  ,  28.47 ] % of the verified credits in the PDD.\n",
      "We find that GS5642 project reports emissions which are 21.26 [ 21.16  ,  21.36 ] times overcredited.\n",
      "Our analysis finds emissions which are 4.7 [ 4.7  ,  4.7 ] % of the verified credits in the PDD.\n",
      "We find that GS1060 project reports emissions which are 21.62 [ 20.97  ,  22.27 ] times overcredited.\n",
      "Our analysis finds emissions which are 4.62 [ 4.62  ,  4.63 ] % of the verified credits in the PDD.\n",
      "We find that GS1060 project reports emissions which are 22.87 [ 22.19  ,  23.56 ] times overcredited.\n",
      "Our analysis finds emissions which are 4.37 [ 4.37  ,  4.37 ] % of the verified credits in the PDD.\n",
      "We find that GS1060 project reports emissions which are 19.66 [ 19.07  ,  20.25 ] times overcredited.\n",
      "Our analysis finds emissions which are 5.09 [ 5.09  ,  5.09 ] % of the verified credits in the PDD.\n",
      "We find that GS1060 project reports emissions which are 21.62 [ 20.97  ,  22.27 ] times overcredited.\n",
      "Our analysis finds emissions which are 4.62 [ 4.62  ,  4.63 ] % of the verified credits in the PDD.\n",
      "We find that GS1060 project reports emissions which are 22.87 [ 22.18  ,  23.55 ] times overcredited.\n",
      "Our analysis finds emissions which are 4.37 [ 4.37  ,  4.37 ] % of the verified credits in the PDD.\n",
      "We find that GS1060 project reports emissions which are 19.66 [ 19.07  ,  20.25 ] times overcredited.\n",
      "Our analysis finds emissions which are 5.09 [ 5.08  ,  5.09 ] % of the verified credits in the PDD.\n",
      "8.173981196497706e-06\n",
      "We find that GS2744 project reports emissions which are 145.58 [ 140.79  ,  150.38 ] times overcredited.\n",
      "Our analysis finds emissions which are 0.69 [ 0.69  ,  0.69 ] % of the verified credits in the PDD.\n",
      "1.8855138721803875e-05\n",
      "We find that GS2744 project reports emissions which are 69.6 [ 67.38  ,  71.82 ] times overcredited.\n",
      "Our analysis finds emissions which are 1.44 [ 1.44  ,  1.44 ] % of the verified credits in the PDD.\n",
      "0.00018093922774734428\n",
      "We find that GS1094 project reports emissions which are 189.7 [ 183.93  ,  195.47 ] times overcredited.\n",
      "Our analysis finds emissions which are 0.53 [ 0.53  ,  0.53 ] % of the verified credits in the PDD.\n",
      "0.00018642155938238238\n",
      "We find that GS1094 project reports emissions which are 183.66 [ 178.11  ,  189.21 ] times overcredited.\n",
      "Our analysis finds emissions which are 0.54 [ 0.54  ,  0.54 ] % of the verified credits in the PDD.\n",
      "0.00018311980981245514\n",
      "We find that GS1094 project reports emissions which are 196.66 [ 190.72  ,  202.61 ] times overcredited.\n",
      "Our analysis finds emissions which are 0.51 [ 0.51  ,  0.51 ] % of the verified credits in the PDD.\n",
      "0.00016188629414325454\n",
      "We find that GS1094 project reports emissions which are 210.86 [ 204.46  ,  217.26 ] times overcredited.\n",
      "Our analysis finds emissions which are 0.47 [ 0.47  ,  0.47 ] % of the verified credits in the PDD.\n",
      "0.00013205953236628844\n",
      "We find that GS1094 project reports emissions which are 238.38 [ 231.16  ,  245.6 ] times overcredited.\n",
      "Our analysis finds emissions which are 0.42 [ 0.42  ,  0.42 ] % of the verified credits in the PDD.\n",
      "0.000141192804598566\n",
      "We find that GS1094 project reports emissions which are 228.77 [ 221.88  ,  235.65 ] times overcredited.\n",
      "Our analysis finds emissions which are 0.44 [ 0.44  ,  0.44 ] % of the verified credits in the PDD.\n",
      "9.326691575716592e-05\n",
      "We find that GS4291 project reports emissions which are 24.46 [ 23.72  ,  25.2 ] times overcredited.\n",
      "Our analysis finds emissions which are 4.09 [ 4.09  ,  4.09 ] % of the verified credits in the PDD.\n",
      "9.300545211051901e-05\n",
      "We find that GS4291 project reports emissions which are 24.53 [ 23.77  ,  25.28 ] times overcredited.\n",
      "Our analysis finds emissions which are 4.08 [ 4.08  ,  4.08 ] % of the verified credits in the PDD.\n",
      "5.5423105665854445e-05\n",
      "We find that GS6604 project reports emissions which are 57.98 [ 56.21  ,  59.74 ] times overcredited.\n",
      "Our analysis finds emissions which are 1.72 [ 1.72  ,  1.72 ] % of the verified credits in the PDD.\n",
      "We find that GS 2896  project reports emissions which are 32.99 [ 31.96  ,  34.02 ] times overcredited.\n",
      "Our analysis finds emissions which are 3.03 [ 3.03  ,  3.03 ] % of the verified credits in the PDD.\n",
      "We find that GS 3422  project reports emissions which are 57.81 [ 56.05  ,  59.56 ] times overcredited.\n",
      "Our analysis finds emissions which are 1.73 [ 1.73  ,  1.73 ] % of the verified credits in the PDD.\n",
      "We find that GS 2456 project reports emissions which are 75.56 [ 73.29  ,  77.84 ] times overcredited.\n",
      "Our analysis finds emissions which are 1.32 [ 1.32  ,  1.32 ] % of the verified credits in the PDD.\n",
      "We find that GS 2456 project reports emissions which are 90.52 [ 87.79  ,  93.24 ] times overcredited.\n",
      "Our analysis finds emissions which are 1.1 [ 1.1  ,  1.1 ] % of the verified credits in the PDD.\n",
      "We find that GS 2456 project reports emissions which are 99.77 [ 96.78  ,  102.76 ] times overcredited.\n",
      "Our analysis finds emissions which are 1.0 [ 1.0  ,  1.0 ] % of the verified credits in the PDD.\n",
      "We find that GS 2456 project reports emissions which are 107.86 [ 104.61  ,  111.11 ] times overcredited.\n",
      "Our analysis finds emissions which are 0.93 [ 0.93  ,  0.93 ] % of the verified credits in the PDD.\n",
      "We find that GS 2456 project reports emissions which are 92.17 [ 89.4  ,  94.93 ] times overcredited.\n",
      "Our analysis finds emissions which are 1.08 [ 1.08  ,  1.09 ] % of the verified credits in the PDD.\n",
      "We find that GS 1028 project reports emissions which are 71.34 [ 69.16  ,  73.52 ] times overcredited.\n",
      "Our analysis finds emissions which are 1.4 [ 1.4  ,  1.4 ] % of the verified credits in the PDD.\n",
      "We find that GS 1028 project reports emissions which are 64.03 [ 62.08  ,  65.97 ] times overcredited.\n",
      "Our analysis finds emissions which are 1.56 [ 1.56  ,  1.56 ] % of the verified credits in the PDD.\n",
      "We find that GS 1028 project reports emissions which are 63.36 [ 61.39  ,  65.33 ] times overcredited.\n",
      "Our analysis finds emissions which are 1.58 [ 1.58  ,  1.58 ] % of the verified credits in the PDD.\n",
      "We find that GS 1028 project reports emissions which are 107.34 [ 104.05  ,  110.63 ] times overcredited.\n",
      "Our analysis finds emissions which are 0.93 [ 0.93  ,  0.93 ] % of the verified credits in the PDD.\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "We find that GS 1028 project reports emissions which are 32.44 [ 31.44  ,  33.44 ] times overcredited.\n",
      "Our analysis finds emissions which are 3.08 [ 3.08  ,  3.08 ] % of the verified credits in the PDD.\n",
      "We find that GS 3018 project reports emissions which are 8.57 [ 8.31  ,  8.83 ] times overcredited.\n",
      "Our analysis finds emissions which are 11.67 [ 11.67  ,  11.67 ] % of the verified credits in the PDD.\n",
      "We find that GS 3018 project reports emissions which are 5.66 [ 5.49  ,  5.83 ] times overcredited.\n",
      "Our analysis finds emissions which are 17.67 [ 17.67  ,  17.68 ] % of the verified credits in the PDD.\n",
      "We find that GS 3018 project reports emissions which are 2.84 [ 2.76  ,  2.93 ] times overcredited.\n",
      "Our analysis finds emissions which are 35.17 [ 35.17  ,  35.17 ] % of the verified credits in the PDD.\n",
      "We find that GS1267 project reports emissions which are 7.24 [ 7.21  ,  7.28 ] times overcredited.\n",
      "Our analysis finds emissions which are 13.81 [ 13.81  ,  13.81 ] % of the verified credits in the PDD.\n",
      "We find that GS1267 project reports emissions which are 7.03 [ 7.0  ,  7.06 ] times overcredited.\n",
      "Our analysis finds emissions which are 14.22 [ 14.22  ,  14.23 ] % of the verified credits in the PDD.\n",
      "We find that GS1267 project reports emissions which are 6.95 [ 6.91  ,  6.98 ] times overcredited.\n",
      "Our analysis finds emissions which are 14.39 [ 14.39  ,  14.4 ] % of the verified credits in the PDD.\n",
      "We find that GS1267 project reports emissions which are 6.46 [ 6.43  ,  6.49 ] times overcredited.\n",
      "Our analysis finds emissions which are 15.47 [ 15.47  ,  15.47 ] % of the verified credits in the PDD.\n",
      "We find that GS1267 project reports emissions which are 6.12 [ 6.09  ,  6.15 ] times overcredited.\n",
      "Our analysis finds emissions which are 16.34 [ 16.34  ,  16.34 ] % of the verified credits in the PDD.\n",
      "We find that GS1267 project reports emissions which are 1.69 [ 1.68  ,  1.7 ] times overcredited.\n",
      "Our analysis finds emissions which are 59.12 [ 59.12  ,  59.12 ] % of the verified credits in the PDD.\n",
      "We find that GS2439 project reports emissions which are 11.01 [ 10.95  ,  11.06 ] times overcredited.\n",
      "Our analysis finds emissions which are 9.09 [ 9.09  ,  9.09 ] % of the verified credits in the PDD.\n",
      "We find that GS2439 project reports emissions which are 10.89 [ 10.84  ,  10.95 ] times overcredited.\n",
      "Our analysis finds emissions which are 9.18 [ 9.18  ,  9.18 ] % of the verified credits in the PDD.\n",
      "We find that GS2439 project reports emissions which are 11.06 [ 11.01  ,  11.12 ] times overcredited.\n",
      "Our analysis finds emissions which are 9.04 [ 9.04  ,  9.04 ] % of the verified credits in the PDD.\n",
      "We find that GS2439 project reports emissions which are 8.94 [ 8.9  ,  8.98 ] times overcredited.\n",
      "Our analysis finds emissions which are 11.19 [ 11.19  ,  11.19 ] % of the verified credits in the PDD.\n",
      "We find that GS2439 project reports emissions which are 8.8 [ 8.76  ,  8.84 ] times overcredited.\n",
      "Our analysis finds emissions which are 11.36 [ 11.36  ,  11.36 ] % of the verified credits in the PDD.\n",
      "We find that GS2439 project reports emissions which are 9.43 [ 9.38  ,  9.47 ] times overcredited.\n",
      "Our analysis finds emissions which are 10.61 [ 10.61  ,  10.61 ] % of the verified credits in the PDD.\n",
      "We find that GS2445 project reports emissions which are 7.12 [ 7.08  ,  7.15 ] times overcredited.\n",
      "Our analysis finds emissions which are 14.05 [ 14.05  ,  14.05 ] % of the verified credits in the PDD.\n",
      "We find that GS7578 project reports emissions which are 8.16 [ 8.1  ,  8.23 ] times overcredited.\n",
      "Our analysis finds emissions which are 12.25 [ 12.25  ,  12.25 ] % of the verified credits in the PDD.\n",
      "We find that GS2513 project reports emissions which are 9.46 [ 9.38  ,  9.53 ] times overcredited.\n",
      "Our analysis finds emissions which are 10.58 [ 10.57  ,  10.58 ] % of the verified credits in the PDD.\n",
      "We find that GS4677 project reports emissions which are 3.76 [ 3.74  ,  3.79 ] times overcredited.\n",
      "Our analysis finds emissions which are 26.56 [ 26.56  ,  26.57 ] % of the verified credits in the PDD.\n",
      "We find that GS2441 project reports emissions which are 10.13 [ 10.05  ,  10.21 ] times overcredited.\n",
      "Our analysis finds emissions which are 9.87 [ 9.87  ,  9.87 ] % of the verified credits in the PDD.\n",
      "We find that GS2441 project reports emissions which are 10.32 [ 10.25  ,  10.4 ] times overcredited.\n",
      "Our analysis finds emissions which are 9.69 [ 9.69  ,  9.69 ] % of the verified credits in the PDD.\n",
      "We find that GS2441 project reports emissions which are 479.82 [ 476.73  ,  482.91 ] times overcredited.\n",
      "Our analysis finds emissions which are 0.21 [ 0.21  ,  0.21 ] % of the verified credits in the PDD.\n",
      "We find that GS2441 project reports emissions which are 480.0 [ 476.93  ,  483.06 ] times overcredited.\n",
      "Our analysis finds emissions which are 0.21 [ 0.21  ,  0.21 ] % of the verified credits in the PDD.\n",
      "We find that GS2441 project reports emissions which are 26.14 [ 25.98  ,  26.31 ] times overcredited.\n",
      "Our analysis finds emissions which are 3.83 [ 3.82  ,  3.83 ] % of the verified credits in the PDD.\n",
      "We find that GS2441 project reports emissions which are 26.13 [ 25.96  ,  26.3 ] times overcredited.\n",
      "Our analysis finds emissions which are 3.83 [ 3.83  ,  3.83 ] % of the verified credits in the PDD.\n",
      "We find that GS2441 project reports emissions which are 26.11 [ 25.94  ,  26.28 ] times overcredited.\n",
      "Our analysis finds emissions which are 3.83 [ 3.83  ,  3.83 ] % of the verified credits in the PDD.\n",
      "We find that GS2441 project reports emissions which are 26.19 [ 26.02  ,  26.36 ] times overcredited.\n",
      "Our analysis finds emissions which are 3.82 [ 3.82  ,  3.82 ] % of the verified credits in the PDD.\n",
      "We find that GS5003 project reports emissions which are 40.91 [ 40.71  ,  41.1 ] times overcredited.\n",
      "Our analysis finds emissions which are 2.44 [ 2.44  ,  2.44 ] % of the verified credits in the PDD.\n",
      "We find that GS5003 project reports emissions which are 5.4 [ 5.38  ,  5.43 ] times overcredited.\n",
      "Our analysis finds emissions which are 18.5 [ 18.5  ,  18.5 ] % of the verified credits in the PDD.\n",
      "We find that GS5003 project reports emissions which are 5.54 [ 5.51  ,  5.56 ] times overcredited.\n",
      "Our analysis finds emissions which are 18.07 [ 18.06  ,  18.07 ] % of the verified credits in the PDD.\n",
      "We find that GS6129 project reports emissions which are 82.74 [ 80.26  ,  85.23 ] times overcredited.\n",
      "Our analysis finds emissions which are 1.21 [ 1.21  ,  1.21 ] % of the verified credits in the PDD.\n",
      "We find that GS6129 project reports emissions which are 73.66 [ 71.46  ,  75.87 ] times overcredited.\n",
      "Our analysis finds emissions which are 1.36 [ 1.36  ,  1.36 ] % of the verified credits in the PDD.\n",
      "We find that GS6129 project reports emissions which are 74.59 [ 72.35  ,  76.83 ] times overcredited.\n",
      "Our analysis finds emissions which are 1.34 [ 1.34  ,  1.34 ] % of the verified credits in the PDD.\n",
      "We find that GS407 project reports emissions which are 10.92 [ 10.87  ,  10.97 ] times overcredited.\n",
      "Our analysis finds emissions which are 9.15 [ 9.15  ,  9.15 ] % of the verified credits in the PDD.\n",
      "We find that GS407 project reports emissions which are 9.88 [ 9.83  ,  9.93 ] times overcredited.\n",
      "Our analysis finds emissions which are 10.12 [ 10.12  ,  10.12 ] % of the verified credits in the PDD.\n",
      "We find that GS407 project reports emissions which are 0.64 [ 0.64  ,  0.64 ] times overcredited.\n",
      "Our analysis finds emissions which are 156.64 [ 156.64  ,  156.64 ] % of the verified credits in the PDD.\n",
      "We find that GS407 project reports emissions which are 0.64 [ 0.64  ,  0.64 ] times overcredited.\n",
      "Our analysis finds emissions which are 156.65 [ 156.65  ,  156.65 ] % of the verified credits in the PDD.\n",
      "We find that GS11509 project reports emissions which are 1.38 [ 1.38  ,  1.38 ] times overcredited.\n",
      "Our analysis finds emissions which are 72.52 [ 72.51  ,  72.52 ] % of the verified credits in the PDD.\n",
      "We find that GS11352 project reports emissions which are 2.0 [ 2.0  ,  2.0 ] times overcredited.\n",
      "Our analysis finds emissions which are 49.91 [ 49.91  ,  49.91 ] % of the verified credits in the PDD.\n",
      "We find that GS11507 project reports emissions which are 1.46 [ 1.46  ,  1.46 ] times overcredited.\n",
      "Our analysis finds emissions which are 68.48 [ 68.48  ,  68.48 ] % of the verified credits in the PDD.\n",
      "We find that GS1146 project reports emissions which are 13.26 [ 7.23  ,  19.28 ] times overcredited.\n",
      "Our analysis finds emissions which are 7.54 [ 7.54  ,  7.54 ] % of the verified credits in the PDD.\n",
      "We find that GS1146 project reports emissions which are 4.4 [ 4.37  ,  4.42 ] times overcredited.\n",
      "Our analysis finds emissions which are 22.75 [ 22.75  ,  22.75 ] % of the verified credits in the PDD.\n",
      "We find that GS1146 project reports emissions which are 9.89 [ 9.8  ,  9.98 ] times overcredited.\n",
      "Our analysis finds emissions which are 10.11 [ 10.11  ,  10.11 ] % of the verified credits in the PDD.\n",
      "We find that GS1146 project reports emissions which are 7.63 [ 1.03  ,  14.23 ] times overcredited.\n",
      "Our analysis finds emissions which are 13.1 [ 13.1  ,  13.1 ] % of the verified credits in the PDD.\n",
      "We find that GS 11330 project reports emissions which are 2.95 [ 2.93  ,  2.96 ] times overcredited.\n",
      "Our analysis finds emissions which are 33.92 [ 33.92  ,  33.92 ] % of the verified credits in the PDD.\n",
      "We find that GS 3071 project reports emissions which are 3.44 [ 3.42  ,  3.45 ] times overcredited.\n",
      "Our analysis finds emissions which are 29.09 [ 29.09  ,  29.09 ] % of the verified credits in the PDD.\n",
      "We find that GS 10777 project reports emissions which are 12.73 [ 12.67  ,  12.79 ] times overcredited.\n",
      "Our analysis finds emissions which are 7.86 [ 7.86  ,  7.86 ] % of the verified credits in the PDD.\n",
      "We find that GS 10886 project reports emissions which are 12.48 [ 12.42  ,  12.55 ] times overcredited.\n",
      "Our analysis finds emissions which are 8.01 [ 8.01  ,  8.01 ] % of the verified credits in the PDD.\n",
      "We find that GS 5107 project reports emissions which are 20.18 [ 20.08  ,  20.27 ] times overcredited.\n",
      "Our analysis finds emissions which are 4.96 [ 4.96  ,  4.96 ] % of the verified credits in the PDD.\n",
      "0.0004882125410621515\n",
      "We find that GS 11195 project reports emissions which are 20.46 [ 19.84  ,  21.07 ] times overcredited.\n",
      "Our analysis finds emissions which are 4.89 [ 4.89  ,  4.89 ] % of the verified credits in the PDD.\n",
      "0.00016341417667722595\n",
      "We find that GS 2077 project reports emissions which are 52.39 [ 50.82  ,  53.96 ] times overcredited.\n",
      "Our analysis finds emissions which are 1.91 [ 1.91  ,  1.91 ] % of the verified credits in the PDD.\n",
      "0.00038394506775862294\n",
      "We find that GS 411 project reports emissions which are 17.57 [ 17.04  ,  18.1 ] times overcredited.\n",
      "Our analysis finds emissions which are 5.69 [ 5.69  ,  5.69 ] % of the verified credits in the PDD.\n",
      "We find that GS 10914 project reports emissions which are 9.35 [ 9.07  ,  9.63 ] times overcredited.\n",
      "Our analysis finds emissions which are 10.69 [ 10.69  ,  10.69 ] % of the verified credits in the PDD.\n",
      "We find that GS 10781 project reports emissions which are 57.08 [ 55.33  ,  58.83 ] times overcredited.\n",
      "Our analysis finds emissions which are 1.75 [ 1.75  ,  1.75 ] % of the verified credits in the PDD.\n",
      "We find that GS 5660 project reports emissions which are 74.47 [ 72.23  ,  76.71 ] times overcredited.\n",
      "Our analysis finds emissions which are 1.34 [ 1.34  ,  1.34 ] % of the verified credits in the PDD.\n"
     ]
    }
   ],
   "source": [
    "i=1 # The row you want to evaluate\n",
    "while i<152:\n",
    "    # first row is units\n",
    "    VERs=VERs_all.iloc[i,1]\n",
    "    protocol_type=df['protocol_type'][i]\n",
    "\n",
    "    protocol_id=df['protocol_id'][i]\n",
    "\n",
    "    monitoring_period=df['monitoring_period'][i]\n",
    "\n",
    "    verified_credits=df['verified_credits'][i]\n",
    "\n",
    "    num_countries=df['num_countries'][i]\n",
    "\n",
    "    prov_project=df['prov_project'][i]\n",
    "\n",
    "    fNRB_project=df['fNRB_project'][i]\n",
    "\n",
    "    mfNRB_bailis=df['mfNRB_bailis'][i]\n",
    "\n",
    "    expfNRB_bailis=df['expfNRB_bailis'][i]\n",
    "\n",
    "    n_devices_domestic=df['n_devices_domestic'][i]\n",
    "\n",
    "    n_devices_commercial=df['n_devices_commercial'][i]\n",
    "\n",
    "    project_days=df['project_days'][i]\n",
    "    \n",
    "    project_days_commercial=df['project_days_commercial'][i]\n",
    "\n",
    "    av_hhsize=df['av_hhsize'][i]\n",
    "\n",
    "    \n",
    "    stove_year_cohort=df['stove_year_cohort'][i]\n",
    "\n",
    "    # Baseline Fuel % What are the baseline fuels in this project location? [put 0% if no use in this location]\n",
    "    prc_firewood_baseline=df['prc_firewood_baseline'][i]\n",
    "\n",
    "    prc_charcoal_baseline=df['prc_charcoal_baseline'][i]\n",
    "\n",
    "    prc_kerosene_baseline=df['prc_kerosene_baseline'][i]\n",
    "\n",
    "    prc_lpg_baseline=df['prc_lpg_baseline'][i]\n",
    "\n",
    "    prc_ethanol_baseline=df['prc_ethanol_baseline'][i]\n",
    "\n",
    "    prc_electric_baseline=df['prc_electric_baseline'][i]\n",
    "\n",
    "    prc_pellets_baseline=df['prc_pellets_baseline'][i]\n",
    "\n",
    "    \n",
    "\n",
    "\n",
    "\n",
    "    # Baseline Fuel % What are the baseline commercial fuels in this project location? [put 0% if no use in this location]\n",
    "    prc_firewood_baseline_commercial=df['prc_firewood_baseline_commercial'][i]\n",
    "\n",
    "    prc_charcoal_baseline_commercial=df['prc_charcoal_baseline_commercial'][i]\n",
    "\n",
    "    prc_kerosene_baseline_commercial=df['prc_kerosene_baseline_commercial'][i]\n",
    "\n",
    "    prc_lpg_baseline_commercial=df['prc_lpg_baseline_commercial'][i]\n",
    "\n",
    "    prc_ethanol_baseline_commercial=df['prc_ethanol_baseline_commercial'][i]\n",
    "\n",
    "    prc_electric_baseline_commercial=df['prc_electric_baseline_commercial'][i]\n",
    "\n",
    "    prc_pellets_baseline_commercial=df['prc_pellets_baseline_commercial'][i]\n",
    "\n",
    "    \n",
    "\n",
    "\n",
    "\n",
    "\n",
    "    # Baseline Consumption: What is the baseline fuel consumption of each baseline fuel in this project location [put 0 if no use in this location]\n",
    "    kg_firewood_baseline=df['kg_firewood_baseline'][i]\n",
    "\n",
    "    kg_charcoal_baseline=df['kg_charcoal_baseline'][i]\n",
    "\n",
    "    kg_kerosene_baseline=df['kg_kerosene_baseline'][i]\n",
    "\n",
    "    kg_lpg_baseline=df['kg_lpg_baseline'][i]\n",
    "\n",
    "    l_ethanol_baseline=df['l_ethanol_baseline'][i]\n",
    "\n",
    "    kg_per_liter_ethanol_baseline=df['kg_per_liter_ethanol_baseline'][i]\n",
    "\n",
    "    kWh_electric_baseline=df['kWh_electric_baseline'][i]\n",
    "\n",
    "    kg_pellets_baseline=df['kg_pellets_baseline'][i]\n",
    "\n",
    "\n",
    "    kg_firewood_baseline_commercial=df['kg_firewood_baseline_commercial'][i]\n",
    "\n",
    "    kg_charcoal_baseline_commercial=df['kg_charcoal_baseline_commercial'][i]\n",
    "\n",
    "    kg_kerosene_baseline_commercial=df['kg_kerosene_baseline_commercial'][i]\n",
    "\n",
    "    kg_lpg_baseline_commercial=df['kg_lpg_baseline_commercial'][i]\n",
    "\n",
    "    l_ethanol_baseline_commercial=df['l_ethanol_baseline_commercial'][i]\n",
    "\n",
    "    kWh_electric_baseline_commercial=df['kWh_electric_baseline_commercial'][i]\n",
    "\n",
    "    kg_pellets_baseline_commercial=df['kg_pellets_baseline_commercial'][i]\n",
    "\n",
    "\n",
    "    # Stove Efficiencies at Baseline\n",
    "\n",
    "\n",
    "    eff_firewood_baseline=df['eff_firewood_baseline'][i]\n",
    "\n",
    "    eff_charcoal_baseline=df['eff_charcoal_baseline'][i]\n",
    "\n",
    "    eff_kerosene_baseline=df['eff_kerosene_baseline'][i]\n",
    "\n",
    "    eff_LPG_baseline=df['eff_LPG_baseline'][i]\n",
    "\n",
    "    eff_ethanol_baseline=df['eff_ethanol_baseline'][i]\n",
    "\n",
    "    eff_electric_baseline=df['eff_electric_baseline'][i]\n",
    "\n",
    "    eff_pellets_baseline=df['eff_pellets_baseline'][i]\n",
    "\n",
    "\n",
    "    # Commercial Baseline Stove Efficiences\n",
    "\n",
    "    eff_firewood_baseline_commercial=df['eff_firewood_baseline_commercial'][i]\n",
    "\n",
    "    eff_charcoal_baseline_commercial=df['eff_charcoal_baseline_commercial'][i]\n",
    "\n",
    "    eff_kerosene_baseline_commercial=df['eff_kerosene_baseline_commercial'][i]\n",
    "\n",
    "    eff_LPG_baseline_commercial=df['eff_LPG_baseline_commercial'][i]\n",
    "\n",
    "    eff_ethanol_baseline_commercial=df['eff_ethanol_baseline_commercial'][i]\n",
    "\n",
    "    eff_electric_baseline_commercial=df['eff_electric_baseline_commercial'][i]\n",
    "\n",
    "    eff_pellet_baseline_commercial=df['eff_pellet_baseline_commercial'][i]\n",
    "\n",
    "\n",
    "    # Caloric Values: What are the assumed calorific values  of the baseline fuels? (TJ/ton) [put 0 if not in use in this location]\n",
    "\n",
    "    ncv_firewood_baseline=df['ncv_firewood_baseline'][i]\n",
    "\n",
    "    ncv_charcoal_baseline=df['ncv_charcoal_baseline'][i]\n",
    "\n",
    "    ncv_kerosene_baseline=df['ncv_kerosene_baseline'][i]\n",
    "\n",
    "    ncv_LPG_baseline=df['ncv_LPG_baseline'][i]\n",
    "\n",
    "    ncv_ethanol_baseline=df['ncv_ethanol_baseline'][i]\n",
    "\n",
    "    ncv_pellets_baseline=df['ncv_pellets_baseline'][i]\n",
    "\n",
    "   \n",
    "    # CO2 Emission Factors: What is the assumed carbon dioxide emission factor of the baseline fuels? \n",
    "    #[if CDM write in the assumed EF of the projected fossil fuel] \n",
    "    #[ put 0 if you do not account for CO2 EF at Baseline or this fuel is not applicable in your location]\n",
    "    co2_ef_firewood=df['co2_ef_firewood'][i]\n",
    "    co2_ef_charcoal=df['co2_ef_charcoal'][i]\n",
    "\n",
    "    co2_ef_kerosene=df['co2_ef_kerosene'][i]\n",
    "\n",
    "    co2_ef_LPG=df['co2_ef_LPG'][i]\n",
    "\n",
    "    co2_ef_ethanol=df['co2_ef_ethanol'][i]\n",
    "\n",
    "    co2_ef_electric=df['co2_ef_electric'][i]\n",
    "\n",
    "    co2_ef_pellets=df['co2_ef_pellets'][i]\n",
    "\n",
    "\n",
    "    #Non-CO2 Baseline Emission Factors \n",
    "\n",
    "    nonco2_ef_firewood=df['nonco2_ef_firewood'][i]\n",
    "    nonco2_ef_charcoal=df['nonco2_ef_charcoal'][i]\n",
    "\n",
    "    nonco2_ef_kerosene=df['nonco2_ef_kerosene'][i]\n",
    "\n",
    "    nonco2_ef_LPG=df['nonco2_ef_LPG'][i]\n",
    "\n",
    "    nonco2_ef_ethanol=df['nonco2_ef_ethanol'][i]\n",
    "\n",
    "    nonco2_ef_electric=df['nonco2_ef_electric'][i]\n",
    "\n",
    "    nonco2_ef_pellets=df['nonco2_ef_pellets'][i]\n",
    "\n",
    "\n",
    "    # Emissions from Stove Production\n",
    "\n",
    "    stoveproduction_ef_co2_firewood_baseline= df['stoveproduction_ef_co2_firewood_baseline'][i]\n",
    "    stoveproduction_ef_ch4_firewood_baseline= df['stoveproduction_ef_ch4_firewood_baseline'][i]\n",
    "    stoveproduction_ef_n2O_firewood_baseline= df['stoveproduction_ef_n2O_firewood_baseline'][i]\n",
    "\n",
    "\n",
    "    stoveproduction_ef_co2_charcoal_baseline=df['stoveproduction_ef_co2_charcoal_baseline'][i]\n",
    "    stoveproduction_ef_ch4_charcoal_baseline=df['stoveproduction_ef_ch4_charcoal_baseline'][i]\n",
    "    stoveproduction_ef_n2O_charcoal_baseline=df['stoveproduction_ef_n2O_charcoal_baseline'][i]\n",
    "\n",
    "\n",
    "    stoveproduction_ef_co2_kerosene_baseline=df['stoveproduction_ef_co2_kerosene_baseline'][i]\n",
    "    stoveproduction_ef_ch4_kerosene_baseline=df['stoveproduction_ef_ch4_kerosene_baseline'][i]\n",
    "    stoveproduction_ef_n2O_kerosene_baseline=df['stoveproduction_ef_n2O_kerosene_baseline'][i]\n",
    "\n",
    "\n",
    "    stoveproduction_ef_co2_LPG_baseline=df['stoveproduction_ef_co2_LPG_baseline'][i]\n",
    "    stoveproduction_ef_ch4_LPG_baseline=df['stoveproduction_ef_ch4_LPG_baseline'][i]\n",
    "    stoveproduction_ef_n2O_LPG_baseline=df['stoveproduction_ef_n2O_LPG_baseline'][i]\n",
    "\n",
    "\n",
    "    stoveproduction_ef_co2_ethanol_baseline=df['stoveproduction_ef_co2_ethanol_baseline'][i]\n",
    "    stoveproduction_ef_ch4_ethanol_baseline=df['stoveproduction_ef_ch4_ethanol_baseline'][i]\n",
    "    stoveproduction_ef_n2O_ethanol_baseline=df['stoveproduction_ef_n2O_ethanol_baseline'][i]\n",
    "\n",
    "\n",
    "    stoveproduction_ef_co2_electric_baseline=df['stoveproduction_ef_co2_electric_baseline'][i]\n",
    "    stoveproduction_ef_ch4_electric_baseline=df['stoveproduction_ef_ch4_electric_baseline'][i]\n",
    "    stoveproduction_ef_n2O_electric_baseline=df['stoveproduction_ef_n2O_electric_baseline'][i]\n",
    "\n",
    "\n",
    "    stoveproduction_ef_co2_pellets_baseline=df['stoveproduction_ef_co2_pellets_baseline'][i]\n",
    "    stoveproduction_ef_ch4_pellets_baseline=df['stoveproduction_ef_ch4_pellets_baseline'][i]\n",
    "    stoveproduction_ef_n2O_pellets_baseline=df['stoveproduction_ef_n2O_pellets_baseline'][i]\n",
    "\n",
    "\n",
    "    feedstockproduction_ef_co2_firewood_baseline=df['feedstockproduction_ef_co2_firewood_baseline'][i]\n",
    "    feedstockproduction_ef_ch4_firewood_baseline=df['feedstockproduction_ef_ch4_firewood_baseline'][i]\n",
    "    feedstockproduction_ef_n2O_firewood_baseline=df['feedstockproduction_ef_n2O_firewood_baseline'][i]\n",
    "\n",
    "\n",
    "    feedstockproduction_ef_co2_charcoal_baseline=df['feedstockproduction_ef_co2_charcoal_baseline'][i]\n",
    "    feedstockproduction_ef_ch4_charcoal_baseline=df['stoveproduction_ef_ch4_charcoal_baseline'][i]\n",
    "    feedstockproduction_ef_n2O_charcoal_baseline=df['feedstockproduction_ef_n2O_charcoal_baseline'][i]\n",
    "\n",
    "\n",
    "    feedstockproduction_ef_co2_kerosene_baseline=df['feedstockproduction_ef_co2_kerosene_baseline'][i]\n",
    "    feedstockproduction_ef_ch4_kerosene_baseline=df['feedstockproduction_ef_ch4_kerosene_baseline'][i]\n",
    "    feedstockproduction_ef_n2O_kerosene_baseline=df['feedstockproduction_ef_n2O_kerosene_baseline'][i]\n",
    "\n",
    "\n",
    "    feedstockproduction_ef_co2_LPG_baseline=df['feedstockproduction_ef_co2_LPG_baseline'][i]\n",
    "    feedstockproduction_ef_ch4_LPG_baseline=df['feedstockproduction_ef_ch4_LPG_baseline'][i]\n",
    "    feedstockproduction_ef_n2O_LPG_baseline=df['feedstockproduction_ef_n2O_LPG_baseline'][i]\n",
    "\n",
    "\n",
    "    feedstockproduction_ef_co2_ethanol_baseline=df['feedstockproduction_ef_co2_firewood_baseline'][i]\n",
    "    feedstockproduction_ef_ch4_ethanol_baseline=df['feedstockproduction_ef_ch4_firewood_baseline'][i]\n",
    "    feedstockproduction_ef_n2O_ethanol_baseline=df['feedstockproduction_ef_n2O_firewood_baseline'][i]\n",
    "\n",
    "\n",
    "    feedstockproduction_ef_co2_electric_baseline=df['feedstockproduction_ef_co2_firewood_baseline'][i]\n",
    "    feedstockproduction_ef_ch4_electric_baseline=df['feedstockproduction_ef_ch4_firewood_baseline'][i]\n",
    "    feedstockproduction_ef_n2O_electric_baseline=df['feedstockproduction_ef_n2O_firewood_baseline'][i]\n",
    "\n",
    "\n",
    "    feedstockproduction_ef_co2_pellets_baseline=df['feedstockproduction_ef_co2_firewood_baseline'][i]\n",
    "    feedstockproduction_ef_ch4_pellets_baseline=df['feedstockproduction_ef_ch4_firewood_baseline'][i]\n",
    "    feedstockproduction_ef_n2O_pellets_baseline=df['feedstockproduction_ef_n2O_firewood_baseline'][i]\n",
    "\n",
    "\n",
    "    distribution_ef_co2_firewood_baseline=df['distribution_ef_co2_firewood_baseline'][i]\n",
    "    distribution_ef_ch4_firewood_baseline=df['distribution_ef_ch4_firewood_baseline'][i]\n",
    "    distribution_ef_n2O_firewood_baseline=df['distribution_ef_n2O_firewood_baseline'][i]\n",
    "\n",
    "\n",
    "    distribution_ef_co2_charcoal_baseline=df['distribution_ef_co2_charcoal_baseline'][i]\n",
    "    distribution_ef_ch4_charcoal_baseline=df['distribution_ef_ch4_charcoal_baseline'][i]\n",
    "    distribution_ef_n2O_charcoal_baseline=df['distribution_ef_n2O_charcoal_baseline'][i]\n",
    "\n",
    "\n",
    "    distribution_ef_co2_kerosene_baseline=df['distribution_ef_co2_kerosene_baseline'][i]\n",
    "    distribution_ef_ch4_kerosene_baseline=df['distribution_ef_ch4_kerosene_baseline'][i]\n",
    "    distribution_ef_n2O_kerosene_baseline=df['distribution_ef_n2O_kerosene_baseline'][i]\n",
    "\n",
    "\n",
    "    distribution_ef_co2_LPG_baseline=df['distribution_ef_co2_LPG_baseline'][i]\n",
    "    distribution_ef_ch4_LPG_baseline=df['distribution_ef_ch4_LPG_baseline'][i]\n",
    "    distribution_ef_n2O_LPG_baseline=df['distribution_ef_n2O_LPG_baseline'][i]\n",
    "\n",
    "\n",
    "    distribution_ef_co2_ethanol_baseline=df['distribution_ef_co2_ethanol_baseline'][i]\n",
    "    distribution_ef_ch4_ethanol_baseline=df['distribution_ef_ch4_ethanol_baseline'][i]\n",
    "    distribution_ef_n2O_ethanol_baseline=df['distribution_ef_n2O_ethanol_baseline'][i]\n",
    "\n",
    "\n",
    "    distribution_ef_co2_electric_baseline=df['distribution_ef_co2_electric_baseline'][i]\n",
    "    distribution_ef_ch4_electric_baseline=df['distribution_ef_ch4_electric_baseline'][i]\n",
    "    distribution_ef_n2O_electric_baseline=df['distribution_ef_n2O_electric_baseline'][i]\n",
    "\n",
    "\n",
    "    distribution_ef_co2_pellets_baseline=df['distribution_ef_co2_pellets_baseline'][i]\n",
    "    distribution_ef_ch4_pellets_baseline=df['distribution_ef_ch4_pellets_baseline'][i]\n",
    "    distribution_ef_n2O_pellets_baseline=df['distribution_ef_n2O_pellets_baseline'][i]\n",
    "\n",
    "\n",
    "    processing_ef_co2_firewood_baseline=df['processing_ef_co2_firewood_baseline'][i]\n",
    "    processing_ef_ch4_firewood_baseline=df['processing_ef_ch4_firewood_baseline'][i]\n",
    "    processing_ef_n2O_firewood_baseline=df['processing_ef_n2O_firewood_baseline'][i]\n",
    "\n",
    "\n",
    "    processing_ef_co2_charcoal_baseline=df['processing_ef_co2_firewood_baseline'][i]\n",
    "    processing_ef_ch4_charcoal_baseline=df['processing_ef_co2_firewood_baseline'][i]\n",
    "    processing_ef_n2O_charcoal_baseline=df['processing_ef_co2_firewood_baseline'][i]\n",
    "\n",
    "\n",
    "    processing_ef_co2_kerosene_baseline=df['processing_ef_co2_kerosene_baseline'][i]\n",
    "    processing_ef_ch4_kerosene_baseline=df['processing_ef_ch4_kerosene_baseline'][i]\n",
    "    processing_ef_n2O_kerosene_baseline=df['processing_ef_n2O_kerosene_baseline'][i]\n",
    "\n",
    "\n",
    "    processing_ef_co2_LPG_baseline=df['processing_ef_co2_LPG_baseline'][i]\n",
    "    processing_ef_ch4_LPG_baseline=df['processing_ef_ch4_LPG_baseline'][i]\n",
    "    processing_ef_n2O_LPG_baseline=df['processing_ef_n2O_LPG_baseline'][i]\n",
    "\n",
    "\n",
    "    processing_ef_co2_ethanol_baseline=df['processing_ef_co2_ethanol_baseline'][i]\n",
    "    processing_ef_ch4_ethanol_baseline=df['processing_ef_ch4_ethanol_baseline'][i]\n",
    "    processing_ef_n2O_ethanol_baseline=df['processing_ef_n2O_ethanol_baseline'][i]\n",
    "\n",
    "\n",
    "    processing_ef_co2_electric_baseline=df['processing_ef_co2_electric_baseline'][i]\n",
    "    processing_ef_ch4_electric_baseline=df['processing_ef_ch4_electric_baseline'][i]\n",
    "    processing_ef_n2O_electric_baseline=df['processing_ef_n2O_electric_baseline'][i]\n",
    "\n",
    "\n",
    "    processing_ef_co2_pellets_baseline=df['processing_ef_co2_firewood_baseline'][i]\n",
    "    processing_ef_ch4_pellets_baseline=df['processing_ef_ch4_firewood_baseline'][i]\n",
    "    processing_ef_n2O_pellets_baseline=df['processing_ef_n2O_firewood_baseline'][i]\n",
    "\n",
    "    # Project Fuel % What are the project domestic fuels in this project location? [put 0% if no use in this location]\n",
    "    prc_firewood_project=df['prc_firewood_project'][i]\n",
    "\n",
    "    prc_charcoal_project=df['prc_charcoal_project'][i]\n",
    "\n",
    "    prc_kerosene_project=df['prc_kerosene_project'][i]\n",
    "\n",
    "    prc_lpg_project=df['prc_lpg_project'][i]\n",
    "\n",
    "    prc_ethanol_project=df['prc_ethanol_project'][i]\n",
    "\n",
    "    prc_electric_project=df['prc_electric_project'][i]\n",
    "\n",
    "    prc_pellets_project=df['prc_pellets_project'][i]\n",
    "\n",
    "    \n",
    "\n",
    "    # Project Fuel % What are the project commercial fuels in this project location? [put 0% if no use in this location]\n",
    "    prc_firewood_project_commercial=df['prc_firewood_project_commercial'][i]\n",
    "\n",
    "    prc_charcoal_project_commercial=df['prc_charcoal_project_commercial'][i]\n",
    "\n",
    "    prc_kerosene_project_commercial=df['prc_kerosene_project_commercial'][i]\n",
    "\n",
    "    prc_lpg_project_commercial=df['prc_lpg_project_commercial'][i]\n",
    "\n",
    "    prc_ethanol_project_commercial=df['prc_ethanol_project_commercial'][i]\n",
    "\n",
    "    prc_electric_project_commercial=df['prc_electric_project_commercial'][i]\n",
    "\n",
    "    prc_pellets_project_commercial=df['prc_pellets_project_commercial'][i]\n",
    "\n",
    "    \n",
    "\n",
    "    # Project Fuel Consumption\n",
    "\n",
    "    kg_firewood_project=df['kg_firewood_project'][i]\n",
    "\n",
    "    kg_charcoal_project=df['kg_charcoal_project'][i]\n",
    "\n",
    "    kg_kerosene_project=df['kg_kerosene_project'][i]\n",
    "\n",
    "    kg_lpg_project=df['kg_lpg_project'][i]\n",
    "\n",
    "    l_ethanol_project=df['l_ethanol_project'][i]\n",
    "\n",
    "\n",
    "    kg_per_liter_ethanol_project=df['kg_per_liter_ethanol_project'][i]\n",
    "\n",
    "    kWh_electric_project=df['kWh_electric_project'][i]\n",
    "\n",
    "    kg_pellets_project=df['kg_pellets_project'][i]\n",
    "\n",
    "\n",
    "    # Project Insitutional Consumption\n",
    "\n",
    "\n",
    "    kg_firewood_project_commercial=df['kg_firewood_project_commercial'][i]\n",
    "\n",
    "    kg_charcoal_project_commercial=df['kg_charcoal_project_commercial'][i]\n",
    "\n",
    "    kg_kerosene_project_commercial=df['kg_kerosene_project_commercial'][i]\n",
    "\n",
    "    kg_lpg_project_commercial=df['kg_lpg_project_commercial'][i]\n",
    "\n",
    "    l_ethanol_project_commercial=df['l_ethanol_project_commercial'][i]\n",
    "\n",
    "    kWh_electric_project_commercial=df['kWh_electric_project_commercial'][i]\n",
    "\n",
    "    kg_pellets_project_commercial=df['kg_pellets_project_commercial'][i]\n",
    "\n",
    "    \n",
    "\n",
    "\n",
    "    # Project Stove Efficiencies\n",
    "\n",
    "    eff_firewood_project=df['eff_firewood_project'][i]\n",
    "\n",
    "    eff_charcoal_project=df['eff_charcoal_project'][i]\n",
    "\n",
    "    eff_kerosene_project=df['eff_kerosene_project'][i]\n",
    "\n",
    "    eff_LPG_project=df['eff_LPG_project'][i]\n",
    "\n",
    "    eff_ethanol_project=df['eff_ethanol_project'][i]\n",
    "\n",
    "    eff_electric_project=df['eff_electric_project'][i]\n",
    "\n",
    "    eff_pellet_project=df['eff_pellet_project'][i]\n",
    "\n",
    "\n",
    "    # Project Instituitonal Stove Efficiencies\n",
    "\n",
    "    eff_firewood_project_commercial=df['eff_firewood_project_commercial'][i]\n",
    "\n",
    "    eff_charcoal_project_commercial=df['eff_charcoal_project_commercial'][i]\n",
    "\n",
    "    eff_kerosene_project_commercial=df['eff_kerosene_project_commercial'][i]\n",
    "\n",
    "    eff_LPG_project_commercial=df['eff_LPG_project_commercial'][i]\n",
    "\n",
    "    eff_ethanol_project_commercial=df['eff_ethanol_project_commercial'][i]\n",
    "\n",
    "    eff_electric_project_commercial=df['eff_electric_project_commercial'][i]\n",
    "\n",
    "    eff_pellet_project_commercial=df['eff_pellet_project_commercial'][i]\n",
    "\n",
    "\n",
    "    #Assumed Caloric Values in the Project Scenario\n",
    "\n",
    "    ncv_firewood_project=df['ncv_firewood_project'][i]\n",
    "\n",
    "    ncv_charcoal_project=df['ncv_charcoal_project'][i]\n",
    "\n",
    "    ncv_kerosene_project=df['ncv_kerosene_project'][i]\n",
    "\n",
    "    ncv_LPG_project=df['ncv_LPG_project'][i]\n",
    "\n",
    "    ncv_ethanol_project=df['ncv_ethanol_project'][i]\n",
    "\n",
    "    ncv_pellets_project=df['ncv_pellets_project'][i]\n",
    "\n",
    "    \n",
    "\n",
    "\n",
    "\n",
    "    \n",
    "\n",
    "\n",
    "    # CO2 Emission Factors in the Project\n",
    "\n",
    "    co2_ef_firewood_project=df['co2_ef_firewood_project'][i]\n",
    "\n",
    "    co2_ef_charcoal_project=df['co2_ef_charcoal_project'][i]\n",
    "\n",
    "    co2_ef_kerosene_project=df['co2_ef_kerosene_project'][i]\n",
    "\n",
    "    co2_ef_LPG_project=df['co2_ef_LPG_project'][i]\n",
    "\n",
    "    co2_ef_ethanol_project=df['co2_ef_ethanol_project'][i]\n",
    "\n",
    "    co2_ef_electric_project=df['co2_ef_electric_project'][i]\n",
    "\n",
    "    co2_ef_pellets_project=df['co2_ef_pellets_project'][i]\n",
    "\n",
    "\n",
    "    # Non-CO2 Emission Factors in the Project\n",
    "\n",
    "    nonco2_ef_firewood_project=df['nonco2_ef_firewood_project'][i]\n",
    "\n",
    "    nonco2_ef_charcoal_project=df['nonco2_ef_charcoal_project'][i]\n",
    "\n",
    "    nonco2_ef_kerosene_project=df['nonco2_ef_kerosene_project'][i]\n",
    "\n",
    "    nonco2_ef_LPG_project=df['nonco2_ef_LPG_project'][i]\n",
    "\n",
    "    nonco2_ef_ethanol_project=df['nonco2_ef_ethanol_project'][i]\n",
    "\n",
    "    nonco2_ef_electric_project=df['nonco2_ef_electric_project'][i]\n",
    "\n",
    "    nonco2_ef_pellets_project=df['nonco2_ef_pellets_project'][i]\n",
    "\n",
    "    # Stove Production Emission Factors\n",
    "\n",
    "\n",
    "    stoveproduction_ef_co2_firewood_project=df['stoveproduction_ef_co2_firewood_project'][i]\n",
    "    stoveproduction_ef_ch4_firewood_project=df['stoveproduction_ef_ch4_firewood_project'][i]\n",
    "    stoveproduction_ef_n2O_firewood_project=df['stoveproduction_ef_n2O_firewood_project'][i]\n",
    "\n",
    "\n",
    "    stoveproduction_ef_co2_charcoal_project=df['stoveproduction_ef_co2_charcoal_project'][i]\n",
    "    stoveproduction_ef_ch4_charcoal_project=df['stoveproduction_ef_ch4_charcoal_project'][i]\n",
    "    stoveproduction_ef_n2O_charcoal_project=df['stoveproduction_ef_n2O_charcoal_project'][i]\n",
    "\n",
    "\n",
    "    stoveproduction_ef_co2_kerosene_project=df['stoveproduction_ef_co2_kerosene_project'][i]\n",
    "    stoveproduction_ef_ch4_kerosene_project=df['stoveproduction_ef_ch4_kerosene_project'][i]\n",
    "    stoveproduction_ef_n2O_kerosene_project=df['stoveproduction_ef_n2O_kerosene_project'][i]\n",
    "\n",
    "\n",
    "    stoveproduction_ef_co2_LPG_project=df['stoveproduction_ef_co2_LPG_project'][i]\n",
    "    stoveproduction_ef_ch4_LPG_project=df['stoveproduction_ef_ch4_LPG_project'][i]\n",
    "    stoveproduction_ef_n2O_LPG_project=df['stoveproduction_ef_n2O_LPG_project'][i]\n",
    "\n",
    "\n",
    "    stoveproduction_ef_co2_ethanol_project=df['stoveproduction_ef_co2_ethanol_project'][i]\n",
    "    stoveproduction_ef_ch4_ethanol_project=df['stoveproduction_ef_ch4_ethanol_project'][i]\n",
    "    stoveproduction_ef_n2O_ethanol_project=df['stoveproduction_ef_n2O_ethanol_project'][i]\n",
    "\n",
    "\n",
    "    stoveproduction_ef_co2_electric_project=df['stoveproduction_ef_co2_electric_project'][i]\n",
    "    stoveproduction_ef_ch4_electric_project=df['stoveproduction_ef_ch4_electric_project'][i]\n",
    "    stoveproduction_ef_n2O_electric_project=df['stoveproduction_ef_n2O_electric_project'][i]\n",
    "\n",
    "\n",
    "    stoveproduction_ef_co2_pellets_project=df['stoveproduction_ef_co2_pellets_project'][i]\n",
    "    stoveproduction_ef_ch4_pellets_project=df['stoveproduction_ef_ch4_pellets_project'][i]\n",
    "    stoveproduction_ef_n2O_pellets_project=df['stoveproduction_ef_n2O_pellets_project'][i]\n",
    "\n",
    "    # Project Feedstock Emission Factors\n",
    "\n",
    "    feedstockproduction_ef_co2_firewood_project=df['feedstockproduction_ef_co2_firewood_project'][i]\n",
    "    feedstockproduction_ef_ch4_firewood_project=df['feedstockproduction_ef_ch4_firewood_project'][i]\n",
    "    feedstockproduction_ef_n2O_firewood_project=df['feedstockproduction_ef_n2O_firewood_project'][i]\n",
    "\n",
    "\n",
    "    feedstockproduction_ef_co2_charcoal_project=df['feedstockproduction_ef_co2_charcoal_project'][i]\n",
    "    feedstockproduction_ef_ch4_charcoal_project=df['feedstockproduction_ef_ch4_charcoal_project'][i]\n",
    "    feedstockproduction_ef_n2O_charcoal_project=df['feedstockproduction_ef_n2O_charcoal_project'][i]\n",
    "\n",
    "\n",
    "    feedstockproduction_ef_co2_kerosene_project=df['feedstockproduction_ef_co2_kerosene_project'][i]\n",
    "    feedstockproduction_ef_ch4_kerosene_project=df['feedstockproduction_ef_ch4_kerosene_project'][i]\n",
    "    feedstockproduction_ef_n2O_kerosene_project=df['feedstockproduction_ef_n2O_kerosene_project'][i]\n",
    "\n",
    "\n",
    "    feedstockproduction_ef_co2_LPG_project=df['feedstockproduction_ef_co2_LPG_project'][i]\n",
    "    feedstockproduction_ef_ch4_LPG_project=df['feedstockproduction_ef_ch4_LPG_project'][i]\n",
    "    feedstockproduction_ef_n2O_LPG_project=df['feedstockproduction_ef_n2O_LPG_project'][i]\n",
    "\n",
    "\n",
    "    feedstockproduction_ef_co2_ethanol_project=df['feedstockproduction_ef_co2_ethanol_project'][i]\n",
    "    feedstockproduction_ef_ch4_ethanol_project=df['feedstockproduction_ef_ch4_ethanol_project'][i]\n",
    "    feedstockproduction_ef_n2O_ethanol_project=df['feedstockproduction_ef_n2O_ethanol_project'][i]\n",
    "\n",
    "\n",
    "    feedstockproduction_ef_co2_electric_project=df['feedstockproduction_ef_co2_electric_project'][i]\n",
    "    feedstockproduction_ef_ch4_electric_project=df['feedstockproduction_ef_ch4_electric_project'][i]\n",
    "    feedstockproduction_ef_n2O_electric_project=df['feedstockproduction_ef_n2O_electric_project'][i]\n",
    "\n",
    "\n",
    "    feedstockproduction_ef_co2_pellets_project=df['feedstockproduction_ef_co2_pellets_project'][i]\n",
    "    feedstockproduction_ef_ch4_pellets_project=df['feedstockproduction_ef_ch4_pellets_project'][i]\n",
    "    feedstockproduction_ef_n2O_pellets_project=df['feedstockproduction_ef_n2O_pellets_project'][i]\n",
    "\n",
    "\n",
    "    # Project Distribution Emission Factors\n",
    "    distribution_ef_co2_firewood_project=df['distribution_ef_co2_firewood_project'][i]\n",
    "    distribution_ef_ch4_firewood_project=df['distribution_ef_ch4_firewood_project'][i]\n",
    "    distribution_ef_n2O_firewood_project=df['distribution_ef_n2O_firewood_project'][i]\n",
    "\n",
    "\n",
    "    distribution_ef_co2_charcoal_project=df['distribution_ef_co2_charcoal_project'][i]\n",
    "    distribution_ef_ch4_charcoal_project=df['distribution_ef_ch4_charcoal_project'][i]\n",
    "    distribution_ef_n2O_charcoal_project=df['distribution_ef_n2O_charcoal_project'][i]\n",
    "\n",
    "\n",
    "    distribution_ef_co2_kerosene_project=df['distribution_ef_co2_kerosene_project'][i]\n",
    "    distribution_ef_ch4_kerosene_project=df['distribution_ef_ch4_kerosene_project'][i]\n",
    "    distribution_ef_n2O_kerosene_project=df['distribution_ef_n2O_kerosene_project'][i]\n",
    "\n",
    "\n",
    "    distribution_ef_co2_LPG_project=df['distribution_ef_co2_LPG_project'][i]\n",
    "    distribution_ef_ch4_LPG_project=df['distribution_ef_ch4_LPG_project'][i]\n",
    "    distribution_ef_n2O_LPG_project=df['distribution_ef_n2O_LPG_project'][i]\n",
    "\n",
    "\n",
    "    distribution_ef_co2_ethanol_project=df['distribution_ef_co2_ethanol_project'][i]\n",
    "    distribution_ef_ch4_ethanol_project=df['distribution_ef_ch4_ethanol_project'][i]\n",
    "    distribution_ef_n2O_ethanol_project=df['distribution_ef_n2O_ethanol_project'][i]\n",
    "\n",
    "\n",
    "    distribution_ef_co2_electric_project=df['distribution_ef_co2_electric_project'][i]\n",
    "    distribution_ef_ch4_electric_project=df['distribution_ef_ch4_electric_project'][i]\n",
    "    distribution_ef_n2O_electric_project=df['distribution_ef_n2O_electric_project'][i]\n",
    "\n",
    "\n",
    "    distribution_ef_co2_pellets_project=df['distribution_ef_co2_pellets_project'][i]\n",
    "    distribution_ef_ch4_pellets_project=df['distribution_ef_ch4_pellets_project'][i]\n",
    "    distribution_ef_n2O_pellets_project=df['distribution_ef_n2O_pellets_project'][i]\n",
    "\n",
    "    # Project Processing Emission Factors\n",
    "    processing_ef_co2_firewood_project=df['processing_ef_co2_firewood_project'][i]\n",
    "    processing_ef_ch4_firewood_project=df['processing_ef_ch4_firewood_project'][i]\n",
    "    processing_ef_n2O_firewood_project=df['processing_ef_n2O_firewood_project'][i]\n",
    "\n",
    "\n",
    "    processing_ef_co2_charcoal_project=df['processing_ef_co2_charcoal_project'][i]\n",
    "    processing_ef_ch4_charcoal_project=df['processing_ef_ch4_charcoal_project'][i]\n",
    "    processing_ef_n2O_charcoal_project=df['processing_ef_n2O_charcoal_project'][i]\n",
    "\n",
    "\n",
    "    processing_ef_co2_kerosene_project=df['processing_ef_co2_kerosene_project'][i]\n",
    "    processing_ef_ch4_kerosene_project=df['processing_ef_ch4_kerosene_project'][i]\n",
    "    processing_ef_n2O_kerosene_project=df['processing_ef_n2O_kerosene_project'][i]\n",
    "\n",
    "\n",
    "    processing_ef_co2_LPG_project=df['processing_ef_co2_LPG_project'][i]\n",
    "    processing_ef_ch4_LPG_project=df['processing_ef_ch4_LPG_project'][i]\n",
    "    processing_ef_n2O_LPG_project=df['processing_ef_n2O_LPG_project'][i]\n",
    "\n",
    "\n",
    "    processing_ef_co2_ethanol_project=df['processing_ef_co2_ethanol_project'][i]\n",
    "    processing_ef_ch4_ethanol_project=df['processing_ef_ch4_ethanol_project'][i]\n",
    "    processing_ef_n2O_ethanol_project=df['processing_ef_n2O_ethanol_project'][i]\n",
    "\n",
    "\n",
    "    processing_ef_co2_electric_project=df['processing_ef_co2_electric_project'][i]\n",
    "    processing_ef_ch4_electric_project=df['processing_ef_ch4_electric_project'][i]\n",
    "    processing_ef_n2O_electric_project=df['processing_ef_n2O_electric_project'][i]\n",
    "\n",
    "\n",
    "    processing_ef_co2_pellets_project=df['processing_ef_co2_pellets_project'][i]\n",
    "    processing_ef_ch4_pellets_project=df['processing_ef_ch4_pellets_project'][i]\n",
    "    processing_ef_n2O_pellets_project=df['processing_ef_n2O_pellets_project'][i]\n",
    "    \n",
    "    \n",
    "    #fNRB\n",
    "    min_fNRB=df['mfNRB_bailis'][i]\n",
    "\n",
    "    exp_fNRB=df['expfNRB_bailis'][i]\n",
    "    max_fNRB=exp_fNRB*1.1\n",
    "\n",
    "    n_fNRB=np.random.triangular(left=min_fNRB, mode=exp_fNRB, right=max_fNRB, size=(N))\n",
    "    \n",
    "    # Adoption, Usage, and Leakage \n",
    "\n",
    "    #adoption_definition\n",
    "\n",
    "    project_adoption=df['project_adoption'][i]\n",
    "\n",
    "\n",
    "    project_usage=df['project_usage'][i]\n",
    "    \n",
    "    kpt_adj_usage=np.random.uniform(0,0.53,N)\n",
    "    \n",
    "    kpt_usage=project_usage*(1-kpt_adj_usage)\n",
    "\n",
    "\n",
    "    project_adoption_commercial=df['project_adoption_commercial'][i]\n",
    "\n",
    "\n",
    "    project_usage_commercial=df['project_usage_commercial'][i]\n",
    "    \n",
    "    project_stacking=df['project_stacking'][i]\n",
    "    \n",
    "    project_stacking_commercial=df['project_stacking_commercial'][i]\n",
    "    \n",
    "    kpt_adj_stacking=np.random.uniform(0,0.29,N)\n",
    "\n",
    "    kpt_stacking=project_stacking*kpt_adj_stacking\n",
    "\n",
    "    leakage_binary=df['leakage_binary'][i]\n",
    "\n",
    "   \n",
    "\n",
    "    leakage_default= df['leakage_default'][i] # put 1 if not being multiplied\n",
    "\n",
    "   \n",
    "\n",
    "    leakage_tCO2e=df['leakage_tCO2e'][i]\n",
    "\n",
    " \n",
    "\n",
    "    more_than_one_stove_adj=df['more_than_one_stove_adj'][i]\n",
    "    \n",
    "    other_tCO2e_reduction=df['other_tCO2e_reduction'][i]\n",
    "    \n",
    "    #rebound=1\n",
    "    \n",
    "    rebound=df['rebound'][i]\n",
    "\n",
    "\n",
    "    \n",
    "    # Baseline Emissions\n",
    "    # Baseline Firewood \n",
    "\n",
    "    be_firewood_gillwiehl=prc_firewood_baseline*(ncv_firewood_agw*kg_firewood_baseline/1000*((n_fNRB*(ef_combustion_co2_firewood\n",
    "                                                                                                      +ef_combustion_ch4_firewood*gwp_ch4_agw\n",
    "                                                                                                      + ef_combustion_n2O_firewood*gwp_n2O_agw) \n",
    "                                                                                              # non-renewable\n",
    "                                                                                              +(1-n_fNRB)*(ef_combustion_ch4_firewood*gwp_ch4_renewable\n",
    "                                                                                                      + ef_combustion_n2O_firewood*gwp_n2O_agw))))\n",
    "                                                                                             # renewable emissions use a different ch4 gwp\n",
    "        \n",
    "    # Baseline Charcoal\n",
    "\n",
    "    be_charcoal_gillwiehl=prc_charcoal_baseline*(ncv_charcoal_agw*kg_charcoal_baseline/1000*((n_fNRB*(ef_combustion_co2_charcoal\n",
    "                                                                                                      +ef_combustion_ch4_charcoal*gwp_ch4_agw\n",
    "                                                                                                      + ef_combustion_n2O_charcoal*gwp_n2O_agw) \n",
    "                                                                                              # non-renewable\n",
    "                                                                                              +(1-n_fNRB)*(ef_combustion_ch4_charcoal*gwp_ch4_renewable\n",
    "                                                                                                      + ef_combustion_n2O_charcoal*gwp_n2O_agw))\n",
    "                                                                                             # renewable emissions use a different ch4 gwp\n",
    "\n",
    "\n",
    "                                                                                             +ef_total_not_combustion_co2_charcoal  \n",
    "                                                                                             +ef_total_not_combustion_ch4_charcoal*gwp_ch4_agw\n",
    "                                                                                             +ef_total_not_combustion_n2O_charcoal*gwp_n2O_agw))\n",
    "    # Baseline LPG\n",
    "\n",
    "    be_LPG_gillwiehl=prc_lpg_baseline*(ncv_LPG_agw*kg_lpg_baseline/1000*(ef_total_co2_LPG +ef_total_ch4_LPG*gwp_ch4_agw\n",
    "                                                                                             +ef_total_n2O_LPG*gwp_n2O_agw))\n",
    "    # Baseline Kerosene\n",
    "\n",
    "    be_kerosene_gillwiehl=prc_kerosene_baseline*(ncv_kerosene_agw*kg_kerosene_baseline/1000*(ef_total_co2_kerosene +ef_total_ch4_kerosene*gwp_ch4_agw\n",
    "                                                                                             +ef_total_n2O_kerosene*gwp_n2O_agw))\n",
    "    # Baseline Ethanol\n",
    "\n",
    "    be_ethanol_gillwiehl=prc_ethanol_baseline*(ncv_ethanol_agw*l_ethanol_baseline_commercial*kg_per_liter_ethanol_baseline/1000*(ef_total_zeroed_co2_ethanol +ef_total_zeroed_ch4_ethanol*gwp_ch4_agw\n",
    "                                                                                             +ef_total_zeroed_n2O_ethanol*gwp_n2O_agw))\n",
    "    # Baseline Pellets\n",
    "\n",
    "    be_pellets_gillwiehl=prc_pellets_baseline*(ncv_pellets_agw*kg_pellets_baseline/1000*(ef_total_zeroed_co2_pellets +ef_total_zeroed_ch4_pellets*gwp_ch4_agw\n",
    "                                                                                             +ef_total_zeroed_n2O_pellets*gwp_n2O_agw))\n",
    "    # Total Baseline\n",
    "    be_total_domestic_gillwiehl = be_firewood_gillwiehl + be_charcoal_gillwiehl + be_LPG_gillwiehl+ be_ethanol_gillwiehl + be_pellets_gillwiehl\n",
    "\n",
    "    be_total_domestic_gillwiehl\n",
    "    \n",
    "   # Baseline Commercial Emissions\n",
    "    # Baseline Commercial Firewood \n",
    "\n",
    "    be_firewood_commercial_gillwiehl=prc_firewood_baseline_commercial*(ncv_firewood_agw*kg_firewood_baseline_commercial/1000*((n_fNRB*(ef_combustion_co2_firewood\n",
    "                                                                                                      +ef_combustion_ch4_firewood*gwp_ch4_agw\n",
    "                                                                                                      + ef_combustion_n2O_firewood*gwp_n2O_agw) \n",
    "                                                                                              # non-renewable\n",
    "                                                                                              +(1-n_fNRB)*(ef_combustion_ch4_firewood*gwp_ch4_renewable\n",
    "                                                                                                      + ef_combustion_n2O_firewood*gwp_n2O_agw))))\n",
    "                                                                                             # renewable emissions use a different ch4 gwp\n",
    "        \n",
    "    # Baseline Charcoal\n",
    "\n",
    "    be_charcoal_commercial_gillwiehl=prc_charcoal_baseline_commercial*(ncv_charcoal_agw*kg_charcoal_baseline_commercial/1000*((n_fNRB*(ef_combustion_co2_charcoal\n",
    "                                                                                                      +ef_combustion_ch4_charcoal*gwp_ch4_agw\n",
    "                                                                                                      + ef_combustion_n2O_charcoal*gwp_n2O_agw) \n",
    "                                                                                              # non-renewable\n",
    "                                                                                              +(1-n_fNRB)*(ef_combustion_ch4_charcoal*gwp_ch4_renewable\n",
    "                                                                                                      + ef_combustion_n2O_charcoal*gwp_n2O_agw))\n",
    "                                                                                             # renewable emissions use a different ch4 gwp\n",
    "\n",
    "\n",
    "                                                                                             +ef_total_not_combustion_co2_charcoal  \n",
    "                                                                                             +ef_total_not_combustion_ch4_charcoal*gwp_ch4_agw\n",
    "                                                                                             +ef_total_not_combustion_n2O_charcoal*gwp_n2O_agw))\n",
    "    # Baseline LPG\n",
    "\n",
    "    be_LPG_commercial_gillwiehl=prc_lpg_baseline_commercial*(ncv_LPG_agw*kg_lpg_baseline_commercial/1000*(ef_total_co2_LPG +ef_total_ch4_LPG*gwp_ch4_agw\n",
    "                                                                                             +ef_total_n2O_LPG*gwp_n2O_agw))\n",
    "    # Baseline Kerosene\n",
    "\n",
    "    be_kerosene_commercial_gillwiehl=prc_kerosene_baseline*(ncv_kerosene_agw*kg_kerosene_baseline_commercial/1000*(ef_total_co2_kerosene +ef_total_ch4_kerosene*gwp_ch4_agw\n",
    "                                                                                             +ef_total_n2O_kerosene*gwp_n2O_agw))\n",
    "    # Baseline Ethanol\n",
    "\n",
    "    be_ethanol_commercial_gillwiehl=prc_ethanol_baseline*(ncv_ethanol_agw*l_ethanol_baseline_commercial*kg_per_liter_ethanol_baseline/1000*(ef_total_zeroed_co2_ethanol +ef_total_zeroed_ch4_ethanol*gwp_ch4_agw\n",
    "                                                                                             +ef_total_zeroed_n2O_ethanol*gwp_n2O_agw))\n",
    "    # Baseline Ethanol\n",
    "\n",
    "    be_pellets_commercial_gillwiehl=prc_pellets_baseline*(ncv_pellets_agw*kg_pellets_baseline/1000*(ef_total_zeroed_co2_pellets +ef_total_zeroed_ch4_pellets*gwp_ch4_agw\n",
    "                                                                                             +ef_total_zeroed_n2O_pellets*gwp_n2O_agw))\n",
    "    # Total Baseline\n",
    "    be_total_commercial_gillwiehl = be_firewood_commercial_gillwiehl + be_charcoal_commercial_gillwiehl + be_LPG_commercial_gillwiehl+ be_ethanol_commercial_gillwiehl + be_pellets_commercial_gillwiehl\n",
    "\n",
    "    be_total_commercial_gillwiehl\n",
    "\n",
    "    # Project Emissions\n",
    "    # Project Firewood \n",
    "\n",
    "    pe_firewood_gillwiehl=prc_firewood_project*(ncv_firewood_agw*kg_firewood_project/1000*((n_fNRB*(ef_combustion_co2_firewood\n",
    "                                                                                                      +ef_combustion_ch4_firewood*gwp_ch4_agw\n",
    "                                                                                                      + ef_combustion_n2O_firewood*gwp_n2O_agw) \n",
    "                                                                                              # non-renewable\n",
    "                                                                                              +(1-n_fNRB)*(ef_combustion_ch4_firewood*gwp_ch4_renewable\n",
    "                                                                                                      + ef_combustion_n2O_firewood*gwp_n2O_agw))))\n",
    "                                                                                             # renewable emissions use a different ch4 gwp\n",
    "        \n",
    "    # Project Charcoal\n",
    "\n",
    "    pe_charcoal_gillwiehl=prc_charcoal_project*(ncv_charcoal_agw*kg_charcoal_project/1000*((n_fNRB*(ef_combustion_co2_charcoal\n",
    "                                                                                                      +ef_combustion_ch4_charcoal*gwp_ch4_agw\n",
    "                                                                                                      + ef_combustion_n2O_charcoal*gwp_n2O_agw) \n",
    "                                                                                              # non-renewable\n",
    "                                                                                              +(1-n_fNRB)*(ef_combustion_ch4_charcoal*gwp_ch4_renewable\n",
    "                                                                                                      + ef_combustion_n2O_charcoal*gwp_n2O_agw))\n",
    "                                                                                             # renewable emissions use a different ch4 gwp\n",
    "\n",
    "\n",
    "                                                                                             +ef_total_not_combustion_co2_charcoal  \n",
    "                                                                                             +ef_total_not_combustion_ch4_charcoal*gwp_ch4_agw\n",
    "                                                                                             +ef_total_not_combustion_n2O_charcoal*gwp_n2O_agw))\n",
    "    # Project LPG\n",
    "\n",
    "    pe_LPG_gillwiehl=prc_lpg_project*(ncv_LPG_agw*kg_lpg_project/1000*(ef_total_co2_LPG +ef_total_ch4_LPG*gwp_ch4_agw\n",
    "                                                                                             +ef_total_n2O_LPG*gwp_n2O_agw))\n",
    "    # Project Kerosene\n",
    "\n",
    "    pe_kerosene_gillwiehl=prc_kerosene_project*(ncv_kerosene_agw*kg_kerosene_project/1000*(ef_total_co2_kerosene +ef_total_ch4_kerosene*gwp_ch4_agw\n",
    "                                                                                             +ef_total_n2O_kerosene*gwp_n2O_agw))\n",
    "    # Project Ethanol\n",
    "\n",
    "    pe_ethanol_gillwiehl=prc_ethanol_project*(ncv_ethanol_agw*l_ethanol_project*kg_per_liter_ethanol_project/1000*(ef_total_zeroed_co2_ethanol +ef_total_zeroed_ch4_ethanol*gwp_ch4_agw\n",
    "                                                                                             +ef_total_zeroed_n2O_ethanol*gwp_n2O_agw))\n",
    "    # Project Pellets\n",
    "\n",
    "    pe_pellets_gillwiehl=prc_pellets_project*(ncv_pellets_agw*kg_pellets_project/1000*(ef_total_zeroed_co2_pellets +ef_total_zeroed_ch4_pellets*gwp_ch4_agw\n",
    "                                                                                             +ef_total_zeroed_n2O_pellets*gwp_n2O_agw))\n",
    "    # Total project\n",
    "    pe_total_domestic_gillwiehl = pe_firewood_gillwiehl + pe_charcoal_gillwiehl + pe_LPG_gillwiehl+ pe_ethanol_gillwiehl + pe_pellets_gillwiehl\n",
    "\n",
    "    pe_total_domestic_gillwiehl\n",
    "    \n",
    "    # Project Commercial Emissions\n",
    "    # Project Commercial Firewood \n",
    "\n",
    "    pe_firewood_commercial_gillwiehl=prc_firewood_project_commercial*(ncv_firewood_agw*kg_firewood_project_commercial/1000*((n_fNRB*(ef_combustion_co2_firewood\n",
    "                                                                                                      +ef_combustion_ch4_firewood*gwp_ch4_agw\n",
    "                                                                                                      + ef_combustion_n2O_firewood*gwp_n2O_agw) \n",
    "                                                                                              # non-renewable\n",
    "                                                                                              +(1-n_fNRB)*(ef_combustion_ch4_firewood*gwp_ch4_renewable\n",
    "                                                                                                      + ef_combustion_n2O_firewood*gwp_n2O_agw))))\n",
    "                                                                                             # renewable emissions use a different ch4 gwp\n",
    "        \n",
    "    # Project Charcoal\n",
    "\n",
    "    pe_charcoal_commercial_gillwiehl=prc_charcoal_project_commercial*(ncv_charcoal_agw*kg_charcoal_project_commercial/1000*((n_fNRB*(ef_combustion_co2_charcoal\n",
    "                                                                                                      +ef_combustion_ch4_charcoal*gwp_ch4_agw\n",
    "                                                                                                      + ef_combustion_n2O_charcoal*gwp_n2O_agw) \n",
    "                                                                                              # non-renewable\n",
    "                                                                                              +(1-n_fNRB)*(ef_combustion_ch4_charcoal*gwp_ch4_renewable\n",
    "                                                                                                      + ef_combustion_n2O_charcoal*gwp_n2O_agw))\n",
    "                                                                                             # renewable emissions use a different ch4 gwp\n",
    "\n",
    "\n",
    "                                                                                             +ef_total_not_combustion_co2_charcoal  \n",
    "                                                                                             +ef_total_not_combustion_ch4_charcoal*gwp_ch4_agw\n",
    "                                                                                             +ef_total_not_combustion_n2O_charcoal*gwp_n2O_agw))\n",
    "    # Project LPG\n",
    "\n",
    "    pe_LPG_commercial_gillwiehl=prc_lpg_project_commercial*(ncv_LPG_agw*kg_lpg_project_commercial/1000*(ef_total_co2_LPG +ef_total_ch4_LPG*gwp_ch4_agw\n",
    "                                                                                             +ef_total_n2O_LPG*gwp_n2O_agw))\n",
    "    # Project Kerosene\n",
    "\n",
    "    pe_kerosene_commercial_gillwiehl=prc_kerosene_project_commercial*(ncv_kerosene_agw*kg_kerosene_project_commercial/1000*(ef_total_co2_kerosene +ef_total_ch4_kerosene*gwp_ch4_agw\n",
    "                                                                                             +ef_total_n2O_kerosene*gwp_n2O_agw))\n",
    "    # Project Ethanol\n",
    "\n",
    "    pe_ethanol_commercial_gillwiehl=prc_ethanol_project_commercial*(ncv_ethanol_agw*l_ethanol_project_commercial*kg_per_liter_ethanol_project/1000*(ef_total_zeroed_co2_ethanol +ef_total_zeroed_ch4_ethanol*gwp_ch4_agw\n",
    "                                                                                             +ef_total_zeroed_n2O_ethanol*gwp_n2O_agw))\n",
    "    # Project Ethanol\n",
    "\n",
    "    pe_pellets_commercial_gillwiehl=prc_pellets_project_commercial*(ncv_pellets_agw*kg_pellets_project_commercial/1000*(ef_total_zeroed_co2_pellets +ef_total_zeroed_ch4_pellets*gwp_ch4_agw\n",
    "                                                                                             +ef_total_zeroed_n2O_pellets*gwp_n2O_agw))\n",
    "    # Total project\n",
    "    pe_total_commercial_gillwiehl = pe_firewood_commercial_gillwiehl + pe_charcoal_commercial_gillwiehl + pe_LPG_commercial_gillwiehl+ pe_ethanol_commercial_gillwiehl + pe_pellets_commercial_gillwiehl\n",
    "\n",
    "    pe_total_commercial_gillwiehl\n",
    "    \n",
    "\n",
    "    \n",
    "    \n",
    "    ## Adjusting adoption, usage, and stacking for the different methodologies \n",
    "    \n",
    "\n",
    "        \n",
    "    be_total_agw=be_total_domestic_gillwiehl\n",
    "\n",
    "    be_total_agw=be_total_agw+be_total_commercial_gillwiehl\n",
    "    \n",
    "    pe_total_agw=pe_total_domestic_gillwiehl\n",
    "    pe_total_agw=pe_total_agw+pe_total_commercial_gillwiehl\n",
    "    \n",
    "    #be_total_agw= be_total_agw*n_adoption*n_usage*n_stacking\n",
    "     \n",
    "    #pe_total_agw= be_total_agw*n_adoption*n_usage*n_stacking\n",
    "    \n",
    "    if df.protocol_type[i]=='GS Metered'and df.stove_year_cohort[i]!='commercial':\n",
    "    \n",
    "        be_total_agw= be_total_agw*project_adoption*project_usage*(1-project_stacking)\n",
    "     \n",
    "        pe_total_agw= pe_total_agw*project_adoption*project_usage*(1-project_stacking)\n",
    "        \n",
    "    if df.protocol_type[i]=='GS TPDDTEC'and df.protocol_id[i]!='GS2758'and df.stove_year_cohort[i]!='commercial':\n",
    "        \n",
    "        be_total_agw= be_total_agw*n_adoption*kpt_usage*(1-kpt_stacking)\n",
    "     \n",
    "        pe_total_agw= pe_total_agw*n_adoption*kpt_usage*(1-kpt_stacking)\n",
    "    \n",
    "    if df.protocol_type[i]=='GS TPDDTEC'and df.stove_year_cohort[i]=='commercial':\n",
    "        \n",
    "        be_total_agw= be_total_agw*project_adoption*project_usage*(1-project_stacking)\n",
    "     \n",
    "        pe_total_agw= pe_total_agw*project_adoption*project_usage*(1-project_stacking)\n",
    "          \n",
    "                                                             \n",
    "    if df.protocol_type[i]=='GS Simplified'and df.stove_year_cohort[i]!='commercial':\n",
    "        \n",
    "        be_total_agw= be_total_agw*n_adoption*n_usage*(1-n_stacking)\n",
    "     \n",
    "        pe_total_agw= pe_total_agw*n_adoption*n_usage*(1-n_stacking)\n",
    "        \n",
    "    if df.protocol_type[i]=='CDM AMS II G'and df.protocol_id[i]!='GS10974'and df.stove_year_cohort[i]!='commercial':\n",
    "        \n",
    "        \n",
    "        be_total_agw= be_total_agw*n_adoption*n_usage*(1-n_stacking)\n",
    "     \n",
    "        pe_total_agw= pe_total_agw*n_adoption*n_usage*(1-n_stacking)\n",
    "        \n",
    "    if df.protocol_type[i]=='CDM AMS I E'and df.stove_year_cohort[i]!='commercial': \n",
    "        \n",
    "        be_total_agw= be_total_agw*n_adoption*n_usage*(1-project_stacking)\n",
    "     \n",
    "        pe_total_agw= pe_total_agw*n_adoption*n_usage*(1-project_stacking)\n",
    "       \n",
    "    if df.protocol_id[i]=='GS2758': \n",
    "       \n",
    "        be_total_agw= be_total_agw*req_adoption*kpt_usage*(1-kpt_stacking)\n",
    "     \n",
    "        pe_total_agw= pe_total_agw*req_adoption*kpt_usage*(1-kpt_stacking)\n",
    "        \n",
    "        \n",
    "    if df.protocol_id[i]=='GS10974'and df.stove_year_cohort[i]!='commercial':\n",
    "    \n",
    "        be_total_agw= be_total_agw*req_adoption*req_usage*(1-req_stacking)\n",
    "     \n",
    "        pe_total_agw= pe_total_agw*req_adoption*req_usage*(1-req_stacking)\n",
    "        \n",
    "        \n",
    "    if df.protocol_type[i]=='CDM AMS I E'and df.stove_year_cohort[i]=='commercial': \n",
    "        \n",
    "        be_total_agw= be_total_agw*project_adoption*project_usage*(1-project_stacking)\n",
    "     \n",
    "        pe_total_agw= pe_total_agw*project_adoption*project_usage*(1-project_stacking)\n",
    "        \n",
    "        \n",
    "    if df.protocol_type[i]=='CDM AMS II G'and df.stove_year_cohort[i]=='commercial': \n",
    "        \n",
    "        be_total_agw= be_total_agw*project_adoption*project_usage*(1-project_stacking)\n",
    "     \n",
    "        pe_total_agw= pe_total_agw*project_adoption*project_usage*(1-project_stacking)\n",
    "        \n",
    "        \n",
    "    # Leakage if not default\n",
    "    le_total_agw=(leakage_tCO2e+other_tCO2e_reduction)/(n_devices_domestic+n_devices_commercial)/(project_days+project_days_commercial) # Project Emission Reductions (PERs) in (tCO2e)\n",
    "    \n",
    "        \n",
    "        \n",
    "    \n",
    "        \n",
    "    # Totally for all\n",
    "    PEs_agw=(pe_total_agw)# Project Emission Reductions (PERs) in (tCO2e)\n",
    "    BEs_agw=(be_total_agw)# Project Emission Reductions (PERs) in (tCO2e)\n",
    "\n",
    "    VERs_agw=rebound*(be_total_agw-pe_total_agw-le_total_agw)*more_than_one_stove_adj*leakage_default\n",
    "    \n",
    "  \n",
    "    # Protocol Analysis: \n",
    "    \n",
    "    \n",
    "    #All \n",
    "    \n",
    "    all_VERS_agw.append(VERs_agw)\n",
    "    \n",
    "    sum_VERS_not_per_stove_per_day_agw=sum_VERS_not_per_stove_per_day_agw+VERs_agw*n_devices_domestic*project_days\n",
    "    \n",
    "   \n",
    "        \n",
    "    #bangladesh\n",
    "    if df.country[i]=='Bangladesh':\n",
    "        bangladesh_N=bangladesh_N+1\n",
    "        bangladesh_VERs_agw.append(VERs_agw)\n",
    "        \n",
    "    #burkina faso\n",
    "    if df.country[i]=='Burkina Faso':\n",
    "        burkinafaso_N=burkinafaso_N+1\n",
    "        burkinafaso_VERs_agw.append(VERs_agw)\n",
    "        \n",
    "    #cameroon\n",
    "    if df.country[i]=='Cameroon':\n",
    "        cameroon_N=cameroon_N+1\n",
    "        cameroon_VERs_agw.append(VERs_agw)\n",
    "        \n",
    "    #china\n",
    "    if df.country[i]=='China':\n",
    "        china_N=china_N+1\n",
    "        china_VERs_agw.append(VERs_agw)\n",
    "        \n",
    "    #ethiopia\n",
    "    if df.country[i]=='Ethiopia':\n",
    "        ethiopia_N=ethiopia_N+1\n",
    "        ethiopia_VERs_agw.append(VERs_agw)\n",
    "        \n",
    "    #ghana\n",
    "    if df.country[i]=='Ghana':\n",
    "        ghana_N=ghana_N+1\n",
    "        ghana_VERs_agw.append(VERs_agw)\n",
    "        \n",
    "    #guatemala\n",
    "    if df.country[i]=='Guatemala':\n",
    "        guatemala_N=guatemala_N+1\n",
    "        guatemala_VERs_agw.append(VERs_agw)\n",
    "        \n",
    "    #haiti\n",
    "    if df.country[i]=='Haiti':\n",
    "        haiti_N=haiti_N+1\n",
    "        haiti_VERs_agw.append(VERs_agw)\n",
    "        \n",
    "    #honduras\n",
    "    if df.country[i]=='Honduras':\n",
    "        honduras_N=honduras_N+1\n",
    "        honduras_VERs_agw.append(VERs_agw)\n",
    "        \n",
    "    #india\n",
    "    if df.country[i]=='India':\n",
    "        india_N=india_N+1\n",
    "        india_VERs_agw.append(VERs_agw)\n",
    "        \n",
    "    #kenya\n",
    "    if df.country[i]=='Kenya':\n",
    "        kenya_N=kenya_N+1\n",
    "        kenya_VERs_agw.append(VERs_agw)\n",
    "        \n",
    "    #lesotho\n",
    "    if df.country[i]=='Lesotho':\n",
    "        lesotho_N=lesotho_N+1\n",
    "        lesotho_VERs_agw.append(VERs_agw)\n",
    "        \n",
    "    #malawi\n",
    "    if df.country[i]=='Malawi':\n",
    "        malawi_N=malawi_N+1\n",
    "        malawi_VERs_agw.append(VERs_agw)\n",
    "        \n",
    "    #mali\n",
    "    if df.country[i]=='Mali':\n",
    "        mali_N=mali_N+1\n",
    "        mali_VERs_agw.append(VERs_agw)\n",
    "        \n",
    "    #mexico\n",
    "    if df.country[i]=='Mexico':\n",
    "        mexico_N=mexico_N+1\n",
    "        mexico_VERs_agw.append(VERs_agw)\n",
    "        \n",
    "    #mozambique\n",
    "    if df.country[i]=='Mozambique':\n",
    "        mozambique_N=mozambique_N+1\n",
    "        mozambique_VERs_agw.append(VERs_agw)\n",
    "        \n",
    "    #myanmar\n",
    "    if df.country[i]=='Myanmar':\n",
    "        myanmar_N=myanmar_N+1\n",
    "        myanmar_VERs_agw.append(VERs_agw)\n",
    "        \n",
    "    #nepal\n",
    "    if df.country[i]=='Nepal':\n",
    "        nepal_N=nepal_N+1\n",
    "        nepal_VERs_agw.append(VERs_agw)\n",
    "        \n",
    "    #nigeria\n",
    "    if df.country[i]=='Nigeria':\n",
    "        nigeria_N=nigeria_N+1\n",
    "        nigeria_VERs_agw.append(VERs_agw)\n",
    "    \n",
    "    #peru\n",
    "    if df.country[i]=='Peru':\n",
    "        peru_N=peru_N+1\n",
    "        peru_VERs_agw.append(VERs_agw)\n",
    "        \n",
    "    #rwanda\n",
    "    if df.country[i]=='Rwanda':\n",
    "        rwanda_N=rwanda_N+1\n",
    "        rwanda_VERs_agw.append(VERs_agw)\n",
    "    \n",
    "    #sudan\n",
    "    if df.country[i]=='Sudan':\n",
    "        sudan_N=sudan_N+1\n",
    "        sudan_VERs_agw.append(VERs_agw)\n",
    "        \n",
    "    #uganda\n",
    "    if df.country[i]=='Uganda':\n",
    "        uganda_N=uganda_N+1\n",
    "        uganda_VERs_agw.append(VERs_agw)\n",
    "    \n",
    "    #zambia\n",
    "    if df.country[i]=='Zambia':\n",
    "        zambia_N=zambia_N+1\n",
    "        zambia_VERs_agw.append(VERs_agw)\n",
    "        \n",
    "    #gs metered pellets\n",
    "    if df.protocol_type[i]=='GS Metered':\n",
    "        GS_Meter_N=GS_Meter_N+1\n",
    "        GS_Meter_VERs_agw.append(VERs_agw)\n",
    "    \n",
    "    #gs_lpg\n",
    "    if df.stove[i]=='gs_lpg':\n",
    "        gs_lpg_N=gs_lpg_N+1\n",
    "        gs_lpg_VERs_agw.append(VERs_agw)\n",
    "        \n",
    "    #gs_imp_firewood\n",
    "    if df.stove[i]=='gs_imp_firewood':\n",
    "        gs_imp_firewood_N=gs_imp_firewood_N+1\n",
    "        gs_imp_firewood_VERs_agw.append(VERs_agw)\n",
    "    \n",
    "    #gs_simple_firewood\n",
    "    if df.stove[i]=='gs_simple_firewood':\n",
    "        gs_simple_firewood_N=gs_simple_firewood_N+1\n",
    "        gs_simple_firewood_VERs_agw.append(VERs_agw)\n",
    "    \n",
    "        \n",
    "    #gs_imp_charcoal\n",
    "    if df.stove[i]=='gs_imp_charcoal':\n",
    "        gs_imp_charcoal_N=gs_imp_charcoal_N+1\n",
    "        gs_imp_charcoal_VERs_agw.append(VERs_agw)\n",
    "        \n",
    "    #cdm_imp_firewood\n",
    "    if df.stove[i]=='cdm_imp_firewood':\n",
    "        cdm_imp_firewood_N=cdm_imp_firewood_N+1\n",
    "        cdm_imp_firewood_VERs_agw.append(VERs_agw)\n",
    "        print(np.mean(VERs_agw))\n",
    "        \n",
    "        \n",
    "    #cdm_imp_charcoal\n",
    "    if df.stove[i]=='cdm_imp_charcoal':\n",
    "        cdm_imp_charcoal_N=cdm_imp_charcoal_N+1\n",
    "        cdm_imp_charcoal_VERs_agw.append(VERs_agw)\n",
    "        \n",
    "    # CDM I E\n",
    "    if df.protocol_type[i]=='CDM AMS I E':\n",
    "        CDM_I_E_N=CDM_I_E_N+1\n",
    "        CDM_I_E_VERs_agw.append(VERs_agw)\n",
    "    \n",
    "    \n",
    "    #Times Over Credited \n",
    "    ratio_analysis_project=VERs_agw/VERs\n",
    "    ratio_analysis_project_mean=np.mean(VERs_agw)/VERs*100\n",
    "    ratio_analysis_project_mean \n",
    "    ratio_analysis_project_mean_compilation.append(ratio_analysis_project_mean)\n",
    "\n",
    "    \n",
    "    ratio_project_analysis=VERs/VERs_agw\n",
    "    ratio_project_analysis_mean=VERs/np.mean(VERs_agw)\n",
    "    ratio_project_analysis_mean\n",
    "    ratio_project_analysis_mean_compilation.append(ratio_project_analysis_mean)\n",
    "    z=1.96 \n",
    "\n",
    "    sd_VERs_agw=np.std(VERs_agw)\n",
    "    sd_ratio_analysis_project=np.std(ratio_analysis_project)\n",
    "    sd_ratio_project_analysis=np.std(ratio_project_analysis)\n",
    "    \n",
    "    \n",
    "    ratio_analysis_project_lower=ratio_analysis_project_mean - z*sd_ratio_analysis_project/np.sqrt(N)\n",
    "    ratio_analysis_project_upper=ratio_analysis_project_mean + z*sd_ratio_analysis_project/np.sqrt(N)\n",
    "\n",
    "    ratio_analysis_project_lower_compilation.append(ratio_analysis_project_lower)\n",
    "    ratio_analysis_project_upper_compilation.append(ratio_analysis_project_upper)\n",
    "    \n",
    "    \n",
    "    ratio_project_analysis_lower=ratio_project_analysis_mean - z*sd_ratio_project_analysis/np.sqrt(N)\n",
    "    ratio_project_analysis_upper=ratio_project_analysis_mean + z*sd_ratio_project_analysis/np.sqrt(N)\n",
    "    \n",
    "    #print(sd_ratio_project_analysis)\n",
    "    \n",
    "    ratio_project_analysis_lower_compilation.append(ratio_project_analysis_lower)\n",
    "    ratio_project_analysis_upper_compilation.append(ratio_project_analysis_upper)\n",
    "    \n",
    "    #if ratio_project_analysis_mean<2:\n",
    "    print('We find that', protocol_id,  'project reports emissions which are',round(ratio_project_analysis_mean,2),'[',round(ratio_project_analysis_lower,2),' , ',round(ratio_project_analysis_upper,2),']','times overcredited.')\n",
    "    print('Our analysis finds emissions which are',round(ratio_analysis_project_mean,2),'[',round(ratio_analysis_project_lower,2),' , ',round(ratio_analysis_project_upper,2),']','% of the verified credits in the PDD.')\n",
    "    \n",
    "    i=i+1\n",
    "    ALL_N=ALL_N+1"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "0    29.95388\n",
      "dtype: float64\n",
      "9.89793418458745\n",
      "9.89793418458745 [ 7.40286761840291 , 12.393000750771991 ]\n"
     ]
    }
   ],
   "source": [
    "ratio_project_analysis_mean_compilation=pd.DataFrame(ratio_project_analysis_mean_compilation)\n",
    "ratio_project_analysis_lower_compilation=pd.DataFrame(ratio_project_analysis_lower_compilation)\n",
    "ratio_project_analysis_upper_compilation=pd.DataFrame(ratio_project_analysis_upper_compilation)\n",
    "ratio_analysis_project_mean_compilation=pd.DataFrame(ratio_analysis_project_mean_compilation)\n",
    "ratio_analysis_project_lower_compilation=pd.DataFrame(ratio_analysis_project_lower_compilation)\n",
    "ratio_analysis_project_upper_compilation=pd.DataFrame(ratio_analysis_project_upper_compilation)\n",
    "ids=pd.DataFrame(df.iloc[:,2])\n",
    "category=pd.DataFrame(df.iloc[:,10])\n",
    "country=pd.DataFrame(df.iloc[:,13])\n",
    "total_credits=pd.DataFrame(df.iloc[:,5])\n",
    "stove_days=pd.DataFrame(df.iloc[:,23])\n",
    "\n",
    "\n",
    "df_info=ids.join(category)\n",
    "df_info=df_info.join(country)\n",
    "df_info=df_info.join(total_credits)\n",
    "df_info=df_info.join(stove_days)\n",
    "df_info=pd.DataFrame(df_info)\n",
    "df_info=pd.concat([df_info[1:].reset_index(drop=True),ratio_project_analysis_lower_compilation], axis=1)\n",
    "df_info= pd.concat([df_info.reset_index(drop=True),ratio_project_analysis_mean_compilation], axis=1)\n",
    "df_info= pd.concat([df_info.reset_index(drop=True),ratio_project_analysis_upper_compilation], axis=1)\n",
    "df_info.columns = ['protocol_id', 'stove', 'country','total_verified_credits','stove_days','lower','mean','upper']\n",
    "df_info['stove_days_mean_overcrediting']=df_info['mean']*df_info['stove_days']\n",
    "df_info\n",
    "df_project_weighted_stove_days_numerator=df_info.groupby('protocol_id')['stove_days_mean_overcrediting'].sum()\n",
    "df_project_weighted_stove_days_numerator=pd.DataFrame(df_project_weighted_stove_days_numerator)\n",
    "\n",
    "\n",
    "df_project_weighted_stove_days_denominator=df_info.groupby('protocol_id')['stove_days'].sum()\n",
    "df_project_weighted_stove_days_denominator=pd.DataFrame(df_project_weighted_stove_days_denominator)\n",
    "df_project_weighted_stove_days=df_project_weighted_stove_days_numerator['stove_days_mean_overcrediting']/df_project_weighted_stove_days_denominator['stove_days']\n",
    "df_project_weighted_stove_days=pd.DataFrame(df_project_weighted_stove_days)\n",
    "df_project_weighted_stove_days.mean() #This is the average project over-crediting\n",
    "print(df_project_weighted_stove_days.mean())\n",
    "\n",
    "df_project_weighted_stove_days\n",
    "protocol_ids=pd.DataFrame(df_info['protocol_id'].unique())\n",
    "protocol_ids.columns=['protocol_id']\n",
    "\n",
    "df_total=df_project_weighted_stove_days.merge(df_info,left_on='protocol_id', right_on='protocol_id',how='left')\n",
    "df_total.columns=['protocol_id','project_over_crediting','category','country','total_verified_credits','stove_days','lower','mean','upper','stove_days_mean_overcrediting']\n",
    "\n",
    "df_total['over_credited_total']=df_total['total_verified_credits']/df_total['project_over_crediting']\n",
    "total_verified_credits=df_total['total_verified_credits'].unique().sum()\n",
    "total_over_crediting =df_total['over_credited_total'].unique().sum()\n",
    "reported_value=total_verified_credits/total_over_crediting\n",
    "print(total_verified_credits/total_over_crediting)\n",
    "sd=df_total['total_verified_credits'].unique().std()/df_total['over_credited_total'].unique().std()\n",
    "ci=sd*1.96/np.sqrt(51)\n",
    "ci\n",
    "ci_lower=reported_value-ci \n",
    "ci_upper=reported_value+ci\n",
    "print(reported_value,'[',ci_lower,',',ci_upper,']')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "7\n"
     ]
    }
   ],
   "source": [
    "df_gs_lpg_VERs_agw=pd.DataFrame(gs_lpg_VERs_agw).transpose()\n",
    "print(df_gs_lpg_VERs_agw.shape[1])\n",
    "gs_lpg_i=0\n",
    "df_gs_lpg_VERs_agw_final=pd.DataFrame(np.zeros((1, 10000)))\n",
    "while gs_lpg_i<df_gs_lpg_VERs_agw.shape[1]:\n",
    "    df_gs_lpg_VERs_agw_final=df_gs_lpg_VERs_agw_final+df_gs_lpg_VERs_agw[gs_lpg_i]\n",
    "    gs_lpg_i=gs_lpg_i+1"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "3\n"
     ]
    }
   ],
   "source": [
    "df_GS_Meter_VERs_agw=pd.DataFrame(GS_Meter_VERs_agw).transpose()\n",
    "print(df_GS_Meter_VERs_agw.shape[1])\n",
    "gs_meter_i=0\n",
    "df_GS_Meter_VERs_agw_final=pd.DataFrame(np.zeros((1, 10000)))\n",
    "while gs_meter_i<df_GS_Meter_VERs_agw.shape[1]:\n",
    "    df_GS_Meter_VERs_agw_final=df_GS_Meter_VERs_agw_final+df_GS_Meter_VERs_agw[gs_meter_i]\n",
    "    gs_meter_i=gs_meter_i+1"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "24\n"
     ]
    }
   ],
   "source": [
    "df_gs_simple_firewood_VERs_agw=pd.DataFrame(gs_simple_firewood_VERs_agw).transpose()\n",
    "print(df_gs_simple_firewood_VERs_agw.shape[1])\n",
    "gs_simple_firewood_i=0\n",
    "df_gs_simple_firewood_VERs_agw_final=pd.DataFrame(np.zeros((1, 10000)))\n",
    "while gs_simple_firewood_i<df_gs_simple_firewood_VERs_agw.shape[1]:\n",
    "    df_gs_simple_firewood_VERs_agw_final=df_gs_simple_firewood_VERs_agw_final+df_gs_simple_firewood_VERs_agw[gs_simple_firewood_i]\n",
    "    gs_simple_firewood_i=gs_simple_firewood_i+1"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "31\n"
     ]
    }
   ],
   "source": [
    "df_gs_imp_firewood_VERs_agw=pd.DataFrame(gs_imp_firewood_VERs_agw).transpose()\n",
    "print(df_gs_imp_firewood_VERs_agw.shape[1])\n",
    "gs_imp_firewood_i=0\n",
    "df_gs_imp_firewood_VERs_agw_final=pd.DataFrame(np.zeros((1, 10000)))\n",
    "while gs_imp_firewood_i<df_gs_imp_firewood_VERs_agw.shape[1]:\n",
    "    df_gs_imp_firewood_VERs_agw_final=df_gs_imp_firewood_VERs_agw_final+df_gs_imp_firewood_VERs_agw[gs_imp_firewood_i]\n",
    "    gs_imp_firewood_i=gs_imp_firewood_i+1"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "15\n"
     ]
    }
   ],
   "source": [
    "df_gs_imp_charcoal_VERs_agw=pd.DataFrame(gs_imp_charcoal_VERs_agw).transpose()\n",
    "print(df_gs_imp_charcoal_VERs_agw.shape[1])\n",
    "gs_imp_charcoal_i=0\n",
    "df_gs_imp_charcoal_VERs_agw_final=pd.DataFrame(np.zeros((1, 10000)))\n",
    "while gs_imp_charcoal_i<df_gs_imp_charcoal_VERs_agw.shape[1]:\n",
    "    df_gs_imp_charcoal_VERs_agw_final=df_gs_imp_charcoal_VERs_agw_final+df_gs_imp_charcoal_VERs_agw[gs_imp_charcoal_i]\n",
    "    gs_imp_charcoal_i=gs_imp_charcoal_i+1"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "19\n"
     ]
    }
   ],
   "source": [
    "df_cdm_imp_charcoal_VERs_agw=pd.DataFrame(cdm_imp_charcoal_VERs_agw).transpose()\n",
    "print(df_cdm_imp_charcoal_VERs_agw.shape[1])\n",
    "cdm_imp_charcoal_i=0\n",
    "df_cdm_imp_charcoal_VERs_agw_final=pd.DataFrame(np.zeros((1, 10000)))\n",
    "while cdm_imp_charcoal_i<df_cdm_imp_charcoal_VERs_agw.shape[1]:\n",
    "    df_cdm_imp_charcoal_VERs_agw_final=df_cdm_imp_charcoal_VERs_agw_final+df_cdm_imp_charcoal_VERs_agw[cdm_imp_charcoal_i]\n",
    "    cdm_imp_charcoal_i=cdm_imp_charcoal_i+1"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "47\n"
     ]
    }
   ],
   "source": [
    "df_cdm_imp_firewood_VERs_agw=pd.DataFrame(cdm_imp_firewood_VERs_agw).transpose()\n",
    "print(df_cdm_imp_firewood_VERs_agw.shape[1])\n",
    "cdm_imp_firewood_i=0\n",
    "df_cdm_imp_firewood_VERs_agw_final=pd.DataFrame(np.zeros((1, 10000)))\n",
    "while cdm_imp_firewood_i<df_cdm_imp_firewood_VERs_agw.shape[1]:\n",
    "    df_cdm_imp_firewood_VERs_agw_final=df_cdm_imp_firewood_VERs_agw_final+df_cdm_imp_firewood_VERs_agw[cdm_imp_firewood_i]\n",
    "    cdm_imp_firewood_i=cdm_imp_firewood_i+1"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "5\n",
      "1\n",
      "2\n",
      "3\n",
      "4\n",
      "5\n"
     ]
    }
   ],
   "source": [
    "df_CDM_I_E_VERs_agw=pd.DataFrame(CDM_I_E_VERs_agw).transpose()\n",
    "print(df_CDM_I_E_VERs_agw.shape[1])\n",
    "cdm_i_e_i=0\n",
    "df_CDM_I_E_VERs_agw_final=pd.DataFrame(np.zeros((1, 10000)))\n",
    "while cdm_i_e_i<df_CDM_I_E_VERs_agw.shape[1]:\n",
    "    df_CDM_I_E_VERs_agw_final=df_CDM_I_E_VERs_agw_final+df_CDM_I_E_VERs_agw[cdm_i_e_i]\n",
    "    cdm_i_e_i=cdm_i_e_i+1\n",
    "    print(cdm_i_e_i)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 21,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "0    0.106229\n",
      "dtype: float64\n",
      "0    0.116393\n",
      "dtype: float64\n",
      "0    0.126557\n",
      "dtype: float64\n",
      "0    9.413584\n",
      "dtype: float64\n",
      "0    8.591577\n",
      "dtype: float64\n",
      "0    7.901599\n",
      "dtype: float64\n",
      "0.12 [ 0.11 , 0.13 ]\n",
      "8.59 [ 7.9 , 9.41 ]\n"
     ]
    }
   ],
   "source": [
    "\n",
    "gs_lpg_ratio_analysis_project=df_gs_lpg_VERs_agw_final/gs_lpg_VERs\n",
    "gs_lpg_ratio_analysis_project_mean=np.mean(df_gs_lpg_VERs_agw_final, axis=1)/gs_lpg_VERs\n",
    "gs_lpg_ratio_project_analysis_mean=gs_lpg_ratio_analysis_project_mean**-1\n",
    "\n",
    "gs_lpg_sd_ratio_analysis_project=np.std(gs_lpg_ratio_analysis_project, axis=1)\n",
    "\n",
    "\n",
    "\n",
    "gs_lpg_ratio_analysis_project_lower=gs_lpg_ratio_analysis_project_mean - z*gs_lpg_sd_ratio_analysis_project/np.sqrt(gs_lpg_N)\n",
    "gs_lpg_ratio_analysis_project_upper=gs_lpg_ratio_analysis_project_mean + z*gs_lpg_sd_ratio_analysis_project/np.sqrt(gs_lpg_N)\n",
    "\n",
    "gs_lpg_ratio_project_analysis_lower=gs_lpg_ratio_analysis_project_lower**-1\n",
    "gs_lpg_ratio_project_analysis_upper=gs_lpg_ratio_analysis_project_upper**-1\n",
    "\n",
    "print(gs_lpg_ratio_analysis_project_lower)\n",
    "print(gs_lpg_ratio_analysis_project_mean)\n",
    "print(gs_lpg_ratio_analysis_project_upper)\n",
    "\n",
    "print(gs_lpg_ratio_project_analysis_lower)\n",
    "print(gs_lpg_ratio_project_analysis_mean)\n",
    "print(gs_lpg_ratio_project_analysis_upper)\n",
    "\n",
    "print(round(gs_lpg_ratio_analysis_project_mean[0],2),'[',round(gs_lpg_ratio_analysis_project_lower[0],2),',',round(gs_lpg_ratio_analysis_project_upper[0],2),']')\n",
    "print(round(gs_lpg_ratio_project_analysis_mean[0],2),'[',round(gs_lpg_ratio_project_analysis_upper[0],2),',',round(gs_lpg_ratio_project_analysis_lower[0],2),']')\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 22,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "0    0.10537\n",
      "dtype: float64\n",
      "0    0.109263\n",
      "dtype: float64\n",
      "0    0.113156\n",
      "dtype: float64\n",
      "0    9.490336\n",
      "dtype: float64\n",
      "0    9.152206\n",
      "dtype: float64\n",
      "0    8.837341\n",
      "dtype: float64\n",
      "0.11 [ 0.11 , 0.11 ]\n",
      "9.15 [ 8.84 , 9.49 ]\n"
     ]
    }
   ],
   "source": [
    "\n",
    "gs_imp_firewood_ratio_analysis_project=df_gs_imp_firewood_VERs_agw_final/gs_imp_firewood_VERs\n",
    "gs_imp_firewood_ratio_analysis_project_mean=np.mean(df_gs_imp_firewood_VERs_agw_final, axis=1)/gs_imp_firewood_VERs\n",
    "gs_imp_firewood_ratio_project_analysis_mean=gs_imp_firewood_ratio_analysis_project_mean**-1\n",
    "\n",
    "gs_imp_firewood_sd_ratio_analysis_project=np.std(gs_imp_firewood_ratio_analysis_project, axis=1)\n",
    "\n",
    "\n",
    "\n",
    "gs_imp_firewood_ratio_analysis_project_lower=gs_imp_firewood_ratio_analysis_project_mean - z*gs_imp_firewood_sd_ratio_analysis_project/np.sqrt(gs_imp_firewood_N)\n",
    "gs_imp_firewood_ratio_analysis_project_upper=gs_imp_firewood_ratio_analysis_project_mean + z*gs_imp_firewood_sd_ratio_analysis_project/np.sqrt(gs_imp_firewood_N)\n",
    "\n",
    "gs_imp_firewood_ratio_project_analysis_lower=gs_imp_firewood_ratio_analysis_project_lower**-1\n",
    "gs_imp_firewood_ratio_project_analysis_upper=gs_imp_firewood_ratio_analysis_project_upper**-1\n",
    "\n",
    "print(gs_imp_firewood_ratio_analysis_project_lower)\n",
    "print(gs_imp_firewood_ratio_analysis_project_mean)\n",
    "print(gs_imp_firewood_ratio_analysis_project_upper)\n",
    "print(gs_imp_firewood_ratio_project_analysis_lower)\n",
    "print(gs_imp_firewood_ratio_project_analysis_mean)\n",
    "print(gs_imp_firewood_ratio_project_analysis_upper)\n",
    "\n",
    "print(round(gs_imp_firewood_ratio_analysis_project_mean[0],2),'[',round(gs_imp_firewood_ratio_analysis_project_lower[0],2),',',round(gs_imp_firewood_ratio_analysis_project_upper[0],2),']')\n",
    "print(round(gs_imp_firewood_ratio_project_analysis_mean[0],2),'[',round(gs_imp_firewood_ratio_project_analysis_upper[0],2),',',round(gs_imp_firewood_ratio_project_analysis_lower[0],2),']')\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 24,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "0    0.020122\n",
      "dtype: float64\n",
      "0    0.025558\n",
      "dtype: float64\n",
      "0    0.030993\n",
      "dtype: float64\n",
      "0    49.696692\n",
      "dtype: float64\n",
      "0    39.127089\n",
      "dtype: float64\n",
      "0    32.264916\n",
      "dtype: float64\n",
      "0.03 [ 0.02 , 0.03 ]\n",
      "39.13 [ 32.26 , 49.7 ]\n"
     ]
    }
   ],
   "source": [
    "gs_simple_firewood_ratio_analysis_project=df_gs_simple_firewood_VERs_agw_final/gs_simple_firewood_VERs\n",
    "gs_simple_firewood_ratio_analysis_project_mean=np.mean(df_gs_simple_firewood_VERs_agw_final, axis=1)/gs_simple_firewood_VERs\n",
    "gs_simple_firewood_ratio_project_analysis_mean=gs_simple_firewood_ratio_analysis_project_mean**-1\n",
    "\n",
    "gs_simple_firewood_sd_ratio_analysis_project=np.std(gs_simple_firewood_ratio_analysis_project, axis=1)\n",
    "\n",
    "\n",
    "\n",
    "gs_simple_firewood_ratio_analysis_project_lower=gs_simple_firewood_ratio_analysis_project_mean - z*gs_simple_firewood_sd_ratio_analysis_project/np.sqrt(gs_simple_firewood_N)\n",
    "gs_simple_firewood_ratio_analysis_project_upper=gs_simple_firewood_ratio_analysis_project_mean + z*gs_simple_firewood_sd_ratio_analysis_project/np.sqrt(gs_simple_firewood_N)\n",
    "\n",
    "gs_simple_firewood_ratio_project_analysis_lower=gs_simple_firewood_ratio_analysis_project_lower**-1\n",
    "gs_simple_firewood_ratio_project_analysis_upper=gs_simple_firewood_ratio_analysis_project_upper**-1\n",
    "\n",
    "print(gs_simple_firewood_ratio_analysis_project_lower)\n",
    "print(gs_simple_firewood_ratio_analysis_project_mean)\n",
    "print(gs_simple_firewood_ratio_analysis_project_upper)\n",
    "print(gs_simple_firewood_ratio_project_analysis_lower)\n",
    "print(gs_simple_firewood_ratio_project_analysis_mean)\n",
    "print(gs_simple_firewood_ratio_project_analysis_upper)\n",
    "\n",
    "print(round(gs_simple_firewood_ratio_analysis_project_mean[0],2),'[',round(gs_simple_firewood_ratio_analysis_project_lower[0],2),',',round(gs_simple_firewood_ratio_analysis_project_upper[0],2),']')\n",
    "print(round(gs_simple_firewood_ratio_project_analysis_mean[0],2),'[',round(gs_simple_firewood_ratio_project_analysis_upper[0],2),',',round(gs_simple_firewood_ratio_project_analysis_lower[0],2),']')\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 25,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "0    0.307373\n",
      "dtype: float64\n",
      "0    0.311342\n",
      "dtype: float64\n",
      "0    0.315311\n",
      "dtype: float64\n",
      "0    3.253375\n",
      "dtype: float64\n",
      "0    3.211901\n",
      "dtype: float64\n",
      "0    3.171472\n",
      "dtype: float64\n",
      "0.31 [ 0.31 , 0.32 ]\n",
      "3.21 [ 3.17 , 3.25 ]\n"
     ]
    }
   ],
   "source": [
    "\n",
    "gs_imp_charcoal_ratio_analysis_project=df_gs_imp_charcoal_VERs_agw_final/gs_imp_charcoal_VERs\n",
    "gs_imp_charcoal_ratio_analysis_project_mean=np.mean(df_gs_imp_charcoal_VERs_agw_final, axis=1)/gs_imp_charcoal_VERs\n",
    "gs_imp_charcoal_ratio_project_analysis_mean=gs_imp_charcoal_ratio_analysis_project_mean**-1\n",
    "\n",
    "gs_imp_charcoal_sd_ratio_analysis_project=np.std(gs_imp_charcoal_ratio_analysis_project, axis=1)\n",
    "\n",
    "\n",
    "\n",
    "gs_imp_charcoal_ratio_analysis_project_lower=gs_imp_charcoal_ratio_analysis_project_mean - z*gs_imp_charcoal_sd_ratio_analysis_project/np.sqrt(gs_imp_charcoal_N)\n",
    "gs_imp_charcoal_ratio_analysis_project_upper=gs_imp_charcoal_ratio_analysis_project_mean + z*gs_imp_charcoal_sd_ratio_analysis_project/np.sqrt(gs_imp_charcoal_N)\n",
    "\n",
    "gs_imp_charcoal_ratio_project_analysis_lower=gs_imp_charcoal_ratio_analysis_project_lower**-1\n",
    "gs_imp_charcoal_ratio_project_analysis_upper=gs_imp_charcoal_ratio_analysis_project_upper**-1\n",
    "\n",
    "print(gs_imp_charcoal_ratio_analysis_project_lower)\n",
    "print(gs_imp_charcoal_ratio_analysis_project_mean)\n",
    "print(gs_imp_charcoal_ratio_analysis_project_upper)\n",
    "print(gs_imp_charcoal_ratio_project_analysis_lower)\n",
    "print(gs_imp_charcoal_ratio_project_analysis_mean)\n",
    "print(gs_imp_charcoal_ratio_project_analysis_upper)\n",
    "\n",
    "print(round(gs_imp_charcoal_ratio_analysis_project_mean[0],2),'[',round(gs_imp_charcoal_ratio_analysis_project_lower[0],2),',',round(gs_imp_charcoal_ratio_analysis_project_upper[0],2),']')\n",
    "print(round(gs_imp_charcoal_ratio_project_analysis_mean[0],2),'[',round(gs_imp_charcoal_ratio_project_analysis_upper[0],2),',',round(gs_imp_charcoal_ratio_project_analysis_lower[0],2),']')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 26,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "0    0.036931\n",
      "dtype: float64\n",
      "0    0.048529\n",
      "dtype: float64\n",
      "0    0.060127\n",
      "dtype: float64\n",
      "0    27.077592\n",
      "dtype: float64\n",
      "0    20.606254\n",
      "dtype: float64\n",
      "0    16.631462\n",
      "dtype: float64\n",
      "0.05 [ 0.04 , 0.06 ]\n",
      "20.61 [ 16.63 , 27.08 ]\n"
     ]
    }
   ],
   "source": [
    "cdm_imp_charcoal_ratio_analysis_project=df_cdm_imp_charcoal_VERs_agw_final/cdm_imp_charcoal_VERs\n",
    "cdm_imp_charcoal_ratio_analysis_project_mean=np.mean(df_cdm_imp_charcoal_VERs_agw_final, axis=1)/cdm_imp_charcoal_VERs\n",
    "cdm_imp_charcoal_ratio_project_analysis_mean=cdm_imp_charcoal_ratio_analysis_project_mean**-1\n",
    "\n",
    "cdm_imp_charcoal_sd_ratio_analysis_project=np.std(cdm_imp_charcoal_ratio_analysis_project, axis=1)\n",
    "\n",
    "\n",
    "\n",
    "cdm_imp_charcoal_ratio_analysis_project_lower=cdm_imp_charcoal_ratio_analysis_project_mean - z*cdm_imp_charcoal_sd_ratio_analysis_project/np.sqrt(cdm_imp_charcoal_N)\n",
    "cdm_imp_charcoal_ratio_analysis_project_upper=cdm_imp_charcoal_ratio_analysis_project_mean + z*cdm_imp_charcoal_sd_ratio_analysis_project/np.sqrt(cdm_imp_charcoal_N)\n",
    "\n",
    "cdm_imp_charcoal_ratio_project_analysis_lower=cdm_imp_charcoal_ratio_analysis_project_lower**-1\n",
    "cdm_imp_charcoal_ratio_project_analysis_upper=cdm_imp_charcoal_ratio_analysis_project_upper**-1\n",
    "\n",
    "print(cdm_imp_charcoal_ratio_analysis_project_lower)\n",
    "print(cdm_imp_charcoal_ratio_analysis_project_mean)\n",
    "print(cdm_imp_charcoal_ratio_analysis_project_upper)\n",
    "print(cdm_imp_charcoal_ratio_project_analysis_lower)\n",
    "print(cdm_imp_charcoal_ratio_project_analysis_mean)\n",
    "print(cdm_imp_charcoal_ratio_project_analysis_upper)\n",
    "\n",
    "print(round(cdm_imp_charcoal_ratio_analysis_project_mean[0],2),'[',round(cdm_imp_charcoal_ratio_analysis_project_lower[0],2),',',round(cdm_imp_charcoal_ratio_analysis_project_upper[0],2),']')\n",
    "print(round(cdm_imp_charcoal_ratio_project_analysis_mean[0],2),'[',round(cdm_imp_charcoal_ratio_project_analysis_upper[0],2),',',round(cdm_imp_charcoal_ratio_project_analysis_lower[0],2),']')\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 28,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "0    0.131344\n",
      "dtype: float64\n",
      "0    0.133859\n",
      "dtype: float64\n",
      "0    0.136375\n",
      "dtype: float64\n",
      "0    7.61362\n",
      "dtype: float64\n",
      "0    7.470544\n",
      "dtype: float64\n",
      "0    7.332746\n",
      "dtype: float64\n",
      "0.13 [ 0.13 , 0.14 ]\n",
      "7.47 [ 7.33 , 7.61 ]\n"
     ]
    }
   ],
   "source": [
    "cdm_imp_firewood_ratio_analysis_project=df_cdm_imp_firewood_VERs_agw_final/cdm_imp_firewood_VERs\n",
    "cdm_imp_firewood_ratio_analysis_project_mean=np.mean(df_cdm_imp_firewood_VERs_agw_final, axis=1)/cdm_imp_firewood_VERs\n",
    "cdm_imp_firewood_ratio_project_analysis_mean=cdm_imp_firewood_ratio_analysis_project_mean**-1\n",
    "\n",
    "cdm_imp_firewood_sd_ratio_analysis_project=np.std(cdm_imp_firewood_ratio_analysis_project, axis=1)\n",
    "\n",
    "\n",
    "\n",
    "cdm_imp_firewood_ratio_analysis_project_lower=cdm_imp_firewood_ratio_analysis_project_mean - z*cdm_imp_firewood_sd_ratio_analysis_project/np.sqrt(cdm_imp_firewood_N)\n",
    "cdm_imp_firewood_ratio_analysis_project_upper=cdm_imp_firewood_ratio_analysis_project_mean + z*cdm_imp_firewood_sd_ratio_analysis_project/np.sqrt(cdm_imp_firewood_N)\n",
    "\n",
    "cdm_imp_firewood_ratio_project_analysis_lower=cdm_imp_firewood_ratio_analysis_project_lower**-1\n",
    "cdm_imp_firewood_ratio_project_analysis_upper=cdm_imp_firewood_ratio_analysis_project_upper**-1\n",
    "\n",
    "print(cdm_imp_firewood_ratio_analysis_project_lower)\n",
    "print(cdm_imp_firewood_ratio_analysis_project_mean)\n",
    "print(cdm_imp_firewood_ratio_analysis_project_upper)\n",
    "print(cdm_imp_firewood_ratio_project_analysis_lower)\n",
    "print(cdm_imp_firewood_ratio_project_analysis_mean)\n",
    "print(cdm_imp_firewood_ratio_project_analysis_upper)\n",
    "\n",
    "print(round(cdm_imp_firewood_ratio_analysis_project_mean[0],2),'[',round(cdm_imp_firewood_ratio_analysis_project_lower[0],2),',',round(cdm_imp_firewood_ratio_analysis_project_upper[0],2),']')\n",
    "print(round(cdm_imp_firewood_ratio_project_analysis_mean[0],2),'[',round(cdm_imp_firewood_ratio_project_analysis_upper[0],2),',',round(cdm_imp_firewood_ratio_project_analysis_lower[0],2),']')\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 30,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "0    0.638156\n",
      "dtype: float64\n",
      "0    0.648721\n",
      "dtype: float64\n",
      "0    0.659286\n",
      "dtype: float64\n",
      "0    1.567016\n",
      "dtype: float64\n",
      "0    1.541495\n",
      "dtype: float64\n",
      "0    1.516792\n",
      "dtype: float64\n",
      "0.65 [ 0.64 , 0.66 ]\n",
      "1.54 [ 1.52 , 1.57 ]\n"
     ]
    }
   ],
   "source": [
    "\n",
    "GS_Meter_ratio_analysis_project=df_GS_Meter_VERs_agw_final/sum(GS_Meter_VERs.iloc[:,1])\n",
    "GS_Meter_ratio_analysis_project_mean=np.mean(df_GS_Meter_VERs_agw_final, axis=1)/sum(GS_Meter_VERs.iloc[:,1])\n",
    "GS_Meter_ratio_project_analysis_mean=GS_Meter_ratio_analysis_project_mean**-1\n",
    "\n",
    "GS_Meter_sd_ratio_analysis_project=np.std(GS_Meter_ratio_analysis_project, axis=1)\n",
    "\n",
    "\n",
    "\n",
    "GS_Meter_ratio_analysis_project_lower=GS_Meter_ratio_analysis_project_mean - z*GS_Meter_sd_ratio_analysis_project/np.sqrt(GS_Meter_N)\n",
    "GS_Meter_ratio_analysis_project_upper=GS_Meter_ratio_analysis_project_mean + z*GS_Meter_sd_ratio_analysis_project/np.sqrt(GS_Meter_N)\n",
    "\n",
    "GS_Meter_ratio_project_analysis_lower=GS_Meter_ratio_analysis_project_lower**-1\n",
    "GS_Meter_ratio_project_analysis_upper=GS_Meter_ratio_analysis_project_upper**-1\n",
    "\n",
    "print(GS_Meter_ratio_analysis_project_lower)\n",
    "print(GS_Meter_ratio_analysis_project_mean)\n",
    "print(GS_Meter_ratio_analysis_project_upper)\n",
    "print(GS_Meter_ratio_project_analysis_lower)\n",
    "print(GS_Meter_ratio_project_analysis_mean)\n",
    "print(GS_Meter_ratio_project_analysis_upper)\n",
    "\n",
    "print(round(GS_Meter_ratio_analysis_project_mean[0],2),'[',round(GS_Meter_ratio_analysis_project_lower[0],2),',',round(GS_Meter_ratio_analysis_project_upper[0],2),']')\n",
    "print(round(GS_Meter_ratio_project_analysis_mean[0],2),'[',round(GS_Meter_ratio_project_analysis_upper[0],2),',',round(GS_Meter_ratio_project_analysis_lower[0],2),']')\n",
    "\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 31,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "0    0.38637\n",
      "dtype: float64\n",
      "0    0.420644\n",
      "dtype: float64\n",
      "0    0.454918\n",
      "dtype: float64\n",
      "0    2.588194\n",
      "dtype: float64\n",
      "0    2.377309\n",
      "dtype: float64\n",
      "0    2.1982\n",
      "dtype: float64\n",
      "0.42 [ 0.39 , 0.45 ]\n",
      "2.38 [ 2.2 , 2.59 ]\n"
     ]
    }
   ],
   "source": [
    "CDM_I_E_ratio_analysis_project=df_CDM_I_E_VERs_agw_final/sum(CDM_I_E_VERs.iloc[:,1])\n",
    "CDM_I_E_ratio_analysis_project_mean=np.mean(df_CDM_I_E_VERs_agw_final, axis=1)/sum(CDM_I_E_VERs.iloc[:,1])\n",
    "\n",
    "CDM_I_E_ratio_project_analysis_mean=CDM_I_E_ratio_analysis_project_mean**-1\n",
    "CDM_I_E_sd_ratio_analysis_project=np.std(CDM_I_E_ratio_analysis_project, axis=1)\n",
    "\n",
    "\n",
    "\n",
    "CDM_I_E_ratio_analysis_project_lower=CDM_I_E_ratio_analysis_project_mean - z*CDM_I_E_sd_ratio_analysis_project/np.sqrt(CDM_I_E_N)\n",
    "CDM_I_E_ratio_analysis_project_upper=CDM_I_E_ratio_analysis_project_mean + z*CDM_I_E_sd_ratio_analysis_project/np.sqrt(CDM_I_E_N)\n",
    "\n",
    "CDM_I_E_ratio_project_analysis_lower=CDM_I_E_ratio_analysis_project_lower**-1\n",
    "CDM_I_E_ratio_project_analysis_upper=CDM_I_E_ratio_analysis_project_upper**-1\n",
    "\n",
    "print(CDM_I_E_ratio_analysis_project_lower)\n",
    "print(CDM_I_E_ratio_analysis_project_mean)\n",
    "print(CDM_I_E_ratio_analysis_project_upper)\n",
    "print(CDM_I_E_ratio_project_analysis_lower)\n",
    "print(CDM_I_E_ratio_project_analysis_mean)\n",
    "print(CDM_I_E_ratio_project_analysis_upper)\n",
    "\n",
    "print(round(CDM_I_E_ratio_analysis_project_mean[0],2),'[',round(CDM_I_E_ratio_analysis_project_lower[0],2),',',round(CDM_I_E_ratio_analysis_project_upper[0],2),']')\n",
    "print(round(CDM_I_E_ratio_project_analysis_mean[0],2),'[',round(CDM_I_E_ratio_project_analysis_upper[0],2),',',round(CDM_I_E_ratio_project_analysis_lower[0],2),']')\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 33,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "8\n",
      "0    0.57986\n",
      "dtype: float64\n",
      "0    0.589386\n",
      "dtype: float64\n",
      "0    0.598912\n",
      "dtype: float64\n",
      "0    1.724555\n",
      "dtype: float64\n",
      "0    1.696682\n",
      "dtype: float64\n",
      "0    1.669695\n",
      "dtype: float64\n",
      "0.59 [ 0.58 , 0.6 ]\n",
      "1.7 [ 1.67 , 1.72 ]\n"
     ]
    }
   ],
   "source": [
    "df_bangladesh_VERs_agw=pd.DataFrame(bangladesh_VERs_agw).transpose()\n",
    "print(df_bangladesh_VERs_agw.shape[1])\n",
    "bangladesh_i=0\n",
    "df_bangladesh_VERs_agw_final=pd.DataFrame(np.zeros((1, 10000)))\n",
    "while bangladesh_i<df_bangladesh_VERs_agw.shape[1]:\n",
    "    df_bangladesh_VERs_agw_final=df_bangladesh_VERs_agw_final+df_bangladesh_VERs_agw[bangladesh_i]\n",
    "    bangladesh_i=bangladesh_i+1\n",
    "\n",
    "bangladesh_ratio_analysis_project=df_bangladesh_VERs_agw_final/bangladesh_VERs\n",
    "bangladesh_ratio_analysis_project_mean=np.mean(df_bangladesh_VERs_agw_final, axis=1)/bangladesh_VERs\n",
    "bangladesh_ratio_project_analysis_mean=bangladesh_ratio_analysis_project_mean**-1\n",
    "\n",
    "bangladesh_sd_ratio_analysis_project=np.std(bangladesh_ratio_analysis_project, axis=1)\n",
    "\n",
    "\n",
    "\n",
    "bangladesh_ratio_analysis_project_lower=bangladesh_ratio_analysis_project_mean - z*bangladesh_sd_ratio_analysis_project/np.sqrt(bangladesh_N)\n",
    "bangladesh_ratio_analysis_project_upper=bangladesh_ratio_analysis_project_mean + z*bangladesh_sd_ratio_analysis_project/np.sqrt(bangladesh_N)\n",
    "\n",
    "bangladesh_ratio_project_analysis_lower=bangladesh_ratio_analysis_project_lower**-1\n",
    "bangladesh_ratio_project_analysis_upper=bangladesh_ratio_analysis_project_upper**-1\n",
    "\n",
    "print(bangladesh_ratio_analysis_project_lower)\n",
    "print(bangladesh_ratio_analysis_project_mean)\n",
    "print(bangladesh_ratio_analysis_project_upper)\n",
    "print(bangladesh_ratio_project_analysis_lower)\n",
    "print(bangladesh_ratio_project_analysis_mean)\n",
    "print(bangladesh_ratio_project_analysis_upper)\n",
    "\n",
    "print(round(bangladesh_ratio_analysis_project_mean[0],2),'[',round(bangladesh_ratio_analysis_project_lower[0],2),',',round(bangladesh_ratio_analysis_project_upper[0],2),']')\n",
    "print(round(bangladesh_ratio_project_analysis_mean[0],2),'[',round(bangladesh_ratio_project_analysis_upper[0],2),',',round(bangladesh_ratio_project_analysis_lower[0],2),']')\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 34,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "0\n",
      "0   NaN\n",
      "dtype: float64\n",
      "0   NaN\n",
      "dtype: float64\n",
      "0   NaN\n",
      "dtype: float64\n",
      "0   NaN\n",
      "dtype: float64\n",
      "0   NaN\n",
      "dtype: float64\n",
      "0   NaN\n",
      "dtype: float64\n",
      "nan [ nan , nan ]\n",
      "nan [ nan , nan ]\n"
     ]
    }
   ],
   "source": [
    "df_burkinafaso_VERs_agw=pd.DataFrame(burkinafaso_VERs_agw).transpose()\n",
    "print(df_burkinafaso_VERs_agw.shape[1])\n",
    "burkinafaso_i=0\n",
    "df_burkinafaso_VERs_agw_final=pd.DataFrame(np.zeros((1, 10000)))\n",
    "while burkinafaso_i<df_burkinafaso_VERs_agw.shape[1]:\n",
    "    df_burkinafaso_VERs_agw_final=df_burkinafaso_VERs_agw_final+df_burkinafaso_VERs_agw[burkinafaso_i]\n",
    "    burkinafaso_i=burkinafaso_i+1\n",
    " \n",
    "burkinafaso_ratio_analysis_project=df_burkinafaso_VERs_agw_final/burkinafaso_VERs\n",
    "burkinafaso_ratio_analysis_project_mean=np.mean(df_burkinafaso_VERs_agw_final, axis=1)/burkinafaso_VERs\n",
    "burkinafaso_ratio_project_analysis_mean=burkinafaso_ratio_analysis_project_mean**-1\n",
    "\n",
    "burkinafaso_sd_ratio_analysis_project=np.std(burkinafaso_ratio_analysis_project, axis=1)\n",
    "\n",
    "\n",
    "\n",
    "burkinafaso_ratio_analysis_project_lower=burkinafaso_ratio_analysis_project_mean - z*burkinafaso_sd_ratio_analysis_project/np.sqrt(burkinafaso_N)\n",
    "burkinafaso_ratio_analysis_project_upper=burkinafaso_ratio_analysis_project_mean + z*burkinafaso_sd_ratio_analysis_project/np.sqrt(burkinafaso_N)\n",
    "\n",
    "burkinafaso_ratio_project_analysis_lower=burkinafaso_ratio_analysis_project_lower**-1\n",
    "burkinafaso_ratio_project_analysis_upper=burkinafaso_ratio_analysis_project_upper**-1\n",
    "\n",
    "print(burkinafaso_ratio_analysis_project_lower)\n",
    "print(burkinafaso_ratio_analysis_project_mean)\n",
    "print(burkinafaso_ratio_analysis_project_upper)\n",
    "print(burkinafaso_ratio_project_analysis_lower)\n",
    "print(burkinafaso_ratio_project_analysis_mean)\n",
    "print(burkinafaso_ratio_project_analysis_upper)\n",
    "\n",
    "print(round(burkinafaso_ratio_analysis_project_mean[0],2),'[',round(burkinafaso_ratio_analysis_project_lower[0],2),',',round(burkinafaso_ratio_analysis_project_upper[0],2),']')\n",
    "print(round(burkinafaso_ratio_project_analysis_mean[0],2),'[',round(burkinafaso_ratio_project_analysis_upper[0],2),',',round(burkinafaso_ratio_project_analysis_lower[0],2),']')\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 35,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "6\n",
      "0    0.002794\n",
      "dtype: float64\n",
      "0    0.004867\n",
      "dtype: float64\n",
      "0    0.006941\n",
      "dtype: float64\n",
      "0    357.957757\n",
      "dtype: float64\n",
      "0    205.445604\n",
      "dtype: float64\n",
      "0    144.065\n",
      "dtype: float64\n",
      "0.0 [ 0.0 , 0.01 ]\n",
      "205.45 [ 144.07 , 357.96 ]\n"
     ]
    }
   ],
   "source": [
    "df_china_VERs_agw=pd.DataFrame(china_VERs_agw).transpose()\n",
    "print(df_china_VERs_agw.shape[1])\n",
    "china_i=0\n",
    "df_china_VERs_agw_final=pd.DataFrame(np.zeros((1, 10000)))\n",
    "while china_i<df_china_VERs_agw.shape[1]:\n",
    "    df_china_VERs_agw_final=df_china_VERs_agw_final+df_china_VERs_agw[china_i]\n",
    "    china_i=china_i+1\n",
    "\n",
    "china_ratio_analysis_project=df_china_VERs_agw_final/china_VERs\n",
    "china_ratio_analysis_project_mean=np.mean(df_china_VERs_agw_final, axis=1)/china_VERs\n",
    "china_ratio_project_analysis_mean=china_ratio_analysis_project_mean**-1\n",
    "\n",
    "china_sd_ratio_analysis_project=np.std(china_ratio_analysis_project, axis=1)\n",
    "\n",
    "\n",
    "\n",
    "china_ratio_analysis_project_lower=china_ratio_analysis_project_mean - z*china_sd_ratio_analysis_project/np.sqrt(china_N)\n",
    "china_ratio_analysis_project_upper=china_ratio_analysis_project_mean + z*china_sd_ratio_analysis_project/np.sqrt(china_N)\n",
    "\n",
    "china_ratio_project_analysis_lower=china_ratio_analysis_project_lower**-1\n",
    "china_ratio_project_analysis_upper=china_ratio_analysis_project_upper**-1\n",
    "\n",
    "print(china_ratio_analysis_project_lower)\n",
    "print(china_ratio_analysis_project_mean)\n",
    "print(china_ratio_analysis_project_upper)\n",
    "print(china_ratio_project_analysis_lower)\n",
    "print(china_ratio_project_analysis_mean)\n",
    "print(china_ratio_project_analysis_upper)\n",
    "\n",
    "print(round(china_ratio_analysis_project_mean[0],2),'[',round(china_ratio_analysis_project_lower[0],2),',',round(china_ratio_analysis_project_upper[0],2),']')\n",
    "print(round(china_ratio_project_analysis_mean[0],2),'[',round(china_ratio_project_analysis_upper[0],2),',',round(china_ratio_project_analysis_lower[0],2),']')\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 36,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "1\n",
      "0    0.168755\n",
      "dtype: float64\n",
      "0    0.290895\n",
      "dtype: float64\n",
      "0    0.413035\n",
      "dtype: float64\n",
      "0    5.925766\n",
      "dtype: float64\n",
      "0    3.437671\n",
      "dtype: float64\n",
      "0    2.421104\n",
      "dtype: float64\n",
      "0.29 [ 0.17 , 0.41 ]\n",
      "3.44 [ 2.42 , 5.93 ]\n"
     ]
    }
   ],
   "source": [
    "df_cameroon_VERs_agw=pd.DataFrame(cameroon_VERs_agw).transpose()\n",
    "print(df_cameroon_VERs_agw.shape[1])\n",
    "cameroon_i=0\n",
    "df_cameroon_VERs_agw_final=pd.DataFrame(np.zeros((1, 10000)))\n",
    "while cameroon_i<df_cameroon_VERs_agw.shape[1]:\n",
    "    df_cameroon_VERs_agw_final=df_cameroon_VERs_agw_final+df_cameroon_VERs_agw[cameroon_i]\n",
    "    cameroon_i=cameroon_i+1\n",
    "\n",
    "cameroon_ratio_analysis_project=df_cameroon_VERs_agw_final/cameroon_VERs\n",
    "cameroon_ratio_analysis_project_mean=np.mean(df_cameroon_VERs_agw_final, axis=1)/cameroon_VERs\n",
    "cameroon_ratio_project_analysis_mean=cameroon_ratio_analysis_project_mean**-1\n",
    "\n",
    "cameroon_sd_ratio_analysis_project=np.std(cameroon_ratio_analysis_project, axis=1)\n",
    "\n",
    "\n",
    "\n",
    "cameroon_ratio_analysis_project_lower=cameroon_ratio_analysis_project_mean - z*cameroon_sd_ratio_analysis_project/np.sqrt(cameroon_N)\n",
    "cameroon_ratio_analysis_project_upper=cameroon_ratio_analysis_project_mean + z*cameroon_sd_ratio_analysis_project/np.sqrt(cameroon_N)\n",
    "\n",
    "cameroon_ratio_project_analysis_lower=cameroon_ratio_analysis_project_lower**-1\n",
    "cameroon_ratio_project_analysis_upper=cameroon_ratio_analysis_project_upper**-1\n",
    "\n",
    "print(cameroon_ratio_analysis_project_lower)\n",
    "print(cameroon_ratio_analysis_project_mean)\n",
    "print(cameroon_ratio_analysis_project_upper)\n",
    "print(cameroon_ratio_project_analysis_lower)\n",
    "print(cameroon_ratio_project_analysis_mean)\n",
    "print(cameroon_ratio_project_analysis_upper)\n",
    "\n",
    "print(round(cameroon_ratio_analysis_project_mean[0],2),'[',round(cameroon_ratio_analysis_project_lower[0],2),',',round(cameroon_ratio_analysis_project_upper[0],2),']')\n",
    "print(round(cameroon_ratio_project_analysis_mean[0],2),'[',round(cameroon_ratio_project_analysis_upper[0],2),',',round(cameroon_ratio_project_analysis_lower[0],2),']')\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 37,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "2\n",
      "0    0.11893\n",
      "dtype: float64\n",
      "0    0.198727\n",
      "dtype: float64\n",
      "0    0.278523\n",
      "dtype: float64\n",
      "0    8.408285\n",
      "dtype: float64\n",
      "0    5.032033\n",
      "dtype: float64\n",
      "0    3.590363\n",
      "dtype: float64\n",
      "0.2 [ 0.12 , 0.28 ]\n",
      "5.03 [ 3.59 , 8.41 ]\n"
     ]
    }
   ],
   "source": [
    "df_ethiopia_VERs_agw=pd.DataFrame(ethiopia_VERs_agw).transpose()\n",
    "print(df_ethiopia_VERs_agw.shape[1])\n",
    "ethiopia_i=0\n",
    "df_ethiopia_VERs_agw_final=pd.DataFrame(np.zeros((1, 10000)))\n",
    "while ethiopia_i<df_ethiopia_VERs_agw.shape[1]:\n",
    "    df_ethiopia_VERs_agw_final=df_ethiopia_VERs_agw_final+df_ethiopia_VERs_agw[ethiopia_i]\n",
    "    ethiopia_i=ethiopia_i+1\n",
    "\n",
    "ethiopia_ratio_analysis_project=df_ethiopia_VERs_agw_final/ethiopia_VERs\n",
    "ethiopia_ratio_analysis_project_mean=np.mean(df_ethiopia_VERs_agw_final, axis=1)/ethiopia_VERs\n",
    "ethiopia_ratio_project_analysis_mean=ethiopia_ratio_analysis_project_mean**-1\n",
    "\n",
    "ethiopia_sd_ratio_analysis_project=np.std(ethiopia_ratio_analysis_project, axis=1)\n",
    "\n",
    "\n",
    "\n",
    "ethiopia_ratio_analysis_project_lower=ethiopia_ratio_analysis_project_mean - z*ethiopia_sd_ratio_analysis_project/np.sqrt(ethiopia_N)\n",
    "ethiopia_ratio_analysis_project_upper=ethiopia_ratio_analysis_project_mean + z*ethiopia_sd_ratio_analysis_project/np.sqrt(ethiopia_N)\n",
    "\n",
    "ethiopia_ratio_project_analysis_lower=ethiopia_ratio_analysis_project_lower**-1\n",
    "ethiopia_ratio_project_analysis_upper=ethiopia_ratio_analysis_project_upper**-1\n",
    "\n",
    "print(ethiopia_ratio_analysis_project_lower)\n",
    "print(ethiopia_ratio_analysis_project_mean)\n",
    "print(ethiopia_ratio_analysis_project_upper)\n",
    "print(ethiopia_ratio_project_analysis_lower)\n",
    "print(ethiopia_ratio_project_analysis_mean)\n",
    "print(ethiopia_ratio_project_analysis_upper)\n",
    "\n",
    "print(round(ethiopia_ratio_analysis_project_mean[0],2),'[',round(ethiopia_ratio_analysis_project_lower[0],2),',',round(ethiopia_ratio_analysis_project_upper[0],2),']')\n",
    "print(round(ethiopia_ratio_project_analysis_mean[0],2),'[',round(ethiopia_ratio_project_analysis_upper[0],2),',',round(ethiopia_ratio_project_analysis_lower[0],2),']')\n",
    "\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 38,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "16\n",
      "0    0.414624\n",
      "dtype: float64\n",
      "0    0.423033\n",
      "dtype: float64\n",
      "0    0.431443\n",
      "dtype: float64\n",
      "0    2.411823\n",
      "dtype: float64\n",
      "0    2.363879\n",
      "dtype: float64\n",
      "0    2.317804\n",
      "dtype: float64\n",
      "0.42 [ 0.41 , 0.43 ]\n",
      "2.36 [ 2.32 , 2.41 ]\n"
     ]
    }
   ],
   "source": [
    "df_ghana_VERs_agw=pd.DataFrame(ghana_VERs_agw).transpose()\n",
    "print(df_ghana_VERs_agw.shape[1])\n",
    "ghana_i=0\n",
    "df_ghana_VERs_agw_final=pd.DataFrame(np.zeros((1, 10000)))\n",
    "while ghana_i<df_ghana_VERs_agw.shape[1]:\n",
    "    df_ghana_VERs_agw_final=df_ghana_VERs_agw_final+df_ghana_VERs_agw[ghana_i]\n",
    "    ghana_i=ghana_i+1\n",
    "\n",
    "ghana_ratio_analysis_project=df_ghana_VERs_agw_final/ghana_VERs\n",
    "ghana_ratio_analysis_project_mean=np.mean(df_ghana_VERs_agw_final, axis=1)/ghana_VERs\n",
    "ghana_ratio_project_analysis_mean=ghana_ratio_analysis_project_mean**-1\n",
    "\n",
    "ghana_sd_ratio_analysis_project=np.std(ghana_ratio_analysis_project, axis=1)\n",
    "\n",
    "\n",
    "\n",
    "ghana_ratio_analysis_project_lower=ghana_ratio_analysis_project_mean - z*ghana_sd_ratio_analysis_project/np.sqrt(ghana_N)\n",
    "ghana_ratio_analysis_project_upper=ghana_ratio_analysis_project_mean + z*ghana_sd_ratio_analysis_project/np.sqrt(ghana_N)\n",
    "\n",
    "ghana_ratio_project_analysis_lower=ghana_ratio_analysis_project_lower**-1\n",
    "ghana_ratio_project_analysis_upper=ghana_ratio_analysis_project_upper**-1\n",
    "\n",
    "print(ghana_ratio_analysis_project_lower)\n",
    "print(ghana_ratio_analysis_project_mean)\n",
    "print(ghana_ratio_analysis_project_upper)\n",
    "print(ghana_ratio_project_analysis_lower)\n",
    "print(ghana_ratio_project_analysis_mean)\n",
    "print(ghana_ratio_project_analysis_upper)\n",
    "\n",
    "print(round(ghana_ratio_analysis_project_mean[0],2),'[',round(ghana_ratio_analysis_project_lower[0],2),',',round(ghana_ratio_analysis_project_upper[0],2),']')\n",
    "print(round(ghana_ratio_project_analysis_mean[0],2),'[',round(ghana_ratio_project_analysis_upper[0],2),',',round(ghana_ratio_project_analysis_lower[0],2),']')\n",
    "\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 39,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "8\n",
      "0    0.084919\n",
      "dtype: float64\n",
      "0    0.091173\n",
      "dtype: float64\n",
      "0    0.097427\n",
      "dtype: float64\n",
      "0    11.775865\n",
      "dtype: float64\n",
      "0    10.968109\n",
      "dtype: float64\n",
      "0    10.264054\n",
      "dtype: float64\n",
      "0.09 [ 0.08 , 0.1 ]\n",
      "10.97 [ 10.26 , 11.78 ]\n"
     ]
    }
   ],
   "source": [
    "df_guatemala_VERs_agw=pd.DataFrame(guatemala_VERs_agw).transpose()\n",
    "print(df_guatemala_VERs_agw.shape[1])\n",
    "guatemala_i=0\n",
    "df_guatemala_VERs_agw_final=pd.DataFrame(np.zeros((1, 10000)))\n",
    "while guatemala_i<df_guatemala_VERs_agw.shape[1]:\n",
    "    df_guatemala_VERs_agw_final=df_guatemala_VERs_agw_final+df_guatemala_VERs_agw[guatemala_i]\n",
    "    guatemala_i=guatemala_i+1\n",
    "\n",
    "guatemala_ratio_analysis_project=df_guatemala_VERs_agw_final/guatemala_VERs\n",
    "guatemala_ratio_analysis_project_mean=np.mean(df_guatemala_VERs_agw_final, axis=1)/guatemala_VERs\n",
    "guatemala_ratio_project_analysis_mean=guatemala_ratio_analysis_project_mean**-1\n",
    "\n",
    "guatemala_sd_ratio_analysis_project=np.std(guatemala_ratio_analysis_project, axis=1)\n",
    "\n",
    "\n",
    "\n",
    "guatemala_ratio_analysis_project_lower=guatemala_ratio_analysis_project_mean - z*guatemala_sd_ratio_analysis_project/np.sqrt(guatemala_N)\n",
    "guatemala_ratio_analysis_project_upper=guatemala_ratio_analysis_project_mean + z*guatemala_sd_ratio_analysis_project/np.sqrt(guatemala_N)\n",
    "\n",
    "guatemala_ratio_project_analysis_lower=guatemala_ratio_analysis_project_lower**-1\n",
    "guatemala_ratio_project_analysis_upper=guatemala_ratio_analysis_project_upper**-1\n",
    "\n",
    "print(guatemala_ratio_analysis_project_lower)\n",
    "print(guatemala_ratio_analysis_project_mean)\n",
    "print(guatemala_ratio_analysis_project_upper)\n",
    "print(guatemala_ratio_project_analysis_lower)\n",
    "print(guatemala_ratio_project_analysis_mean)\n",
    "print(guatemala_ratio_project_analysis_upper)\n",
    "\n",
    "print(round(guatemala_ratio_analysis_project_mean[0],2),'[',round(guatemala_ratio_analysis_project_lower[0],2),',',round(guatemala_ratio_analysis_project_upper[0],2),']')\n",
    "print(round(guatemala_ratio_project_analysis_mean[0],2),'[',round(guatemala_ratio_project_analysis_upper[0],2),',',round(guatemala_ratio_project_analysis_lower[0],2),']')\n",
    "\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 40,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "1\n",
      "0    0.047825\n",
      "dtype: float64\n",
      "0    0.08283\n",
      "dtype: float64\n",
      "0    0.117835\n",
      "dtype: float64\n",
      "0    20.909575\n",
      "dtype: float64\n",
      "0    12.07289\n",
      "dtype: float64\n",
      "0    8.486412\n",
      "dtype: float64\n",
      "0.08 [ 0.05 , 0.12 ]\n",
      "12.07 [ 8.49 , 20.91 ]\n"
     ]
    }
   ],
   "source": [
    "df_haiti_VERs_agw=pd.DataFrame(haiti_VERs_agw).transpose()\n",
    "print(df_haiti_VERs_agw.shape[1])\n",
    "haiti_i=0\n",
    "df_haiti_VERs_agw_final=pd.DataFrame(np.zeros((1, 10000)))\n",
    "while haiti_i<df_haiti_VERs_agw.shape[1]:\n",
    "    df_haiti_VERs_agw_final=df_haiti_VERs_agw_final+df_haiti_VERs_agw[haiti_i]\n",
    "    haiti_i=haiti_i+1\n",
    "\n",
    "haiti_ratio_analysis_project=df_haiti_VERs_agw_final/haiti_VERs\n",
    "haiti_ratio_analysis_project_mean=np.mean(df_haiti_VERs_agw_final, axis=1)/haiti_VERs\n",
    "haiti_ratio_project_analysis_mean=haiti_ratio_analysis_project_mean**-1\n",
    "\n",
    "haiti_sd_ratio_analysis_project=np.std(haiti_ratio_analysis_project, axis=1)\n",
    "\n",
    "\n",
    "\n",
    "haiti_ratio_analysis_project_lower=haiti_ratio_analysis_project_mean - z*haiti_sd_ratio_analysis_project/np.sqrt(haiti_N)\n",
    "haiti_ratio_analysis_project_upper=haiti_ratio_analysis_project_mean + z*haiti_sd_ratio_analysis_project/np.sqrt(haiti_N)\n",
    "\n",
    "haiti_ratio_project_analysis_lower=haiti_ratio_analysis_project_lower**-1\n",
    "haiti_ratio_project_analysis_upper=haiti_ratio_analysis_project_upper**-1\n",
    "\n",
    "print(haiti_ratio_analysis_project_lower)\n",
    "print(haiti_ratio_analysis_project_mean)\n",
    "print(haiti_ratio_analysis_project_upper)\n",
    "print(haiti_ratio_project_analysis_lower)\n",
    "print(haiti_ratio_project_analysis_mean)\n",
    "print(haiti_ratio_project_analysis_upper)\n",
    "\n",
    "print(round(haiti_ratio_analysis_project_mean[0],2),'[',round(haiti_ratio_analysis_project_lower[0],2),',',round(haiti_ratio_analysis_project_upper[0],2),']')\n",
    "print(round(haiti_ratio_project_analysis_mean[0],2),'[',round(haiti_ratio_project_analysis_upper[0],2),',',round(haiti_ratio_project_analysis_lower[0],2),']')\n",
    "\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 41,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "6\n",
      "0    0.067106\n",
      "dtype: float64\n",
      "0    0.126718\n",
      "dtype: float64\n",
      "0    0.186329\n",
      "dtype: float64\n",
      "0    14.901738\n",
      "dtype: float64\n",
      "0    7.891559\n",
      "dtype: float64\n",
      "0    5.366849\n",
      "dtype: float64\n",
      "0.13 [ 0.07 , 0.19 ]\n",
      "7.89 [ 5.37 , 14.9 ]\n"
     ]
    }
   ],
   "source": [
    "df_honduras_VERs_agw=pd.DataFrame(honduras_VERs_agw).transpose()\n",
    "print(df_honduras_VERs_agw.shape[1])\n",
    "honduras_i=0\n",
    "df_honduras_VERs_agw_final=pd.DataFrame(np.zeros((1, 10000)))\n",
    "while honduras_i<df_honduras_VERs_agw.shape[1]:\n",
    "    df_honduras_VERs_agw_final=df_honduras_VERs_agw_final+df_honduras_VERs_agw[honduras_i]\n",
    "    honduras_i=honduras_i+1\n",
    "\n",
    "honduras_ratio_analysis_project=df_honduras_VERs_agw_final/honduras_VERs\n",
    "honduras_ratio_analysis_project_mean=np.mean(df_honduras_VERs_agw_final, axis=1)/honduras_VERs\n",
    "honduras_ratio_project_analysis_mean=honduras_ratio_analysis_project_mean**-1\n",
    "\n",
    "honduras_sd_ratio_analysis_project=np.std(honduras_ratio_analysis_project, axis=1)\n",
    "\n",
    "\n",
    "\n",
    "honduras_ratio_analysis_project_lower=honduras_ratio_analysis_project_mean - z*honduras_sd_ratio_analysis_project/np.sqrt(honduras_N)\n",
    "honduras_ratio_analysis_project_upper=honduras_ratio_analysis_project_mean + z*honduras_sd_ratio_analysis_project/np.sqrt(honduras_N)\n",
    "\n",
    "honduras_ratio_project_analysis_lower=honduras_ratio_analysis_project_lower**-1\n",
    "honduras_ratio_project_analysis_upper=honduras_ratio_analysis_project_upper**-1\n",
    "\n",
    "print(honduras_ratio_analysis_project_lower)\n",
    "print(honduras_ratio_analysis_project_mean)\n",
    "print(honduras_ratio_analysis_project_upper)\n",
    "print(honduras_ratio_project_analysis_lower)\n",
    "print(honduras_ratio_project_analysis_mean)\n",
    "print(honduras_ratio_project_analysis_upper)\n",
    "\n",
    "print(round(honduras_ratio_analysis_project_mean[0],2),'[',round(honduras_ratio_analysis_project_lower[0],2),',',round(honduras_ratio_analysis_project_upper[0],2),']')\n",
    "print(round(honduras_ratio_project_analysis_mean[0],2),'[',round(honduras_ratio_project_analysis_upper[0],2),',',round(honduras_ratio_project_analysis_lower[0],2),']')\n",
    "\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 42,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "2\n",
      "0    0.010439\n",
      "dtype: float64\n",
      "0    0.040831\n",
      "dtype: float64\n",
      "0    0.071223\n",
      "dtype: float64\n",
      "0    95.796856\n",
      "dtype: float64\n",
      "0    24.491309\n",
      "dtype: float64\n",
      "0    14.040435\n",
      "dtype: float64\n",
      "0.04 [ 0.01 , 0.07 ]\n",
      "24.49 [ 14.04 , 95.8 ]\n"
     ]
    }
   ],
   "source": [
    "df_india_VERs_agw=pd.DataFrame(india_VERs_agw).transpose()\n",
    "print(df_india_VERs_agw.shape[1])\n",
    "india_i=0\n",
    "df_india_VERs_agw_final=pd.DataFrame(np.zeros((1, 10000)))\n",
    "while india_i<df_india_VERs_agw.shape[1]:\n",
    "    df_india_VERs_agw_final=df_india_VERs_agw_final+df_india_VERs_agw[india_i]\n",
    "    india_i=india_i+1\n",
    "\n",
    "india_ratio_analysis_project=df_india_VERs_agw_final/india_VERs\n",
    "india_ratio_analysis_project_mean=np.mean(df_india_VERs_agw_final, axis=1)/india_VERs\n",
    "india_ratio_project_analysis_mean=india_ratio_analysis_project_mean**-1\n",
    "\n",
    "india_sd_ratio_analysis_project=np.std(india_ratio_analysis_project, axis=1)\n",
    "\n",
    "\n",
    "\n",
    "india_ratio_analysis_project_lower=india_ratio_analysis_project_mean - z*india_sd_ratio_analysis_project/np.sqrt(india_N)\n",
    "india_ratio_analysis_project_upper=india_ratio_analysis_project_mean + z*india_sd_ratio_analysis_project/np.sqrt(india_N)\n",
    "\n",
    "india_ratio_project_analysis_lower=india_ratio_analysis_project_lower**-1\n",
    "india_ratio_project_analysis_upper=india_ratio_analysis_project_upper**-1\n",
    "\n",
    "print(india_ratio_analysis_project_lower)\n",
    "print(india_ratio_analysis_project_mean)\n",
    "print(india_ratio_analysis_project_upper)\n",
    "print(india_ratio_project_analysis_lower)\n",
    "print(india_ratio_project_analysis_mean)\n",
    "print(india_ratio_project_analysis_upper)\n",
    "\n",
    "print(round(india_ratio_analysis_project_mean[0],2),'[',round(india_ratio_analysis_project_lower[0],2),',',round(india_ratio_analysis_project_upper[0],2),']')\n",
    "print(round(india_ratio_project_analysis_mean[0],2),'[',round(india_ratio_project_analysis_upper[0],2),',',round(india_ratio_project_analysis_lower[0],2),']')\n",
    "\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 43,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "7\n",
      "0    0.46996\n",
      "dtype: float64\n",
      "0    0.485762\n",
      "dtype: float64\n",
      "0    0.501564\n",
      "dtype: float64\n",
      "0    2.127841\n",
      "dtype: float64\n",
      "0    2.058621\n",
      "dtype: float64\n",
      "0    1.993763\n",
      "dtype: float64\n",
      "0.49 [ 0.47 , 0.5 ]\n",
      "2.06 [ 1.99 , 2.13 ]\n"
     ]
    }
   ],
   "source": [
    "df_kenya_VERs_agw=pd.DataFrame(kenya_VERs_agw).transpose()\n",
    "print(df_kenya_VERs_agw.shape[1])\n",
    "kenya_i=0\n",
    "df_kenya_VERs_agw_final=pd.DataFrame(np.zeros((1, 10000)))\n",
    "while kenya_i<df_kenya_VERs_agw.shape[1]:\n",
    "    df_kenya_VERs_agw_final=df_kenya_VERs_agw_final+df_kenya_VERs_agw[kenya_i]\n",
    "    kenya_i=kenya_i+1\n",
    "kenya_ratio_analysis_project=df_kenya_VERs_agw_final/kenya_VERs\n",
    "kenya_ratio_analysis_project_mean=np.mean(df_kenya_VERs_agw_final, axis=1)/kenya_VERs\n",
    "kenya_ratio_project_analysis_mean=kenya_ratio_analysis_project_mean**-1\n",
    "\n",
    "kenya_sd_ratio_analysis_project=np.std(kenya_ratio_analysis_project, axis=1)\n",
    "\n",
    "\n",
    "\n",
    "kenya_ratio_analysis_project_lower=kenya_ratio_analysis_project_mean - z*kenya_sd_ratio_analysis_project/np.sqrt(kenya_N)\n",
    "kenya_ratio_analysis_project_upper=kenya_ratio_analysis_project_mean + z*kenya_sd_ratio_analysis_project/np.sqrt(kenya_N)\n",
    "\n",
    "kenya_ratio_project_analysis_lower=kenya_ratio_analysis_project_lower**-1\n",
    "kenya_ratio_project_analysis_upper=kenya_ratio_analysis_project_upper**-1\n",
    "\n",
    "print(kenya_ratio_analysis_project_lower)\n",
    "print(kenya_ratio_analysis_project_mean)\n",
    "print(kenya_ratio_analysis_project_upper)\n",
    "print(kenya_ratio_project_analysis_lower)\n",
    "print(kenya_ratio_project_analysis_mean)\n",
    "print(kenya_ratio_project_analysis_upper)\n",
    "\n",
    "print(round(kenya_ratio_analysis_project_mean[0],2),'[',round(kenya_ratio_analysis_project_lower[0],2),',',round(kenya_ratio_analysis_project_upper[0],2),']')\n",
    "print(round(kenya_ratio_project_analysis_mean[0],2),'[',round(kenya_ratio_project_analysis_upper[0],2),',',round(kenya_ratio_project_analysis_lower[0],2),']')\n",
    "\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 44,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "1\n",
      "0   -0.002349\n",
      "dtype: float64\n",
      "0    0.047214\n",
      "dtype: float64\n",
      "0    0.096777\n",
      "dtype: float64\n",
      "0   -425.748178\n",
      "dtype: float64\n",
      "0    21.180103\n",
      "dtype: float64\n",
      "0    10.333028\n",
      "dtype: float64\n",
      "0.05 [ -0.0 , 0.1 ]\n",
      "21.18 [ 10.33 , -425.75 ]\n"
     ]
    }
   ],
   "source": [
    "df_lesotho_VERs_agw=pd.DataFrame(lesotho_VERs_agw).transpose()\n",
    "print(df_lesotho_VERs_agw.shape[1])\n",
    "lesotho_i=0\n",
    "df_lesotho_VERs_agw_final=pd.DataFrame(np.zeros((1, 10000)))\n",
    "while lesotho_i<df_lesotho_VERs_agw.shape[1]:\n",
    "    df_lesotho_VERs_agw_final=df_lesotho_VERs_agw_final+df_lesotho_VERs_agw[lesotho_i]\n",
    "    lesotho_i=lesotho_i+1\n",
    "lesotho_ratio_analysis_project=df_lesotho_VERs_agw_final/lesotho_VERs\n",
    "lesotho_ratio_analysis_project_mean=np.mean(df_lesotho_VERs_agw_final, axis=1)/lesotho_VERs\n",
    "lesotho_ratio_project_analysis_mean=lesotho_ratio_analysis_project_mean**-1\n",
    "\n",
    "lesotho_sd_ratio_analysis_project=np.std(lesotho_ratio_analysis_project, axis=1)\n",
    "\n",
    "\n",
    "\n",
    "lesotho_ratio_analysis_project_lower=lesotho_ratio_analysis_project_mean - z*lesotho_sd_ratio_analysis_project/np.sqrt(lesotho_N)\n",
    "lesotho_ratio_analysis_project_upper=lesotho_ratio_analysis_project_mean + z*lesotho_sd_ratio_analysis_project/np.sqrt(lesotho_N)\n",
    "\n",
    "lesotho_ratio_project_analysis_lower=lesotho_ratio_analysis_project_lower**-1\n",
    "lesotho_ratio_project_analysis_upper=lesotho_ratio_analysis_project_upper**-1\n",
    "\n",
    "print(lesotho_ratio_analysis_project_lower)\n",
    "print(lesotho_ratio_analysis_project_mean)\n",
    "print(lesotho_ratio_analysis_project_upper)\n",
    "print(lesotho_ratio_project_analysis_lower)\n",
    "print(lesotho_ratio_project_analysis_mean)\n",
    "print(lesotho_ratio_project_analysis_upper)\n",
    "\n",
    "print(round(lesotho_ratio_analysis_project_mean[0],2),'[',round(lesotho_ratio_analysis_project_lower[0],2),',',round(lesotho_ratio_analysis_project_upper[0],2),']')\n",
    "print(round(lesotho_ratio_project_analysis_mean[0],2),'[',round(lesotho_ratio_project_analysis_upper[0],2),',',round(lesotho_ratio_project_analysis_lower[0],2),']')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 45,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "5\n",
      "0    0.093181\n",
      "dtype: float64\n",
      "0    0.104579\n",
      "dtype: float64\n",
      "0    0.115976\n",
      "dtype: float64\n",
      "0    10.731764\n",
      "dtype: float64\n",
      "0    9.562166\n",
      "dtype: float64\n",
      "0    8.622451\n",
      "dtype: float64\n",
      "0.1 [ 0.09 , 0.12 ]\n",
      "9.56 [ 8.62 , 10.73 ]\n"
     ]
    }
   ],
   "source": [
    "df_malawi_VERs_agw=pd.DataFrame(malawi_VERs_agw).transpose()\n",
    "print(df_malawi_VERs_agw.shape[1])\n",
    "malawi_i=0\n",
    "df_malawi_VERs_agw_final=pd.DataFrame(np.zeros((1, 10000)))\n",
    "while malawi_i<df_malawi_VERs_agw.shape[1]:\n",
    "    df_malawi_VERs_agw_final=df_malawi_VERs_agw_final+df_malawi_VERs_agw[malawi_i]\n",
    "    malawi_i=malawi_i+1\n",
    "\n",
    "malawi_ratio_analysis_project=df_malawi_VERs_agw_final/malawi_VERs\n",
    "malawi_ratio_analysis_project_mean=np.mean(df_malawi_VERs_agw_final, axis=1)/malawi_VERs\n",
    "malawi_ratio_project_analysis_mean=malawi_ratio_analysis_project_mean**-1\n",
    "\n",
    "malawi_sd_ratio_analysis_project=np.std(malawi_ratio_analysis_project, axis=1)\n",
    "\n",
    "\n",
    "\n",
    "malawi_ratio_analysis_project_lower=malawi_ratio_analysis_project_mean - z*malawi_sd_ratio_analysis_project/np.sqrt(malawi_N)\n",
    "malawi_ratio_analysis_project_upper=malawi_ratio_analysis_project_mean + z*malawi_sd_ratio_analysis_project/np.sqrt(malawi_N)\n",
    "\n",
    "malawi_ratio_project_analysis_lower=malawi_ratio_analysis_project_lower**-1\n",
    "malawi_ratio_project_analysis_upper=malawi_ratio_analysis_project_upper**-1\n",
    "\n",
    "print(malawi_ratio_analysis_project_lower)\n",
    "print(malawi_ratio_analysis_project_mean)\n",
    "print(malawi_ratio_analysis_project_upper)\n",
    "print(malawi_ratio_project_analysis_lower)\n",
    "print(malawi_ratio_project_analysis_mean)\n",
    "print(malawi_ratio_project_analysis_upper)\n",
    "\n",
    "print(round(malawi_ratio_analysis_project_mean[0],2),'[',round(malawi_ratio_analysis_project_lower[0],2),',',round(malawi_ratio_analysis_project_upper[0],2),']')\n",
    "print(round(malawi_ratio_project_analysis_mean[0],2),'[',round(malawi_ratio_project_analysis_upper[0],2),',',round(malawi_ratio_project_analysis_lower[0],2),']')\n",
    "\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 46,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "2\n",
      "0    0.002562\n",
      "dtype: float64\n",
      "0    0.010802\n",
      "dtype: float64\n",
      "0    0.019042\n",
      "dtype: float64\n",
      "0    390.388989\n",
      "dtype: float64\n",
      "0    92.577745\n",
      "dtype: float64\n",
      "0    52.515722\n",
      "dtype: float64\n",
      "0.01 [ 0.0 , 0.02 ]\n",
      "92.58 [ 52.52 , 390.39 ]\n"
     ]
    }
   ],
   "source": [
    "df_mali_VERs_agw=pd.DataFrame(mali_VERs_agw).transpose()\n",
    "print(df_mali_VERs_agw.shape[1])\n",
    "mali_i=0\n",
    "df_mali_VERs_agw_final=pd.DataFrame(np.zeros((1, 10000)))\n",
    "while mali_i<df_mali_VERs_agw.shape[1]:\n",
    "    df_mali_VERs_agw_final=df_mali_VERs_agw_final+df_mali_VERs_agw[mali_i]\n",
    "    mali_i=mali_i+1\n",
    "\n",
    "mali_ratio_analysis_project=df_mali_VERs_agw_final/mali_VERs\n",
    "mali_ratio_analysis_project_mean=np.mean(df_mali_VERs_agw_final, axis=1)/mali_VERs\n",
    "mali_ratio_project_analysis_mean=mali_ratio_analysis_project_mean**-1\n",
    "\n",
    "mali_sd_ratio_analysis_project=np.std(mali_ratio_analysis_project, axis=1)\n",
    "\n",
    "\n",
    "\n",
    "mali_ratio_analysis_project_lower=mali_ratio_analysis_project_mean - z*mali_sd_ratio_analysis_project/np.sqrt(mali_N)\n",
    "mali_ratio_analysis_project_upper=mali_ratio_analysis_project_mean + z*mali_sd_ratio_analysis_project/np.sqrt(mali_N)\n",
    "\n",
    "mali_ratio_project_analysis_lower=mali_ratio_analysis_project_lower**-1\n",
    "mali_ratio_project_analysis_upper=mali_ratio_analysis_project_upper**-1\n",
    "\n",
    "print(mali_ratio_analysis_project_lower)\n",
    "print(mali_ratio_analysis_project_mean)\n",
    "print(mali_ratio_analysis_project_upper)\n",
    "print(mali_ratio_project_analysis_lower)\n",
    "print(mali_ratio_project_analysis_mean)\n",
    "print(mali_ratio_project_analysis_upper)\n",
    "\n",
    "print(round(mali_ratio_analysis_project_mean[0],2),'[',round(mali_ratio_analysis_project_lower[0],2),',',round(mali_ratio_analysis_project_upper[0],2),']')\n",
    "print(round(mali_ratio_project_analysis_mean[0],2),'[',round(mali_ratio_project_analysis_upper[0],2),',',round(mali_ratio_project_analysis_lower[0],2),']')\n",
    "\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 47,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "32\n",
      "0    0.029389\n",
      "dtype: float64\n",
      "0    0.032539\n",
      "dtype: float64\n",
      "0    0.035688\n",
      "dtype: float64\n",
      "0    34.02576\n",
      "dtype: float64\n",
      "0    30.732374\n",
      "dtype: float64\n",
      "0    28.020266\n",
      "dtype: float64\n",
      "0.03 [ 0.03 , 0.04 ]\n",
      "30.73 [ 28.02 , 34.03 ]\n"
     ]
    }
   ],
   "source": [
    "df_mexico_VERs_agw=pd.DataFrame(mexico_VERs_agw).transpose()\n",
    "print(df_mexico_VERs_agw.shape[1])\n",
    "mexico_i=0\n",
    "df_mexico_VERs_agw_final=pd.DataFrame(np.zeros((1, 10000)))\n",
    "while mexico_i<df_mexico_VERs_agw.shape[1]:\n",
    "    df_mexico_VERs_agw_final=df_mexico_VERs_agw_final+df_mexico_VERs_agw[mexico_i]\n",
    "    mexico_i=mexico_i+1\n",
    "mexico_ratio_analysis_project=df_mexico_VERs_agw_final/mexico_VERs\n",
    "mexico_ratio_analysis_project_mean=np.mean(df_mexico_VERs_agw_final, axis=1)/mexico_VERs\n",
    "mexico_ratio_project_analysis_mean=mexico_ratio_analysis_project_mean**-1\n",
    "\n",
    "mexico_sd_ratio_analysis_project=np.std(mexico_ratio_analysis_project, axis=1)\n",
    "\n",
    "\n",
    "\n",
    "mexico_ratio_analysis_project_lower=mexico_ratio_analysis_project_mean - z*mexico_sd_ratio_analysis_project/np.sqrt(mexico_N)\n",
    "mexico_ratio_analysis_project_upper=mexico_ratio_analysis_project_mean + z*mexico_sd_ratio_analysis_project/np.sqrt(mexico_N)\n",
    "\n",
    "mexico_ratio_project_analysis_lower=mexico_ratio_analysis_project_lower**-1\n",
    "mexico_ratio_project_analysis_upper=mexico_ratio_analysis_project_upper**-1\n",
    "\n",
    "print(mexico_ratio_analysis_project_lower)\n",
    "print(mexico_ratio_analysis_project_mean)\n",
    "print(mexico_ratio_analysis_project_upper)\n",
    "print(mexico_ratio_project_analysis_lower)\n",
    "print(mexico_ratio_project_analysis_mean)\n",
    "print(mexico_ratio_project_analysis_upper)\n",
    "\n",
    "print(round(mexico_ratio_analysis_project_mean[0],2),'[',round(mexico_ratio_analysis_project_lower[0],2),',',round(mexico_ratio_analysis_project_upper[0],2),']')\n",
    "print(round(mexico_ratio_project_analysis_mean[0],2),'[',round(mexico_ratio_project_analysis_upper[0],2),',',round(mexico_ratio_project_analysis_lower[0],2),']')\n",
    "\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 48,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "3\n",
      "0    0.075192\n",
      "dtype: float64\n",
      "0    0.10126\n",
      "dtype: float64\n",
      "0    0.127327\n",
      "dtype: float64\n",
      "0    13.299287\n",
      "dtype: float64\n",
      "0    9.87561\n",
      "dtype: float64\n",
      "0    7.853786\n",
      "dtype: float64\n",
      "0.1 [ 0.08 , 0.13 ]\n",
      "9.88 [ 7.85 , 13.3 ]\n"
     ]
    }
   ],
   "source": [
    "df_mozambique_VERs_agw=pd.DataFrame(mozambique_VERs_agw).transpose()\n",
    "print(df_mozambique_VERs_agw.shape[1])\n",
    "mozambique_i=0\n",
    "df_mozambique_VERs_agw_final=pd.DataFrame(np.zeros((1, 10000)))\n",
    "while mozambique_i<df_mozambique_VERs_agw.shape[1]:\n",
    "    df_mozambique_VERs_agw_final=df_mozambique_VERs_agw_final+df_mozambique_VERs_agw[mozambique_i]\n",
    "    mozambique_i=mozambique_i+1\n",
    "\n",
    "mozambique_ratio_analysis_project=df_mozambique_VERs_agw_final/mozambique_VERs\n",
    "mozambique_ratio_analysis_project_mean=np.mean(df_mozambique_VERs_agw_final, axis=1)/mozambique_VERs\n",
    "mozambique_ratio_project_analysis_mean=mozambique_ratio_analysis_project_mean**-1\n",
    "\n",
    "mozambique_sd_ratio_analysis_project=np.std(mozambique_ratio_analysis_project, axis=1)\n",
    "\n",
    "\n",
    "\n",
    "mozambique_ratio_analysis_project_lower=mozambique_ratio_analysis_project_mean - z*mozambique_sd_ratio_analysis_project/np.sqrt(mozambique_N)\n",
    "mozambique_ratio_analysis_project_upper=mozambique_ratio_analysis_project_mean + z*mozambique_sd_ratio_analysis_project/np.sqrt(mozambique_N)\n",
    "\n",
    "mozambique_ratio_project_analysis_lower=mozambique_ratio_analysis_project_lower**-1\n",
    "mozambique_ratio_project_analysis_upper=mozambique_ratio_analysis_project_upper**-1\n",
    "\n",
    "print(mozambique_ratio_analysis_project_lower)\n",
    "print(mozambique_ratio_analysis_project_mean)\n",
    "print(mozambique_ratio_analysis_project_upper)\n",
    "print(mozambique_ratio_project_analysis_lower)\n",
    "print(mozambique_ratio_project_analysis_mean)\n",
    "print(mozambique_ratio_project_analysis_upper)\n",
    "\n",
    "print(round(mozambique_ratio_analysis_project_mean[0],2),'[',round(mozambique_ratio_analysis_project_lower[0],2),',',round(mozambique_ratio_analysis_project_upper[0],2),']')\n",
    "print(round(mozambique_ratio_project_analysis_mean[0],2),'[',round(mozambique_ratio_project_analysis_upper[0],2),',',round(mozambique_ratio_project_analysis_lower[0],2),']')\n",
    "\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 49,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "4\n",
      "0    0.006273\n",
      "dtype: float64\n",
      "0    0.013093\n",
      "dtype: float64\n",
      "0    0.019912\n",
      "dtype: float64\n",
      "0    159.402599\n",
      "dtype: float64\n",
      "0    76.377174\n",
      "dtype: float64\n",
      "0    50.219936\n",
      "dtype: float64\n",
      "0.01 [ 0.01 , 0.02 ]\n",
      "76.38 [ 50.22 , 159.4 ]\n"
     ]
    }
   ],
   "source": [
    "df_myanmar_VERs_agw=pd.DataFrame(myanmar_VERs_agw).transpose()\n",
    "print(df_myanmar_VERs_agw.shape[1])\n",
    "myanmar_i=0\n",
    "df_myanmar_VERs_agw_final=pd.DataFrame(np.zeros((1, 10000)))\n",
    "while myanmar_i<df_myanmar_VERs_agw.shape[1]:\n",
    "    df_myanmar_VERs_agw_final=df_myanmar_VERs_agw_final+df_myanmar_VERs_agw[myanmar_i]\n",
    "    myanmar_i=myanmar_i+1\n",
    "\n",
    "myanmar_ratio_analysis_project=df_myanmar_VERs_agw_final/myanmar_VERs\n",
    "myanmar_ratio_analysis_project_mean=np.mean(df_myanmar_VERs_agw_final, axis=1)/myanmar_VERs\n",
    "myanmar_ratio_project_analysis_mean=myanmar_ratio_analysis_project_mean**-1\n",
    "\n",
    "myanmar_sd_ratio_analysis_project=np.std(myanmar_ratio_analysis_project, axis=1)\n",
    "\n",
    "\n",
    "\n",
    "myanmar_ratio_analysis_project_lower=myanmar_ratio_analysis_project_mean - z*myanmar_sd_ratio_analysis_project/np.sqrt(myanmar_N)\n",
    "myanmar_ratio_analysis_project_upper=myanmar_ratio_analysis_project_mean + z*myanmar_sd_ratio_analysis_project/np.sqrt(myanmar_N)\n",
    "\n",
    "myanmar_ratio_project_analysis_lower=myanmar_ratio_analysis_project_lower**-1\n",
    "myanmar_ratio_project_analysis_upper=myanmar_ratio_analysis_project_upper**-1\n",
    "\n",
    "print(myanmar_ratio_analysis_project_lower)\n",
    "print(myanmar_ratio_analysis_project_mean)\n",
    "print(myanmar_ratio_analysis_project_upper)\n",
    "print(myanmar_ratio_project_analysis_lower)\n",
    "print(myanmar_ratio_project_analysis_mean)\n",
    "print(myanmar_ratio_project_analysis_upper)\n",
    "\n",
    "print(round(myanmar_ratio_analysis_project_mean[0],2),'[',round(myanmar_ratio_analysis_project_lower[0],2),',',round(myanmar_ratio_analysis_project_upper[0],2),']')\n",
    "print(round(myanmar_ratio_project_analysis_mean[0],2),'[',round(myanmar_ratio_project_analysis_upper[0],2),',',round(myanmar_ratio_project_analysis_lower[0],2),']')\n",
    "\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 50,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "5\n",
      "0    0.062441\n",
      "dtype: float64\n",
      "0    0.117005\n",
      "dtype: float64\n",
      "0    0.171569\n",
      "dtype: float64\n",
      "0    16.015242\n",
      "dtype: float64\n",
      "0    8.546655\n",
      "dtype: float64\n",
      "0    5.828553\n",
      "dtype: float64\n",
      "0.12 [ 0.06 , 0.17 ]\n",
      "8.55 [ 5.83 , 16.02 ]\n"
     ]
    }
   ],
   "source": [
    "df_nepal_VERs_agw=pd.DataFrame(nepal_VERs_agw).transpose()\n",
    "print(df_nepal_VERs_agw.shape[1])\n",
    "nepal_i=0\n",
    "df_nepal_VERs_agw_final=pd.DataFrame(np.zeros((1, 10000)))\n",
    "while nepal_i<df_nepal_VERs_agw.shape[1]:\n",
    "    df_nepal_VERs_agw_final=df_nepal_VERs_agw_final+df_nepal_VERs_agw[nepal_i]\n",
    "    nepal_i=nepal_i+1\n",
    "    \n",
    "#df_nepal_VERs_agw_final.to_excel('df_nepal_VERs_agw_final.xlsx')\n",
    "#df_nepal_VERs_agw.to_excel('df_nepal_VERs_agw.xlsx')\n",
    "\n",
    "nepal_ratio_analysis_project=df_nepal_VERs_agw_final/nepal_VERs\n",
    "nepal_ratio_analysis_project_mean=np.mean(df_nepal_VERs_agw_final, axis=1)/nepal_VERs\n",
    "nepal_ratio_project_analysis_mean=nepal_ratio_analysis_project_mean**-1\n",
    "\n",
    "nepal_sd_ratio_analysis_project=np.std(nepal_ratio_analysis_project, axis=1)\n",
    "\n",
    "\n",
    "\n",
    "nepal_ratio_analysis_project_lower=nepal_ratio_analysis_project_mean - z*nepal_sd_ratio_analysis_project/np.sqrt(nepal_N)\n",
    "nepal_ratio_analysis_project_upper=nepal_ratio_analysis_project_mean + z*nepal_sd_ratio_analysis_project/np.sqrt(nepal_N)\n",
    "\n",
    "nepal_ratio_project_analysis_lower=nepal_ratio_analysis_project_lower**-1\n",
    "nepal_ratio_project_analysis_upper=nepal_ratio_analysis_project_upper**-1\n",
    "\n",
    "print(nepal_ratio_analysis_project_lower)\n",
    "print(nepal_ratio_analysis_project_mean)\n",
    "print(nepal_ratio_analysis_project_upper)\n",
    "print(nepal_ratio_project_analysis_lower)\n",
    "print(nepal_ratio_project_analysis_mean)\n",
    "print(nepal_ratio_project_analysis_upper)\n",
    "\n",
    "print(round(nepal_ratio_analysis_project_mean[0],2),'[',round(nepal_ratio_analysis_project_lower[0],2),',',round(nepal_ratio_analysis_project_upper[0],2),']')\n",
    "print(round(nepal_ratio_project_analysis_mean[0],2),'[',round(nepal_ratio_project_analysis_upper[0],2),',',round(nepal_ratio_project_analysis_lower[0],2),']')\n",
    "\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 51,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "3\n",
      "0    0.087368\n",
      "dtype: float64\n",
      "0    0.105354\n",
      "dtype: float64\n",
      "0    0.123341\n",
      "dtype: float64\n",
      "0    11.445813\n",
      "dtype: float64\n",
      "0    9.491773\n",
      "dtype: float64\n",
      "0    8.107631\n",
      "dtype: float64\n",
      "0.11 [ 0.09 , 0.12 ]\n",
      "9.49 [ 8.11 , 11.45 ]\n"
     ]
    }
   ],
   "source": [
    "df_nigeria_VERs_agw=pd.DataFrame(nigeria_VERs_agw).transpose()\n",
    "print(df_nigeria_VERs_agw.shape[1])\n",
    "nigeria_i=0\n",
    "df_nigeria_VERs_agw_final=pd.DataFrame(np.zeros((1, 10000)))\n",
    "while nigeria_i<df_nigeria_VERs_agw.shape[1]:\n",
    "    df_nigeria_VERs_agw_final=df_nigeria_VERs_agw_final+df_nigeria_VERs_agw[nigeria_i]\n",
    "    nigeria_i=nigeria_i+1\n",
    "\n",
    "nigeria_ratio_analysis_project=df_nigeria_VERs_agw_final/nigeria_VERs\n",
    "nigeria_ratio_analysis_project_mean=np.mean(df_nigeria_VERs_agw_final, axis=1)/nigeria_VERs\n",
    "nigeria_ratio_project_analysis_mean=nigeria_ratio_analysis_project_mean**-1\n",
    "\n",
    "nigeria_sd_ratio_analysis_project=np.std(nigeria_ratio_analysis_project, axis=1)\n",
    "\n",
    "\n",
    "\n",
    "nigeria_ratio_analysis_project_lower=nigeria_ratio_analysis_project_mean - z*nigeria_sd_ratio_analysis_project/np.sqrt(nigeria_N)\n",
    "nigeria_ratio_analysis_project_upper=nigeria_ratio_analysis_project_mean + z*nigeria_sd_ratio_analysis_project/np.sqrt(nigeria_N)\n",
    "\n",
    "nigeria_ratio_project_analysis_lower=nigeria_ratio_analysis_project_lower**-1\n",
    "nigeria_ratio_project_analysis_upper=nigeria_ratio_analysis_project_upper**-1\n",
    "\n",
    "print(nigeria_ratio_analysis_project_lower)\n",
    "print(nigeria_ratio_analysis_project_mean)\n",
    "print(nigeria_ratio_analysis_project_upper)\n",
    "print(nigeria_ratio_project_analysis_lower)\n",
    "print(nigeria_ratio_project_analysis_mean)\n",
    "print(nigeria_ratio_project_analysis_upper)\n",
    "\n",
    "print(round(nigeria_ratio_analysis_project_mean[0],2),'[',round(nigeria_ratio_analysis_project_lower[0],2),',',round(nigeria_ratio_analysis_project_upper[0],2),']')\n",
    "print(round(nigeria_ratio_project_analysis_mean[0],2),'[',round(nigeria_ratio_project_analysis_upper[0],2),',',round(nigeria_ratio_project_analysis_lower[0],2),']')\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 52,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "1\n",
      "0    0.028626\n",
      "dtype: float64\n",
      "0    0.049559\n",
      "dtype: float64\n",
      "0    0.070493\n",
      "dtype: float64\n",
      "0    34.933125\n",
      "dtype: float64\n",
      "0    20.177815\n",
      "dtype: float64\n",
      "0    14.185879\n",
      "dtype: float64\n",
      "0.05 [ 0.03 , 0.07 ]\n",
      "20.18 [ 14.19 , 34.93 ]\n"
     ]
    }
   ],
   "source": [
    "df_peru_VERs_agw=pd.DataFrame(peru_VERs_agw).transpose()\n",
    "print(df_peru_VERs_agw.shape[1])\n",
    "peru_i=0\n",
    "df_peru_VERs_agw_final=pd.DataFrame(np.zeros((1, 10000)))\n",
    "while peru_i<df_peru_VERs_agw.shape[1]:\n",
    "    df_peru_VERs_agw_final=df_peru_VERs_agw_final+df_peru_VERs_agw[peru_i]\n",
    "    peru_i=peru_i+1\n",
    "\n",
    "peru_ratio_analysis_project=df_peru_VERs_agw_final/peru_VERs\n",
    "peru_ratio_analysis_project_mean=np.mean(df_peru_VERs_agw_final, axis=1)/peru_VERs\n",
    "peru_ratio_project_analysis_mean=peru_ratio_analysis_project_mean**-1\n",
    "\n",
    "peru_sd_ratio_analysis_project=np.std(peru_ratio_analysis_project, axis=1)\n",
    "\n",
    "\n",
    "\n",
    "peru_ratio_analysis_project_lower=peru_ratio_analysis_project_mean - z*peru_sd_ratio_analysis_project/np.sqrt(peru_N)\n",
    "peru_ratio_analysis_project_upper=peru_ratio_analysis_project_mean + z*peru_sd_ratio_analysis_project/np.sqrt(peru_N)\n",
    "\n",
    "peru_ratio_project_analysis_lower=peru_ratio_analysis_project_lower**-1\n",
    "peru_ratio_project_analysis_upper=peru_ratio_analysis_project_upper**-1\n",
    "\n",
    "print(peru_ratio_analysis_project_lower)\n",
    "print(peru_ratio_analysis_project_mean)\n",
    "print(peru_ratio_analysis_project_upper)\n",
    "print(peru_ratio_project_analysis_lower)\n",
    "print(peru_ratio_project_analysis_mean)\n",
    "print(peru_ratio_project_analysis_upper)\n",
    "\n",
    "print(round(peru_ratio_analysis_project_mean[0],2),'[',round(peru_ratio_analysis_project_lower[0],2),',',round(peru_ratio_analysis_project_upper[0],2),']')\n",
    "print(round(peru_ratio_project_analysis_mean[0],2),'[',round(peru_ratio_project_analysis_upper[0],2),',',round(peru_ratio_project_analysis_lower[0],2),']')\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 53,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "13\n",
      "0    0.14208\n",
      "dtype: float64\n",
      "0    0.14947\n",
      "dtype: float64\n",
      "0    0.156861\n",
      "dtype: float64\n",
      "0    7.038306\n",
      "dtype: float64\n",
      "0    6.690291\n",
      "dtype: float64\n",
      "0    6.37507\n",
      "dtype: float64\n",
      "0.15 [ 0.14 , 0.16 ]\n",
      "6.69 [ 6.38 , 7.04 ]\n"
     ]
    }
   ],
   "source": [
    "df_rwanda_VERs_agw=pd.DataFrame(rwanda_VERs_agw).transpose()\n",
    "print(df_rwanda_VERs_agw.shape[1])\n",
    "rwanda_i=0\n",
    "df_rwanda_VERs_agw_final=pd.DataFrame(np.zeros((1, 10000)))\n",
    "while rwanda_i<df_rwanda_VERs_agw.shape[1]:\n",
    "    df_rwanda_VERs_agw_final=df_rwanda_VERs_agw_final+df_rwanda_VERs_agw[rwanda_i]\n",
    "    rwanda_i=rwanda_i+1\n",
    "\n",
    "    rwanda_ratio_analysis_project=df_rwanda_VERs_agw_final/rwanda_VERs\n",
    "rwanda_ratio_analysis_project_mean=np.mean(df_rwanda_VERs_agw_final, axis=1)/rwanda_VERs\n",
    "rwanda_ratio_project_analysis_mean=rwanda_ratio_analysis_project_mean**-1\n",
    "\n",
    "rwanda_sd_ratio_analysis_project=np.std(rwanda_ratio_analysis_project, axis=1)\n",
    "\n",
    "\n",
    "\n",
    "rwanda_ratio_analysis_project_lower=rwanda_ratio_analysis_project_mean - z*rwanda_sd_ratio_analysis_project/np.sqrt(rwanda_N)\n",
    "rwanda_ratio_analysis_project_upper=rwanda_ratio_analysis_project_mean + z*rwanda_sd_ratio_analysis_project/np.sqrt(rwanda_N)\n",
    "\n",
    "rwanda_ratio_project_analysis_lower=rwanda_ratio_analysis_project_lower**-1\n",
    "rwanda_ratio_project_analysis_upper=rwanda_ratio_analysis_project_upper**-1\n",
    "\n",
    "print(rwanda_ratio_analysis_project_lower)\n",
    "print(rwanda_ratio_analysis_project_mean)\n",
    "print(rwanda_ratio_analysis_project_upper)\n",
    "print(rwanda_ratio_project_analysis_lower)\n",
    "print(rwanda_ratio_project_analysis_mean)\n",
    "print(rwanda_ratio_project_analysis_upper)\n",
    "\n",
    "print(round(rwanda_ratio_analysis_project_mean[0],2),'[',round(rwanda_ratio_analysis_project_lower[0],2),',',round(rwanda_ratio_analysis_project_upper[0],2),']')\n",
    "print(round(rwanda_ratio_project_analysis_mean[0],2),'[',round(rwanda_ratio_project_analysis_upper[0],2),',',round(rwanda_ratio_project_analysis_lower[0],2),']')\n",
    "\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 54,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "2\n",
      "0    0.014892\n",
      "dtype: float64\n",
      "0    0.028456\n",
      "dtype: float64\n",
      "0    0.042021\n",
      "dtype: float64\n",
      "0    67.149981\n",
      "dtype: float64\n",
      "0    35.14161\n",
      "dtype: float64\n",
      "0    23.797879\n",
      "dtype: float64\n",
      "0.03 [ 0.01 , 0.04 ]\n",
      "35.14 [ 23.8 , 67.15 ]\n"
     ]
    }
   ],
   "source": [
    "df_sudan_VERs_agw=pd.DataFrame(sudan_VERs_agw).transpose()\n",
    "print(df_sudan_VERs_agw.shape[1])\n",
    "sudan_i=0\n",
    "df_sudan_VERs_agw_final=pd.DataFrame(np.zeros((1, 10000)))\n",
    "while sudan_i<df_sudan_VERs_agw.shape[1]:\n",
    "    df_sudan_VERs_agw_final=df_sudan_VERs_agw_final+df_sudan_VERs_agw[sudan_i]\n",
    "    sudan_i=sudan_i+1\n",
    "\n",
    "sudan_ratio_analysis_project=df_sudan_VERs_agw_final/sudan_VERs\n",
    "sudan_ratio_analysis_project_mean=np.mean(df_sudan_VERs_agw_final, axis=1)/sudan_VERs\n",
    "sudan_ratio_project_analysis_mean=sudan_ratio_analysis_project_mean**-1\n",
    "\n",
    "sudan_sd_ratio_analysis_project=np.std(sudan_ratio_analysis_project, axis=1)\n",
    "\n",
    "\n",
    "\n",
    "sudan_ratio_analysis_project_lower=sudan_ratio_analysis_project_mean - z*sudan_sd_ratio_analysis_project/np.sqrt(sudan_N)\n",
    "sudan_ratio_analysis_project_upper=sudan_ratio_analysis_project_mean + z*sudan_sd_ratio_analysis_project/np.sqrt(sudan_N)\n",
    "\n",
    "sudan_ratio_project_analysis_lower=sudan_ratio_analysis_project_lower**-1\n",
    "sudan_ratio_project_analysis_upper=sudan_ratio_analysis_project_upper**-1\n",
    "\n",
    "print(sudan_ratio_analysis_project_lower)\n",
    "print(sudan_ratio_analysis_project_mean)\n",
    "print(sudan_ratio_analysis_project_upper)\n",
    "print(sudan_ratio_project_analysis_lower)\n",
    "print(sudan_ratio_project_analysis_mean)\n",
    "print(sudan_ratio_project_analysis_upper)\n",
    "\n",
    "print(round(sudan_ratio_analysis_project_mean[0],2),'[',round(sudan_ratio_analysis_project_lower[0],2),',',round(sudan_ratio_analysis_project_upper[0],2),']')\n",
    "print(round(sudan_ratio_project_analysis_mean[0],2),'[',round(sudan_ratio_project_analysis_upper[0],2),',',round(sudan_ratio_project_analysis_lower[0],2),']')\n",
    "\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 55,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "6\n",
      "0    0.073567\n",
      "dtype: float64\n",
      "0    0.084269\n",
      "dtype: float64\n",
      "0    0.094972\n",
      "dtype: float64\n",
      "0    13.593054\n",
      "dtype: float64\n",
      "0    11.866706\n",
      "dtype: float64\n",
      "0    10.529443\n",
      "dtype: float64\n",
      "0.08 [ 0.07 , 0.09 ]\n",
      "11.87 [ 10.53 , 13.59 ]\n"
     ]
    }
   ],
   "source": [
    "df_uganda_VERs_agw=pd.DataFrame(uganda_VERs_agw).transpose()\n",
    "print(df_uganda_VERs_agw.shape[1])\n",
    "uganda_i=0\n",
    "df_uganda_VERs_agw_final=pd.DataFrame(np.zeros((1, 10000)))\n",
    "while uganda_i<df_uganda_VERs_agw.shape[1]:\n",
    "    df_uganda_VERs_agw_final=df_uganda_VERs_agw_final+df_uganda_VERs_agw[uganda_i]\n",
    "    uganda_i=uganda_i+1\n",
    "uganda_ratio_analysis_project=df_uganda_VERs_agw_final/uganda_VERs\n",
    "uganda_ratio_analysis_project_mean=np.mean(df_uganda_VERs_agw_final, axis=1)/uganda_VERs\n",
    "uganda_ratio_project_analysis_mean=uganda_ratio_analysis_project_mean**-1\n",
    "\n",
    "uganda_sd_ratio_analysis_project=np.std(uganda_ratio_analysis_project, axis=1)\n",
    "\n",
    "\n",
    "\n",
    "uganda_ratio_analysis_project_lower=uganda_ratio_analysis_project_mean - z*uganda_sd_ratio_analysis_project/np.sqrt(uganda_N)\n",
    "uganda_ratio_analysis_project_upper=uganda_ratio_analysis_project_mean + z*uganda_sd_ratio_analysis_project/np.sqrt(uganda_N)\n",
    "\n",
    "uganda_ratio_project_analysis_lower=uganda_ratio_analysis_project_lower**-1\n",
    "uganda_ratio_project_analysis_upper=uganda_ratio_analysis_project_upper**-1\n",
    "\n",
    "print(uganda_ratio_analysis_project_lower)\n",
    "print(uganda_ratio_analysis_project_mean)\n",
    "print(uganda_ratio_analysis_project_upper)\n",
    "print(uganda_ratio_project_analysis_lower)\n",
    "print(uganda_ratio_project_analysis_mean)\n",
    "print(uganda_ratio_project_analysis_upper)\n",
    "\n",
    "print(round(uganda_ratio_analysis_project_mean[0],2),'[',round(uganda_ratio_analysis_project_lower[0],2),',',round(uganda_ratio_analysis_project_upper[0],2),']')\n",
    "print(round(uganda_ratio_project_analysis_mean[0],2),'[',round(uganda_ratio_project_analysis_upper[0],2),',',round(uganda_ratio_project_analysis_lower[0],2),']')\n",
    "\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 56,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "2\n",
      "0    0.275308\n",
      "dtype: float64\n",
      "0    0.295176\n",
      "dtype: float64\n",
      "0    0.315044\n",
      "dtype: float64\n",
      "0    3.63229\n",
      "dtype: float64\n",
      "0    3.387809\n",
      "dtype: float64\n",
      "0    3.174163\n",
      "dtype: float64\n",
      "0.3 [ 0.28 , 0.32 ]\n",
      "3.39 [ 3.17 , 3.63 ]\n"
     ]
    }
   ],
   "source": [
    "df_zambia_VERs_agw=pd.DataFrame(zambia_VERs_agw).transpose()\n",
    "print(df_zambia_VERs_agw.shape[1])\n",
    "zambia_i=0\n",
    "df_zambia_VERs_agw_final=pd.DataFrame(np.zeros((1, 10000)))\n",
    "while zambia_i<df_zambia_VERs_agw.shape[1]:\n",
    "    df_zambia_VERs_agw_final=df_zambia_VERs_agw_final+df_zambia_VERs_agw[zambia_i]\n",
    "    zambia_i=zambia_i+1\n",
    "\n",
    "zambia_ratio_analysis_project=df_zambia_VERs_agw_final/zambia_VERs\n",
    "zambia_ratio_analysis_project_mean=np.mean(df_zambia_VERs_agw_final, axis=1)/zambia_VERs\n",
    "zambia_ratio_project_analysis_mean=zambia_ratio_analysis_project_mean**-1\n",
    "\n",
    "zambia_sd_ratio_analysis_project=np.std(zambia_ratio_analysis_project, axis=1)\n",
    "\n",
    "\n",
    "\n",
    "zambia_ratio_analysis_project_lower=zambia_ratio_analysis_project_mean - z*zambia_sd_ratio_analysis_project/np.sqrt(zambia_N)\n",
    "zambia_ratio_analysis_project_upper=zambia_ratio_analysis_project_mean + z*zambia_sd_ratio_analysis_project/np.sqrt(zambia_N)\n",
    "\n",
    "zambia_ratio_project_analysis_lower=zambia_ratio_analysis_project_lower**-1\n",
    "zambia_ratio_project_analysis_upper=zambia_ratio_analysis_project_upper**-1\n",
    "\n",
    "print(zambia_ratio_analysis_project_lower)\n",
    "print(zambia_ratio_analysis_project_mean)\n",
    "print(zambia_ratio_analysis_project_upper)\n",
    "print(zambia_ratio_project_analysis_lower)\n",
    "print(zambia_ratio_project_analysis_mean)\n",
    "print(zambia_ratio_project_analysis_upper)\n",
    "\n",
    "print(round(zambia_ratio_analysis_project_mean[0],2),'[',round(zambia_ratio_analysis_project_lower[0],2),',',round(zambia_ratio_analysis_project_upper[0],2),']')\n",
    "print(round(zambia_ratio_project_analysis_mean[0],2),'[',round(zambia_ratio_project_analysis_upper[0],2),',',round(zambia_ratio_project_analysis_lower[0],2),']')\n",
    "\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 57,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "151\n",
      "1\n",
      "2\n",
      "3\n",
      "4\n",
      "5\n",
      "6\n",
      "7\n",
      "8\n",
      "9\n",
      "10\n",
      "11\n",
      "12\n",
      "13\n",
      "14\n",
      "15\n",
      "16\n",
      "17\n",
      "18\n",
      "19\n",
      "20\n",
      "21\n",
      "22\n",
      "23\n",
      "24\n",
      "25\n",
      "26\n",
      "27\n",
      "28\n",
      "29\n",
      "30\n",
      "31\n",
      "32\n",
      "33\n",
      "34\n",
      "35\n",
      "36\n",
      "37\n",
      "38\n",
      "39\n",
      "40\n",
      "41\n",
      "42\n",
      "43\n",
      "44\n",
      "45\n",
      "46\n",
      "47\n",
      "48\n",
      "49\n",
      "50\n",
      "51\n",
      "52\n",
      "53\n",
      "54\n",
      "55\n",
      "56\n",
      "57\n",
      "58\n",
      "59\n",
      "60\n",
      "61\n",
      "62\n",
      "63\n",
      "64\n",
      "65\n",
      "66\n",
      "67\n",
      "68\n",
      "69\n",
      "70\n",
      "71\n",
      "72\n",
      "73\n",
      "74\n",
      "75\n",
      "76\n",
      "77\n",
      "78\n",
      "79\n",
      "80\n",
      "81\n",
      "82\n",
      "83\n",
      "84\n",
      "85\n",
      "86\n",
      "87\n",
      "88\n",
      "89\n",
      "90\n",
      "91\n",
      "92\n",
      "93\n",
      "94\n",
      "95\n",
      "96\n",
      "97\n",
      "98\n",
      "99\n",
      "100\n",
      "101\n",
      "102\n",
      "103\n",
      "104\n",
      "105\n",
      "106\n",
      "107\n",
      "108\n",
      "109\n",
      "110\n",
      "111\n",
      "112\n",
      "113\n",
      "114\n",
      "115\n",
      "116\n",
      "117\n",
      "118\n",
      "119\n",
      "120\n",
      "121\n",
      "122\n",
      "123\n",
      "124\n",
      "125\n",
      "126\n",
      "127\n",
      "128\n",
      "129\n",
      "130\n",
      "131\n",
      "132\n",
      "133\n",
      "134\n",
      "135\n",
      "136\n",
      "137\n",
      "138\n",
      "139\n",
      "140\n",
      "141\n",
      "142\n",
      "143\n",
      "144\n",
      "145\n",
      "146\n",
      "147\n",
      "148\n",
      "149\n",
      "150\n",
      "151\n",
      "0    0.145494\n",
      "dtype: float64\n",
      "0    0.146864\n",
      "dtype: float64\n",
      "0    0.148234\n",
      "dtype: float64\n",
      "0    6.873124\n",
      "dtype: float64\n",
      "0    6.809009\n",
      "dtype: float64\n",
      "0    6.746079\n",
      "dtype: float64\n",
      "0.146864255 [ 0.1454942465 , 0.1482342635 ]\n",
      "6.8090087685 [ 6.7460786485 , 6.8731240186 ]\n"
     ]
    }
   ],
   "source": [
    "\n",
    "df_ALL_VERs_agw=pd.DataFrame(all_VERS_agw).transpose()\n",
    "print(df_ALL_VERs_agw.shape[1])\n",
    "ALL_i=0\n",
    "df_ALL_VERs_agw_final=pd.DataFrame(np.zeros((1, 10000)))\n",
    "while ALL_i<df_ALL_VERs_agw.shape[1]:\n",
    "    df_ALL_VERs_agw_final=df_ALL_VERs_agw_final+df_ALL_VERs_agw[ALL_i]\n",
    "    ALL_i=ALL_i+1\n",
    "    print(ALL_i)\n",
    "\n",
    "\n",
    "\n",
    "\n",
    "ALL_VERs=VERs_all\n",
    "\n",
    "ALL_ratio_analysis_project=df_ALL_VERs_agw_final/sum(ALL_VERs.iloc[:,1])\n",
    "ALL_ratio_analysis_project_mean=np.mean(df_ALL_VERs_agw_final, axis=1)/sum(ALL_VERs.iloc[:,1])\n",
    "ALL_ratio_project_analysis_mean=ALL_ratio_analysis_project_mean**-1\n",
    "\n",
    "ALL_sd_ratio_analysis_project=np.std(ALL_ratio_analysis_project, axis=1)\n",
    "\n",
    "\n",
    "\n",
    "ALL_ratio_analysis_project_lower=ALL_ratio_analysis_project_mean - z*ALL_sd_ratio_analysis_project/np.sqrt(ALL_N)\n",
    "ALL_ratio_analysis_project_upper=ALL_ratio_analysis_project_mean + z*ALL_sd_ratio_analysis_project/np.sqrt(ALL_N)\n",
    "\n",
    "ALL_ratio_project_analysis_lower=ALL_ratio_analysis_project_lower**-1\n",
    "ALL_ratio_project_analysis_upper=ALL_ratio_analysis_project_upper**-1\n",
    "\n",
    "print(ALL_ratio_analysis_project_lower)\n",
    "print(ALL_ratio_analysis_project_mean)\n",
    "print(ALL_ratio_analysis_project_upper)\n",
    "print(ALL_ratio_project_analysis_lower)\n",
    "print(ALL_ratio_project_analysis_mean)\n",
    "print(ALL_ratio_project_analysis_upper)\n",
    "\n",
    "print(round(ALL_ratio_analysis_project_mean[0],10),'[',round(ALL_ratio_analysis_project_lower[0],10),',',round(ALL_ratio_analysis_project_upper[0],10),']')\n",
    "print(round(ALL_ratio_project_analysis_mean[0],10),'[',round(ALL_ratio_project_analysis_upper[0],10),',',round(ALL_ratio_project_analysis_lower[0],10),']')\n"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "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.9.13"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}