Functional additive model

<h2 id="overview">Overview</h2>
<h3>Functional Additive Model</h3>
<p>In these models, functional predictors (
  
    
      
        X
      
    
    {\displaystyle X}
  
) are paired with responses (
  
    
      
        Y
      
    
    {\displaystyle Y}
  
) that can be either scalar or functional. The response can follow a continuous or discrete distribution and this distribution may be in the <a href="/facts/Exponential_family/1LkkqEIf">exponential family</a>. In the latter case, there would be a canonical link that connects predictors and responses. Functional predictors (or responses) can be viewed as random trajectories generated by a square-integrable stochastic process. Using <a href="/facts/Functional_principal_component_analysis/LN8Vsqmg">functional principal component analysis</a> and the <a href="/facts/Karhunen-Lo%25C3%25A8ve_expansion/HX3cIv6g">Karhunen-Loève expansion</a>, these processes can be equivalently expressed as a countable sequence of their functional principal component scores (FPCs) and eigenfunctions. In the FAM<a class="footnote-ref" id="fnref:1" href="#fn:1"><sup>1</sup></a> the responses (scalar or functional) conditional on the predictor functions are modeled as function of the functional principal component scores of the predictor function in an additive structure. This model can be categorized as a Frequency Additive Model since it is additive in the predictor FPC scores.
</p>
<h3>Continuously Additive Model</h3>
<p>The Continuously Additive Model (CAM)<a class="footnote-ref" id="fnref:2" href="#fn:2"><sup>2</sup></a> assumes additivity in the time domain. The functional predictors are assumed to be smooth across the time domain since the times contained in an interval domain are an uncountable set, an unrestricted time-additive model is not feasible. This motivates to approximate sums of additive functions by integrals so that the traditional vector additive model be replaced by a smooth additive surface. CAM can handle generalized responses paired with multiple functional predictors.
</p>
<h3>Functional Generalized Additive Model</h3>
<p>The Functional Generalized Additive Model (FGAM)<a class="footnote-ref" id="fnref:3" href="#fn:3"><sup>3</sup></a> is an extension of <a href="/facts/Generalized_additive_model/qgKRK09U">generalized additive model</a> with a scalar response and a functional predictor. This model can also deal with multiple functional predictors.
The CAM and the FGAM are essentially equivalent apart from implementation details and therefore can be covered under one description. They can be categorized as Time-Additive Models.
</p>
<h2 id="functional-additive-model">Functional Additive Model</h2>
<h3>Model</h3>
<p>Functional Additive Model for scalar and functional responses respectively, are given by
</p>

E
        (
        Y
        ∣
        X
        )
        =
        
          μ
          
            Y
          
        
        +
        
          ∑
          
            k
            =
            1
          
          
            ∞
          
        
        
          f
          
            k
          
        
        (
        
          ξ
          
            k
          
        
        )
      
    
    {\displaystyle E(Y\mid X)=\mu _{Y}+\sum _{k=1}^{\infty }f_{k}(\xi _{k})}

E
        (
        Y
        (
        t
        )
        ∣
        X
        )
        =
        
          μ
          
            Y
          
        
        (
        t
        )
        +
        
          ∑
          
            k
            =
            1
          
          
            ∞
          
        
        
          ∑
          
            m
            =
            1
          
          
            ∞
          
        
        
          f
          
            k
            m
          
        
        (
        
          ξ
          
            k
          
        
        )
        
          ψ
          
            m
          
        
        (
        t
        )
        ,
      
    
    {\displaystyle E(Y(t)\mid X)=\mu _{Y}(t)+\sum _{k=1}^{\infty }\sum _{m=1}^{\infty }f_{km}(\xi _{k})\psi _{m}(t),}

<p>where 
  
    
      
        
          ξ
          
            k
          
        
      
    
    {\displaystyle \xi _{k}}
  
 and 
  
    
      
        
          ζ
          
            m
          
        
      
    
    {\displaystyle \zeta _{m}}
  
 are FPC scores of the processes 
  
    
      
        X
      
    
    {\displaystyle X}
  
 and 
  
    
      
        Y
      
    
    {\displaystyle Y}
  
 respectively,

ϕ
          
            k
          
        
      
    
    {\displaystyle \phi _{k}}
  
 and 
  
    
      
        
          ψ
          
            m
          
        
      
    
    {\displaystyle \psi _{m}}
  
 are the eigenfunctions of processes 
  
    
      
        X
      
    
    {\displaystyle X}
  
 and 
  
    
      
        Y
      
    
    {\displaystyle Y}
  
 respectively, and 
  
    
      
        
          f
          
            k
          
        
      
    
    {\displaystyle f_{k}}
  
 and 
  
    
      
        
          f
          
            k
            m
          
        
      
    
    {\displaystyle f_{km}}
  
 are arbitrary smooth functions.
