# Load packages
# Core
library(tidyverse)
library(tidyquant)
Collect individual returns into a portfolio by assigning a weight to each stock
Choose your stocks.
from 2012-12-31 to 2017-12-31
symbols <- c ("GS", "JPM", "MS")
prices <- tq_get(x = symbols,
get = "stockprices",
from = "2012-12-31",
to = "2017-12-31")
asset_returns_tbl <- prices %>%
group_by(symbol) %>%
tq_transmute(select = adjusted,
mutate_fun = periodReturn,
period = "quarterly",
type = "log") %>%
slice(-1) %>%
ungroup() %>%
set_names(c("asset", "date", "returns"))
# symbols
symbols <- asset_returns_tbl %>% distinct(asset) %>% pull()
symbols
## [1] "GS" "JPM" "MS"
# weights
weights <- c(0.33, 0.33, 0.33)
weights
## [1] 0.33 0.33 0.33
w_tbl <- tibble(symbols, weights)
w_tbl
## # A tibble: 3 × 2
## symbols weights
## <chr> <dbl>
## 1 GS 0.33
## 2 JPM 0.33
## 3 MS 0.33
# ?tq_portfolio
portfolio_returns_tbl <- asset_returns_tbl %>%
tq_portfolio(assets_col = asset,
returns_col = returns,
weights = w_tbl,
rebalance_on = "quarters")
portfolio_returns_tbl
## # A tibble: 20 × 2
## date portfolio.returns
## <date> <dbl>
## 1 2013-03-28 0.122
## 2 2013-06-28 0.0829
## 3 2013-09-30 0.0444
## 4 2013-12-31 0.132
## 5 2014-03-31 -0.0118
## 6 2014-06-30 0.00625
## 7 2014-09-30 0.0715
## 8 2014-12-31 0.0728
## 9 2015-03-31 -0.0443
## 10 2015-06-30 0.104
## 11 2015-09-30 -0.160
## 12 2015-12-31 0.0466
## 13 2016-03-31 -0.155
## 14 2016-06-30 0.0159
## 15 2016-09-30 0.125
## 16 2016-12-30 0.312
## 17 2017-03-31 0.000947
## 18 2017-06-30 0.0192
## 19 2017-09-29 0.0669
## 20 2017-12-29 0.0937
portfolio_returns_tbl %>%
ggplot(mapping = aes(x = date, y = portfolio.returns)) +
geom_point(color = "cornflowerblue") +
# Formatting
scale_x_date(date_breaks = "1 year",
date_labels = "%Y") +
# Labeling
labs(y = "monthly returns",
x = NULL,
title = "Portfolio Returns Scatter")
portfolio_returns_tbl %>%
ggplot(mapping = aes(x = portfolio.returns)) +
geom_histogram(fill = "cornflowerblue", binwidth = 0.05) +
labs(x = "returns",
title = "Portfolio Returns Distribution")
portfolio_returns_tbl %>%
ggplot(mapping = aes(x = portfolio.returns)) +
geom_histogram(fill = "cornflowerblue", binwidth = 0.01) +
geom_density() +
# Formatting
scale_x_continuous(labels = scales::percent_format()) +
labs(x = "returns",
y = "distribution",
title = "Portfolio Histogram & Density")