Get seasonal forcing

Convert Fourier function parameters into a smooth vector of daily seasonal forcing. Used within set_seasonality function

get_seasonal_forcing(season_params, n_days, floor = 0.001)

Arguments

season_params

List of Fourier function parameters describing seasonality

n_days

Numerical value of output vector length

floor

Minimum permitted value of output

Examples

season_params <- list(g0 = 0.28,
                          g = c(-0.30, -0.03, 0.17),
                          h = c(-0.35, 0.32, -0.08),
                          theta_c = 0.28)
get_seasonal_forcing(season_params,
                      n_days = 1000,
                      floor = 0.005)
#>    [1] 0.431058045 0.432371381 0.432516563 0.431502665 0.429342627 0.426053169
#>    [7] 0.421654688 0.416171150 0.409629971 0.402061886 0.393500817 0.383983725
#>   [13] 0.373550456 0.362243583 0.350108237 0.337191938 0.323544409 0.309217399
#>   [19] 0.294264488 0.278740900 0.262703301 0.246209605 0.229318767 0.212090587
#>   [25] 0.194585498 0.176864367 0.158988287 0.141018375 0.123015571 0.105040433
#>   [31] 0.087152943 0.069412309 0.051876777 0.034603441 0.017648063 0.001064895
#>   [37] 0.005000000 0.005000000 0.005000000 0.005000000 0.005000000 0.005000000
#>   [43] 0.005000000 0.005000000 0.005000000 0.005000000 0.005000000 0.005000000
#>   [49] 0.005000000 0.005000000 0.005000000 0.005000000 0.005000000 0.005000000
#>   [55] 0.005000000 0.005000000 0.005000000 0.005000000 0.005000000 0.005000000
#>   [61] 0.005000000 0.005000000 0.005000000 0.005000000 0.005000000 0.005000000
#>   [67] 0.005000000 0.005000000 0.005000000 0.005000000 0.005000000 0.005000000
#>   [73] 0.005000000 0.005000000 0.005000000 0.005000000 0.005000000 0.005000000
#>   [79] 0.005000000 0.003651030 0.016892749 0.030127458 0.043301221 0.056360448
#>   [85] 0.069252084 0.081923809 0.094324230 0.106403075 0.118111380 0.129401677
#>   [91] 0.140228178 0.150546948 0.160316085 0.169495882 0.178048993 0.185940583
#>   [97] 0.193138486 0.199613335 0.205338702 0.210291223 0.214450709 0.217800259
#>  [103] 0.220326354 0.222018943 0.222871526 0.222881218 0.222048803 0.220378785
#>  [109] 0.217879420 0.214562736 0.210444551 0.205544472 0.199885879 0.193495908
#>  [115] 0.186405416 0.178648933 0.170264606 0.161294133 0.151782677 0.141778783
#>  [121] 0.131334270 0.120504121 0.109346360 0.097921922 0.086294501 0.074530411
#>  [127] 0.062698411 0.050869545 0.039116960 0.027515716 0.016142601 0.005075920
#>  [133] 0.005000000 0.005000000 0.005000000 0.005000000 0.005000000 0.005000000
#>  [139] 0.005000000 0.005000000 0.005000000 0.005000000 0.005000000 0.005000000
#>  [145] 0.005000000 0.005000000 0.005000000 0.005000000 0.005000000 0.005000000
#>  [151] 0.005000000 0.005000000 0.005000000 0.002305742 0.020056647 0.039863399
#>  [157] 0.061768847 0.085811660 0.112026161 0.140442181 0.171084914 0.203974793
#>  [163] 0.239127363 0.276553182 0.316257716 0.358241264 0.402498880 0.449020315
#>  [169] 0.497789975 0.548786882 0.601984659 0.657351518 0.714850272 0.774438352
#>  [175] 0.836067842 0.899685523 0.965232938 1.032646463 1.101857394 1.172792048
#>  [181] 1.245371876 1.319513586 1.395129287 1.472126633 1.550408987 1.629875597
#>  [187] 1.710421777 1.791939107 1.874315632 1.957436083 2.041182100 2.125432464
#>  [193] 2.210063342 2.294948533 2.379959727 2.464966768 2.549837921 2.634440148
#>  [199] 2.718639389 2.802300839 2.885289241 2.967469170 3.048705325 3.128862822
#>  [205] 3.207807485 3.285406138 3.361526898 3.436039463 3.508815399 3.579728424
#>  [211] 3.648654692 3.715473064 3.780065381 3.842316730 3.902115702 3.959354644
#>  [217] 4.013929902 4.065742056 4.114696151 4.160701909 4.203673939 4.243531934
#>  [223] 4.280200856 4.313611110 4.343698709 4.370405421 4.393678908 4.413472850
#>  [229] 4.429747058 4.442467567 4.451606726 4.457143261 4.459062334 4.457355583
#>  [235] 4.452021148 4.443063683 4.430494353 4.414330814 4.394597186 4.371323999
#>  [241] 4.344548139 4.314312766 4.280667230 4.243666962 4.203373361 4.159853660
#>  [247] 4.113180785 4.063433198 4.010694724 3.955054376 3.896606162 3.835448875
#>  [253] 3.771685888 3.705424925 3.636777827 3.565860312 3.492791724 3.417694773
#>  [259] 3.340695270 3.261921855 3.181505721 3.099580328 3.016281115 2.931745212
#>  [265] 2.846111146 2.759518538 2.672107814 2.584019896 2.495395912 2.406376889
#>  [271] 2.317103459 2.227715566 2.138352166 2.049150943 1.960248022 1.871777684
#>  [277] 1.783872097 1.696661041 1.610271647 1.524828142 1.440451599 1.357259699
#>  [283] 1.275366500 1.194882217 1.115913011 1.038560787 0.962923005 0.889092502
#>  [289] 0.817157330 0.747200592 0.679300311 0.613529294 0.549955018 0.488639527
#>  [295] 0.429639345 0.373005397 0.318782948 0.267011557 0.217725040 0.170951451
#>  [301] 0.126713075 0.085026437 0.045902320 0.009345805 0.005000000 0.005000000
#>  [307] 0.005000000 0.005000000 0.005000000 0.005000000 0.005000000 0.005000000
#>  [313] 0.005000000 0.005000000 0.005000000 0.005000000 0.005000000 0.005000000
#>  [319] 0.005000000 0.005000000 0.005000000 0.005000000 0.005000000 0.005000000
#>  [325] 0.005000000 0.005000000 0.005000000 0.005000000 0.005000000 0.005000000
#>  [331] 0.005000000 0.005000000 0.005000000 0.005000000 0.005000000 0.005000000
#>  [337] 0.005000000 0.010435620 0.031898806 0.053516104 0.075198158 0.096857394
#>  [343] 0.118408244 0.139767342 0.160853723 0.181589012 0.201897600 0.221706820
#>  [349] 0.240947100 0.259552119 0.277458946 0.294608169 0.310944013 0.326414449
#>  [355] 0.340971290 0.354570275 0.367171143 0.378737691 0.389237827 0.398643606
#>  [361] 0.406931253 0.414081177 0.420077976 0.424910419 0.428571429 0.431058045
#>  [367] 0.432371381 0.432516563 0.431502665 0.429342627 0.426053169 0.421654688
#>  [373] 0.416171150 0.409629971 0.402061886 0.393500817 0.383983725 0.373550456
#>  [379] 0.362243583 0.350108237 0.337191938 0.323544409 0.309217399 0.294264488
#>  [385] 0.278740900 0.262703301 0.246209605 0.229318767 0.212090587 0.194585498
#>  [391] 0.176864367 0.158988287 0.141018375 0.123015571 0.105040433 0.087152943
#>  [397] 0.069412309 0.051876777 0.034603441 0.017648063 0.001064895 0.005000000
#>  [403] 0.005000000 0.005000000 0.005000000 0.005000000 0.005000000 0.005000000
#>  [409] 0.005000000 0.005000000 0.005000000 0.005000000 0.005000000 0.005000000
#>  [415] 0.005000000 0.005000000 0.005000000 0.005000000 0.005000000 0.005000000
#>  [421] 0.005000000 0.005000000 0.005000000 0.005000000 0.005000000 0.005000000
#>  [427] 0.005000000 0.005000000 0.005000000 0.005000000 0.005000000 0.005000000
#>  [433] 0.005000000 0.005000000 0.005000000 0.005000000 0.005000000 0.005000000
#>  [439] 0.005000000 0.005000000 0.005000000 0.005000000 0.005000000 0.005000000
#>  [445] 0.003651030 0.016892749 0.030127458 0.043301221 0.056360448 0.069252084
#>  [451] 0.081923809 0.094324230 0.106403075 0.118111380 0.129401677 0.140228178
#>  [457] 0.150546948 0.160316085 0.169495882 0.178048993 0.185940583 0.193138486
#>  [463] 0.199613335 0.205338702 0.210291223 0.214450709 0.217800259 0.220326354
#>  [469] 0.222018943 0.222871526 0.222881218 0.222048803 0.220378785 0.217879420
#>  [475] 0.214562736 0.210444551 0.205544472 0.199885879 0.193495908 0.186405416
#>  [481] 0.178648933 0.170264606 0.161294133 0.151782677 0.141778783 0.131334270
#>  [487] 0.120504121 0.109346360 0.097921922 0.086294501 0.074530411 0.062698411
#>  [493] 0.050869545 0.039116960 0.027515716 0.016142601 0.005075920 0.005000000
#>  [499] 0.005000000 0.005000000 0.005000000 0.005000000 0.005000000 0.005000000
#>  [505] 0.005000000 0.005000000 0.005000000 0.005000000 0.005000000 0.005000000
#>  [511] 0.005000000 0.005000000 0.005000000 0.005000000 0.005000000 0.005000000
#>  [517] 0.005000000 0.005000000 0.002305742 0.020056647 0.039863399 0.061768847
#>  [523] 0.085811660 0.112026161 0.140442181 0.171084914 0.203974793 0.239127363
#>  [529] 0.276553182 0.316257716 0.358241264 0.402498880 0.449020315 0.497789975
#>  [535] 0.548786882 0.601984659 0.657351518 0.714850272 0.774438352 0.836067842
#>  [541] 0.899685523 0.965232938 1.032646463 1.101857394 1.172792048 1.245371876
#>  [547] 1.319513586 1.395129287 1.472126633 1.550408987 1.629875597 1.710421777
#>  [553] 1.791939107 1.874315632 1.957436083 2.041182100 2.125432464 2.210063342
#>  [559] 2.294948533 2.379959727 2.464966768 2.549837921 2.634440148 2.718639389
#>  [565] 2.802300839 2.885289241 2.967469170 3.048705325 3.128862822 3.207807485
#>  [571] 3.285406138 3.361526898 3.436039463 3.508815399 3.579728424 3.648654692
#>  [577] 3.715473064 3.780065381 3.842316730 3.902115702 3.959354644 4.013929902
#>  [583] 4.065742056 4.114696151 4.160701909 4.203673939 4.243531934 4.280200856
#>  [589] 4.313611110 4.343698709 4.370405421 4.393678908 4.413472850 4.429747058
#>  [595] 4.442467567 4.451606726 4.457143261 4.459062334 4.457355583 4.452021148
#>  [601] 4.443063683 4.430494353 4.414330814 4.394597186 4.371323999 4.344548139
#>  [607] 4.314312766 4.280667230 4.243666962 4.203373361 4.159853660 4.113180785
#>  [613] 4.063433198 4.010694724 3.955054376 3.896606162 3.835448875 3.771685888
#>  [619] 3.705424925 3.636777827 3.565860312 3.492791724 3.417694773 3.340695270
#>  [625] 3.261921855 3.181505721 3.099580328 3.016281115 2.931745212 2.846111146
#>  [631] 2.759518538 2.672107814 2.584019896 2.495395912 2.406376889 2.317103459
#>  [637] 2.227715566 2.138352166 2.049150943 1.960248022 1.871777684 1.783872097
#>  [643] 1.696661041 1.610271647 1.524828142 1.440451599 1.357259699 1.275366500
#>  [649] 1.194882217 1.115913011 1.038560787 0.962923005 0.889092502 0.817157330
#>  [655] 0.747200592 0.679300311 0.613529294 0.549955018 0.488639527 0.429639345
#>  [661] 0.373005397 0.318782948 0.267011557 0.217725040 0.170951451 0.126713075
#>  [667] 0.085026437 0.045902320 0.009345805 0.005000000 0.005000000 0.005000000
#>  [673] 0.005000000 0.005000000 0.005000000 0.005000000 0.005000000 0.005000000
#>  [679] 0.005000000 0.005000000 0.005000000 0.005000000 0.005000000 0.005000000
#>  [685] 0.005000000 0.005000000 0.005000000 0.005000000 0.005000000 0.005000000
#>  [691] 0.005000000 0.005000000 0.005000000 0.005000000 0.005000000 0.005000000
#>  [697] 0.005000000 0.005000000 0.005000000 0.005000000 0.005000000 0.005000000
#>  [703] 0.010435620 0.031898806 0.053516104 0.075198158 0.096857394 0.118408244
#>  [709] 0.139767342 0.160853723 0.181589012 0.201897600 0.221706820 0.240947100
#>  [715] 0.259552119 0.277458946 0.294608169 0.310944013 0.326414449 0.340971290
#>  [721] 0.354570275 0.367171143 0.378737691 0.389237827 0.398643606 0.406931253
#>  [727] 0.414081177 0.420077976 0.424910419 0.428571429 0.431058045 0.432371381
#>  [733] 0.432516563 0.431502665 0.429342627 0.426053169 0.421654688 0.416171150
#>  [739] 0.409629971 0.402061886 0.393500817 0.383983725 0.373550456 0.362243583
#>  [745] 0.350108237 0.337191938 0.323544409 0.309217399 0.294264488 0.278740900
#>  [751] 0.262703301 0.246209605 0.229318767 0.212090587 0.194585498 0.176864367
#>  [757] 0.158988287 0.141018375 0.123015571 0.105040433 0.087152943 0.069412309
#>  [763] 0.051876777 0.034603441 0.017648063 0.001064895 0.005000000 0.005000000
#>  [769] 0.005000000 0.005000000 0.005000000 0.005000000 0.005000000 0.005000000
#>  [775] 0.005000000 0.005000000 0.005000000 0.005000000 0.005000000 0.005000000
#>  [781] 0.005000000 0.005000000 0.005000000 0.005000000 0.005000000 0.005000000
#>  [787] 0.005000000 0.005000000 0.005000000 0.005000000 0.005000000 0.005000000
#>  [793] 0.005000000 0.005000000 0.005000000 0.005000000 0.005000000 0.005000000
#>  [799] 0.005000000 0.005000000 0.005000000 0.005000000 0.005000000 0.005000000
#>  [805] 0.005000000 0.005000000 0.005000000 0.005000000 0.005000000 0.003651030
#>  [811] 0.016892749 0.030127458 0.043301221 0.056360448 0.069252084 0.081923809
#>  [817] 0.094324230 0.106403075 0.118111380 0.129401677 0.140228178 0.150546948
#>  [823] 0.160316085 0.169495882 0.178048993 0.185940583 0.193138486 0.199613335
#>  [829] 0.205338702 0.210291223 0.214450709 0.217800259 0.220326354 0.222018943
#>  [835] 0.222871526 0.222881218 0.222048803 0.220378785 0.217879420 0.214562736
#>  [841] 0.210444551 0.205544472 0.199885879 0.193495908 0.186405416 0.178648933
#>  [847] 0.170264606 0.161294133 0.151782677 0.141778783 0.131334270 0.120504121
#>  [853] 0.109346360 0.097921922 0.086294501 0.074530411 0.062698411 0.050869545
#>  [859] 0.039116960 0.027515716 0.016142601 0.005075920 0.005000000 0.005000000
#>  [865] 0.005000000 0.005000000 0.005000000 0.005000000 0.005000000 0.005000000
#>  [871] 0.005000000 0.005000000 0.005000000 0.005000000 0.005000000 0.005000000
#>  [877] 0.005000000 0.005000000 0.005000000 0.005000000 0.005000000 0.005000000
#>  [883] 0.005000000 0.002305742 0.020056647 0.039863399 0.061768847 0.085811660
#>  [889] 0.112026161 0.140442181 0.171084914 0.203974793 0.239127363 0.276553182
#>  [895] 0.316257716 0.358241264 0.402498880 0.449020315 0.497789975 0.548786882
#>  [901] 0.601984659 0.657351518 0.714850272 0.774438352 0.836067842 0.899685523
#>  [907] 0.965232938 1.032646463 1.101857394 1.172792048 1.245371876 1.319513586
#>  [913] 1.395129287 1.472126633 1.550408987 1.629875597 1.710421777 1.791939107
#>  [919] 1.874315632 1.957436083 2.041182100 2.125432464 2.210063342 2.294948533
#>  [925] 2.379959727 2.464966768 2.549837921 2.634440148 2.718639389 2.802300839
#>  [931] 2.885289241 2.967469170 3.048705325 3.128862822 3.207807485 3.285406138
#>  [937] 3.361526898 3.436039463 3.508815399 3.579728424 3.648654692 3.715473064
#>  [943] 3.780065381 3.842316730 3.902115702 3.959354644 4.013929902 4.065742056
#>  [949] 4.114696151 4.160701909 4.203673939 4.243531934 4.280200856 4.313611110
#>  [955] 4.343698709 4.370405421 4.393678908 4.413472850 4.429747058 4.442467567
#>  [961] 4.451606726 4.457143261 4.459062334 4.457355583 4.452021148 4.443063683
#>  [967] 4.430494353 4.414330814 4.394597186 4.371323999 4.344548139 4.314312766
#>  [973] 4.280667230 4.243666962 4.203373361 4.159853660 4.113180785 4.063433198
#>  [979] 4.010694724 3.955054376 3.896606162 3.835448875 3.771685888 3.705424925
#>  [985] 3.636777827 3.565860312 3.492791724 3.417694773 3.340695270 3.261921855
#>  [991] 3.181505721 3.099580328 3.016281115 2.931745212 2.846111146 2.759518538
#>  [997] 2.672107814 2.584019896 2.495395912 2.406376889