2026-11-06 | | Total: 54
Pocock and Simon's minimization method is a widely used approach for covariate-adaptive randomization in clinical trials. However, it can exhibit the shift problem under unequal allocation and variance inflation of the imbalances in additional covariates. In this paper, we propose a new minimization method that preserves the conventional biased-coin allocation probability while ensuring marginal imbalances bounded in probability. Meanwhile, the proposed procedure guarantees that neither the shift problem nor variance inflation occurs. Simulation and empirical studies demonstrate the excellent performance of the proposed procedure.
At the edge of stability, the iterates of gradient descent oscillate along the sharp directions of the Hessian. This poses a significant challenge for continuous-time models as the usual continuous-time approximation -- the gradient flow -- admits no oscillations. Several continuous-time models have been proposed to capture the behavior at the edge of stability, among them central flow and rod flow. Although both models can be interpreted as performing a coarse-graining in time, the connection between them remains poorly understood. Here we show that both flows can be derived from a third object, the edge coupling: a functional on consecutive iterates whose critical points are the fixed points and period-two orbits of gradient descent. It is revealed that the two flows correspond to coarse-grainings at different time scales. Rod flow coarse-grains over only the fastest time scale, while central flow coarse-grains over both the fast and intermediate time scales. Furthermore, using the edge coupling, we derive rod flows and central flows for momentum methods such as heavy ball and Adam. We verify the derived flows experimentally on representative neural network architectures, where they track the discrete iterates throughout all phases of training -- including the transition from the stable regime to the edge of stability.
We consider the problem of cross-validation (CV) in time series forecasting, with a focus on global forecasting models based on deep learning. The prevailing CV strategy in forecasting is to partition data across time with the last part of each series used for testing. With global models trained on collections of time series, however, a fundamentally different approach becomes possible: partitioning across series and time. We formalize series-wise CV for global forecasting models and benchmark several instantiations, including K-fold CV, bootstrapping, and Monte Carlo CV, on 11 datasets with 10 neural network architectures. We evaluate CV procedures along five dimensions: performance estimation, overoptimism rate, model selection accuracy, regret, and downstream forecasting performance. We withhold a random fraction of series from training so that they appear only in the test set. This protocol simulates the emergence of new series not seen during training, a common scenario in global forecasting. Under this protocol, our results reveal a nuanced trade-off: series-wise CV yields better downstream forecasting accuracy by guiding model selection toward configurations that generalize across both time and series, while time-wise CV retains a slight advantage for estimating out-of-sample error.
Two-phase designs provide a cost-effective framework for biomarker studies and have been widely used to reduce the burden of collecting expensive biomarker measurements. Most methodological research has focused on efficient Phase II sampling and estimation of biomarker--outcome associations. A natural subsequent step, however, is to incorporate identified biomarkers into risk prediction models, where measures of predictive performance become important inferential targets. Even for binary outcomes, methodological developments for two-phase designs that jointly address association parameters and predictive performance remain limited. Building on classical case-control and matched case-control designs developed primarily for association inference with binary outcomes, we propose a risk-matched case-control (RMCC) design that uses a scalar disease-risk score constructed from Phase I variables to define matching strata and guide Phase II sampling. To address the limitations of existing likelihood-based approaches, which typically require correct specification of the nuisance distribution of biomarkers, we develop a unified design-based inference framework based on inverse probability weighting for both association parameters and predictive performance measures, including positive predictive value and the area under the receiver operating characteristic curve. We further introduce target-specific weight scaling and calibration procedures and derive variance estimators that account for sampling variability in both phases. Simulation studies demonstrate the efficiency and robustness of the proposed design and inference procedures, and an application to a Taiwanese oral cancer cohort illustrates their potential utility.
Functional clustering identifies coordinated molecular profiles, but repeated numerical covariance optimization can make mixture-model fitting computationally demanding. We develop a fast functional clustering procedure for joint analysis across conditions, combining linear mean profiles on a standardized log niche coordinate with first-order structured antedependence (SAD(1)) covariance. Weighted residual statistics yield closed-form conditional covariance updates under four parameter-sharing specifications, while an innovation representation simplifies likelihood evaluation. In same-R comparisons with numerical optimization of identical conditional objectives, the closed-form solver was 91.5--201.4 times as fast, with objective agreement in all 60 targets. Under matched complete mixture fits, it achieved a 1.87--2.37-fold speedup with identical cluster assignments. Recovery simulations showed higher module classification accuracy and lower mean-curve error than $K$-means and a Gaussian mixture with shared full covariance when the mean structure was correctly specified. SAD(1) provided an additional recovery gain under strong residual correlation, whereas curved-mean controls identified the importance of an appropriate mean family. Applications to single-condition gut microbiota, paired leaf--root microbiota and three-visit human gut multi-omics illustrate the framework across different numbers of conditions. Closed-form covariance estimation makes joint functional clustering substantially faster while retaining the fitted model's likelihood and assignments.
Attention-based time series forecasters commonly group consecutive observations into patches. Patch length is known to have a substantial impact on forecasting performance, yet its selection remains largely empirical. We develop a theoretical framework that explains why patch length matters for forecasting accuracy and provides a principled basis for its selection. For linear state-space models, we derive an exact factorization of the forecasting coefficient matrix that separates forecast dynamics from filtering dynamics. We show that dominance of a single forecast component at long horizons does not generally imply a rank-one representation. This explains why longer horizons alone need not reduce representational complexity. We combine the factorization with estimation bounds to develop a patch-selection criterion that balances representation and estimation errors under temporal dependence. Our analysis establishes conditions under which balanced patching is justified and clarifies how predictive structure, head count, training sample size, and long-run dependence enter the selection criterion.
We propose a novel class of Bounded GARCH (B-GARCH) processes designed to characterise the conditional heteroskedasticity of continuous time series constrained within a compact interval. While standard GARCH specifications are tailored for R-valued data, many empirical series, such as realised correlation, unemployment rates, expenditure shares, and relative humidity, are bounded. Fitting a traditional GARCH model to such data creates a severe support mismatch, frequently generating predictions and density forecasts that breach the natural boundaries of the data (boundary violation). Furthermore, conventional approaches utilising non-linear link transformations onto the real line might induce interpretability issues and estimation biases. To address these limitations, our proposed B-GARCH framework ensures a proper predictive distribution and an equation for the dynamic conditional variance with a parameter constraint that strictly preserves the bounded support of the entire process almost surely. We establish the fundamental probabilistic properties of the process, proving strict stationarity and ergodicity, and address parameter estimation via conditional maximum likelihood. The asymptotic consistency and normality of the estimators are rigorously derived, and their finite-sample performance is validated through Monte Carlo simulations. An empirical application to the monthly realised correlation between the Vanguard Real Estate ETF (VNQ) and the iShares 20+ Year Treasury Bond ETF (TLT) demonstrates that B-GARCH models effectively capture volatility clustering while completely eliminating boundary violation issue that plagues traditional GARCH models.
State-space models (SSMs) recover low-dimensional latent dynamics from high-dimensional observations, but latents are generally not identifiable from observational data alone. Interventional state-space models (iSSMs) use perturbations to break observational symmetries and identify parameters up to permutation, scaling, and shift. However, prior work restricts identifiability to linear latent dynamics and piecewise-linear emissions, excluding many behaviors characteristic of neural circuits. We extend iSSMs to a substantially broader model class: we prove that (i) affine identifiability holds under general injective emission functions and Gaussian latent transitions, and (ii) with interventions, identifiability up to permutation, scaling, and shift is preserved under general nonlinear latent dynamics. We validate our theory on simulated FitzHugh-Nagumo dynamics, and then demonstrate two applications: adjudicating between competing models of decision-making; and predicting responses to unseen micro-stimulation in macaque and human neural recordings, with proof-of-concept generalization to a held-out stimulation site in the latter. Finally, simulations yield practical guidance for designing stimulation protocols that improve identifiability.
We propose a doubly-robust method to estimate the causal effects of multi-valued treatment regimes on patient outcomes. A multi-valued treatment regime assigns personalized treatments from a set of competing options based on a patient's history. Existing methods for optimal treatment regimes are largely limited to binary or continuous treatments. We extend the longitudinal targeted minimum loss-based estimator (LTMLE) to accommodate multi-valued treatments by modeling treatment assignments with a multinomial model. We further integrate recurrent neural networks (RNNs) into the LTMLE framework as learners for the outcome, treatment, and censoring models, and compare this approach to ensemble machine learning. Simulation studies show that LTMLE with either ensemble machine learning or RNNs is approximately unbiased, with close to nominal coverage, whereas sequential g-computation with ensemble machine learning is substantially biased under a dynamic rule. A multinomial treatment model performs similarly to a series of binomial regressions with ensemble machine learning, but yields more stable inference when the treatment model is an RNN.
Let $\mathbf X_1,\ldots,\mathbf X_n$ be independent random vectors in $\mathbb R^d$ with common mean $\boldsymbolμ_n$, where $d\ge2$ is a fixed integer. Let $\bar{\mathbf X}_n=n^{-1}\sum_{i=1}^n\mathbf X_i$ denote the sample mean and let \[ \widehatΣ_n=\frac1n\sum_{i=1}^n (\mathbf X_i-\bar{\mathbf X}_n) (\mathbf X_i-\bar{\mathbf X}_n)^\top \] denote the sample covariance matrix. Assume that the average covariance matrix is positive definite and that $\widehatΣ_n$ is nonsingular almost surely. For testing $H_0:\boldsymbolμ_n=\boldsymbolμ_0$, the Hotelling statistic is \[ T^2_{H,n} =n(\bar{\mathbf X}_n-\boldsymbolμ_0)^\top \widehatΣ_n^{-1} (\bar{\mathbf X}_n-\boldsymbolμ_0). \] Under $H_0$, we establish finite-sample Cramér-type moderate deviation results and uniform Berry--Esseen bounds under finite $p$th moments for $2<p\le3$. The moderate deviation error order is optimal over independent triangular arrays, as shown by explicit constructions. For identically distributed observations with a fixed distribution, we obtain \[ \sup_{0\le x\le b_n} \left| \frac{\mathbb P(T^2_{H,n}\ge x^2)} {\mathbb P(χ_d^2\ge x^2)}-1 \right| = O\!\left((1+b_n)^p n^{-(p-2)/2}\right) \] for every nonnegative sequence $b_n=o(n^{(p-2)/(2p)})$. The uniform Berry--Esseen bound is $O(n^{-(p-2)/2})$. The proof uses a radial Stein identity and a recursive error estimate. We apply the results to chi-squared calibration in simultaneous testing.
Extreme value theory provides a natural framework for modelling rare events and extrapolating tail probabilities, but its direct application to survival analysis is complicated by right censoring, heterogeneous follow-up, and the need for coherent inference beyond the range of the observed data. We propose a Bayesian survival model that integrates extreme value theory with a hazard-based formulation built upon the extended generalized Pareto distribution. Rather than separating the bulk and tail of the distribution through an explicit threshold, the proposed model represents the entire hazard function within a unified parametric framework, allowing moderate and extreme survival times to be analysed simultaneously while preserving the asymptotic tail behaviour implied by extreme value theory. Inference is performed through a two-stage procedure. We further develop bias-correction procedures based on latent representative locations and adaptive partition refinement, and establish theoretical results showing the asymptotic consistency of the proposed approximation under standard regularity conditions. Simulation studies demonstrate that the methodology provides reliable uncertainty quantification across a wide range of tail regimes, censoring proportions, and sample sizes, while remaining robust under moderate model misspecification. The practical performance of the approach is illustrated through an analysis of a heart failure cohort, where the proposed model enables inference on long-term survival probabilities and tail behaviour beyond the observed follow-up period. The resulting framework offers a computationally feasible and theoretically grounded alternative for modelling censored survival data when inference on extreme survival times is of primary scientific interest.
Spearman's rho treats all ranks alike, yet in many comparisons of two rankings an agreement at the top matters more than one at the bottom. Blest's rank correlation builds this preference into Spearman's rho, and its symmetrized version $η$, introduced by Genest and Plante, does so without favoring either of the two rankings. How far can such a top-weighted coefficient stray from other rank correlations? We answer this question for Spearman's rho, Blomqvist's beta and Chatterjee's rank correlation $ξ$ by determining the exact set of values that $η$ and each of these coefficients can attain jointly. Given Spearman's rho, $η$ can differ from it by at most $27/128$, and the extreme values of $η$ are attained only when the two rankings agree perfectly on the top ranks and are completely reversed on the bottom ranks, or vice versa. Given Blomqvist's beta, $η$ obeys exactly the same sharp bounds as Spearman's rho. For Chatterjee's rank correlation, both boundary curves are explicit and are attained by copulas with densities. Here symmetrization has a surprising effect. Compared with Blest's original coefficient, $η$ has a narrower range at fixed Spearman's rho, but near independence a wider range at fixed Chatterjee's rank correlation, where $|η|$ can come close to $\sqrt{11ξ/5}$ while Blest's coefficient stays below $\sqrt{32ξ/15}$. The proofs combine explicit Kantorovich potentials from optimal transport with a convex relaxation that turns out to be exact.
We record two observations concerning \cite{HLW24}. First, the proof of its Theorem~1 yields the population linearization (13)--(14) at \textit{any critical point} of the Fréchet function $F_μ$, under its remaining geometric and measure-theoretic assumptions; neither isolation nor a nonsingular Hessian of $F_μ$ is needed. Second, its Theorem~2 is given a shorter proof by applying the M-estimator theorem \cite[Theorem~5.23, pp.~53--54]{vdV98} to negative squared distance in normal coordinates, as opposed to the Donsker argument used there. We verify the pointwise differentiability, local square-integrable Lipschitz bound, second order expansion of population loss, consistency, and empirical maximization hypotheses. The resulting covariance retains the cut-locus correction. The same argument applies to suitably selected local sample Fréchet minima when the population Hessian is positive definite.
Functional data analysis provides a framework for studying time-varying associations between covariates and high-frequency biological signals. In digital health, participants often contribute repeated functional trajectories, yet existing multilevel regression methods largely focus on the conditional mean, overlooking variability and its role in biological heterogeneity. We propose a heteroskedastic functional mixed-effects regression framework that jointly models the conditional mean and variance of repeated trajectories, with participant-specific random effects in both components. We also develop a participant-level bootstrap to quantify uncertainty in estimated functional parameters. We apply the method to postprandial glucose trajectories from the AEGIS study. Premeal glucose, carbohydrate intake, and meal type show distinct time-dependent associations with glucose variability, whereas HbA1c is more strongly associated with the conditional mean than with variability. Substantial between-participant differences in residual variability remain even after accounting for observed covariates. These findings show that modeling conditional variability alongside the mean provides insights into glucose regulation and individual heterogeneity that mean-focused analyses may miss. More broadly, the proposed framework enables joint characterization of average responses, variability, and participant-specific heterogeneity in repeated biological signals.
Continuous diagnostic tests are often used to discriminate between diseased and healthy populations. Therefore, it is important to look for measures to assess the effectiveness of diagnostic markers in discriminating between these two states. The overlap coefficient (OVL), which is defined as the area of intersection of the graphs of two probability density functions, is considered as a summary measure of diagnostic accuracy. Often additional information is available on some of the covariates which are known to influence the accuracy of such measures. Consequently, in this work we incorporate covariate information using linear regression models with normally distributed errors, enabling the construction of the OVL estimator based on the ordinary least squares estimates of the regression coefficients. We further develop a method for constructing generalized confidence intervals for the covariate-specific overlap coefficient. A simulation study is conducted to assess the finite-sample performance of the proposed generalized confidence interval method and to compare it with other procedures. Finally, an illustration to assess the age-specific accuracy of glucose as a biomarker of diabetes is provided.
A provider holding a watermark key may audit an external document and choose passages after reading it. A bound computed for a fixed passage is not valid for a passage chosen after seeing the evidence. One remedy splits the evidence between selection and validation. We take the other route and make the statement hold for every interval at once. CountScan gives finite-sample lower bounds on the number of retained historical watermark matches, simultaneously over all text intervals and count thresholds. Its validity rests on one condition: given the auditor's information, the pivots of positions that are not history matches are conditionally super-uniform. This condition holds when generation uses the key only through the Gumbel-max rule and editors have no access to the key, so the guarantee covers edits that depend on the generated text when the calibrators are nonincreasing. Simultaneity raises the threshold by the logarithm of the document length; for a long interval this costs about that logarithm divided by the critical evidence value, in certified tokens. We give a finite-length expansion of the certified proportion, finite-sample rates for two-level rules with lower bounds of the same order within a pivot-only class, and an upper bound over a calibrator class. On pre-registered holdouts with OPT and Qwen base models at temperature one, CountScan certifies substantial counts; under heavy random edits it remains valid but loses most of its power.
Estimating causal effects from observational network data requires accounting for both individual treatment and neighborhood exposure. Existing outcome models commonly incorporate these intervention variables into a joint predictor and estimate causal effects through prediction contrasts, leaving the corresponding response components implicit. We propose the $\underline{N}$etwork $\underline{E}$ffect $\underline{D}$ecomposition $\underline{Net}$work (NEDNet), which parameterizes the conditional outcome response as a query-invariant pre-treatment baseline plus separate treatment- and exposure-response pathways. Path-specific neighborhood contextualization supplies local pre-treatment information to the two pathways. We characterize the fitted response through a unique reference-normalized decomposition, which is equivalent to a zero mixed treatment--exposure contrast. Experiments on semi-synthetic benchmarks derived from BlogCatalog and Flickr show that NEDNet achieves the lowest PEHE among the compared methods for individual, peer, and total effects across homogeneous and covariate-heterogeneous response settings.
This study presents a methodology for detecting change point in precipitation series using a smooth version of segmented regression model. We examine several existing change point detection methods and identify their limitations in capturing gradual changes commonly observed in climate data. To address these limitations, we discuss an extension of the traditional segmented regression model enabling the detection of gradual shifts in trend. The developed methodology is further extended to accommodate dependent error structures, which are common in time series data. Extensive simulation studies are conducted to evaluate the performance and computational efficiency of the proposed approach relative to existing methods. The methodology is subsequently applied to 60 years of monsoon rainfall data from two districts of West Bengal, India. This approach offers improved insights into evolving rainfall patterns and serves as a useful tool for climate data analysis and long-term forecasting.
Electric robotaxis are rapidly expanding from pilots to scaled commercial deployment; unlike traditional human-driven ride-hailing vehicles, they require centralized fleet control and careful management of vehicle state of charge (SoC). We study how a fleet manager should jointly make service, empty repositioning, idling, and charging decisions in such systems. We develop a continuous-time fluid model in which vehicles evolve across location-charge states, and identify two forms of charge-aware control that are absent from classic repositioning models: early return-to-depot (RTD) and charge-aware repositioning. In particular, the optimal solution admits simple threshold structures. Each location has a battery threshold such that vehicles should return to the depot to recharge when their SoC is below it and should not do so when their SoC is above it. Vehicles should not reposition below the same threshold either; above it, the choice of destination depends on whether a vehicle's SoC matches the destination's battery sensitivity, a notion we introduce in the paper. We then quantify the value of charge awareness by comparing fully charge-aware policies with charge-oblivious policies that ignore SoC in RTD and repositioning decisions. We show that the performance gap can be unbounded, driven by a dilemma between under-utilizing high-charge vehicles and over-sending low-charge vehicles to locations from which charging is costly. A numerical case study calibrated using 2025 NYC Taxi and Limousine Commission trip records illustrates these thresholds and shows that charge-aware control can generate substantial gains in a realistic urban network. Importantly, simple threshold-based heuristics recover most of the value of the fully charge-aware solution.
Wastewater surveillance provides population-level evidence of infectious disease activity. As surveillance expands across multiple sites, repeated measurements must be analyzed jointly while accounting for differences among sites and the public health regions they represent. We develop a Bayesian spatio-temporal framework that decomposes wastewater measurements into a shared epidemic curve, spatially correlated regional deviations, persistent site effects, and within-site temporal dependence. This hierarchical decomposition accommodates unequal observation windows and sampling schedules, allowing regional trajectories to be compared while preserving site-level heterogeneity. A direct marginal spatial covariance specification further separates the magnitude of regional variation from the strength of association between neighboring regions, providing distinct interpretations of these two sources of spatial variation. We also introduce region-specific time shifts to distinguish regional leads and lags in epidemic timing from differences in signal levels and local fluctuations. Simulation studies demonstrate the recovery of these components and show that omitting true time shifts can lead to increased estimates of residual temporal variation. Application to Ontario SARS-CoV-2 wastewater surveillance data reveals persistent site heterogeneity, within-site temporal dependence, and spatial association among neighboring regions. Regional timing differences are more pronounced during the first year of Omicron than in the XBB-and-later period.
We study the cost of certifying a scalar target driven by a product of centered Bernoulli means when an unobserved terminal gate leaves irreducible uncertainty. Each factor can be queried independently at a price. We derive the exact identified interval and determine when a requested accuracy is feasible. At every fixed strictly feasible interior instance, the optimal high-confidence expected cost is $(1+o(1))\log(1/δ)/D^*$, where $D^*$ is an answer-and-allocation information game; a globally honest procedure attains this benchmark. At product cancellation, we solve the game through a one-dimensional inverse-KL equation and explicit cost weights. For two tickets, quartic projections and rational bounds give checkable global optimality gaps. Joint confidence--precision limits expose a failure of the analyzed fixed pilot: it can cost asymptotically at least twice the optimum. A charged, precision-aware pilot with checkpoint certification restores optimal expected cost in exact arithmetic when $\log\{e+\log(1/g)\}=o(\log(1/δ))$ for effective precision $g$, including exponential and certain double-exponential precision scales. Pooled rescue preserves validity and cannot increase observation cost on coupled paths. Implemented auditors verify returns by rational inequalities; reproducible experiments show both savings and finite-confidence overhead. These guarantees concern the specified independent-reference model and its scalar target.
The joint expectation of functions of potential outcomes $(Y_1,Y_0)$ is a fundamental quantity in causal inference, particularly when investigating the variation or heterogeneity of causal effects. This paper provides novel assumptions for their identification: comonotonicity and countercomonotonicity. These assumptions offer a unified framework for their identification across discrete, continuous, and mixed outcome variables. In the absence of these assumptions, we derive sharp bounds for two broad classes of functions, yielding new results for bounding moments of individual causal effects, which are central to measuring effect heterogeneity. We present corresponding estimation methods and illustrate them in simulation studies and real-world datasets.
An experiment can be informative about the parameters of a working model and uninformative about whether that model is wrong. Adaptive design can reinforce this mismatch by repeatedly selecting settings where an incorrect model fits. We bound the power of any valid adequacy test in terms of the experimenter's own sampling distribution, then distinguish information absent from the selected data from information lost by the detector. A computable diagnostic measures unresolved error relative to a declared expansion of the working model. It motivates a training-free rule that reserves part of the design budget for falsification. The empirical record is mixed: the independently specified primary study fails detection eligibility, and a later numerator diagnostic fails its improvement criteria. After the frozen rule fails on fresh continuous-null panels, a disclosed secondary horizon-calibrated revision detects 35.8% of misspecified campaigns versus 31.7% for a class-separation control at T = 240, but has no resolved advantage at T = 480. In a corrected 289-unit measured replay it improves on random design and trails class-discrimination EIG. A conditional futility bound explains why remaining time matters without proving the gate constant or its power. These results support separating parameter learning, model checking and detector performance; they do not establish a universally superior design policy or prospective scientific discovery.
Renewal processes model inter-event times between recurrent events using a parametric distribution, but selecting a single model can ignore model uncertainty and affect inference, particularly in small samples. This paper introduces marp, an R package for model-averaged inference in parametric renewal processes. The package fits six commonly used distributions, compares them using Akaike and Bayesian information criteria, and reports selected-model and model-averaged estimates of the mean inter-event time, event probability and hazard function, with percentile and Studentised bootstrap confidence intervals. A simulation study under three data-generating models showed that model averaging did not improve point estimation but often improved interval coverage relative to the AIC-selected model. Applications to seizure-recurrence and engineering failure-time data illustrate the proposed methods. The marp package provides a reproducible framework for inference that accounts for model uncertainty in parametric renewal processes.
In this article, we study deep nonparametric regression trained via stochastic gradient Langevin dynamics (SGLD), a variant of stochastic gradient descent that injects carefully calibrated Gaussian noise during optimization. Using a PAC-Bayesian analysis tailored to deep ReLU networks, we establish an excess-risk bound that simultaneously accounts for network approximation error, statistical complexity, and the algorithmic cost accumulated along the SGLD trajectory. Under a Bernstein condition and within a stable SGLD regime, the excess-risk bound contains a stochastic term whose coefficient is inversely proportional to the variance of injected noise, while the estimator retains the optimal rate for nonparametric regression under the Hölder smoothness condition up to logarithmic factors. This result provides an approximation-error-aware, rate-optimal statistical perspective on the widely observed generalization benefits of moderate SGLD noise. For per-sample clipped SGLD, we further establish an exact f-DP guarantee and the corresponding asymptotic $μ$-GDP guarantee for the proposed estimator. We establish statistical convergence rates and differential privacy guarantees for DNN estimators trained using SGLD, while explicitly accounting for both network approximation error and the effect of injected noise in the excess-risk bound. Extensive simulation studies and analyses of two real datasets examine the finite-sample predictive performance and privacy-utility behavior of the proposed procedure.