Tolerance bands for exponential family functional data

Galappaththige S. R. de Silva,Pankaj K. Choudhary
DOI: https://doi.org/10.1002/cjs.11808
2024-06-24
Canadian Journal of Statistics
Abstract:A tolerance band for a functional response provides a region that is expected to contain a given fraction of observations from the sampled population at each point in the domain. This band is a functional analogue of the tolerance interval for a univariate response. Although the problem of constructing functional tolerance bands has been considered for a Gaussian response, it has not been investigated for non‐Gaussian responses, which are common in biomedical applications. We describe a methodology for constructing tolerance bands for two non‐Gaussian members of the exponential family: binomial and Poisson. The approach is to first model the data using the framework of generalized functional principal components analysis. Then, a parameter is identified in which the marginal distribution of the response is stochastically monotone. We show that the tolerance limits can be readily obtained from confidence limits for this parameter, which in turn can be computed using large‐sample theory and bootstrapping. Our proposed methodology works for both dense and sparse functional data. We report the results of simulation studies designed to evaluate its performance and get recommendations for practical applications. We illustrate our proposed method using two actual biomedical studies, and also provide computer source code that implements our method. Résumé Une bande de tolérance pour une réponse fonctionnelle fournit une région censée contenir une fraction donnée d'observations de la population échantillonnée à chaque point du domaine. Cette bande est un analogue fonctionnel de l'intervalle de tolérance pour une réponse univariée. Bien que la construction de bandes de tolérance fonctionnelles ait été étudiée pour les réponses gaussiennes, peu de travaux se sont intéressés au cas non gaussien, pourtant fréquent en applications biomédicales. Dans ce travail, les auteurs décrivent une approche pour construire des bandes de tolérance pour deux lois de la famille exponentielle non gaussiennes: binomiale et Poisson. Leur méthodologie consiste premièrement à modéliser les données fonctionnelles par une analyse par composantes principales généralisée. Un paramètre, dont la loi marginale de la réponse est stochastiquement monotone, est ensuite identifié. Les auteurs montrent que les limites de tolérance peuvent alors être déduites des limites de confiance de ce paramètre, limites qui peuvent à leur tour être calculées à l'aide de la théorie asymptotique et du bootstrap. L'approche proposée s'applique aussi bien aux données fonctionnelles denses qu'éparses. Ses performances ont été évaluées au moyen d'études de simulation, ce qui a permis aux auteurs de formuler des recommandations pour des applications pratiques. Ils ont également fourni le code informatique implémentant la méthode et l'ont illustrée par deux applications biomédicales réelles.
statistics & probability
What problem does this paper attempt to address?