Skip to content
Snippets Groups Projects
Select Git revision
  • 97aba0d265abdd15736367c8f9989b2131983fc4
  • master default protected
  • feature/quantileforecast
  • develop
  • add_kseq
5 results

onlineforecast-package.R

Blame
  • Code owners
    Assign users and groups as approvers for specific file changes. Learn more.
    fs.Rd 1.24 KiB
    % Generated by roxygen2: do not edit by hand
    % Please edit documentation in R/fs.R
    \name{fs}
    \alias{fs}
    \title{Generation of Fourrier series.}
    \usage{
    fs(X, nharmonics)
    }
    \arguments{
    \item{X}{must be a dataframe with columns k1,k2,..., . One period is from 0 to 1
    (so for example if X is hour of day, then divide X by 24 to obtain a daily period).}
    
    \item{nharmonics}{the number of harmonics, so creates double as many inputs! i.e. one sine and one cos for each harmonic.}
    }
    \value{
    Returns a list of dataframes (two for each i in \code{1:nharmonics}) with same number of columns as X.
    }
    \description{
    Function for generating Fourrier series as a function of x E.g. use for
    harmonic functions for modelling the diurnal patterns or for basis functions.
    }
    \examples{
    # Make a data.frame with time of day in hours for different horizons
    tday <- make_tday(seq(ct("2019-01-01"), ct("2019-01-04"), by=3600), kseq=1:5)
    # See whats in it
    str(tday)
    head(tday)
    
    # Now use the function to generate Fourier series
    L <- fs(tday/24, nharmonics=2)
    # See what is in it
    str(L)
    
    # Make a plot to see the harmonics
    par(mfrow=c(2,1))
    # The first harmonic
    plot(L$sin1$k1, type="l")
    lines(L$cos1$k1, type="l")
    # The second harmonic
    plot(L$sin2$k1, type="l")
    lines(L$cos2$k1, type="l")
    
    
    }