{"id":"779639b9-85de-4a92-b33f-2fab615fa3e8","arxiv_id":"2411.18416","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"A Bayesian mixed model using norm-preserving phase actions can recover the size-and-shape, rather than the exact values, of a population-level function from noisy, misaligned functional data.","lead":"This paper proposes a Bayesian functional mixed model in which random phase distortions act as rotations of the function space, and shows that the size-and-shape of the average function can be reliably estimated even when the full function cannot. The approach reframes a known identifiability problem in functional data analysis by targeting a geometric invariant, and reports gains over the warpMix baseline on simulated and real growth and ECG data.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The claim that the size-and-shape of mu is recoverable rests on an unproved identifiability of phase versus amplitude in Eq. (6); without a formal check, the numerical demonstrations do not establish the claim.","rationale":"Reading the paper in good faith, the central claim is not that every square-integrable mu can be pointwise recovered, but that its norm-preserving orbit [mu] can be recovered. For this to be true, the marginal likelihood in Eq. (6) must separate the common amplitude coefficients a from the per-observation phase functions gamma_i and from the random-effect coefficients c_i. The paper asserts in Section 3.2 that the phase priors 'guard against confounding', and Section 5 explicitly says theoretical support is lacking. That admission is decisive: the only evidence for the central claim is a small set of single-chain MCMC experiments, without standard errors or code, plus comparisons that select the best of four model variants. Appendix H does show robustness to over-specification of Bf and to misspecification of Br and theta_gamma, which is genuine supporting evidence, but it also shows that under-specifying Bf destroys recovery, which is exactly the kind of truncation/prior dependence that needs a formal identifiability statement. The proposed Fisher-information check is a concrete way to decide whether the phase-amplitude trade-off is locally flat; if it is flat, the posterior mean cannot be trusted to concentrate on [mu], and the abstract's claim would need to be qualified. Because the reader's conditional verdict already reflects this missing support, the verdict is unchanged, but the test identifies precisely what would strengthen or overturn it.","tokens_in":21299,"tokens_out":11055,"duration_ms":114132,"concrete_test":"Derive the expected Fisher information matrix of the marginal likelihood in Eq. (6) with respect to theta = (a, sigma^2, sigma_c^2, gamma_1,...,gamma_n) at the true parameter values used in Example 1 (Bf = Br = 6, n = 30, PM2 with T_gamma = 7), and compute its null space after removing the orbit direction generated by the norm-preserving action. If any non-orbit direction in (a, gamma_i) has zero information, the likelihood is locally flat along a phase-amplitude trade-off and the claimed recovery of size-and-shape is not identifiable; if the null space is exactly the orbit direction, local identifiability of [mu] holds and the concern is resolved.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The weakest load-bearing point is the unproved identifiability of the size-and-shape of mu from the marginal likelihood in Eq. (6). The parameters theta = (a, sigma^2, sigma_c^2, gamma_1,...,gamma_n) enter the likelihood through both the rotated mean Phi_i a and the rotated covariance sigma^2 diag(gamma_dot_i) + sigma_c^2 Phi_tilde_i Phi_tilde_i^T. Section 3.2 states that the priors are chosen 'to guard against confounding of inference between {gamma_i} and mu', but it supplies no theorem that the map from theta to the sampling distribution is injective on the quotient by the norm-preserving action, nor that the posterior contracts to the true orbit [mu]. Section 5 explicitly concedes: 'What is lacking is theoretical support for the same.' This gap is load-bearing because the central claim is exactly that size-and-shape recovery is possible. The numerical evidence consists of single MCMC runs on data generated either from the same model (Example 1) or from a different warping model (Example 2), with no standard errors and no released code, and Appendix H shows that under-specifying Bf visibly destroys recovery of mu. If the prior/truncation combination does not fully break the phase-amplitude trade-off, the centered posterior mean can be biased toward a prior-preferred representative unrelated to the true size-and-shape, so the abstract's claim is not established for general square-integrable mu.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a Bayesian functional mixed model for a fixed effect function μ in the presence of object-level phase and amplitude variability. The observation model is fi = D_{γ_i}(μ + v_i + ε_i), where D_γ is the norm-preserving (square-root slope) action of the phase group Γ. The likelihood is obtained by marginalizing the random-effect coefficients, leading to the multivariate normal model in Eq. (6). Two prior models for phase functions are proposed, one one-parameter and one Dirichlet-based, and posterior inference is carried out by MCMC with Metropolis-Hastings updates. The main claim is that the size-and-shape of μ, i.e., its equivalence class under the norm-preserving action, can be recovered from the posterior, even though pointwise recovery is impossible. Evidence is provided through two simulation studies and several real-data examples, with comparisons to the warpMix estimator and a Bayesian registration approach.","tokens_in":21665,"tokens_out":2783,"duration_ms":29438,"significance":"If the central claim is established, the paper would make a useful contribution: it offers a Bayesian functional mixed model with unrestricted phase variation, avoids finite-rank covariance assumptions on the error process, and provides a principled way to summarize what can be learned about μ under a strong symmetry group. The simulations show competitive or better performance than warpMix in several examples, and the real-data analyses are visually plausible. The paper also ships detailed derivations of the marginal likelihood and Metropolis-Hastings ratios (Appendices C and D), which are valuable. However, the central claim is currently supported only by numerical demonstrations, and the paper explicitly concedes in Section 5 that theoretical support is lacking. The numerical evidence consists of single MCMC runs without repeated-experiment error bars, and the evaluation metric is partially tied to the chosen centering procedure. The significance of the work therefore depends on whether the identifiability and posterior-concentration issue can be resolved or at least sharply characterized.","major_comments":[{"comment":"The abstract claims to 'demonstrate that it is possible to recover the size-and-shape of a square-integrable μ', but the paper explicitly states in Section 5: 'What is lacking is theoretical support for the same.' This is a load-bearing gap: the model in Eq. (6) involves parameters θ = (a, σ², σ_c², γ_1,...,γ_n) where the phase functions enter both the mean Φ_i a and the covariance σ²diag(γ̇_i) + σ_c² Φ̃_i Φ̃_i^T. No theorem is given that the map from θ to the sampling distribution is identifiable on the quotient by the norm-preserving action, nor that the posterior contracts to the true orbit [μ]. Without such a result, the numerical demonstrations do not establish the general claim for square-integrable μ. A concrete remedy would be to prove identifiability of the orbit under the stated priors, or to state the claim as a conjecture supported by a systematic simulation study with varied μ, n, signal-to-noise ratio, and repeated seeds.","section":"Section 5 and Eq. (6)"},{"comment":"The numerical evidence consists of single MCMC runs: Example 1 draws data from the model itself, and Example 2 draws from the warpMix model, but no repeated simulations, Monte Carlo standard errors, or error bars are reported for the entries in Table 1. Consequently, it is unclear whether the reported improvements over warpMix are systematic or within Monte Carlo noise. I recommend reporting repeated-seed summaries (means and standard errors) for the error criterion Δμ, and ideally reporting an orbit-level error that is invariant to the norm-preserving action, rather than only a pointwise criterion applied after the posterior centering step.","section":"Section 4.1 and Table 1"},{"comment":"The recovery of μ is evaluated after centering posterior samples by the average posterior phase γ̄, i.e., (μ̂_j ∘ γ̄)√(γ̄̇). This centering is reasonable for visualization, but it means the 'recovered' representative is partly determined by the posterior phase estimates. If the posterior for γ_i is biased, the centered mean may lie in the wrong equivalence class while still looking well-centered. The paper does not provide an orbit-invariant measure of recovery error, so the numerical demonstrations do not fully separate phase recovery from amplitude recovery. A simple additional check would be to report distances between orbits, e.g., inf_{γ∈Γ} ‖μ̂_centered - D_γ(μ_true)‖, or to assess recovery of γ_i and μ jointly under a known generative model.","section":"Section 4 and Appendix H"},{"comment":"Appendix H shows that under-specifying Bf, the number of basis functions for μ, visibly destroys recovery of μ, while the abstract promises recovery for general square-integrable μ. Since the model represents μ in a finite-dimensional basis and no data-driven selection of Bf is provided, the central claim is sensitive to a tuning parameter whose correct value is generally unknown. The paper should either restrict the claim to functions well-approximated by the chosen basis, or provide a practical procedure for choosing Bf and evidence that recovery is robust to reasonable misspecification.","section":"Appendix H"}],"minor_comments":[{"comment":"The caption contains a typo: 'warpMix esimate' should be 'warpMix estimate'.","section":"Section 4, Figure 2 caption"},{"comment":"The one-parameter phase family γ(t) = t + α t(t−1) with α ∈ (−1,1) is stated to have γ̇(t) ∈ (0,2); please verify and state the exact range of γ̇, since the boundary behavior affects the prior support.","section":"Eq. (7)"},{"comment":"The notation γ̇^{-1}_{cur}(·) in the Jacobian expression is confusing; it would be clearer to write the derivative of γ_cur^{-1} evaluated at the relevant point.","section":"Appendix D, Eq. (20)"},{"comment":"The statement that the measure of an equivalence class under the norm-preserving action is 'larger' than under the value-preserving action is heuristic and is not used in the subsequent theory or experiments; consider clarifying or removing it to avoid overclaiming.","section":"Section 2, item (iii)"}],"recommendation":"major_revision","confidential_remarks":"The paper has a clear and interesting idea, and the algorithmic derivations appear sound. My main concern is that the central claim—recovery of the size-and-shape of μ—is explicitly conceded in Section 5 to lack theoretical support, and the numerical evidence is not strong enough to fully compensate (single runs, no repeated-seed error bars, orbit-level error not reported). This is fixable within the manuscript's scope by adding an identifiability or posterior-concentration analysis, or by carefully restating the claim and providing a more systematic simulation study. I therefore recommend major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Reading this made me think the editors should send it out, not desk reject it. The model itself is new: a Bayesian functional mixed model with phase variation through the norm-preserving action, and with the target of inference explicitly the size-and-shape equivalence class of the fixed effect mu. That framing is genuinely different from Claeskens et al., who use value-preserving warping, and from Cheng et al. The derivations in Appendices C and D look correct; the MH ratios are carefully worked out. The warpMix-based simulation is the strongest part: data generated under a value-preserving model, yet the proposed estimator beat warpMix on two of three targets. That is an independent check, not a circle. The real data examples are plausible and the sensitivity analysis in Appendix H is honest about what breaks. The soft spot is exactly where the reader and the stress-test note point: there is no identifiability theorem. Section 5 says 'What is lacking is theoretical support for the same.' That is an honest sentence, but it means the central claim that recovery of the size-and-shape is possible is a numerical demonstration, not an established property. The phase prior is supposed to guard against confounding, but there is no proof that the map from parameters to the sampling distribution is injective on the quotient by the norm-preserving action. If the prior or basis truncation does not fully separate phase from amplitude, the centered posterior mean could be a prior-biased representative rather than the true orbit. Also, the simulations are single runs without repeated-experiment error bars, and the comparison highlights the best of four model variants. No code is released. That said, the paper does not hide its limitations; Appendix H shows that under-specified Bf visibly destroys recovery, which is useful and honest. Who is this for? Statisticians working on functional mixed models or elastic functional data analysis. It is serious work, not a hack. But it should not be accepted on the strength of the abstract's claim as stated. My recommendation: send to peer review, with a referee who knows identifiability theory. The authors have the algebra; they need either a contraction or identifiability result for the quotient, or a weakening of the claim to 'numerical evidence suggests.' A conditional accept is plausible after that.","headline":"A genuinely new Bayesian functional mixed model with sound algebra and one solid independent simulation, but the headline recovery claim is only numerically supported and the missing identifiability theory should shape the review.","tokens_in":22184,"tokens_out":2095,"would_cite":false,"duration_ms":20755,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["62F15","62R10","62G08"],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper demonstrates that a Bayesian functional mixed model can recover the size-and-shape of a square-integrable fixed effect function from noisy replicates with phase variation, although pointwise recovery is impossible.","keywords":["functional mixed models","size-and-shape analysis","phase variability","norm-preserving action","Bayesian inference","unitary transformations","function alignment"],"falsifier":"Simulate $n$ curves from $f_i=D_{\\gamma_i}(\\mu+v_i+\\epsilon_i)$ with phase functions drawn from a flexible family outside the two priors and with a true $\\mu$ that is not exactly representable in the chosen basis; compute $d=\\inf_{\\gamma\\in\\Gamma}\\|\\hat\\mu_{\\mathrm{center}}-D_\\gamma\\mu\\|$ for the centered posterior mean. If $d$ does not decrease as $n$ grows, the claim that the size-and-shape is recoverable is falsified.","tokens_in":21115,"feed_emoji":"📐","tokens_out":10236,"duration_ms":85064,"temperature":0.7,"pith_summary":"In functional mixed models, observations are the sum of a population average, individual random deviations, and noise, and when individual time-warps (phase variation) are present, the average function cannot be recovered pointwise. This paper identifies a coarser target that can be recovered: the size-and-shape of the population function, defined as its equivalence class under the norm-preserving action $D_\\gamma(f)=(f\\circ\\gamma)\\sqrt{\\dot\\gamma}$, where $\\gamma$ ranges over the group of smooth increasing time-warps of $[0,1]$. The authors build a Bayesian mixed model in $L^2[0,1]$ in which each observation is $f_i=D_{\\gamma_i}(\\mu+v_i+\\epsilon_i)$, so $\\gamma_i$ acts as a rotation of the function space rather than a simple time-shift. They report that phase-centering the posterior samples of $\\mu$ gives a posterior mean that tracks the true size-and-shape in simulations and real data, and that it improves on the current state-of-the-art estimator. The practical payoff is that meaningful uncertainty quantification survives for a geometric summary of the population function even when the full function is unidentifiable.","feed_headline":"Bayesian model recovers a function's size-and-shape, not its values","feed_subtitle":"A size-and-shape target makes an average function meaningful when individual timing differences hide it.","key_machinery":"The load-bearing object is the norm-preserving action $D_\\gamma(f)=(f\\circ\\gamma)\\sqrt{\\dot\\gamma}$, where $\\gamma$ belongs to the group $\\Gamma$ of increasing diffeomorphisms of $[0,1]$; because it preserves the $L^2$ norm, each $D_\\gamma$ is a unitary transformation, i.e., an infinite-dimensional rotation of the Hilbert space. The paper defines the size-and-shape of $f$ as the equivalence class $\\{D_\\gamma(f):\\gamma\\in\\Gamma\\}$, treats the random phase $\\gamma_i$ as a rotation of the coordinate system in which $\\mu$ and $v_i$ are expressed, and carries out inference on the phase-rotated basis coefficients. Two prior families on $\\Gamma$---a one-parameter quadratic family and a Dirichlet-increment prior on discretized warps---regularize the rotation, and the marginal likelihood (integrating out $v_i$) is multivariate normal with mean $\\Phi_i a$ and covariance $\\sigma^2\\mathrm{diag}(\\dot\\gamma_i)+\\sigma_c^2\\tilde\\Phi_i\\tilde\\Phi_i^T$. Phase centering of posterior draws, using the average posterior phase, converts the posterior into summaries of the size-and-shape orbit.","core_discovery":"The paper's central claim is that the size-and-shape of a square-integrable fixed effect $\\mu$---its $D_\\gamma$-orbit in $L^2[0,1]$---is reliably recoverable under the Bayesian functional mixed model $f_i=D_{\\gamma_i}(\\mu+v_i+\\epsilon_i)$, despite the well-known non-identifiability of $\\mu$ itself. The norm-preserving operator $D_\\gamma(f)=(f\\circ\\gamma)\\sqrt{\\dot\\gamma}$ is a unitary rotation of the Hilbert space, so the model separates size-and-shape preserving deviations (the random phase $\\gamma_i$) from size-and-shape altering ones (the random function $v_i$ and measurement error $\\epsilon_i$). The authors estimate the posterior of $\\mu$, then center every posterior draw by the average posterior phase $\\bar\\gamma$: $(\\hat\\mu^j\\circ\\bar\\gamma)\\sqrt{\\dot{\\bar\\gamma}}$; the centered draws and their pointwise mean are taken as summaries of the size-and-shape class. In simulations the centered posterior mean recovers the true orbit even when the data were generated with value-preserving warping (the mechanism assumed by the state-of-the-art), and in real data examples it preserves topological features such as growth spurts and ECG complexes that the state-of-the-art estimate misses. The paper is explicit that theoretical support for this recovery is not yet provided; the claim is established through posterior sampling and numerical comparisons.","pith_inferences":["A likely formal consequence of the numerical findings is a posterior contraction theorem for the orbit $[\\mu]$ under the norm-preserving action; the natural next step is to prove that the phase-centered posterior accumulates at the true equivalence class as $n$ grows, with a rate depending on the support of the phase prior.","If the equivalence class under the norm-preserving action is genuinely larger than under the value-preserving action, the method should exhibit faster orbit recovery in high phase-variability regimes; this is a testable prediction that could be checked by comparing centered posterior distances across simulated phase distributions with increasing entropy.","The paper's own loss function for comparing estimators (squared pointwise error on centered representatives) implicitly rewards pointwise fidelity within the orbit; a stricter evaluation would report the distance to the orbit itself, which may change the apparent margin over the state of the art.","The framework suggests a general recipe for mixed effects models with nuisance symmetries: choose a group action whose orbits define the identifiable quantity, put a prior on the group element, and summarize the fixed effect by an orbit-centered posterior; applying this recipe with rotation groups rather than time-warps would yield Bayesian size-and-shape analysis for landmark or curve data."],"forward_implications":["If the claim holds, practitioners can report a meaningful average function for populations whose members differ in timing, because the centered posterior mean targets the size-and-shape orbit rather than a misaligned pointwise mean.","The model produces posterior credible intervals for that orbit, so uncertainty quantification for the geometric average is available without finite-rank covariance assumptions on the error process.","Because the phase functions rotate the chosen basis, a modest number of fixed Fourier or B-spline basis functions can represent complex fixed effects; the paper reports that this beats an empirical FPCA basis fitted to the same data.","The size-and-shape target remains estimable when the true data generator uses value-preserving warping, which is the assumption of the state-of-the-art approach, so the two modeling traditions are not incompatible at the level of the recovered object.","The authors point out that the same construction extends to sparsely observed or fragmented functions and to curves and surfaces in higher dimensions through the same norm-preserving action."],"supporting_citations":[{"why":"Defines the warpMix functional mixed model with value-preserving phase variation; it is the state-of-the-art baseline and the data-generating mechanism in Example 2.","marker":"Claeskens et al. [2021]"},{"why":"Supplies the group $\\Gamma$ of phase functions, the norm-preserving action, and the one-parameter family used as Prior Model 1.","marker":"Srivastava and Klassen [2016]"},{"why":"Provides the Dirichlet-process-based prior on phase functions used as Prior Model 2.","marker":"Bharath and Kurtek [2020]"},{"why":"Introduces the Fisher-Rao registration and orbit-centering approach that the paper's phase centering follows, and supplies the Berkeley growth-rate data.","marker":"Srivastava et al. [2011b]"},{"why":"Gives the landmark size-and-shape identifiability argument that motivates treating only the orbit of $\\mu$ as recoverable.","marker":"Lele and McCulloch [2002]"},{"why":"Establishes identifiability and rank conditions for registration with value-preserving warping, justifying the claim that pointwise recovery is impossible in general.","marker":"Chakraborty and Panaretos [2021]"},{"why":"Provides the Bayesian registration of curves under the square-root velocity transform that serves as the alternative Bayesian comparison (BRFC).","marker":"Cheng et al. [2016]"},{"why":"Details the Dirichlet-increment phase prior in Bayesian functional data settings and supports the extension to sparse and fragmented observations.","marker":"Matuk et al. [2022]"}],"fun_headline_variants":["Bayesian method nails a function's size-and-shape","Shape beats values when functions are misaligned","Timing warps no obstacle for Bayesian shape recovery","Function shape recovered despite timing differences","Bayesian twist: recover shape, skip the values"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole result depends on the assumption that the prior choices and the finite basis used to represent the mean keep time-warping from being confused with amplitude differences or noise; the paper demonstrates this by simulation but provides no proof.","fun_headline_variants_meta":{"raw":{"variants":["Bayesian method nails a function's size-and-shape","Shape beats values when functions are misaligned","Timing warps no obstacle for Bayesian shape recovery","Function shape recovered despite timing differences","Bayesian twist: recover shape, skip the values"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000303,"raw_usage":{"total_tokens":1810,"prompt_tokens":1080,"completion_tokens":730,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":696,"completion_tokens_details":{"reasoning_tokens":658}},"tokens_in":696,"tokens_out":730,"duration_ms":6842,"temperature":1.0,"reasoning_tokens":658,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T11:13:31.593695+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Simulate $n$ curves from $f_i=D_{\\gamma_i}(\\mu+v_i+\\epsilon_i)$ with phase functions drawn from a flexible family outside the two priors and with a true $\\mu$ that is not exactly representable in the chosen basis; compute $d=\\inf_{\\gamma\\in\\Gamma}\\|\\hat\\mu_{\\mathrm{center}}-D_\\gamma\\mu\\|$ for the centered posterior mean. If $d$ does not decrease as $n$ grows, the claim that the size-and-shape is recoverable is falsified.","supporting_citations":[],"review_version":1}