REVIEW 4 major objections 7 minor 49 references
FLAME: A Model for Duration-Dependent Risk Accumulation in Episodic Temporal Exposures
T0 review · 4 major / 7 minor · reviewed 2026-08-06 · deepseek-v4-flash
Pith's one-line read A single 60-minute hypotensive episode is associated with roughly 0.09 higher AKI probability than sixty 1-minute episodes in a cardiac-surgery cohort.
desk verdict A genuinely useful episode-level GAM for duration-dependent risk, with honest simulations, but the applied headline contrast is overstated and the underlying exchangeability assumption is unexamined. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The central object is the risk accumulation function $f(z)$: the amount of logit-scale risk contributed by a single exposure episode of duration $z$. It is modeled with cubic B-spline bases, equally spaced knots, a second-order difference penalty, and a half-Cauchy prior on the smoothing variance, with the constraint $f(0) \approx 0$ encoded through a tight prior on the first spline coefficient. The sum over episodes makes total duration a special case and turns the scientific question 'is one long episode riskier than many short ones?' into a posterior comparison between $f(60)$ and $60 f(1)$.
What would settle it
Re-fit FLAME on the same cohort with an added term for surgical phase or episode order; if the estimated risk curve differs by phase, or if the 0.09 contrast between one 60-minute episode and sixty 1-minute episodes shrinks to zero, the key claim is not supported.
Extended reading notes
Core claim
FLAME claims that the risk contribution of an exposure process is additive over episodes: $g(\mathbb{E}[y_i]) = X_i\beta + \sum_j f(z_{ij})$, where $z_{ij}$ is the duration of episode $j$ and $f$ is a smooth risk accumulation function estimated with Bayesian penalized B-splines. Under this model, total duration as a single covariate is the special linear case $f(z) = z^\gamma$. The paper reports that the estimated $f$ is approximately linear below 20 minutes and superlinear above, so the same total hypotension time carries different risk depending on how it is packaged: the posterior mean AKI probability rises from 0.21 with no hypotension to 0.23 for sixty one-minute episodes and to 0.32 for one sixty-minute episode, a difference of 0.09 on the probability scale with 95% credible interval [0, 0.21].
Load-bearing premise
The model assumes that every episode of low blood pressure adds its risk separately, with only its length mattering, so the timing, order, and depth of episodes have no effect.
Editorial extensions
If this is right
- If risk accumulation is superlinear beyond 20 minutes, then intraoperative intervention targets should be stated as a maximum continuous duration of hypotension, not just a total-minutes budget.
- Because total duration is nested as the linear special case, a clearly nonlinear estimated $f$ is evidence against the usual scalar-summary model.
- The posterior contrast $f(60) - 60 f(1)$ provides a single interpretable quantity for communicating risk, and its credible interval shows the strength of evidence for the episodic effect.
- Similar estimated AKI probabilities across $K = 20$, $30$, and $40$ B-spline bases indicate the headline conclusion is not an artifact of the spline basis count.
Reading between the lines
- If the additive duration-only structure holds, a natural extension is a phase-specific FLAME with a separate risk accumulation function before, during, and after cardiopulmonary bypass; material differences would refine the clinical trigger.
- The same framework could test dose packaging in other settings, such as whether short bursts of air pollution separated by clean periods carry less risk than one sustained exposure of equal average concentration.
- Because the paper excludes the eight patients with episodes over 75 minutes, the superlinear tail beyond one hour rests on sparse data; a larger cohort with longer episodes would test whether the bend at 20 minutes continues.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper proposes FLAME, a Bayesian semiparametric model for outcomes in which each exposure episode contributes a term f(duration) to the linear predictor, with f estimated by penalized B-splines under a nonnegativity constraint. The authors demonstrate through simulations that the spline estimator recovers linear, piecewise-linear, logarithmic, and sigmoid risk accumulation functions, and they apply the model to 1,188 cardiac surgery patients to estimate how the duration of intraoperative hypotension episodes relates to postoperative acute kidney injury. The application reports that a single 60-minute hypotensive episode is associated with an AKI probability of about 0.32 versus 0.23 for 60 one-minute episodes, a difference of 0.09 with a 95% credible interval that includes zero. The paper also announces an R package for the method.
Significance. If the additive episode-level model is correct, FLAME is a useful and conceptually simple addition to the toolkit for high-resolution exposure data, since it avoids collapsing episodes into total duration and directly targets contrasts such as f(60) versus 60f(1). The strengths are the systematic simulation study, the transparent Bayesian implementation in Stan, the sensitivity analysis with respect to the number of spline bases, and the provision of software. However, the application's headline contrast is not statistically significant, is estimated from very sparse long-duration data, and is only interpretable under an exchangeability/additivity assumption that the paper does not test; these issues currently prevent the paper's conclusions from being accepted at face value.
major comments (4)
- [Section 3.1, Eq. (1)] The central claim that FLAME recovers the duration-response and that f(60) versus 60f(1) is a valid contrast depends entirely on the assumption that episodes contribute additively and exchangeably on the logit scale, ignoring timing, order, episode-count interactions, and surgical phase; the model uses a single RAF over the entire surgery, as acknowledged in Section 6. The simulations in Section 4 generate outcomes from exactly this model, so they cannot validate the assumption, and the real-data analysis in Section 5 never compares FLAME with a model using total hypotensive time plus number of episodes or with phase- or order-dependent specifications. If long episodes occur disproportionately during cardiopulmonary bypass or late in surgery, or if repeated episodes sensitize the kidney, the estimated f and the headline contrast are biased summaries of duration-specific risk; please add sensitivity analyses that relax the additivity/exchangeability assumption or clearly label the application as conditional on this untested assumption.
- [Section 5.2, Table 3] The application's headline finding is not statistically significant: the 95% credible interval for the probability difference between one 60-minute episode and 60 one-minute episodes is [0, 0.21] for K=30 and [-0.01, 0.20] for K=40, both of which include zero or have a boundary at zero. The relevant exposure range is extremely sparse, since only 2% of episodes last 20 minutes or longer and only 23 patients have an episode of at least 60 minutes (Table 2), so the estimated acceleration beyond 20 minutes and the value f(60) are driven by very few observations. The abstract's '38% increase' and 'reveals' overstate the evidence; please report the posterior probability that f(60) exceeds 60f(1), discuss the implications of the sparse support, and temper the conclusion accordingly.
- [Section 4 and Table 1] The simulation study demonstrates that the spline estimator recovers f when the data are generated from the same additive episode model, but it does not test robustness to any plausible misspecification of Eq. (1), such as phase-specific effects, order dependence, or an additional effect of episode count. The simulations also only cover durations up to 30 minutes, so the estimator's behavior at f(60), the key application contrast, is never evaluated. Consequently, the statement in Section 6 that simulations 'confirm FLAME's ability to recover the true RAF' is too strong; at minimum, the claim should be limited to the correctly specified case, and ideally the authors would add a misspecification scenario, such as generating from a phase-specific or total-duration-plus-count model and examining bias in the contrast.
- [Section 3.2, Table 3] Because the half-normal priors enforce f(z) ≥ 0, the posterior of the contrast f(60) − 60f(1) is supported on [0, ∞), so the interval [0, 0.21] is a one-sided interval whose lower endpoint is a boundary value; it does not by itself convey the strength of evidence for a positive contrast. The paper should report the posterior probability P(f(60) > 60f(1)) and should also present an unconstrained fit, or otherwise justify the nonnegativity constraint, so that readers can assess how much of the estimated effect is due to the constraint.
minor comments (7)
- [Abstract and Section 5.2] The abstract and the body report inconsistent numbers: the abstract states probabilities 0.24 and 0.33 and a 38% increase, while the body reports 0.23 and 0.32 and a 37% increase; harmonize these values.
- [Section 4.1] The text refers to the third simulation row as 'exponential,' but the four simulation shapes are linear, piecewise linear, logarithm, and sigmoid; correct this wording.
- [Section 3.1] The notation f(z)=zγ is ambiguous: if γ is a multiplier, write γz, and if it is an exponent, the displayed identity sum_j z_j^γ = (sum_j z_j)γ is algebraically incorrect.
- [Abstract and Section 5.2] The R package is called 'flameRisk' in the abstract and 'flame' in the body; make the name consistent and provide a repository or link.
- [Sections 4 and 5] No MCMC convergence diagnostics, such as R-hat or effective sample size, are reported; add them for both the simulation and application analyses.
- [Section 2] The acronym 'FOSR' is introduced without being defined; spell out the term at first use.
- [Table 3] The predicted probabilities in Table 3 do not state the covariate values at which they are evaluated; specify whether these are marginal predictions or predictions at mean covariates.
Circularity Check
No significant circularity: the estimated RAF and the f(60) versus 60 f(1) contrast are posterior summaries of outcome data under a stated model, not inputs or self-citations.
full rationale
The central derivation is Eq. (1), which defines the outcome model; f is estimated by Bayesian P-splines from the outcome and exposure data. The clinical contrast expit[f(60)] - expit[60 f(1)] is a functional of the posterior of f, so it is an output of the fit, not a fitted constant relabeled as a prediction. The simulation section generates binary outcomes from the same additive model used for estimation; this is an internal consistency check of the spline estimator and does not make the real-data estimate circular, though it does not validate the additivity/exchangeability assumption. The priors (half-Normal on spline coefficients, tiny variance on gamma_1 to encode f(0) near 0) are stated modeling assumptions, not hidden identities that define the estimand. Citations to Goeddel et al. (2024) supply the dataset, and citations to Eilers and Marx, Lang and Brezger, and others supply standard P-spline machinery; none is a load-bearing self-citation or an imported uniqueness theorem. The conclusion is conditional on the additive model, but conditionality is not circularity.
Assumptions & free parameters
free parameters (4)
- Spline coefficients gamma_k =
Estimated via MCMC; values not tabulated
- Smoothing variance tau_gamma^2 =
Not reported
- Number of B-spline bases K =
30 (sensitivity: 20, 40)
- Prior hyperparameters for gamma_1 and gamma_2 =
half-Normal(0, 1e-6) and half-Normal(0, 1)
assumptions (4)
- domain assumption Episodes are exchangeable and contribute additively on the linear predictor scale, with timing and order ignored.
- domain assumption A single RAF f applies to all patients and all phases of surgery.
- domain assumption Bayesian penalized splines yield valid posterior inference for this non-Gaussian additive model.
- domain assumption f(0)=0 and non-negativity near zero, enforced by half-Normal priors.
Cite this review
Pith. "Pith review of FLAME: A Model for Duration-Dependent Risk Accumulation in Episodic Temporal Exposures." pith.science (2026). https://pith.science/paper/24ZEXILU
@misc{pith2026250622805,
author = {Pith},
title = {Pith review of: FLAME: A Model for Duration-Dependent Risk Accumulation in Episodic Temporal Exposures},
year = {2026},
howpublished = {\url{https://pith.science/paper/24ZEXILU}},
note = {Machine review of arXiv:2506.22805}
}
read the original abstract
Emerging technologies enable continuous monitoring of temporal exposures to disease risk factors, leading to complex exposure processes characterized by subject-specific numbers and durations of exposure episodes. A key scientific question is how the number and duration of such episodes influence disease risk. Existing methods typically rely on scalar summaries or time-indexed representations and are not naturally suited to model duration-dependent risk accumulation at the episode level. We introduce the FLexible Accumulation ModEl (FLAME), a semiparametric model for risk accumulation at the level of individual exposure episodes, with duration as the primary driver of risk. FLAME is motivated by and applied to quantifying the association between the duration of intraoperative hypotension and acute kidney injury (AKI) following cardiac surgery. The estimated risk accumulation function reveals that, although 60 one-minute hypotensive episodes are associated with an AKI probability of 0.24, a single sustained 60-minute episode increases that probability to 0.33, representing a 38% increase despite identical total duration. These findings provide actionable insights for intraoperative hemodynamic management and demonstrate the importance of accounting for episodic exposure patterns. While motivated by cardiac surgery, FLAME is broadly applicable to other settings involving high-resolution temporal exposures. An R package, flameRisk, is provided to facilitate application of the method in practice.
Figures
Reference graph
Works this paper leans on
-
[1]
bach2025posterior APACrefauthors Bach, P. \ Klein, N. APACrefauthors \ 2025 . Posterior Concentration Rates for Bayesian Penalized Splines Posterior concentration rates for bayesian penalized splines . Bayesian Analysis 1 1 1--21
work page 2025
-
[2]
Bai2022spikeGAM APACrefauthors Bai, R. , Moran, G E. , Antonelli, J L. , Chen, Y. \ Boland, M R. APACrefauthors \ 2022 . Spike-and-slab group lassos for grouped regression and sparse generalized additive models Spike-and-slab group lassos for grouped regression and sparse generalized additive models . Journal of the American Statistical Association 117 53...
work page 2022
-
[3]
berry1979asbestosis APACrefauthors Berry, G. , Gilson, J. , Holmes, S. , Lewinsohn, H. \ Roach, S. APACrefauthors \ 1979 . Asbestosis: a study of dose-response relationships in an asbestos textile factory Asbestosis: a study of dose-response relationships in an asbestos textile factory . Occupational and Environmental Medicine 36 2 98--112
work page 1979
-
[4]
breslow1983multiplicative APACrefauthors Breslow, N E. , Lubin, J H. , Marek, P. \ Langholz, B. APACrefauthors \ 1983 . Multiplicative models and cohort analysis Multiplicative models and cohort analysis . Journal of the American Statistical Association 78 381 1--12
work page 1983
-
[5]
claeskens2009asymptotic APACrefauthors Claeskens, G. , Krivobokova, T. \ Opsomer, J D. APACrefauthors \ 2009 . Asymptotic properties of penalized spline estimators Asymptotic properties of penalized spline estimators . Biometrika 96 3 529--544
work page 2009
-
[6]
crainiceanu2023book APACrefauthors Crainiceanu, C M. , Goldsmith, J. , Leroux, A. \ Cui, E. APACrefauthors \ 2023 . Functional Data Analysis with R Functional Data Analysis with R . Springer New York, NY, USA
work page 2023
-
[7]
Crainiceanu2005WinBUGS APACrefauthors Crainiceanu, C M. , Ruppert, D. \ Wand, M P. APACrefauthors \ 2005 . Bayesian Analysis for Penalized Spline Regression Using WinBUGS Bayesian analysis for penalized spline regression using WinBUGS . Journal of Statistical Software 14 14 1–24
work page 2005
-
[8]
de2012adaptive APACrefauthors de Jonge, R. \ van Zanten, J. APACrefauthors \ 2012 . Adaptive estimation of multivariate functions using conditionally Gaussian tensor-product spline priors Adaptive estimation of multivariate functions using conditionally Gaussian tensor-product spline priors . Electronic Journal of Statistics 6 none 1984 -- 2001 . APACrefU...
Show all 49 references
-
[9]
\ Marx, B D
Eilers1996psplines APACrefauthors Eilers, P H C. \ Marx, B D. APACrefauthors \ 1996 . Flexible smoothing with B-splines and penalties Flexible smoothing with B-splines and penalties . Statistical Science 11 2 89 -- 121 . APACrefURL https://doi.org/10.1214/ss/1038425655 APACref...
1996
-
[10]
, Colantuoni, E
gellar2014 APACrefauthors Gellar, J. , Colantuoni, E. , Needham, D. \ Crainiceanu, C. APACrefauthors \ 2014 . Variable-Domain Functional Regression for Modeling ICU Data Variable-domain functional regression for modeling ICU data . Journal of American Statistical Association 1...
2014
-
[11]
APACrefauthors \ 2006
Gelman2006VariancePrior APACrefauthors Gelman, A. APACrefauthors \ 2006 . Prior distributions for variance parameters in hierarchical models (comment on article by Browne and Draper) Prior distributions for variance parameters in hierarchical models (comment on article by Brow...
2006 doi
-
[12]
\ van der Vaart, A
ghosal2007convergence APACrefauthors Ghosal, S. \ van der Vaart, A. APACrefauthors \ 2007 . Convergence rates of posterior distributions for noniid observations Convergence rates of posterior distributions for noniid observations . The Annals of Statistics 35 1 192 -- 223 . AP...
2007 doi
-
[13]
, Hernandez, M
goeddel2024mapcvp APACrefauthors Goeddel, L A. , Hernandez, M. , Koffman, L. , Murphy, Z. , Khanna, A K. , Robich, M. Faraday, N. APACrefauthors \ 2024 . Fine-Mapping the Association of Acute Kidney Injury with Mean Arterial and Central Venous Pressures during Coronary Artery ...
2024
-
[14]
, Ruppert, D
HarezlakRuppertWand2018Book APACrefauthors Harezlak, J. , Ruppert, D. \ Wand, M. APACrefauthors \ 2018 . Semiparametric Regression with R Semiparametric Regression with R . Springer New York, NY, USA
2018
-
[15]
\ Tibshirani, R
hastie1986gam APACrefauthors Hastie, T. \ Tibshirani, R. APACrefauthors \ 1986 . Generalized additive models Generalized additive models . Statistical Science 1 3 297 -- 310 . APACrefURL https://doi.org/10.1214/ss/1177013604 APACrefURL APACrefDOI doi:10.1214/ss/1177013604 APACrefDOI
1986
-
[16]
, Wellmann, J
hauptmann2000analysis APACrefauthors Hauptmann, M. , Wellmann, J. , Lubin, J H. , Rosenberg, P S. \ Kreienbrock, L. APACrefauthors \ 2000 . Analysis of exposure-time-response relationships using a spline weight function Analysis of exposure-time-response relationships using a ...
2000
-
[17]
, Liu, X
Hu2023Exposome APACrefauthors Hu, H. , Liu, X. , Zheng, Y. , He, X. , Hart, J. , James, P. Bian, J. APACrefauthors \ 2023 . Methodological challenges in spatial and contextual exposome-health studies Methodological challenges in spatial and contextual exposome-health studies ....
2023
-
[18]
Johns, V Z., Ciprian Crainiceanu \ Gellar, J
Johns02102019 APACrefauthors Jordan T. Johns, V Z., Ciprian Crainiceanu \ Gellar, J. APACrefauthors \ 2019 . Variable-Domain Functional Principal Component Analysis Variable-domain functional principal component analysis . Journal of Computational and Graphical Statistics 28 4...
2019
-
[19]
APACrefauthors \ 2012
khwaja2012kdigo APACrefauthors Khwaja, A. APACrefauthors \ 2012 . KDIGO clinical practice guidelines for acute kidney injury KDIGO clinical practice guidelines for acute kidney injury . Nephron Clinical Practice 120 4 c179--c184
2012
-
[20]
\ Brezger, A
lang2004bayesianpsplines APACrefauthors Lang, S. \ Brezger, A. APACrefauthors \ 2004 . Bayesian P-splines Bayesian p-splines . Journal of computational and graphical statistics 13 1 183--212
2004
-
[21]
\ Ruppert, D
li2008asymptotics APACrefauthors Li, Y. \ Ruppert, D. APACrefauthors \ 2008 . On the asymptotics of penalized splines On the asymptotics of penalized splines . Biometrika 95 2 415--436
2008
-
[22]
\ Zhang, H H
lin2006component APACrefauthors Lin, Y. \ Zhang, H H. APACrefauthors \ 2006 . Component selection and smoothing in multivariate nonparametric regression Component selection and smoothing in multivariate nonparametric regression . The Annals of Statistics 34 5 2272 -- 2297 . AP...
2006 doi
-
[23]
\ Wood, S N
marra2011practical APACrefauthors Marra, G. \ Wood, S N. APACrefauthors \ 2011 . Practical variable selection for generalized additive models Practical variable selection for generalized additive models . Computational Statistics & Data Analysis 55 7 2372--2387
2011
-
[24]
, Hooker, G
mclean2014functional APACrefauthors McLean, M. , Hooker, G. , Staicu, A M. , Scheipl, F. \ Ruppert, D. APACrefauthors \ 2014 . Functional generalized additive models Functional generalized additive models . Journal of Computational and Graphical Statistics 23 1 249--269
2014
-
[25]
\ Yao, F
muller2008fam APACrefauthors M \"u ller, H G. \ Yao, F. APACrefauthors \ 2008 . Functional additive models Functional additive models . Journal of the American Statistical Association 103 484 1534--1544
2008
-
[26]
APACrefauthors \ 2003
Nieuwenhuijsen2003Exposure APACrefauthors Nieuwenhuijsen, M J. APACrefauthors \ 2003 . Exposure Assessment in Occupational and Environmental Epidemiology Exposure assessment in occupational and environmental epidemiology . Oxford University Press . APACrefURL https://doi.org/1...
2003
-
[27]
APACrefauthors \ 1986
osullivan1986 APACrefauthors O’Sullivan, F. APACrefauthors \ 1986 . A statistical perspective on ill-posed inverse problems (with discussion) A statistical perspective on ill-posed inverse problems (with discussion) . Statistical Science 1 4 505–527
1986
-
[28]
\ Smith, M
Panagiotelis2008FunctionSelection APACrefauthors Panagiotelis, A. \ Smith, M. APACrefauthors \ 2008 . Bayesian identification, selection and estimation of semiparametric functions in high-dimensional additive models Bayesian identification, selection and estimation of semipara...
2008
-
[29]
\ Silverman, B
ramsaysilv2005 APACrefauthors Ramsay, J. \ Silverman, B. APACrefauthors \ 2005 . Functional Data Analysis Functional Data Analysis . Springer New York, NY, USA
2005
-
[30]
, Lafferty, J
ravikumar2009SAM APACrefauthors Ravikumar, P. , Lafferty, J. , Liu, H. \ Wasserman, L. APACrefauthors \ 2009 . Sparse additive models Sparse additive models . Journal of the Royal Statistical Society Series B: Statistical Methodology 71 5 1009--1030
2009
-
[31]
, Goldsmith, J
reiss2017sofr APACrefauthors Reiss, P T. , Goldsmith, J. , Shang, H L. \ Ogden, R T. APACrefauthors \ 2017 . Methods for Scalar-on-Function Regression Methods for scalar-on-function regression . International Statistical Review 85 2 228-249
2017
-
[32]
, Greenland, S
rothman2008modernEpi APACrefauthors Rothman, K. , Greenland, S. \ Lash, T. APACrefauthors \ 2008 . Modern Epidemiology Modern epidemiology . Wolters Kluwer Health/Lippincott Williams & Wilkins . APACrefURL https://books.google.com/books?id=Z3vjT9ALxHUC APACrefURL
2008
-
[33]
, Wand, M
ruppert2003semiparametric APACrefauthors Ruppert, D. , Wand, M. \ Carroll, R. APACrefauthors \ 2003 . Semiparametric Regression Semiparametric regression . Cambridge University Press . APACrefURL https://books.google.com/books?id=Y4uEvXFP2voC APACrefURL
2003
-
[34]
, Fahrmeir, L
Scheipl2012STAR APACrefauthors Scheipl, F. , Fahrmeir, L. \ Kneib, T. APACrefauthors \ 2012 . Spike-and-Slab Priors for Function Selection in Structured Additive Regression Models Spike-and-slab priors for function selection in structured additive regression models . Journal o...
2012
-
[35]
\ Krivobokova, T
schwarz2016unified APACrefauthors Schwarz, K. \ Krivobokova, T. APACrefauthors \ 2016 . A unified framework for spline estimators A unified framework for spline estimators . Biometrika 103 1 121--131
2016
-
[36]
\ APSF Hemodynamic Instability Writing Group
Scott2024perioperative APACrefauthors Scott, M J. \ APSF Hemodynamic Instability Writing Group . APACrefauthors \ 2024 . Perioperative patients with hemodynamic instability: consensus recommendations of the anesthesia patient safety foundation Perioperative patients with hemod...
2024
-
[37]
\ Abrahamowicz, M
sylvestre2009flexible APACrefauthors Sylvestre, M P. \ Abrahamowicz, M. APACrefauthors \ 2009 . Flexible modeling of the cumulative effects of time-dependent exposures on the hazard Flexible modeling of the cumulative effects of time-dependent exposures on the hazard . Statist...
2009
-
[38]
APACrefauthors \ 1988
thomas1988models APACrefauthors Thomas, D C. APACrefauthors \ 1988 . Models for exposure-time-response relationships with applications to cancer epidemiology Models for exposure-time-response relationships with applications to cancer epidemiology . Annual Review of Public Heal...
1988
-
[39]
APACrefauthors \ 1963
tikhonov1963regularization APACrefauthors Tikhonov, A. APACrefauthors \ 1963 . Solution of incorrectly formulated problems and the regularization method Solution of incorrectly formulated problems and the regularization method . Soviet Mathematics Doklady
1963
-
[40]
, Klein, N
umlauf2018bamlss APACrefauthors Umlauf, N. , Klein, N. \ Zeileis, A. APACrefauthors \ 2018 . BAMLSS : Bayesian additive models for location, scale, and shape (and beyond) BAMLSS : Bayesian additive models for location, scale, and shape (and beyond) . Journal of Computational a...
2018
-
[41]
APACrefauthors \ 1997
vacek1997assessing APACrefauthors Vacek, P M. APACrefauthors \ 1997 . Assessing the effect of intensity when exposure varies over time Assessing the effect of intensity when exposure varies over time . Statistics in Medicine 16 5 505--513
1997
-
[42]
, Gelman, A
Vehtari2017loo APACrefauthors Vehtari, A. , Gelman, A. \ Gabry, J. APACrefauthors \ 2017 . Practical Bayesian model evaluation using leave-one-out cross-validation and WAIC Practical bayesian model evaluation using leave-one-out cross-validation and waic . Statistics and Compu...
2017 doi
-
[43]
APACrefauthors \ 1990
wahba1990spline APACrefauthors Wahba, G. APACrefauthors \ 1990 . Spline models for observational data Spline models for observational data . Society for Industrial and Applied Mathematics
1990
-
[44]
\ Ormerod, J
Wand2008OSul APACrefauthors Wand, M. \ Ormerod, J. APACrefauthors \ 2008 06 . On semiparametric regression with O'Sullivan penalized splines On semiparametric regression with O'Sullivan penalized splines . Australian & New Zealand Journal of Statistics 50 179 - 198 . APACrefDO...
2008 arXiv
-
[45]
APACrefauthors \ 2017
wood2017gam APACrefauthors Wood, S. APACrefauthors \ 2017 . Generalized Additive Models: An Introduction with R Generalized Additive Models: An Introduction with R \ ( 2 \ ). Chapman and Hall/CRC
2017
-
[46]
APACrefauthors \ 2019
xiao2019theory APACrefauthors Xiao, L. APACrefauthors \ 2019 . Asymptotic theory of penalized splines Asymptotic theory of penalized splines . Electronic Journal of Statistics 13 1 747 -- 794 . APACrefURL https://doi.org/10.1214/19-EJS1541 APACrefURL APACrefDOI doi:10.1214/19-...
2019 doi
-
[47]
xiao2012local APACrefauthors Xiao, L. , Li, Y. , Apanasovich, T V. \ Ruppert, D. APACrefauthors \ 2012 . Local asymptotics of P-splines Local asymptotics of p-splines . arXiv preprint arXiv:1201.0708
2012 arXiv
-
[48]
, Abrahamowicz, M
xiao2014flexible APACrefauthors Xiao, Y. , Abrahamowicz, M. , Moodie, E E. , Weber, R. \ Young, J. APACrefauthors \ 2014 . Flexible marginal structural models for estimating the cumulative effect of a time-dependent treatment on the hazard: reassessing the cardiovascular risks...
2014
-
[49]
\ Ghosal, S
yoo2016supremum APACrefauthors Yoo, W W. \ Ghosal, S. APACrefauthors \ 2016 . SUPREMUM NORM POSTERIOR CONTRACTION AND CREDIBLE SETS FOR NONPARAMETRIC MULTIVARIATE REGRESSION Supremum norm posterior contraction and credible sets for nonparametric multivariate regression . The A...
2016
Reviewed August 6, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.