Stock Portfolio Analyses Part 1
Python for Finance: Stock Portfolio Analyses
My two most recent blog posts were about Scaling Analytical Insights with Python; part 1 can be found here and part 2 can be found here. It has been several months since I wrote those, largely due to the fact that I relocated my family to Seattle to join Amazon in November; Iâve spent most of the time on my primary project determining our global rollout plan and related business intelligence roadmap.
Prior to my departure at my former company, FloSports, we were in the process of overhauling our analytics reporting across the organization (data, marketing, product et al), and part of this overhaul included our financial reporting. While I left early on in that implementation, over the past several months Iâve continued using Python extensively for financial analyses, particularly pandas
. In this post, I will share how I leveraged some very helpful online resources, the Yahoo Finance API
(requires a work around and may require a future data source replacement), and Jupyter notebook
to largely automate the tracking and benchmarking of a stock portfolioâs performance.
Overview of PME and benchmarking individual stock performance
As a quick background, I have been investing in my own stock portfolio since 2002 and developed a financial model for my portfolio a number of years ago. For years, I would download historical prices and load the data into the financial model â while online brokers calculate realized and unrealized returns, as well as income and dividends, I like to have historical data in the model as I conduct my own analyses to evaluate positions. One view / report which Iâve never found from online brokers and services is a âPublic Market Equivalentâ-like analysis. In short, the Public Market Equivalent (PME) is a set of analyses used in the private equity industry to compare the performance of a private equity fund relative to an industry benchmark. Much more detail here.
Related, the vast majority of equity portfolio managers are unable to select a portfolio of stocks which outperforms the broader market, e.g., S&P 500, over the long-term (~1 in 20 actively managed domestic funds beat index funds). Even when some individual stocks outperform, the underperformance of others often outweighs the better performing stocks, meaning overall an investor is worse off than simply investing in an index fund. During business school I learned about PME, and I incorporated a conceptually similar analysis into the evaluation of my current public equity holdings. To do this properly, you should measure the timing of investment inflows specific to each portfolio position (holding periods) relative to an S&P 500 equivalent dollar investment over the identical holding period. As an example, if you bought a stock on 6/1/2016 and you still own it, you would want to compare the stockâs return over that period to the return of an equal dollar investment on 6/1/2016 in the S&P 500 (our benchmark example). Among other things, you may find that even if a stock has done relatively well it may still trail the S&P 500âs return over the same time period.
In the past, I downloaded historical price data from Yahoo Finance and used INDEX and MATCH functions in excel to calculate the relative holding period performance of each position versus the S&P 500. While this is an OK way to accomplish this goal, conducting the same using pandas
in Jupyter notebook is more scalable and extensible. Whenever you download new data and load into excel, you inevitably need to modify some formulas and validate for errors. Using pandas
, adding new calculations, such as a cumulative ROI multiple (which Iâll cover), takes almost no time to implement. And the visualizations, for which I use Plotly
, are highly reproducible and much more useful in generating insights.
Disclosure: Nothing in this post should be considered investment advice. Past performance is not necessarily indicative of future returns. These are general examples about how to import data using pandas for a small sample of stocks across different time intervals and to benchmark their individual performance against an index. You should direct all investment related questions that you have to your financial advisor.
In addition to contributing this tutorial, Iâm continuing to revise and build upon this approach, and I outline some considerations for further development at the end of this post. I believe this post will be helpful for novice to intermediate-level data science oriented finance professionals, especially since this should extend to many other types of financial analyses. This approach is âPME-likeâ in the sense thatâs itâs measuring investment inflows over equal holding periods. As public market investments are much more liquid than private equity, and presuming you follow a trailing stop approach, from my perspective itâs more important to focus on active holdings â itâs generally advisable to divest holdings which underperform a benchmark or which you no longer want to own for various reasons, while I take a long-term view and am happy to own outperforming stocks for as long as theyâll have me.
Resources:
I am a current DataCamp subscriber (future post forthcoming on DataCamp) and this community tutorial on Python for Finance is great.
I have created a repo for this post including the Python notebook here, and the excel file here.
If you want to see the full interactive version (because Jupyter <â>> GitHub integration is awesome), you can view using nbviewer here.
Outline of what we want to accomplish:
Import S&P 500 and sample ticker data, using the Yahoo Finance API
Create a merged portfolio âmasterâ file which combines the sample portfolio dataframe with the historical ticker and historical S&P 500 data
Determine what the S&P 500 close was on the date of acquisition of each investment, which allows us to calculate the S&P 500 equivalent share position with the same dollars invested
Calculate the relative % and dollar value returns for the portfolio positions versus S&P 500 returns over that time
Calculate cumulative portfolio returns and ROI multiple, in order to assess how well this example portfolio compared to a market index
One of the more important items: dynamically calculate how each position is doing relative to a trailing stop, e.g., if a position closes 25% below its closing high, consider selling the position on the next trading day.
Visualizations
Total Return Comparisons â % return of each position relative to index benchmark
Cumulative Returns Over Time â $ Gain / (Loss) of each position relative to benchmark
Cumulative Investments Over Time â given the above, how do the overall investment returns compare to the equal weighting and time period of S&P 500 investments?
Adjusted Close % off of High Comparison â what is each positionâs most recent close relative to its adjusted closing high since purchased?
Investment Portfolio Python Notebook
Data Import and Dataframe Manipulation
You will begin by importing the necessary Python libraries, import the Plotly
offline module, and read in our sample portfolio dataframe.
Now that you have read in the sample portfolio file, youâll create a few variables which capture the date ranges for the S&P 500 and all of the portfolioâs tickers. Note that this is one of the few aspects of this notebook which requires an update each week (adjust the date range to include the most recent trading week â here, we are running this off of prices through 3/9/2018).
As mentioned in the Python Finance training post, the pandas-datareader
package enables us to read in data from sources like Google, Yahoo! Finance and the World Bank. Here Iâll focus on Yahoo! Finance, although Iâve worked very preliminarily with Quantopian and have also begun looking into quandl
as a data source. As also mentioned in the DataCamp post, the Yahoo API endpoint recently changed and this requires the installation of a temporary fix in order for Yahoo! Finance to work. Iâve made this needed slight adjustment in the code below. I have noticed some minor data issues where the data does not always read in as expected, or the last trading day is sometimes missing. While these issues have been relatively infrequent, Iâm continuing to monitor whether or not Yahoo! Finance will be the best and most reliable data source going forward.
If youâre following along with your own notebook, you should see something like the below once youâve successfully read in the data from Yahooâs API:
After loading in the S&P 500 data, youâll see that I inspect the head and tail of the dataframe, as well as condense the dataframe to only include the Adj Close
column. The difference between the Adjusted Close
and the Close
columns is that an adjusted close reflects dividends (see future areas for development below). When a company issues a dividend, the share price is reduced by the size of the dividend per share, as the company is distributing a portion of the companyâs earnings. For purposes of this analysis, you will only need to analyze this column. I also create a dataframe which only includes the S&Pâs adjusted close on the last day of 2017 (start of 2018); this is in order to run YTD comparisons of individual tickers relative to the S&P 500âs performance.
In the below code, you create an array of all of the tickers in our sample portfolio dataframe. You then write a function to read in all of the tickers and their relevant data into a new dataframe, which is essentially the same approach you took for the S&P500 but applied to all of the portfolioâs tickers.
As with the S&P 500 dataframe, youâll create an adj_close
dataframe which only has the Adj Close
column for all of your stock tickers. If you look at the notebook in the repo I link to above, this code is chunked out in more code blocks than shown below. For purposes of describing this here, Iâve included below all of the code which leads up to our initial merged_portfolio
dataframe.
Depending on your level of familiarity with pandas
, this will be very straightforward to slightly overwhelming. Below, Iâll unpack what these lines are doing:
The overall approach you are taking is an example of split-apply-combine (note this downloads a PDF).
The
all_data[['Adj Close']]
line creates a new dataframe with only the columns provided in the list; hereAdj Close
is the only item provided in the list.Using this line of code,
adj_close[adj_close['Date']==end_of_last_year]
, you are filtering theadj_close
dataframe to only the row where the dataâsDate
column equals the date which you earlier specified in theend_of_last_year
variable (2017, 12, 29).You also set the index of the
adj_close_latest
andportfolio_df
dataframes. I did this because this is how youâll merge the two dataframes. Themerge
function, very similar to SQL joins, is an extremely useful function which I use very often.Within the
merge
function, you specify the left dataframe (portfolio_df
) and our right dataframe (adj_close_latest
). By specifyingleft_index
andright_index
equal True, you are stating that the two dataframes share a common index and you will join both on this.Last, you create a new column called
'ticker return'
. This calculates the percent return for each stock position by dividing theAdj Close
by theUnit Cost
(initial purchase price for stock) and subtracting 1. This is similar to calculating a formula in excel and carrying it down, but inpandas
this is accomplished with one-line of code.
You have taken the individual dataframes for the S&P 500 and individual stocks, and you are beginning to develop a âmasterâ dataframe which weâll use for calculations, visualizations and any further analysis. Next, you continue to build on this âmasterâ dataframe with further use of pandas merge
function. Below, you reset the current dataframeâs index and begin joining your smaller dataframes with the master one. Once again, the below code block is broken out further in the Jupyter
notebook; here I take a similar approach to before where Iâll share the code below and then break down the key callouts below the code block.
You use
reset_index
on themerged_portfolio
in order to flatten the master dataframe and join on the smaller dataframesâ relevant columns.In the
merged_portfolio_sp
line, you merge the current master dataframe (merged_portfolio) with thesp_500_adj_close
; you do this in order to have the S&Pâs closing price on each positionâs purchase date â this allows you to track the S&P performance over the same time period that each position is held (from acquisition date to most recent market close date).The merge here is slightly different than before, in that we join on the left dataframeâs
Acquisition Date
column and on the right dataframeâsDate
column.After completing this merge, you will have extra columns which you do not need â since our master dataframe will eventually have a considerable number of columns for analysis, it is important to prune duplicative and unnecessary columns along the way.
There are several ways to remove unnecessary columns and perform various column name cleanups; for simplicity, I use
python
del
and then rename a few columns with pandasrename
method, clarifying the tickerâsAdj Close
column by renaming toTicker Adj Close
; and you distinguish the S&Pâs initial adjusted close withSP 500 Initial Close
.When you calculate
merged_portfolio_sp['Equiv SP Shares']
, you do so in order to be able to calculate the S&P 500âs equivalent value for the close on the date you acquired each ticker position: if you spend $5,000 on a new stock position, you could have spent $5,000 on the S&P 500; continuing the example, if the S&P 500 was trading at $2,500 per share at the time of purchase, you would have been able to purchase 2 shares. Later, if the S&P 500 is trading for $3,000 per share, your stake would be worth $6,000 (2 equivalent shares * $3,000 per share) and you would have $1,000 in paper profits over this comparable time period.In the rest of the code block, you next perform a similar merge, this time joining on the S&P 500âs latest close â this provides the second piece needed to calculate the S&Pâs comparable return relative to each positionâs holding period: the S&P 500 price on each tickerâs acquisition day and the S&P 500âs latest market close.
You have now further developed your âmasterâ dataframe with the following:
Each portfolio positionâs price, shares and value on the position acquisition day, as well as the latest marketâs closing price.
An equivalent S&P 500 price, shares and value on the equivalent position acquisition day for each ticker, as well as the latest S&P 500 closing price.
Given the above, you will next perform the requisite calculations in order to compare each positionâs performance, as well as the overall performance of this strategy / basket of stocks, relative to comparable dollar investment and holding times of the S&P 500.
Below is a summary of the new columns which you are adding to the âmasterâ dataframe.
In the first column,
['SP Return']
, you create a column which calculates the absolute percent return of the S&P 500 over the holding period of each position (note, this is an absolute return and is not an annualized return). In the second column (['Abs. Return Compare']
), you compare the['ticker return']
(each positionâs return) relative to the['SP Return']
over the same time period.In the next three columns,
['Ticker Share Value']
,['SP 500 Value']
and['Abs Value Compare']
, we calculate the dollar value (market value) equivalent based on the shares we hold multiplied by the latest adjusted close price (and subtract the S&P return from the ticker to calculate over / (under) performance).Last, the
['Stock Gain / (Loss)']
and['SP 500 Gain / (Loss)']
columns calculate our unrealized dollar gain / loss on each position and comparable S&P 500 gain / loss; this allows us to compare the value impact of each position versus simply investing those dollars in the S&P 500.You now have what you need in order to compare your portfolioâs performance to a portfolio equally invested in the S&P 500. The next two code block sections allow you to i) compare YTD performance of each position relative to the S&P 500 (a measure of momentum and how your positions are pacing) and ii) compare the most recent closing price for each portfolio position relative to its most recent closing high (this allows you to assess if a position has triggered a trailing stop, e.g., closed 25% below closing high).
Below, Iâll start with the YTD performance code block and provide details regarding the code further below.
When creating the
merged_portfolio_sp_latest_YTD
dataframe, you are now merging the âmasterâ dataframe with theadj_close_start
dataframe; as a quick reminder, you created this dataframe by filtering on theadj_close
dataframe where the'Date'
column equaled the variableend_of_last_year
; you do this because itâs how YTD (year-to-date) stock and index performances are measured; last yearâs ending close is the following yearâs starting price.From here, we once again use
del
to remove unnecessary columns and therename
method to clarify the âmasterâ dataframeâs newly added columns.Last, we take each Ticker (in the
['Ticker Adj Close']
column) and calculate the YTD return for each (we also have an S&P 500 equivalent value for each value in the'SP 500 Latest Close'
column).
In the below code block, you use the sort_values
method to re-sort our âmasterâ dataframe and then you calculate cumulative portfolio investments (sum of your position acquisition costs), as well the cumulative value of portfolio positions and the cumulative value of the theoretical S&P 500 investments. This allows you to be able to see how your total portfolio, with investments in positions made at different times across the entire period, compares overall to a strategy where you had simply invested in an index. Later on, youâll use the ['Cum Ticker ROI Mult']
to help you visualize how much each investment contributed to or decreased your overall return on investment (ROI).
You are now nearing the home stretch and almost ready to start visualizing your data and assessing the strengths and weaknesses of your portfolioâs individual ticker and overall strategy performance.
As before, Iâve included the main code block for determining where positions are trading relative to their recent closing high; Iâll then unpack the code further below.
To begin, you merge the
adj_close
dataframe with theportfolio_df
dataframe; this is the third time that youâve leveraged thisadj_close
dataframe in order to conduct an isolated analysis which youâll then combine with the overall âmasterâ dataframe.This initial merge is not particularly useful, as you have dates and adjusted close prices which pre-date your acquisition date for each position; as a result, weâll subset the data post our acquisition date, and then find the
max
closing price since that time.Once again, I used
del
to delete the merged dataframeâs unneeded columns; this is code I should refactor, as creating a list, e.g.,cols_to_keep
, and then filtering the dataframe with this would be a better approach â as an FYI, running thedel
code block more than once will throw an error and you would need to re-initialize your dataframe then run thedel
code block again.After removing the unnecessary columns, you then use the
sort_values
method and sort the values by the'Ticker'
,'Acquisition Date'
, and'Date'
columns (all ascending); you do this to make sure all of the ticker rows are sorted together, and we sort by Acquisition Date (in case weâve purchased the same stock more than once) and Date ascending in order to filter out the dates prior to your positionsâ acquisition dates. In other words, you are only concerned with the closing high since youâve held the position.In order to filter our dataframe, you create a new column
['Date Delta']
which is calculated by the difference between the Date and Acquisition Date columns.You then convert this column into a numeric column, and you create a new dataframe called
adj_close_acq_date_modified
where the['Date Delta']
is >= 0. This ensures that you are only evaluating closing highs since the date that you purchased each position.Now that you have the
adj_close_acq_date_modified
dataframe, weâll use a very powerful pandas function calledpivot_table
. If youâre familiar with pivot tables in Excel, this function is similar in that you can pivot data based on a single or multi-index, specify values to calculate and columns to pivot on, and also useagg functions
(which leveragenumpy
).Using the
pivot_table
function, we pivot on Ticker and Acquisition Date and specify that we would like to find the maximum (np.max
)Adj Close
for each position; this allows you to compare the recent Adjusted Close for each position relative to this High Adjusted Close.Now you have an
adj_close_pivot
dataframe, and you reset the index and join this once again on theadj_close
dataframe. This creates theadj_close_pivot_merged
dataframe, which tells you when you purchased each position and the date on which it hit its closing high since acquisition.Finally, we will combine our âmasterâ dataframe with this last smaller dataframe,
adj_close_pivot_merged
.After doing so, you are now able to calculate the final column needed,
['Pct off High']
. You take the['Ticker Adj Close']
and divide it by the['Closing High Adj Close']
and subtract 1. Note, that this percentage will always be negative, unless the stock happened to have its highest close (in this case it will be zero) on the most recent trading day evaluated (this is generally a very good sign if itâs the case).
This has been a pretty significant lift, and itâs now time for our long-awaited visualizations. If youâve continued to follow along in your own notebook, you now have a very rich dataframe with a number of calculated portfolio metrics, as shown in the below:
Total Return and Cumulative Return Visualizations
For all of these visualizations youâll use Plotly
, which allows you to make D3 charts entirely without code. While I also use Matplotlib
and Seaborn
, I really value the interactivity of Plotly
; and once you are used to it, the syntax becomes fairly straightforward and dynamic charts are easily attainable.
Your first chart below compares each individual positionâs total return relative to the S&P 500 (same holding periods for the position and hypothetical investment in the S&P 500). In the below, youâll see that over their distinct holding periods, 6 of the 8 positions outperformed the S&P. The last two, Twitter (which actually has had a negative return) and Walmart underperformed an equal timed investment in the S&P 500.
As each of these visualizations are relatively similar, Iâll explain the code required to generate the above Plotly visualization, and for the remaining ones Iâll only summarize observations from each visualization.
When using
Plotly
, you createtraces
which will plot the x and y data you specify. Here, you specify in trace1 that you want to plot a bar chart, with each Ticker on the x-axis and each tickerâs return on the y-axis.In trace2, to break up the data a bit, weâll use a Scatter line chart for the Ticker on the x-axis and the S&P Return on the y-axis.
Where the bar is above the line, the individual ticker (6 of 8 times) has outperformed the S&P 500.
You then create a data object with these traces, and then you provide a layout for the chart; in this case you specify a title, barmode, and the position of the legend; you also pass in a title and tick format (percent format to two decimal places) for the y-axis series.
You then create a figure object using
go.Figure
, specifying the data and layout objects, which you previously nameddata
andlayout
.
The next chart below shows the gain / (loss) dollar amount for each position, relative to the S&P 500, as well as shows the Ticker Total Return %. While it is generally recommended that you allocate an equal position size to your positions (or potentially determine positition sizing based on implied volatility), this may not always be the case. For a less volatile investment, you may invest more than in a riskier position (or you may have other position sizing rules). Given this, this visualization shows both each positionâs return and the dollar value contribution to your overall portfolioâs return.
Here, you can see that although you invested slightly less in Facebook (FB) than other positions, this stock has returned an ~$20k in this mock portfolio, greater than a 4x return relative to an equivalent S&P 500 investment over the same holding period.
The next chart below leverages the cumulative columns which you created: 'Cum Invst'
, 'Cum SP Returns'
, 'Cum Ticker Returns'
, and 'Cum Ticker ROI Mult'
.
Across the x-axis you have sorted the portfolio alphabetically. Each position shows the initial investment and total value (investment plus returns or less losses) for that position, combined with the positions preceding it.
To explain further, based on the ~$8k investment in AAPL, this grew to ~$22.5k (>$14k in gains), versus $15k in total value for the S&P. This is a 2.75x return over the initial investment in AAPL ($22.5k value from $8k investment is ~2.75x ROI).
Continuing to FB, you have invested ~$16k in aggregate ($8k in both positions), and this has grown to over $50k, a greater than 3x total return â this means that FB expanded your overall portfolio ROI.
Further down the x-axis, you see that both TWTR and WMT have reduced the overall portfolio ROI â this is obvious, as both have underperformed the S&P, but I believe that the magnitude of the contribution is clearer with this visualization.
As a caveat, this cumulative approach, given the different holding periods, is a bit of an apples and oranges combination for some positions based on when they were acquired. However, you can always isolate this analysis by sub-setting into smaller dataframes and separately compare positions which have more consistent holding periods. For example, you could compare your 2H 2016 and 1H 2017 purchases separate of one another.
Adjusted Close % off of High Comparison
Your final chart compares how far off each positionâs latest close price is from its adjusted closing high since the position was purchased. This is generally an important visualization to consider:
When a stock closes at higher prices, itâs generally recommended to adjust your trailing stop up as well. To illustrate, hereâs an example:
A position is acquired at $10 and doubles to $20 â using a 25% trailing stop, you would want to consider selling this position the next day if it closed at $15 ($15 / $20â1 = (25%)).
If the position increased to $25, you would want to consider moving your trailing stop up to $18.75 ($18.75 / $25â1 = (25%)).
As mentioned early on, nothing in here is intended to be financial advice; different trading systems have different rules for trailing stops, and this is an illustrative example.
Trailing stops are meant to help preserve gains and are generally important in mitigating the emotions of investing; while itâs easy to see your positionâs current return, what tends to be manual (or somewhat expensive if you use a trailing stop service) is calculating how close your positions are to your trailing stops.
This final visualization makes this easy to evaluate for any date you are reviewing; in the chart, we see that AAPL, MTCH, and NFLX all closed on 3/9/2018 at their closing highs (typically a very good sign).
However, TWTR is greater than 25% below its highest close (33% below as of 3/9/2018) and WMT is ~20% below its highest close.
In this case, you might want to sell TWTR and continue to keep a close eye on the performance of WMT.
Limitations to Approach and Closing Summary
Now you have a relatively extensible Jupyter notebook and portfolio dataset, which you are able to use to evaluate your stock portfolio, as well as add in new metrics and visualizations as you see fit.
Please note that while this notebook provides a fairly thorough review of a portfolio, the below have not yet been taken into consideration, would have an impact on the overall comparison, and likely present great areas for future development:
As noted initially, this notebook focuses on active holdings â ideally, we would evaluate all positions, both exited and active, in order to have a truly holistic view on oneâs investment strategy relative to alternatives, such as an index comparison.
The approach in here does not factor in dividends; while we evaluate adjusted close prices (which reflect dividends), total shareholder return combines share price appreciation and dividends to show a stockâs total return; while this is more difficult to do, it is something Iâm evaluating to include in the future.
On a related note, investors can also reinvest dividends in a position, rather than take a cash distribution; this is arguably even more complicated than accounting for dividends, as the acquisition costs are low and spread out, and over several years of holding a position you could have four (or more) acquisition dates each year for stocks where you reinvest dividends.
With those future areas in mind, we accomplished a lot here; this includes importing S&P 500 and ticker data using Yahoo! Financeâs API and creating a master dataframe which combines your portfolio with historical ticker and comparative S&P 500 prices. In doing this, you are able to calculate the absolute percent and dollar value returns for each position (and as compared to equally timed S&P 500 investments), as well as the cumulative impact of each position on your overall portfolioâs performance. You can also dynamically monitor your trailing stops, based on your own trading rules. And you have created visualizations which allow you to have much better insight into your master dataframe, focusing on the different metrics and each positionâs contribution to each.
I hope that you found this tutorial useful, and I welcome any feedback in the comments. Feel free to also reach out to me on twitter, @kevinboller, and my personal blog can be found here.
Last updated