166
|
1 |
// Advanvced Part 3 about a really dumb investment strategy
|
130
|
2 |
//==========================================================
|
|
3 |
|
|
4 |
object CW6c {
|
|
5 |
|
|
6 |
|
|
7 |
//two test portfolios
|
|
8 |
|
160
|
9 |
val blchip_portfolio = List("GOOG", "AAPL", "MSFT", "IBM", "FB", "AMZN", "BIDU")
|
166
|
10 |
val rstate_portfolio = List("PLD", "PSA", "AMT", "AIV", "AVB", "BXP", "CCI",
|
|
11 |
"DLR", "EQIX", "EQR", "ESS", "EXR", "FRT", "GGP", "HCP")
|
130
|
12 |
|
166
|
13 |
// (1) The function below should obtain the first trading price
|
|
14 |
// for a stock symbol by using the query
|
|
15 |
//
|
|
16 |
// http://ichart.yahoo.com/table.csv?s=<<symbol>>&a=0&b=1&c=<<year>>&d=1&e=1&f=<<year>>
|
|
17 |
//
|
|
18 |
// and extracting the first January Adjusted Close price in a year.
|
|
19 |
|
130
|
20 |
|
|
21 |
import io.Source
|
|
22 |
import scala.util._
|
|
23 |
|
166
|
24 |
def get_january_data(symbol: String, year: Int) : List[String] =
|
|
25 |
Source.fromFile(symbol ++ ".csv")("ISO-8859-1").getLines.toList.filter(_.startsWith(year.toString))
|
130
|
26 |
|
|
27 |
|
|
28 |
def get_first_price(symbol: String, year: Int) : Option[Double] = {
|
166
|
29 |
val data = Try(Some(get_january_data(symbol, year).head)) getOrElse None
|
|
30 |
data.map(_.split(",").toList(1).toDouble)
|
130
|
31 |
}
|
|
32 |
|
166
|
33 |
get_first_price("GOOG", 1980)
|
|
34 |
get_first_price("GOOG", 2010)
|
|
35 |
get_first_price("FB", 2014)
|
130
|
36 |
|
166
|
37 |
|
|
38 |
// Complete the function below that obtains all first prices
|
|
39 |
// for the stock symbols from a portfolio for the given
|
|
40 |
// range of years
|
|
41 |
|
|
42 |
def get_prices(portfolio: List[String], years: Range): List[List[Option[Double]]] =
|
|
43 |
for (year <- years.toList) yield
|
|
44 |
for (symbol <- portfolio) yield get_first_price(symbol, year)
|
130
|
45 |
|
|
46 |
|
166
|
47 |
// test case
|
|
48 |
val prices_fb = get_prices(List("FB"), 2012 to 2014)
|
|
49 |
val prices_goap = get_prices(List("GOOG", "AAPL"), 2010 to 2014)
|
|
50 |
val prices_bidu = get_prices(List("BIDU"), 2004 to 2008)
|
|
51 |
val prices_goapfb = get_prices(List("GOOG", "AAPL", "FB"), 2010 to 2016)
|
|
52 |
|
|
53 |
// (2) The first function below calculates the change factor (delta) between
|
|
54 |
// a price in year n and a price in year n+1. The second function calculates
|
|
55 |
// all change factors for all prices (from a portfolio).
|
|
56 |
|
|
57 |
def get_delta(price_old: Option[Double], price_new: Option[Double]) : Option[Double] = {
|
|
58 |
(price_old, price_new) match {
|
|
59 |
case (Some(x), Some(y)) => Some((y - x) / x)
|
|
60 |
case _ => None
|
|
61 |
}
|
161
|
62 |
}
|
|
63 |
|
166
|
64 |
def get_deltas(data: List[List[Option[Double]]]): List[List[Option[Double]]] =
|
|
65 |
for (i <- (0 until (data.length - 1)).toList) yield
|
|
66 |
for (j <- (0 until (data(0).length)).toList) yield get_delta(data(i)(j), data(i + 1)(j))
|
130
|
67 |
|
|
68 |
|
166
|
69 |
// test case using the prices calculated above
|
|
70 |
val deltas_fb = get_deltas(prices_fb)
|
|
71 |
val deltas_goap = get_deltas(prices_goap)
|
|
72 |
val deltas_bidu = get_deltas(prices_bidu)
|
|
73 |
val deltas_goapfb = get_deltas(prices_goapfb)
|
130
|
74 |
|
|
75 |
|
|
76 |
// (3) Write a function that given change factors, a starting balance and a year
|
166
|
77 |
// calculates the yearly yield, i.e. new balanace, according to our dump investment
|
|
78 |
// strategy. Another function calculates given the same data calculates the
|
|
79 |
// compound yield up to a given year. Finally a function combines all
|
|
80 |
// calculations by taking a portfolio, a range of years and a start balance
|
|
81 |
// as arguments.
|
130
|
82 |
|
|
83 |
|
166
|
84 |
def yearly_yield(data: List[List[Option[Double]]], balance: Long, year: Int): Long = {
|
|
85 |
val somes = data(year).flatten
|
|
86 |
val somes_length = somes.length
|
|
87 |
if (somes_length == 0) balance
|
|
88 |
else {
|
|
89 |
val portion: Double = balance.toDouble / somes_length.toDouble
|
|
90 |
(for (x <- somes) yield ((1.0 + x) * portion)).sum.toLong
|
|
91 |
//balance + (for (x <- somes) yield (x * portion)).sum.toLong
|
|
92 |
}
|
|
93 |
}
|
|
94 |
|
|
95 |
def yearly_yield_double(data: List[List[Option[Double]]], balance: Double, year: Int): Double = {
|
|
96 |
val somes = data(year).flatten
|
|
97 |
val somes_length = somes.length
|
|
98 |
if (somes_length == 0) balance
|
|
99 |
else {
|
|
100 |
val portion: Double = balance / somes_length.toDouble
|
|
101 |
balance + (for (x <- somes) yield (x * portion)).sum
|
|
102 |
}
|
|
103 |
}
|
|
104 |
|
|
105 |
def compound_yield(data: List[List[Option[Double]]], balance: Long, year: Int): Long = {
|
|
106 |
println(balance)
|
|
107 |
if (year >= data.length) balance else {
|
|
108 |
val new_balance = yearly_yield(data, balance, year)
|
|
109 |
compound_yield(data, new_balance, year + 1)
|
|
110 |
}
|
|
111 |
}
|
|
112 |
|
|
113 |
def compound_yield_double(data: List[List[Option[Double]]], balance: Double, year: Int): Long = {
|
|
114 |
if (year >= data.length) balance.toLong else {
|
|
115 |
val new_balance = yearly_yield_double(data, balance.toDouble, year)
|
|
116 |
compound_yield_double(data, new_balance, year + 1)
|
|
117 |
}
|
130
|
118 |
}
|
|
119 |
|
166
|
120 |
val prices_fb = get_prices(List("FB"), 2012 to 2014)
|
|
121 |
val prices_goap = get_prices(List("GOOG", "AAPL"), 2010 to 2014)
|
|
122 |
val prices_bidu = get_prices(List("BIDU"), 2004 to 2008)
|
|
123 |
val prices_goapfb = get_prices(List("GOOG", "AAPL", "FB"), 2010 to 2016)
|
|
124 |
|
|
125 |
val deltas_fb = get_deltas(prices_fb)
|
|
126 |
val deltas_goap = get_deltas(prices_goap)
|
|
127 |
val deltas_bidu = get_deltas(prices_bidu)
|
|
128 |
val deltas_goapfb = get_deltas(prices_goapfb)
|
|
129 |
|
|
130 |
|
|
131 |
yearly_yield(deltas_bidu, 100, 0)
|
|
132 |
yearly_yield(deltas_bidu, 100, 1)
|
|
133 |
yearly_yield(deltas_bidu, 100, 2)
|
|
134 |
yearly_yield(deltas_bidu, 192, 3)
|
|
135 |
//598
|
|
136 |
|
|
137 |
yearly_yield(deltas_fb, 100, 0)
|
|
138 |
yearly_yield(deltas_fb, 100, 1)
|
|
139 |
//yearly_yield(deltas_fb, 195, 2)
|
|
140 |
|
|
141 |
yearly_yield(deltas_goap, 100, 0)
|
|
142 |
yearly_yield(deltas_goap, 125, 1)
|
|
143 |
yearly_yield(deltas_goap, 146, 2)
|
|
144 |
yearly_yield(deltas_goap, 177, 3)
|
|
145 |
//yearly_yield(deltas_goap, 227, 4)
|
|
146 |
|
|
147 |
yearly_yield(deltas_goapfb, 100, 0)
|
|
148 |
yearly_yield(deltas_goapfb, 125, 1)
|
|
149 |
yearly_yield(deltas_goapfb, 146, 2)
|
|
150 |
yearly_yield(deltas_goapfb, 177, 3)
|
|
151 |
yearly_yield(deltas_goapfb, 267, 4)
|
|
152 |
yearly_yield(deltas_goapfb, 337, 5)
|
|
153 |
//416
|
|
154 |
|
|
155 |
compound_yield(deltas_goapfb.take(1), 100, 0)
|
|
156 |
compound_yield(deltas_goapfb.take(2), 100, 0)
|
|
157 |
compound_yield(deltas_goapfb.take(3), 100, 0)
|
|
158 |
compound_yield(deltas_goapfb.take(4), 100, 0)
|
|
159 |
compound_yield(deltas_goapfb.take(5), 100, 0)
|
|
160 |
compound_yield(deltas_goapfb.take(6), 100, 0)
|
|
161 |
|
|
162 |
compound_yield_double(deltas_goapfb.take(1), 100.0, 0)
|
|
163 |
compound_yield_double(deltas_goapfb.take(2), 100.0, 0)
|
|
164 |
compound_yield_double(deltas_goapfb.take(3), 100.0, 0)
|
|
165 |
compound_yield_double(deltas_goapfb.take(4), 100.0, 0)
|
|
166 |
compound_yield_double(deltas_goapfb.take(5), 100.0, 0)
|
|
167 |
compound_yield_double(deltas_goapfb.take(6), 100.0, 0)
|
|
168 |
|
|
169 |
//yearly_yield(d, 100, 0)
|
|
170 |
//compound_yield(d.take(6), 100, 0)
|
|
171 |
|
|
172 |
//test case
|
|
173 |
//yearly_yield(d, 100, 0)
|
|
174 |
//yearly_yield(d, 225, 1)
|
|
175 |
//yearly_yield(d, 246, 2)
|
|
176 |
//yearly_yield(d, 466, 3)
|
|
177 |
//yearly_yield(d, 218, 4)
|
|
178 |
|
|
179 |
//yearly_yield(d, 100, 0)
|
|
180 |
//yearly_yield(d, 125, 1)
|
|
181 |
|
|
182 |
def investment(portfolio: List[String], years: Range, start_balance: Long): Long = {
|
|
183 |
compound_yield(get_deltas(get_prices(portfolio, years)), start_balance, 0)
|
130
|
184 |
}
|
|
185 |
|
166
|
186 |
def investment_double(portfolio: List[String], years: Range, start_balance: Long): Long = {
|
|
187 |
compound_yield_double(get_deltas(get_prices(portfolio, years)), start_balance.toDouble, 0)
|
|
188 |
}
|
130
|
189 |
|
|
190 |
|
161
|
191 |
investment(rstate_portfolio, 1978 to 2017, 100)
|
|
192 |
investment(blchip_portfolio, 1978 to 2017, 100)
|
130
|
193 |
|
166
|
194 |
investment_double(rstate_portfolio, 1978 to 2017, 100)
|
|
195 |
investment_double(blchip_portfolio, 1978 to 2017, 100)
|
|
196 |
|
|
197 |
|
|
198 |
/*
|
|
199 |
val q1 = get_deltas(get_prices(List("GOOG", "AAPL", "BIDU"), 2000 to 2017))
|
|
200 |
yearly_yield(q1, 100, 0)
|
|
201 |
yearly_yield(q1, 100, 1)
|
|
202 |
yearly_yield(q1, 100, 2)
|
|
203 |
yearly_yield(q1, 100, 3)
|
|
204 |
yearly_yield(q1, 100, 4)
|
|
205 |
yearly_yield(q1, 100, 5)
|
|
206 |
yearly_yield(q1, 100, 6)
|
|
207 |
|
|
208 |
investment(List("GOOG", "AAPL", "BIDU"), 2004 to 2017, 100)
|
|
209 |
val one = get_deltas(get_prices(rstate_portfolio, 1978 to 1984))
|
|
210 |
val two = get_deltas(get_prices(blchip_portfolio, 1978 to 1984))
|
|
211 |
|
|
212 |
val one_full = get_deltas(get_prices(rstate_portfolio, 1978 to 2017))
|
|
213 |
val two_full = get_deltas(get_prices(blchip_portfolio, 1978 to 2017))
|
|
214 |
|
|
215 |
one_full.map(_.flatten).map(_.sum).sum
|
|
216 |
two_full.map(_.flatten).map(_.sum).sum
|
|
217 |
|
|
218 |
//test cases for the two portfolios given above
|
|
219 |
|
|
220 |
//println("Real data: " + investment(rstate_portfolio, 1978 to 2017, 100))
|
|
221 |
//println("Blue data: " + investment(blchip_portfolio, 1978 to 2017, 100))
|
|
222 |
|
|
223 |
for (i <- 2000 to 2017) {
|
|
224 |
println("Year " + i)
|
|
225 |
//println("Real data: " + investment(rstate_portfolio, 1978 to i, 100))
|
|
226 |
//println("Blue data: " + investment(blchip_portfolio, 1978 to i, 100))
|
|
227 |
println("test: " + investment(List("GOOG", "AAPL", "BIDU"), 2000 to i, 100))
|
144
|
228 |
}
|
166
|
229 |
|
|
230 |
|
|
231 |
*/
|
|
232 |
//1984
|
|
233 |
//1992
|
|
234 |
}
|