</p><p>To ensure identifiability one may require, 
  
    
      
        E
        
          f
          
            k
          
        
        (
        
          ξ
          
            k
          
        
        )
        =
        0
        ,
        
        k
        =
        1
        ,
        2
        ,
        …
        E
        
          f
          
            k
            m
          
        
        (
        
          ξ
          
            k
          
        
        )
        =
        0
        ,
        k
        =
        1
        ,
        2
        ,
        …
        m
        =
        1
        ,
        2
        ,
        …
      
    
    {\displaystyle Ef_{k}(\xi _{k})=0,\quad k=1,2,\ldots Ef_{km}(\xi _{k})=0,k=1,2,\ldots m=1,2,\ldots }

</p>
<h3>Implementation</h3>
<p>The above model is considered under the assumption that the true FPC scores 
  
    
      
        
          ξ
          
            k
          
        
      
    
    {\displaystyle \xi _{k}}
  
 for predictor processes are known. 
In general, estimation in the <a href="/facts/Generalized_additive_model/qgKRK09U">generalized additive model</a> requires <a href="/facts/Backfitting_algorithm/fpy26cXT">backfitting algorithm</a> or smooth backfitting to account for the dependencies between predictors. Now FPCs are always uncorrelated and if the predictor processes are assumed to be <a href="/facts/Gaussian/PcIz5lFw">gaussian</a> then the FPCs are independent. Then 
</p>

E
        (
        Y
        −
        
          μ
          
            Y
          
        
        
          |
        
        
          ξ
          
            k
          
        
        )
        =
        E
        {
        E
        (
        Y
        −
        
          μ
          
            Y
          
        
        
          |
        
        X
        )
        
          |
        
        
          ξ
          
            k
          
        
        }
        =
        E
        {
        
          ∑
          
            j
            =
            1
          
          
            ∞
          
        
        
          f
          
            j
          
        
        (
        
          ξ
          
            j
          
        
        )
        
          |
        
        
          ξ
          
            k
          
        
        }
        =
        
          f
          
            k
          
        
        (
        
          ξ
          
            k
          
        
        )
        ,
      
    
    {\displaystyle E(Y-\mu _{Y}|\xi _{k})=E\{E(Y-\mu _{Y}|X)|\xi _{k}\}=E\{\sum _{j=1}^{\infty }f_{j}(\xi _{j})|\xi _{k}\}=f_{k}(\xi _{k}),}

<p>similarly for functional responses 
</p>

E
        (
        
          ζ
          
            m
          
        
        
          |
        
        
          ξ
          
            k
          
        
        )
        =
        
          f
          
            k
            m
          
        
        (
        
          ξ
          
            k
          
        
        )
        ,
      
    
    {\displaystyle E(\zeta _{m}|\xi _{k})=f_{km}(\xi _{k}),}

<p>This simplifies the estimation and requires only one-dimensional <a href="/facts/Smoothing/2svFjsbI">smoothing</a> of responses against individual predictor scores and will yield consistent estimates of 
  
    
      
        
          f
          
            j
          
        
        .
      
    
    {\displaystyle f_{j}.}

In data analysis one needs to estimate 
  
    
      
        
          ξ
          
            k
          
        
      
    
    {\displaystyle \xi _{k}}
  
 before proceeding to infer the functions 
  
    
      
        
          f
          
            k
          
        
      
    
    {\displaystyle f_{k}}
  
 and 
  
    
      
        
          f
          
            k
            m
          
        
      
    
    {\displaystyle f_{km}}
  
