Justine Dushimirimana , Isaac Kipchirchir Chumba , Lydia Musiga , Joseph Nzabanita , Ronald Waliaula Wanyonyi
{"title":"Test for a general trilinear hypothesis in the generalized growth curve model","authors":"Justine Dushimirimana , Isaac Kipchirchir Chumba , Lydia Musiga , Joseph Nzabanita , Ronald Waliaula Wanyonyi","doi":"10.1016/j.jmva.2025.105470","DOIUrl":"10.1016/j.jmva.2025.105470","url":null,"abstract":"<div><div>In this paper, we consider the problem of testing a general trilinear hypothesis in the generalized growth curve model. The general trilinear hypothesis was formulated to test for example the significance of the generalized growth curves or the equality of the trilinear mean between groups in the two dimensions. The null hypothesis considered is of the form <span><math><mrow><mi>ℬ</mi><mspace></mspace><mo>×</mo><mrow><mo>{</mo><mi>L</mi><mo>,</mo><mi>M</mi><mo>,</mo><mi>N</mi><mo>}</mo></mrow><mo>=</mo><mi>O</mi></mrow></math></span>, where <span><math><mrow><mi>L</mi><mo>,</mo><mi>M</mi></mrow></math></span> and <span><math><mi>N</mi></math></span> are known matrices, <span><math><mi>ℬ</mi></math></span> is unknown parameter tensor and <span><math><mi>O</mi></math></span> is a tensor of zeros. The estimators of the parameters were obtained using a flip-flop algorithm under the null and alternative hypotheses. The likelihood ratio test for testing the general trilinear hypothesis was discussed. The proposed test is an extension of the likelihood ratio test for the general linear hypothesis under the growth curve model. A simulation study was performed to evaluate the performance of the proposed test and a real dataset was used for an illustrative example.</div></div>","PeriodicalId":16431,"journal":{"name":"Journal of Multivariate Analysis","volume":"210 ","pages":"Article 105470"},"PeriodicalIF":1.4,"publicationDate":"2025-07-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144662186","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Jiazhen Xu, Janice L. Scealy, Andrew T.A. Wood, Tao Zou
{"title":"Generalized score matching","authors":"Jiazhen Xu, Janice L. Scealy, Andrew T.A. Wood, Tao Zou","doi":"10.1016/j.jmva.2025.105473","DOIUrl":"10.1016/j.jmva.2025.105473","url":null,"abstract":"<div><div>Score matching is an estimation procedure that has been developed for statistical models whose probability density function or probability mass function is known up to proportionality but whose normalizing constant is intractable, so that maximum likelihood is difficult or impossible to implement. To date, applications of score matching have focused more on continuous IID models. Motivated by various data modeling problems, this article proposes a unified asymptotic theory of generalized score matching developed under the independence assumption, covering both continuous and discrete response data, thereby giving a sound basis for score-matching-based inference. Real data analyses and simulation studies provide convincing evidence of strong practical performance of the proposed methods.</div></div>","PeriodicalId":16431,"journal":{"name":"Journal of Multivariate Analysis","volume":"210 ","pages":"Article 105473"},"PeriodicalIF":1.4,"publicationDate":"2025-07-05","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144597334","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Moment-type estimators for the Dirichlet and the multivariate gamma distributions","authors":"Ioannis Oikonomidis, Samis Trevezas","doi":"10.1016/j.jmva.2025.105471","DOIUrl":"10.1016/j.jmva.2025.105471","url":null,"abstract":"<div><div>This study presents new closed-form estimators for the Dirichlet and the multivariate gamma distribution families, whose maximum likelihood estimator cannot be explicitly derived. The methodology builds upon the score-adjusted estimators for the beta and gamma distributions, extending their applicability to the Dirichlet and multivariate gamma distributions. Expressions for the asymptotic variance–covariance matrices are provided, demonstrating the superior performance of score-adjusted estimators over the traditional moment ones. Leveraging well-established connections between the Dirichlet and multivariate gamma distributions, a novel class of estimators for the latter is introduced, referred to as “Dirichlet-based moment-type estimators”. The general asymptotic variance–covariance matrix form for this estimator class is derived. To facilitate the application of these innovative estimators, an <span>R</span> package called <span>joker</span> is developed and made publicly available.</div></div>","PeriodicalId":16431,"journal":{"name":"Journal of Multivariate Analysis","volume":"210 ","pages":"Article 105471"},"PeriodicalIF":1.4,"publicationDate":"2025-07-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144571661","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Robust functional inverse regression","authors":"Haoyang Cheng , Jianjun Xu , Qian Huang","doi":"10.1016/j.jmva.2025.105472","DOIUrl":"10.1016/j.jmva.2025.105472","url":null,"abstract":"<div><div>Functional sufficient dimension reduction (FSDR) is a popular approach for supervised dimensionality reduction in regression settings, as it allows for the reduction of functional predictors to a lower-dimensional subspace without loss of information. However, most existing FSDR methods are vulnerable to heavy-tailedness or outliers, which are common in many real-world applications. To address this limitation, we propose a robust FSDR method that utilizes a functional pairwise spatial sign (PASS) operator. This approach is suitable for both completely observed functional data and sparsely observed longitudinal data. Our method provides a more robust approach to FSDR, by taking into account the spatial information of the data and assigning greater weights to the less outlier-prone observations. We also provide a convergence rate analysis of the estimator, demonstrating that our method yields a consistent estimate of the dimension reduction directions. The effectiveness of our proposed method is demonstrated through extensive simulations and real data analysis. Our method outperforms existing methods in terms of robustness and accuracy, making it a valuable tool for analyzing functional data across various applications.</div></div>","PeriodicalId":16431,"journal":{"name":"Journal of Multivariate Analysis","volume":"210 ","pages":"Article 105472"},"PeriodicalIF":1.4,"publicationDate":"2025-07-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144562999","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Robust signal recovery in Hadamard spaces","authors":"Georg Köstenberger, Thomas Stark","doi":"10.1016/j.jmva.2025.105469","DOIUrl":"10.1016/j.jmva.2025.105469","url":null,"abstract":"<div><div>We analyze the stability of (strong) laws of large numbers in Hadamard spaces with respect to distributional perturbations. For the inductive means of a sequence of independent but not necessarily identically distributed random variables, we provide a concentration inequality in quadratic mean and a strong law of large numbers, generalizing a classical result of K.-T. Sturm. For the Fréchet mean, we generalize H. Ziezold’s law of large numbers in Hadamard spaces. In this case, we neither require our data to be independent nor identically distributed; reasonably mild conditions on the first two moments of our sample are enough. Additionally, we look at data contamination via a model inspired by Huber’s <span><math><mi>ɛ</mi></math></span>-contamination model, in which we replace a random portion of the data with noise. In the most general setup, we neither require the data nor the noise to be i.i.d., nor do we require the noise to be independent of the data. A resampling scheme is introduced to analyze the stability of the (non-symmetric) inductive mean with respect to data loss, data permutation, and noise, and sufficient conditions for its convergence are provided. These results suggest that means in Hadamard spaces are as robust as those in Euclidean spaces. This is underlined by a small simulation study in which we compare the robustness of means on the manifold of positive definite matrices with means on open books.</div></div>","PeriodicalId":16431,"journal":{"name":"Journal of Multivariate Analysis","volume":"210 ","pages":"Article 105469"},"PeriodicalIF":1.4,"publicationDate":"2025-06-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144549931","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Robust factorization for high-dimensional matrix-variate observations","authors":"Yalin Wang , Long Yu","doi":"10.1016/j.jmva.2025.105467","DOIUrl":"10.1016/j.jmva.2025.105467","url":null,"abstract":"<div><div>Large-dimensional matrix-variate observations have been ubiquitous in the big data era, while unsupervised low-rank approximation technique would help reveal their hidden patterns and structures. In this paper, we study a hierarchical Canonical Polyadic (CP) product matrix factor model under the elliptical framework, which essentially assumes that the matrix-variate observations are from a matrix elliptical distribution. The proposed model not only incorporates the row-wise and column-wise interrelated information, but also adapts to the tail properties of the matrix-variate observations. We resort to the matrix Kendall’s tau introduced in the recent literature to recover the loading spaces, and minimize the square loss function to estimate the factor scores. We also propose an eigenvalue-ratio method to estimate the pair of factor numbers. Thorough theories for the model estimation, including statistical consistency and rates of convergence, are established under regular conditions. It is worth highlighting that the proposed method exhibits superior performance compared to other methods for estimating the signal part, particularly in the heavy-tailed cases. This superiority has been thoroughly validated through extensive simulations. The effectiveness in matrix reconstruction of the proposed method is demonstrated by applying it to a macroeconomic dataset of China.</div></div>","PeriodicalId":16431,"journal":{"name":"Journal of Multivariate Analysis","volume":"210 ","pages":"Article 105467"},"PeriodicalIF":1.4,"publicationDate":"2025-06-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144501685","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Error analysis for a statistical finite element method","authors":"Toni Karvonen , Fehmi Cirak , Mark Girolami","doi":"10.1016/j.jmva.2025.105468","DOIUrl":"10.1016/j.jmva.2025.105468","url":null,"abstract":"<div><div>The recently proposed statistical finite element (statFEM) approach synthesises measurement data with finite element models and allows for making predictions about the unknown true system response. We provide a probabilistic error analysis for a prototypical statFEM setup based on a Gaussian process prior under the assumption that the noisy measurement data are generated by a deterministic true system response function that satisfies a second-order elliptic partial differential equation for an unknown true source term. In certain cases, properties such as the smoothness of the source term may be misspecified by the Gaussian process model. The error estimates we derive are for the expectation with respect to the measurement noise of the <span><math><msup><mrow><mi>L</mi></mrow><mrow><mn>2</mn></mrow></msup></math></span>-norm of the difference between the true system response and the mean of the statFEM posterior. The estimates imply polynomial rates of convergence in the numbers of measurement points and finite element basis functions and depend on the Sobolev smoothness of the true source term and the Gaussian process model. A numerical example for Poisson’s equation is used to illustrate these theoretical results.</div></div>","PeriodicalId":16431,"journal":{"name":"Journal of Multivariate Analysis","volume":"210 ","pages":"Article 105468"},"PeriodicalIF":1.4,"publicationDate":"2025-06-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144492150","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Random projection-based response best-subset selector for ultra-high dimensional multivariate data","authors":"Jianhua Hu , Tao Li , Xiaoqian Liu , Xu Liu","doi":"10.1016/j.jmva.2025.105465","DOIUrl":"10.1016/j.jmva.2025.105465","url":null,"abstract":"<div><div>In this article, we propose a random projection-based response best-subset selector to perform response variable selection in ultra-high dimensional multivariate data, where both the dimensions of response and predictor variables are substantially greater than the sample size. This method is developed by integrating the response best-subset selector and random projection technique which is applied to reduce dimensionality of predictors. Under a multivariate tail eigenvalue condition, such a random projection-based dimensionality reduction of predictors only leads to an ignorable error between the original and dimension-reduced models. A computational procedure is presented. The proposed method exhibits model consistency under some certain conditions. The efficiency and merit of the proposed method are strongly supported by extensive finite-sample simulation studies. A real breast cancer dataset spanning 22 chromosomes are analyzed to demonstrate the proposed method.</div></div>","PeriodicalId":16431,"journal":{"name":"Journal of Multivariate Analysis","volume":"210 ","pages":"Article 105465"},"PeriodicalIF":1.4,"publicationDate":"2025-06-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144514351","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Additive regression for Riemannian functional responses","authors":"Jeong Min Jeon, Germain Van Bever","doi":"10.1016/j.jmva.2025.105466","DOIUrl":"10.1016/j.jmva.2025.105466","url":null,"abstract":"<div><div>We explore additive regression for a functional response whose values lie on a general Riemannian manifold. Due to the absence of a vector space structure on the manifold, we transform the response into a vector space. The transformation involves an unknown quantity that needs to be estimated. An additive model on the vector space is estimated using the smooth backfitting method, which effectively avoids the curse of dimensionality. We derive its rates of convergence and demonstrate its practical performance through a simulation study and a real data application.</div></div>","PeriodicalId":16431,"journal":{"name":"Journal of Multivariate Analysis","volume":"210 ","pages":"Article 105466"},"PeriodicalIF":1.4,"publicationDate":"2025-06-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144335660","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Deconvolution density estimation on Lie groups without auxiliary data","authors":"Jeong Min Jeon","doi":"10.1016/j.jmva.2025.105464","DOIUrl":"10.1016/j.jmva.2025.105464","url":null,"abstract":"<div><div>In this paper, we study density estimation on a general Lie group when data contain measurement errors and the distribution of measurement error is unknown. We estimate the target density without additional observations, such as an observable random sample from the measurement error distribution or repeated measurements. To achieve this, we take a semiparametric approach assuming that the measurement error distribution belongs to a parametric family. We also discuss maximum likelihood estimation for the case where the target density is also parametric. We establish the identifiability of a measurement error model and derive various asymptotic properties for our estimators. The performance of our estimators is demonstrated via simulation studies.</div></div>","PeriodicalId":16431,"journal":{"name":"Journal of Multivariate Analysis","volume":"209 ","pages":"Article 105464"},"PeriodicalIF":1.4,"publicationDate":"2025-06-11","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144290886","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}