Demo of CCC for PSL Meetup#
This notebook provides a demonstration of the Cost-of-Capital-Calculator (CCC) for the PSL Meetup on April 29, 2019.
To run this notebook on your machine, you will need to follow the instructions to install CCC as described in the CCC README here. In particular, you need to:
Install the Anaconda distribution of Python
Install the CCC package by typing
conda install -c conda-forge ccc(orpip install cost-of-captial-calculator) in the command prompt.
Once you follow the above, you will be ready to work with this Jupyter Notebook.
First things first, import necessary packages#
# To install ccc package (if not already):
import sys
if 'ccc' not in sys.modules:
!pip install cost-of-capital-calculator
Requirement already satisfied: cost-of-capital-calculator in /home/runner/miniconda3/envs/ccc-dev/lib/python3.12/site-packages (2.0.0)
Requirement already satisfied: taxcalc in /home/runner/miniconda3/envs/ccc-dev/lib/python3.12/site-packages (from cost-of-capital-calculator) (4.4.0)
Requirement already satisfied: pandas in /home/runner/miniconda3/envs/ccc-dev/lib/python3.12/site-packages (from cost-of-capital-calculator) (2.2.3)
Requirement already satisfied: bokeh in /home/runner/miniconda3/envs/ccc-dev/lib/python3.12/site-packages (from cost-of-capital-calculator) (3.6.2)
Requirement already satisfied: numpy in /home/runner/miniconda3/envs/ccc-dev/lib/python3.12/site-packages (from cost-of-capital-calculator) (2.0.2)
Requirement already satisfied: paramtools in /home/runner/miniconda3/envs/ccc-dev/lib/python3.12/site-packages (from cost-of-capital-calculator) (0.19.0)
Requirement already satisfied: Jinja2>=2.9 in /home/runner/miniconda3/envs/ccc-dev/lib/python3.12/site-packages (from bokeh->cost-of-capital-calculator) (3.1.5)
Requirement already satisfied: contourpy>=1.2 in /home/runner/miniconda3/envs/ccc-dev/lib/python3.12/site-packages (from bokeh->cost-of-capital-calculator) (1.3.1)
Requirement already satisfied: packaging>=16.8 in /home/runner/miniconda3/envs/ccc-dev/lib/python3.12/site-packages (from bokeh->cost-of-capital-calculator) (24.2)
Requirement already satisfied: pillow>=7.1.0 in /home/runner/miniconda3/envs/ccc-dev/lib/python3.12/site-packages (from bokeh->cost-of-capital-calculator) (11.1.0)
Requirement already satisfied: PyYAML>=3.10 in /home/runner/miniconda3/envs/ccc-dev/lib/python3.12/site-packages (from bokeh->cost-of-capital-calculator) (6.0.2)
Requirement already satisfied: tornado>=6.2 in /home/runner/miniconda3/envs/ccc-dev/lib/python3.12/site-packages (from bokeh->cost-of-capital-calculator) (6.4.2)
Requirement already satisfied: xyzservices>=2021.09.1 in /home/runner/miniconda3/envs/ccc-dev/lib/python3.12/site-packages (from bokeh->cost-of-capital-calculator) (2024.9.0)
Requirement already satisfied: python-dateutil>=2.8.2 in /home/runner/miniconda3/envs/ccc-dev/lib/python3.12/site-packages (from pandas->cost-of-capital-calculator) (2.9.0.post0)
Requirement already satisfied: pytz>=2020.1 in /home/runner/miniconda3/envs/ccc-dev/lib/python3.12/site-packages (from pandas->cost-of-capital-calculator) (2024.1)
Requirement already satisfied: tzdata>=2022.7 in /home/runner/miniconda3/envs/ccc-dev/lib/python3.12/site-packages (from pandas->cost-of-capital-calculator) (2024.2)
Requirement already satisfied: marshmallow>=3.0.0 in /home/runner/miniconda3/envs/ccc-dev/lib/python3.12/site-packages (from paramtools->cost-of-capital-calculator) (3.24.0)
Requirement already satisfied: fsspec in /home/runner/miniconda3/envs/ccc-dev/lib/python3.12/site-packages (from paramtools->cost-of-capital-calculator) (2024.12.0)
Requirement already satisfied: sortedcontainers in /home/runner/miniconda3/envs/ccc-dev/lib/python3.12/site-packages (from paramtools->cost-of-capital-calculator) (2.4.0)
Requirement already satisfied: numba in /home/runner/miniconda3/envs/ccc-dev/lib/python3.12/site-packages (from taxcalc->cost-of-capital-calculator) (0.60.0)
Requirement already satisfied: MarkupSafe>=2.0 in /home/runner/miniconda3/envs/ccc-dev/lib/python3.12/site-packages (from Jinja2>=2.9->bokeh->cost-of-capital-calculator) (3.0.2)
Requirement already satisfied: six>=1.5 in /home/runner/miniconda3/envs/ccc-dev/lib/python3.12/site-packages (from python-dateutil>=2.8.2->pandas->cost-of-capital-calculator) (1.17.0)
Requirement already satisfied: llvmlite<0.44,>=0.43.0dev0 in /home/runner/miniconda3/envs/ccc-dev/lib/python3.12/site-packages (from numba->taxcalc->cost-of-capital-calculator) (0.43.0)
# import packages
import pandas as pd
import numpy as np
from bokeh.plotting import figure, show
from bokeh.io import output_notebook
# import CCC classes that we'll work with
from ccc.data import Assets
from ccc.parameters import Specification, DepreciationParams
from ccc.calculator import Calculator
# to print bokeh plots inline
output_notebook()
Create an instance of the Assets class#
This is the class that contains the data that underlie CCC. The basic object is a Pandas DataFrame where each row represents a combination of a specific type of asset, industry (approximately 6-digit NAICS) and tax treatment. The columns represent names and codes for the asset and industry classifications, the tax depreciation rules used for that asset, and the rate of economic depreciation (delta).
assets = Assets()
assets.df.head(n=5)
| Unnamed: 0 | tax_treat | assets | bea_asset_code | bea_ind_code | Industry | minor_code_alt | major_asset_group | minor_asset_group | major_industry | asset_name | delta | |
|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 0 | 0 | corporate | 0.0 | RD32 | 110C | Farms | 111 | Intellectual Property | Intellectual Property | Agriculture, forestry, fishing, and hunting | Aerospace products and parts manufacturing | 0.22 |
| 1 | 1 | corporate | 0.0 | RD32 | 113F | Forestry, fishing, and related activities | 113 | Intellectual Property | Intellectual Property | Agriculture, forestry, fishing, and hunting | Aerospace products and parts manufacturing | 0.22 |
| 2 | 2 | corporate | 0.0 | RD32 | 113F | Forestry, fishing, and related activities | 114 | Intellectual Property | Intellectual Property | Agriculture, forestry, fishing, and hunting | Aerospace products and parts manufacturing | 0.22 |
| 3 | 3 | corporate | 0.0 | RD32 | 2110 | Oil and gas extraction | 211110 | Intellectual Property | Intellectual Property | Mining | Aerospace products and parts manufacturing | 0.22 |
| 4 | 4 | corporate | 0.0 | RD32 | 2120 | Mining, except oil and gas | 212110 | Intellectual Property | Intellectual Property | Mining | Aerospace products and parts manufacturing | 0.22 |
Create instances of the two parameters classes#
The Specification class contains many of the model parameters, although depreciation system parameters are contained in the DepreciationParams class. Both are required arguments for the Calculator object.
A Specification object has methods that load the data from a file that contains the default parameter values (default_parameters.json) and then stores them as attributes of the Specification class object.
p.u looks into the instance of this class named p and executing the cell below will show you that p.u is dictionary containtin marginal tax rates on corporate and pass-through income.
Note that the Specification class has some arguments, such as call_tc, which defaults to False, but if set to True will call the Tax-Calculator to estimate marginal tax rates on individual filers’ income.
Similar to Specification, DepreciationParams loads default parameters from a JSON file, in this case the file tax_depreciation_rules.json.
# Create an instance of the Specification class
# p = Specification()
p = Specification(call_tc=True)
# Look at attributes
p.u
Using CPS
Calculator initial year = 2014
Calculator year = 2024
year: 2024
{'tau_pt': array([0.30514124]), 'tau_div': array([0.17460724]), 'tau_int': array([0.31309077]), 'tau_scg': array([0.28526829]), 'tau_lcg': array([0.18438764]), 'tau_td': array([0.20617536]), 'tau_h': array([0.03500787])}
{'c': array([0.21]), 'pt': array([0.30514124])}
# Create an instance of the DepreciationParams class
dp = DepreciationParams()
# Look at an attribute of the dp object
dp.EP1A
[OrderedDict([('value', {'system': 'GDS', 'life': 5.0, 'method': 'DB 200%'}),
('year', 2013)])]
Create an instance of the Calculator class#
This class does the calculations on the data. It takes as arguments the data object (named assets here) and the parameters object (named p here).
Excuting the cell below creates an instance of the calculator class with these data and parameters, but does not yet excecute any calculations.
# Create an instance of the Calculator class
calc1 = Calculator(p, dp, assets)
With an instance of the Calculator class created, we can start doing some calculations with these data.
For instance, we can compute a table of the share of corporate vs non-corporate assets across each industry.
# Look at shares of assets across industry
calc1.asset_share_table()
| Industry | Corporate | Pass-Through | |
|---|---|---|---|
| 0 | Agriculture, forestry, fishing, and hunting | 0.309688 | 0.690312 |
| 1 | Mining | 0.826839 | 0.173161 |
| 2 | Utilities | 0.948540 | 0.051460 |
| 3 | Construction | 0.286840 | 0.713160 |
| 4 | Manufacturing | 0.859014 | 0.140986 |
| 5 | Wholesale trade | 0.715166 | 0.284834 |
| 6 | Retail trade | 0.728628 | 0.271372 |
| 7 | Transportation and warehousing | 0.789652 | 0.210348 |
| 8 | Information | 0.918722 | 0.081278 |
| 9 | Finance and insurance | 0.786056 | 0.213944 |
| 10 | Real estate and rental and leasing | 0.212052 | 0.787948 |
| 11 | Professional, scientific, and technical services | 0.501699 | 0.498301 |
| 12 | Management of companies and enterprises | 0.824397 | 0.175603 |
| 13 | Administrative and waste management services | 0.569110 | 0.430890 |
| 14 | Educational services | 0.590168 | 0.409832 |
| 15 | Health care and social assistance | 0.510425 | 0.489575 |
| 16 | Arts, entertainment, and recreation | 0.490811 | 0.509189 |
| 17 | Accommodation and food services | 0.482236 | 0.517764 |
| 18 | Other services, except government | 0.343467 | 0.656533 |
Specifying a reform policy#
To see some more interesting results, we will want to create another Calculator object with a change in policy or economic assumptions. We can do this in a way analogous to our original instantiation of the Calculator object above. In particular, we’ll need to create a new Specification object (we can use the same underlying data, which was in the object we named assets).
In the code below, we’ll specify our “reform” as current law tax policy for 2026 (the baseline parameters above came from the default model year, 2019). In addition, we’ll increase the corporate income tax rate from 21 to 25%.
# Create another policy
p2 = Specification(year=2026)
p2.update_specification({'CIT_rate': 0.35})
calc2 = Calculator(p2, dp, assets)
Tabular output#
Now with two Calculator objects named calc1 and calc2 (representing the baseline and reform policies), we are ready to compute some of the changes in effective tax rates, cost of capital, or other variables measured in this model.
We start with an overall summary table showing the marginal effective total tax rates (METTRs) for all investments, corporate investments, and pass-through investments under varying financing assumptions. This is done through the summary_table function. It takes a calculator object as an argument.
# Look at differences in METTRs between the two policies
calc1.summary_table(calc2) # calc1 is the baseline, calc2 the reform
| Marginal Effective Total Tax Rate Under Baseline Policy | Marginal Effective Total Tax Rate Under Reform Policy | Change from Baseline (pp) | ||
|---|---|---|---|---|
| 0 | Overall | 21.526876 | 36.559180 | 15.032303 |
| 1 | Corporations | 19.289825 | 42.290696 | 23.000871 |
| 2 | Equity Financed | 22.125764 | 46.019861 | 23.894097 |
| 3 | Debt Financed | 9.098922 | 26.810908 | 17.711986 |
| 4 | Pass-Through Entities | 25.445118 | 29.374410 | 3.929292 |
| 5 | Equity Financed | 27.142495 | 32.437430 | 5.294935 |
| 6 | Debt Financed | 17.726241 | 13.580154 | -4.146087 |
the Calculator.summary_table() method defaults to showing the results for the METTR, but there is a keyword argument that would allow you to view the output for other variables computed in CCC. We can use this to see changes in the cost of capital (denoted by \(\rho\) in the model):
# Look at how the cost of capital changed
calc1.summary_table(calc2, output_variable='rho')
| Cost of Capital Under Baseline Policy | Cost of Capital Under Reform Policy | Change from Baseline (pp) | ||
|---|---|---|---|---|
| 0 | Overall | 5.991266 | 7.410897 | 1.419631 |
| 1 | Corporations | 5.825206 | 8.146925 | 2.321719 |
| 2 | Equity Financed | 6.721989 | 9.697451 | 2.975462 |
| 3 | Debt Financed | 3.925764 | 4.875811 | 0.950047 |
| 4 | Pass-Through Entities | 6.189211 | 6.533551 | 0.344340 |
| 5 | Equity Financed | 7.184843 | 7.747925 | 0.563082 |
| 6 | Debt Financed | 3.762559 | 3.582046 | -0.180513 |
One can also save results to disk by specifying an output type (‘excel’, ‘json’, ‘csv’, ‘tex’) and a file path:
# Save these results to disk
calc1.summary_table(calc2, output_variable='rho', output_type='excel', path='cc_table.xlsx')
There are also Calculator methods to compute summary tables by asset type or industry. These are computed in the next two cells.
# Summary by asset type
calc1.asset_summary_table(calc2)
| Category | Marginal Effective Total Tax Rate Under Baseline Policy | Marginal Effective Total Tax Rate Under Reform Policy | Change from Baseline (pp) | |
|---|---|---|---|---|
| 0 | Overall | 21.526876 | 36.559180 | 15.032303 |
| 1 | Corporate | 19.289825 | 42.290696 | 23.000871 |
| 2 | Equipment | 12.268577 | 51.221226 | 38.952649 |
| 3 | Structures | 19.432117 | 37.101159 | 17.669042 |
| 4 | Intellectual Property | 9.379703 | 52.586404 | 43.206701 |
| 5 | Inventories | 28.801548 | 39.348000 | 10.546452 |
| 6 | Land | 26.056496 | 35.232345 | 9.175849 |
| 7 | Pass-through | 25.445118 | 29.374410 | 3.929292 |
| 8 | Equipment | 7.946231 | 18.018260 | 10.072029 |
| 9 | Structures | 23.314464 | 26.872234 | 3.557770 |
| 10 | Intellectual Property | 2.192332 | 5.866596 | 3.674263 |
| 11 | Inventories | 33.053519 | 36.968629 | 3.915110 |
| 12 | Land | 28.857306 | 32.236065 | 3.378758 |
# Summary by industry
calc1.industry_summary_table(calc2)
| Category | Marginal Effective Total Tax Rate Under Baseline Policy | Marginal Effective Total Tax Rate Under Reform Policy | Change from Baseline (pp) | |
|---|---|---|---|---|
| 0 | Overall | 21.526876 | 36.559180 | 15.032303 |
| 1 | Corporate | 19.289825 | 42.290696 | 23.000871 |
| 2 | Agriculture, forestry, fishing, and hunting | 20.534322 | 39.825970 | 19.291648 |
| 3 | Mining | 12.219031 | 40.334556 | 28.115526 |
| 4 | Utilities | 12.014254 | 35.617403 | 23.603149 |
| 5 | Construction | 20.688305 | 44.152531 | 23.464226 |
| 6 | Manufacturing | 18.982439 | 44.975737 | 25.993298 |
| 7 | Wholesale trade | 24.082118 | 44.175721 | 20.093603 |
| 8 | Retail trade | 24.739980 | 41.380972 | 16.640992 |
| 9 | Transportation and warehousing | 13.656523 | 39.123823 | 25.467300 |
| 10 | Information | 15.067806 | 45.626685 | 30.558879 |
| 11 | Finance and insurance | 20.429720 | 48.923895 | 28.494175 |
| 12 | Real estate and rental and leasing | 24.425292 | 37.289023 | 12.863731 |
| 13 | Professional, scientific, and technical ser... | 16.882275 | 52.644895 | 35.762620 |
| 14 | Management of companies and enterprises | 22.944526 | 45.750770 | 22.806244 |
| 15 | Administrative and waste management services | 17.130679 | 48.198467 | 31.067788 |
| 16 | Educational services | 22.748061 | 40.858827 | 18.110766 |
| 17 | Health care and social assistance | 21.572828 | 42.546117 | 20.973289 |
| 18 | Arts, entertainment, and recreation | 17.693016 | 40.552789 | 22.859773 |
| 19 | Accommodation and food services | 24.357793 | 40.406482 | 16.048689 |
| 20 | Other services, except government | 21.339215 | 39.923188 | 18.583973 |
| 21 | Pass-through | 25.445118 | 29.374410 | 3.929292 |
| 22 | Agriculture, forestry, fishing, and hunting | 22.890823 | 28.258745 | 5.367922 |
| 23 | Mining | 8.178225 | 13.282775 | 5.104550 |
| 24 | Utilities | 7.674243 | 15.386904 | 7.712661 |
| 25 | Construction | 22.783365 | 27.715992 | 4.932628 |
| 26 | Manufacturing | 19.656924 | 23.876224 | 4.219300 |
| 27 | Wholesale trade | 28.582124 | 32.982385 | 4.400262 |
| 28 | Retail trade | 28.669117 | 32.818548 | 4.149431 |
| 29 | Transportation and warehousing | 10.988363 | 16.814706 | 5.826343 |
| 30 | Information | 13.080848 | 22.380169 | 9.299321 |
| 31 | Finance and insurance | 25.182230 | 29.686517 | 4.504287 |
| 32 | Real estate and rental and leasing | 27.634351 | 31.037967 | 3.403616 |
| 33 | Professional, scientific, and technical ser... | 15.409318 | 21.332319 | 5.923000 |
| 34 | Management of companies and enterprises | 24.491387 | 28.707654 | 4.216267 |
| 35 | Administrative and waste management services | 15.727368 | 22.533481 | 6.806113 |
| 36 | Educational services | 24.105494 | 27.924088 | 3.818595 |
| 37 | Health care and social assistance | 22.439307 | 27.127672 | 4.688364 |
| 38 | Arts, entertainment, and recreation | 16.351020 | 22.227996 | 5.876976 |
| 39 | Accommodation and food services | 26.723881 | 30.833800 | 4.109919 |
| 40 | Other services, except government | 21.759689 | 25.712610 | 3.952921 |
Visualizations#
If one wants to visualize the effects of changes in tax policy, the Calculator class has a few methods for this.
We can use the grouped_bar method to show differential effects across assets (the default):
# Visualizing changes by asset type
aplot = calc1.grouped_bar(calc2)
show(aplot)
or we can change from the default asset category split to a split by industry by changing the value of the group_by_asset keyword argument:
# Visualizing changes by industry
iplot = calc1.grouped_bar(calc2, group_by_asset=False)
show(iplot)
There’s also a plot that illustrates the range of the effects of taxes on investments across asset types, by showing the min, max, and mean values:
# Plot to show variation in METTRs across assets
rplot = calc1.range_plot(calc2, output_variable='metr')
show(rplot)
changing to show the effects on pass-through businesses and for a different output variable:
show(rplot)
# Plot to show variation in METTRs across assets
output_notebook()
rplot2 = calc1.range_plot(calc2, corporate=False, output_variable='mettr')
show(rplot2)
There’s a bubble plot too, though it doesn’t yet allow for many options (currently only plots METTRs for corporate entities):
bplot = calc1.asset_bubble(calc2)
show(bplot)
WARNING:bokeh.core.validation.check:W-1000 (MISSING_RENDERERS): Plot has no renderers: figure(id='p1470', ...)
WARNING:bokeh.core.validation.check:W-1000 (MISSING_RENDERERS): Plot has no renderers: figure(id='p1409', ...)
Summary#
This notebook provides a brief example of how one would work with CCC. Please explore the source code to see additional flexiblity in the functions. And please leave any questions or suggestions in the CCC repo at PSLmodels/Cost-of-Capital-Calculator#issues.