11 page 363
sim.time.series = function(start = 100, n = 10000)
{
time.series = c()
for (i in 1:n)
{
time.series[i] = start + sum(rnorm(364, 0, 0.25))
}
return(time.series)
}
sim.series = sim.time.series()
p.gt.100 = length(sim.series[sim.series >=100])/ length(sim.series)
p.gt.110 = length(sim.series[sim.series >=110])/ length(sim.series)
p.gt.120 = length(sim.series[sim.series >=120])/ length(sim.series)
print(p.gt.100)
## [1] 0.4918
## [1] 0.0165
## [1] 0
knitr::include_graphics("https://github.com/mlforsachid/MSDSQ2Data605/blob/master/Week9/HW-9/HW-9.png?raw=true")