, so there are errors in the predictors. <a href="/facts/Functional_principal_component_analysis/LN8Vsqmg">functional principal component analysis</a> generates estimates 
  
    
      
        
          
            
              
                ξ
                
                  k
                
              
              ^
            
          
        
      
    
    {\displaystyle {\hat {\xi _{k}}}}
  
 of 
  
    
      
        
          ξ
          
            k
          
        
      
    
    {\displaystyle \xi _{k}}
  
 for individual predictor trajectories along with estimates for eigenfunctions, eigenvalues, mean functions and covariance functions. Different <a href="/facts/Smoothing/2svFjsbI">smoothing</a> methods can be applied to the data 
  
    
      
        {
        
          
            
              
                ξ
                ^
              
            
          
          
            i
            k
          
        
        ,
        
          Y
          
            i
          
        
        
          }
          
            i
            =
            1
            ,
            .
            .
            .
            ,
            n
          
        
      
    
    {\displaystyle \{{\hat {\xi }}_{ik},Y_{i}\}_{i=1,...,n}}
  
 and 
  
    
      
        {
        
          
            
              
                ξ
                ^
              
            
          
          
            i
            k
          
        
        ,
        
          
            
              
                ζ
                ^
              
            
          
          
            i
            m
          
        
        
          }
          
            i
            =
            1
            ,
            .
            .
            .
            ,
            n
          
        
      
    
    {\displaystyle \{{\hat {\xi }}_{ik},{\hat {\zeta }}_{im}\}_{i=1,...,n}}
  
 to estimate 
  
    
      
        
          f
          
            k
          
        
      
    
    {\displaystyle f_{k}}
  
 and 
  
    
      
        
          f
          
            k
            m
          
        
      
    
    {\displaystyle f_{km}}
  
 respectively.
</p><p>The fitted Functional Additive Model for scalar response is given by
</p>

E
              ^
            
          
        
        (
        Y
        
          |
        
        X
        )
        =
        
          
            
              Y
              ¯
            
          
        
        +
        
          ∑
          
            k
            =
            1
          
          
            K
          
        
        
          
            
              
                f
                
                  k
                
              
              ^
            
          
        
        (
        
          ξ
          
            k
          
        
        )
        ,
      
    
    {\displaystyle {\hat {E}}(Y|X)={\bar {Y}}+\sum _{k=1}^{K}{\hat {f_{k}}}(\xi _{k}),}
  
 and the fitted Functional Additive Model for functional responses is by

E
              ^
            
          
        
        (
        Y
        (
        t
        )
        
          |
        
        X
        )
        =
        
          
            
              
                μ
                ^
              
            
          
          
            Y
          
        
        (
        t
        )
        +
        
          ∑
          
            m
            =
            1
          
          
            M
          
        
        
          ∑
          
            k
            =
            1
          
          
            K
          
        
        
          
            
              
                f
                ^
              
            
          
          
            k
            m
          
        
        (
        
          ξ
          
            k
          
        
        )
        
          
            
              
                ψ
                ^
              
            
          
          
            m
          
        
        (
        t
        )
        ,
        t
        ∈
        
          T
        
      
    
    {\displaystyle {\hat {E}}(Y(t)|X)={\hat {\mu }}_{Y}(t)+\sum _{m=1}^{M}\sum _{k=1}^{K}{\hat {f}}_{km}(\xi _{k}){\hat {\psi }}_{m}(t),t\in {T}}

<p>Note: The truncation points 
  
    
      
        K
      
    
    {\displaystyle K}
  
 and 
  
    
      
        M
      
    
    {\displaystyle M}
  
 need to be chosen data-adaptively. Possible methods include pseudo-<a href="/facts/Akaike_information_criterion/RjKKgcnF">AIC</a>, fraction of variance explained or minimization of prediction error or <a href="/facts/Cross-validation_(statistics)/pDWg8s0e">cross-validation</a>.
</p>
<h3>Extensions</h3>
<p>For the case of multiple functional predictors with a scalar response, the Functional Additive Model can be extended by fitting a functional regression which is additive in the FPCs of each of the predictor processes 
  
    
      
        
          X
          
            j
          
        
        ,
        j
        =
        1
        ,
        .
        .
        .
        ,
        d
      
    
    {\displaystyle X_{j},j=1,...,d}
  
. The model considered here is Additive Functional Score Model (AFSM) given by
</p>

E
        (
        Y
        
          |
        
        
          X
          
            1
          
        
        ,
        
          X
          
            2
          
        
        ,
        .
        .
        .
        ,
        
          X
          
            d
          
        
        )
        =
        
          ∑
          
            j
            =
            1
          
          
            d
          
        
        
          ∑
          
            k
          
        
        
          f
          
            j
            k
          
        
        (
        
          ξ
          
            j
            k
          
        
        )
      
    
    {\displaystyle E(Y|X_{1},X_{2},...,X_{d})=\sum _{j=1}^{d}\sum _{k}f_{jk}(\xi _{jk})}

<p>In case of multiple predictors the FPCs of different predictors are in general correlated and a smooth backfitting<a class="footnote-ref" id="fnref:4" href="#fn:4"><sup>4</sup></a> technique has been developed to obtain consistent estimates of the component functions 
  
    
      
        
          f
          
            j
            k
          
        
      
    
    {\displaystyle f_{jk}}
  
 when the predictors are observed with errors having unknown distribution.
