We propose an omnibus nonparametric test for structural changes in the multivariate volatility matrix. The test aggregates bounded generalized quantile scores over a range of quantile levels and has a weighted leave-$q$-out $U$-statistic representation. Deleting nearby index pairs renders the centering effect induced by serial dependence asymptotically negligible. All quantities required for implementation, including the variance estimator used for standardization, are…
Sokbae Lee · Yuan Liao · Myung Hwan Seo · Youngki Shin·首次提交 2026-08-26 · 修订 2026-08-26 · v1·arXiv预印本
SAUSS:有限因变量模型的无偏模拟得分随机逼近AI 解读
Multinomial choice models allow flexible substitution patterns but become computationally demanding with many alternatives or observations. With a fixed per-observation simulation budget, simulated maximum likelihood introduces simulation bias, while each optimization step requires a full-sample likelihood evaluation. We propose Stochastic Approximation with Unbiased Simulated Scores (SAUSS), an averaged stochastic approximation based on conditionally unbiased mini-batch…
The graphical lasso is a widely used algorithm for fitting undirected Gaussian graphical models. However, for inference on functionals of edge values in the learned graph, standard tools lack formal statistical guarantees, such as control of the type I error rate. We introduce a selective inference method for asymptotically valid inference after graphical lasso selection with added randomization. We obtain a selective likelihood, conditional on the event of selection,…
This article surveys two deeply connected literatures that approach the same fundamental problem from different disciplinary traditions: dynamic discrete choice (DDC) in structural econometrics and inverse reinforcement learning (IRL) in machine learning. Both seek to infer the preferences of decision makers from observed sequential behavior, assuming that individuals act to maximize an expected reward function within a dynamic, uncertain environment formalized as a Markov…
We study a dynamic spatial panel model with observed regressors, interactive effects, and contemporaneous and lagged dependence in a large-$N$, fixed-$T$ framework. The spatial model constitutes an $N$-dimensional simultaneous-equations system. In this $N$-equation view, the interactive effects introduce $N$ unit-specific loading vectors. Estimating them individually when $T$ is fixed creates the type of incidental-parameters problem underlying Nickell bias. We use the $N$-eq…
David McCoy · Yi Li·首次提交 2026-08-20 · 修订 2026-08-20 · v1·arXiv预印本
靶向深度生存对比:神经网络下处理特定生存收益的有效推断AI 解读
Neural survival models are increasingly asked to support counterfactual claims---how much a treatment would change survival in a population---rather than only prognostic risk scores. Answering such questions from observational data requires valid inference for treatment-specific survival contrasts under confounding and covariate-dependent censoring, targets for which standard deep survival estimators are biased and provide no honest uncertainty. We propose Targeted Deep…
Jack Storror Carter·首次提交 2025-05-26 · 修订 2026-08-20 · v2·arXiv预印本
高斯精度矩阵的惩罚似然估计存在性与可分离先验的后验适定性AI 解读
Penalised likelihoods are often used for sparse estimation of a Gaussian precision matrix. In high dimensional settings where the matrix dimension is larger than the sample size, the sample covariance matrix $S$ is not of full rank and the maximum likelihood estimate of the precision matrix does not exist. An additional advantage of some penalised likelihood estimates, for example the graphical lasso, is that it can exist even in such high dimensional settings. This paper…
Wenxin Zhang · Rachael Phillips · Mark van der Laan·首次提交 2026-08-19 · 修订 2026-08-19 · v1·arXiv预印本
诊断测试与人工智能医疗器械的因果评估:效应修饰与信息增强决策AI 解读
Diagnostic medical tests and devices provide useful information for evaluating the potential benefits and risks of therapeutic treatments. However, unlike treatments, their impact on health outcomes is generally indirect because measuring diagnostic information typically does not itself affect patient outcomes, which complicates evaluation of their effectiveness. In this work, we develop a causal inference approach for evaluating diagnostic tests by distinguishing…
We propose a family of control variate estimators for variance reduction in design-based survey sampling and causal inference, with and without interference. In these settings, inverse probability weighting (IPW) estimators are widely used, but may have large variance when sampling, treatment, or exposure probabilities are small. Building on the observation that several common estimators, including the Hajek estimator, the normalized estimator, the augmented inverse…
We study finite-sample parameter estimation in logistic regression with Gaussian design, where the goal is to estimate $\mathbfθ^*\in \mathbb{R}^d$ with $R=\|\mathbfθ^*\|_2\ge 1$ from i.i.d. samples $\{(\mathbf{x}_i,y_i)\}_{i=1}^n,$ $\mathbf{x}_i \sim N(0,\mathbf{I}_d)$, $y_i\mid \mathbf{x}_i \sim \mathrm{Bernoulli}((1+\exp(-\mathbf{x}_i^\top \mathbfθ^*))^{-1})$. In this paper, we provide the first minimax optimal estimator, and improve on the best known finite-sample error…
Runjia Li · Victor B. Talisa · Chung-Chou H. Chang·首次提交 2024-07-25 · 修订 2026-08-12 · v2·arXiv预印本
竞争风险生存结局的亚组平均处理效应双重稳健目标估计AI 解读
In recent years, precision treatment strategies have gained significant attention in medical research, particularly for patient care. We propose a novel framework for estimating subgroup average treatment effects (ATE) in time-to-event data with competing risks, using ICU patients with sepsis as an illustrative example. Our approach, based on cumulative incidence functions and targeted maximum likelihood estimation (TMLE), achieves both asymptotic efficiency and double…
Existing panel beta regression estimators do not accommodate higher-order autoregressive dependence within a joint likelihood framework. We develop such an estimator: a joint conditional maximum-likelihood beta-AR(k) model extended from a single time series to panel data, and extend it to panel inference through two new covariance estimators - a generalization of Beck-Katz panel-corrected standard errors (PCSE) and a Driscoll-Kraay-style alternative (DK) - to account for the…
Classical estimators for ARIMA parameters (MLE, CSS, OLS) assume Gaussian innovations, an assumption frequently violated in financial and economic data exhibiting asymmetric distributions with heavy tails. We develop and validate the second-order polynomial maximization method (PMM2) for estimating ARIMA$(p,d,q)$ models with non-Gaussian innovations. PMM2 is a semiparametric technique that exploits higher-order moments and cumulants without requiring full distributional…
M. Ehsan Karim·首次提交 2026-08-08 · 修订 2026-08-08 · v1·arXiv预印本
通向同一答案的三条路径:协调 AIPW、TMLE 与双重机器学习AI 解读
Augmented inverse-probability weighting (AIPW), targeted maximum likelihood estimation (TMLE), and double/debiased machine learning (DML) are three routes to the same efficient influence function for the average treatment effect --- settled theory we treat as background. This tutorial's contribution is its worked, shared-nuisance reconciliation on real data: what a practitioner must actually match for the routes, and the software packages, to agree. Working the effect of…
This paper introduces Mixtures of Geodesic Factor Analyzers (MGFA) on Riemannian homogeneous spaces. MGFA uses a geodesic factor model within each mixture component, providing greater expressiveness than mixtures of Riemannian radial distributions and enabling clustering of manifold-valued data with anisotropic subpopulations. We establish root-$n$ consistency for the MGFA maximum likelihood estimator (MLE), thereby filling a theoretical gap for mixtures of Riemannian radial…
A central problem in the theory of empirical Bayes is to control the regret (excess risk) of a learned Bayes rule by the Hellinger distance between the estimated and true marginal densities. In the normal means model, the classical result of Jiang and Zhang (2009) achieves this only after regularizing the Bayes rule and incurs an extraneous cubic logarithmic factor through a delicate recursive argument. This paper introduces a new technique, based on polynomial approximation…
Daniel Paulin · Victor Elvira·首次提交 2026-08-06 · 修订 2026-08-06 · v1·arXiv预印本
可扩展的 VARMA 模型估计AI 解读
Vector autoregressive moving-average (VARMA) models have long been considered impractical beyond moderate dimensions: the likelihood is non-convex, the parametrization is identified only up to equivalence, and every evaluation costs a pass over the entire series. Yet their moving-average term captures with a few parameters what a pure autoregression matches only with many lags. We introduce an estimation framework that removes this computational barrier: each optimization…
I develop an estimation and inference framework for distribution regression in dyadic network settings with two-way fixed effects that vary across thresholds of the outcome. I show that identification of the structural parameters is achieved through binarization of the outcome at each threshold, and estimate the model by conditional maximum likelihood, which "differences out" the fixed effects and circumvents the incidental parameter problem. The estimator remains…
We developed a unified covariate-adjusted causal inference framework for estimating the desirability of outcome ranking (DOOR) probability for benefit-risk evaluation in randomized trials and observational studies. The framework expresses the DOOR probability as a bilinear functional of the marginal ordinal outcome distributions under the two treatment strategies, estimates conditional ordinal distributions through sequential risk-set hazards, and derives the efficient…
Martin Chak · Zoraida F. Rico · Giacomo Zanella·首次提交 2026-05-29 · 修订 2026-08-04 · v2·arXiv预印本
Probit 回归中的强对数凹性AI 解读
We show that strong log-concavity emerges in probit regression likelihoods without ridge penalization (i.e. Gaussian priors), unlike for the logistic case. Specifically, we provide: (a) a non-asymptotic characterization of strong log-concavity for fixed designs, similar to that for the existence of the maximum likelihood estimator (MLE) and (b) an asymptotic analysis for random designs, in the proportional regime when the sample size $n$ and the number of covariates $d$ grow…
We consider statistical inference for a class of continuous semimartingale regression models based on high-frequency observations subject to contamination by finite-activity jumps and spike noise. By employing density-power weighting and Hölder-inequality-based normalization, we propose easy-to-implement, robustified versions of the conventional Gaussian quasi-maximum-likelihood estimator that require only a single tuning parameter. We prove their asymptotic mixed normality…
Empirical Bayes provides a powerful framework for learning and adapting to latent structure in data. In sequence models where an independent observation is associated to each latent parameter, theory and methods around empirical Bayes are well-developed. However, in models where latent parameters and observed data interact through more complex designs, many statistical and algorithmic questions remain unanswered. In this work, we study a canonical setting of empirical Bayes…
The long-term relationship between radiative forcing and surface temperature is central to predicting the impacts of climate change. This study employs multicointegration to characterize this relationship and provides the first application of the Transformed and Augmented Ordinary Least Squares (TAOLS) estimator to estimate the multicointegration model. The main objective is to estimate the Equilibrium Climate Sensitivity (ECS), defined as the global mean surface temperature…
Sohyeon Kim · Luo Xiao · Xinming An · Wenbin Lu·首次提交 2026-08-03 · 修订 2026-08-03 · v1·arXiv预印本
相关分类结果的多变量多项Logit模型与ANOVA分解AI 解读
In medical research, patients often have multiple interdependent outcomes, such as posttraumatic stress disorder (PTSD), depression, and pain among trauma survivors. Most existing research uses multinomial regression to analyze these interdependent outcomes separately, which ignores correlations between concurrent conditions. This omission may lead to loss of information and reduced predictive accuracy. Accounting for correlations between multiple categorical outcomes…
We propose an omnibus nonparametric test for structural changes in the multivariate volatility matrix. The test aggregates bounded generalized quantile scores over a range of quantile levels and has a weighted leave-q-out U-statistic representation. Deleting nearby index pairs renders the centering effect induced by serial dependence asymptotically negligible. All quantities required for implementation, including the variance estimator used for standardization, are…
Fixed-T Dynamic Spatial Panel Model with Common Shocks
We study a dynamic spatial panel model with observed regressors, interactive effects, and contemporaneous and lagged dependence in a large-N, fixed-T framework. The spatial model constitutes an N-dimensional simultaneous-equations system. In this N-equation view, the interactive effects introduce N unit-specific loading vectors. Estimating them individually when T is fixed creates the type of incidental-parameters problem underlying Nickell bias. We use the N-eq…
Penalised likelihoods are often used for sparse estimation of a Gaussian precision matrix. In high dimensional settings where the matrix dimension is larger than the sample size, the sample covariance matrix S is not of full rank and the maximum likelihood estimate of the precision matrix does not exist. An additional advantage of some penalised likelihood estimates, for example the graphical lasso, is that it can exist even in such high dimensional settings. This paper…
We study finite-sample parameter estimation in logistic regression with Gaussian design, where the goal is to estimate θ∗∈Rd with R=∥θ∗∥2≥1 from i.i.d. samples {(xi,yi)}i=1n,xi∼N(0,Id), yi∣xi∼Bernoulli((1+exp(−xi⊤θ∗))−1). In this paper, we provide the first minimax optimal estimator, and improve on the best known finite-sample error…
Classical estimators for ARIMA parameters (MLE, CSS, OLS) assume Gaussian innovations, an assumption frequently violated in financial and economic data exhibiting asymmetric distributions with heavy tails. We develop and validate the second-order polynomial maximization method (PMM2) for estimating ARIMA(p,d,q) models with non-Gaussian innovations. PMM2 is a semiparametric technique that exploits higher-order moments and cumulants without requiring full distributional…
This paper introduces Mixtures of Geodesic Factor Analyzers (MGFA) on Riemannian homogeneous spaces. MGFA uses a geodesic factor model within each mixture component, providing greater expressiveness than mixtures of Riemannian radial distributions and enabling clustering of manifold-valued data with anisotropic subpopulations. We establish root-n consistency for the MGFA maximum likelihood estimator (MLE), thereby filling a theoretical gap for mixtures of Riemannian radial…
We show that strong log-concavity emerges in probit regression likelihoods without ridge penalization (i.e. Gaussian priors), unlike for the logistic case. Specifically, we provide: (a) a non-asymptotic characterization of strong log-concavity for fixed designs, similar to that for the existence of the maximum likelihood estimator (MLE) and (b) an asymptotic analysis for random designs, in the proportional regime when the sample size n and the number of covariates d grow…