</p>
<h2 id="continuously-additive-model">Continuously Additive Model</h2>
<h3>Model</h3>
<p>Since the number of time points on an interval domain is uncountable, an unrestricted time-additive model 
  
    
      
        E
        (
        Y
        
          |
        
        X
        )
        =
        
          ∑
          
            {
            t
            ∈
            [
            0
            ,
            T
            ]
            }
          
        
        
          f
          
            t
          
        
        (
        X
        (
        t
        )
        )
      
    
    {\displaystyle E(Y|X)=\sum _{\{t\in [0,T]\}}f_{t}(X(t))}
  
 is not feasible. Thus a sequence of time-additive models is considered on an increasingly dense finite time grid 
  
    
      
        
          t
          
            1
          
        
        ,
        
          t
          
            2
          
        
        ,
        .
        .
        .
        ,
        
          t
          
            m
          
        
      
    
    {\displaystyle t_{1},t_{2},...,t_{m}}
  
 in 
  
    
      
        T
      
    
    {\displaystyle T}
  
 leading to
</p>

E
        (
        Y
        
          |
        
        X
        (
        
          t
          
            1
          
        
        )
        ,
        .
        .
        .
        ,
        X
        (
        
          t
          
            m
          
        
        )
        )
        =
        E
        (
        Y
        )
        +
        
          ∑
          
            j
            =
            1
          
          
            m
          
        
        
          f
          
            j
          
        
        (
        
          X
          
            
              t
              
                j
              
            
          
        
        )
      
    
    {\displaystyle E(Y|X(t_{1}),...,X(t_{m}))=E(Y)+\sum _{j=1}^{m}f_{j}(X_{t_{j}})}
  
 where 
  
    
      
        
          f
          
            j
          
        
        (
        ⋅
        )
        =
        g
        (
        
          t
          
            j
          
        
        ,
        ⋅
        )
      
    
    {\displaystyle f_{j}(\cdot )=g(t_{j},\cdot )}
  
 for a smooth bivariate function 
  
    
      
        g
      
    
    {\displaystyle g}
  
 with 
  
    
      
        E
        (
        {
        g
        (
        
          t
          
            j
          
        
        ,
        X
        (
        
          t
          
            j
          
        
        )
        }
        )
        =
        0
      
    
    {\displaystyle E(\{g(t_{j},X(t_{j})\})=0}
  
 (to ensure identifiability). In the limit 
  
    
      
        m
        →
        ∞
      
    
    {\displaystyle m\rightarrow \infty }
  
 this becomes the continuously additive model

E
        (
        Y
        
          |
        
        X
        )
        =
        E
        (
        Y
        )
        +
        
          lim
          
            m
            →
            ∞
          
        
        
          
            1
            m
          
        
        
          ∑
          
            j
            =
            1
          
          
            m
          
        
        g
        {
        
          t
          
            j
          
        
        ,
        X
        (
        
          t
          
            j
          
        
        )
        }
        =
        E
        (
        Y
        )
        +
        
          ∫
          
            T
          
        
        g
        {
        t
        ,
        X
        (
        t
        )
        }
        d
        t
        .
      
    
    {\displaystyle E(Y|X)=E(Y)+\lim _{m\to \infty }{\frac {1}{m}}\sum _{j=1}^{m}g\{t_{j},X(t_{j})\}=E(Y)+\int _{T}g\{t,X(t)\}dt.}

<h3>Special Cases</h3>
<h4>Generalized Functional Linear Model</h4>
<p>For 
  
    
      
        g
        {
        t
        ,
        X
        (
        t
        )
        }
        =
        β
        (
        t
        )
        {
        X
        (
        t
        )
        −
        E
        X
        (
        t
        )
        }
      
    
    {\displaystyle g\{t,X(t)\}=\beta (t)\{X(t)-EX(t)\}}
  
 the model reduces to <a href="/facts/Generalized_functional_linear_model/YYb5w8ch">generalized functional linear model</a>
</p>
<h4>Functional Transformation Model</h4>
<p>For non-Gaussian predictor process, 
  
    
      
        g
        {
        t
        ,
        X
        (
        t
        )
        }
        =
        β
        (
        t
        )
        [
        ζ
        {
        X
        (
        t
        )
        }
        −
        E
        ζ
        {
        X
        (
        t
        )
        }
        ,
      
    
    {\displaystyle g\{t,X(t)\}=\beta (t)[\zeta \{X(t)\}-E\zeta \{X(t)\},}
  
 where 
  
    
      
        ζ
      
    
    {\displaystyle \zeta }
  
 is a smooth transformation of 
  
    
      
        X
        (
        t
        )
      
    
    {\displaystyle X(t)}
  
 reduces CAM to a Functional Transformation model.
</p>
<h3>Extensions</h3>
<p>This model has also been introduced with a different notation under the name Functional Generalized Additive Model (FGAM). Adding a <a href="/facts/Link_function/Sf1htm5W">link function</a> 
  
    
      
        h
      
    
    {\displaystyle h}
  
 to the mean-response and applying a probability transformation 
  
    
      
        
          G
          
            t
          
        
      
    
    {\displaystyle G_{t}}
  
 to 
  
    
      
        X
        (
        t
        )
      
    
    {\displaystyle X(t)}
  
 yields the FGAM given by 
</p>

h
        (
        E
        (
        Y
        
          |
        
        X
        )
        =
        
          θ
          
            0
          
        
        +
        
          ∫
          
            T
          
        
        F
        [
        
          G
          
            t
          
        
        {
        
          X
          
            i
          
        
        (
        t
        )
        }
        ,
        t
        ]
        d
        t
        ,
      
    
    {\displaystyle h(E(Y|X)=\theta _{0}+\int _{T}F[G_{t}\{X_{i}(t)\},t]dt,}

<p>where 
  
    
      
        
          θ
          
            0
          
        
      
    
    {\displaystyle \theta _{0}}
  
 is the intercept.
Note: For estimation and implementation see<a class="footnote-ref" id="fnref:5" href="#fn:5"><sup>5</sup></a><a class="footnote-ref" id="fnref:6" href="#fn:6"><sup>6</sup></a>
</p>

<h2 id="references">References</h2>

<ol>
<li id="fn:1"><p>Müller and Yao (2008). "Functional Additive Models". Journal of the American Statistical Association. 103 (484): 1534–1544. doi:10.1198/016214508000000751. S2CID 1927777. <a href="/wiki/Doi_(identifier)" target="_blank">/wiki/Doi_(identifier)</a> <a href="#fnref:1" class="footnote-back-ref">↩</a></p></li>
<li id="fn:2"><p>Müller, Wu and Yao (2013). "Continuously Additive models for nonlinear functional regression". Biometrika. 100 (3): 607–622. CiteSeerX 10.1.1.698.4344. doi:10.1093/biomet/ast004. <a href="/wiki/CiteSeerX_(identifier)" target="_blank">/wiki/CiteSeerX_(identifier)</a> <a href="#fnref:2" class="footnote-back-ref">↩</a></p></li>
<li id="fn:3"><p>McLean; et al. (2014). "Functional Generalized additive models". Journal of Computational and Graphical Statistics. 23 (1): 249–269. doi:10.1080/10618600.2012.729985. PMC 3982924. PMID 24729671. <a href="https://www.ncbi.nlm.nih.gov/pmc/articles/PMC3982924" target="_blank">https://www.ncbi.nlm.nih.gov/pmc/articles/PMC3982924</a> <a href="#fnref:3" class="footnote-back-ref">↩</a></p></li>
<li id="fn:4"><p>Han, Müller and Park (2017). "Smooth Backfitting for Additive Modeling with Small Errors-in-Variables, with an Application to Additive Functional Regression for Multiple Predictor Functions". Bernoulli. <a href="#fnref:4" class="footnote-back-ref">↩</a></p></li>
<li id="fn:5"><p>Müller, Wu and Yao (2013). "Continuously Additive models for nonlinear functional regression". Biometrika. 100 (3): 607–622. CiteSeerX 10.1.1.698.4344. doi:10.1093/biomet/ast004. <a href="/wiki/CiteSeerX_(identifier)" target="_blank">/wiki/CiteSeerX_(identifier)</a> <a href="#fnref:5" class="footnote-back-ref">↩</a></p></li>
<li id="fn:6"><p>McLean; et al. (2014). "Functional Generalized additive models". Journal of Computational and Graphical Statistics. 23 (1): 249–269. doi:10.1080/10618600.2012.729985. PMC 3982924. PMID 24729671. <a href="https://www.ncbi.nlm.nih.gov/pmc/articles/PMC3982924" target="_blank">https://www.ncbi.nlm.nih.gov/pmc/articles/PMC3982924</a> <a href="#fnref:6" class="footnote-back-ref">↩</a></p></li>
</ol>

Functional additive model open-in-new

Functional additive model