The heteroscedastic multivariate linear model with multivariate normal error distribution has been considered. Using the structural relation of the model, the prediction distribution of future responses of the model h...
详细信息
The heteroscedastic multivariate linear model with multivariate normal error distribution has been considered. Using the structural relation of the model, the prediction distribution of future responses of the model has been derived. It is observed that for known covariance parameters the prediction distribution of the model has a product of m multivariate Student t distribution. It is to be noted that the prediction distribution for the Student terror also has a product of m multivariate Student t distribution. Some special cases have been discussed.
In this paper, we propose a novel variable selection approach in the framework of multivariate linear models taking into account the dependence that may exist between the responses. It consists in estimating beforehan...
详细信息
In this paper, we propose a novel variable selection approach in the framework of multivariate linear models taking into account the dependence that may exist between the responses. It consists in estimating beforehand the covariance matrix. of the responses and to plug this estimator in a Lasso criterion, in order to obtain a sparse estimator of the coefficient matrix. The properties of our approach are investigated both from a theoretical and a numerical point of view. More precisely, we give general conditions that the estimators of the covariance matrix and its inverse have to satisfy in order to recover the positions of the null and non null entries of the coefficient matrix when the size of Z is not fixed and can tend to infinity. We prove that these conditions are satisfied in the particular case of some Toeplitz matrices. Our approach is implemented in the R package MultiVarSel available from the Comprehensive R Archive Network (CRAN) and is very attractive since it benefits from a low computational load. We also assess the performance of our methodology using synthetic data and compare it with alternative approaches. Our numerical experiments show that including the estimation of the covariance matrix in the Lasso criterion dramatically improves the variable selection performance in many cases. (C) 2018 Elsevier Inc. All rights reserved.
Although there are many imaging studies on traditional ROI-based amygdala volumetry, there are very few studies on modeling amygdala shape variations This paper presents a unified computational and statistical framewo...
详细信息
Although there are many imaging studies on traditional ROI-based amygdala volumetry, there are very few studies on modeling amygdala shape variations This paper presents a unified computational and statistical framework for modeling amygdala shape variations in a clinical population The weighted spherical harmonic representation is used to parameterize, smooth out, and normalize amygdala surfaces The representation is subsequently used as an input for multivariate linear models accounting for nuisance covariates such as age and brain size difference using the SurfStat package that completely avoids the complexity of specifying design matrices. The methodology has been applied for quantifying abnormal local amygdala shape variations in 22 high functioning autistic subjects (C) 2010 Elsevier Inc All rights reserved
Hypothesis-error (or "HE") plots, introduced by Friendly (J Stat Softw 17(6):1-42, 2006a;J Comput Graph Stat 16:421-444, 2006b), permit the visualization of hypothesis tests in multivariate linear models by ...
详细信息
Hypothesis-error (or "HE") plots, introduced by Friendly (J Stat Softw 17(6):1-42, 2006a;J Comput Graph Stat 16:421-444, 2006b), permit the visualization of hypothesis tests in multivariate linear models by representing hypothesis and error matrices of sums of squares and cross-products as ellipses. This paper describes the implementation of these methods in the heplots package for R, as well as their extension, for example from two to three dimensions and by scaling hypothesis ellipses and ellipsoids in a natural manner relative to error.
It is common knowledge that Akaike's information criterion (AIC) is not a consistent model selection criterion, and Bayesian information criterion (BIC) is. These have been confirmed from an asymptotic selection p...
详细信息
It is common knowledge that Akaike's information criterion (AIC) is not a consistent model selection criterion, and Bayesian information criterion (BIC) is. These have been confirmed from an asymptotic selection probability evaluated from a large-sample framework. However, when a high-dimensional asymptotic framework, such that the dimension of the response variables and the sample size are approaching oo, is used for evaluating the selection probability, there are cases that the AIC for selecting variables in multivariate linear models is consistent, but the BIC is not. The AIC and BIC are included in a family of information criteria defined by adding a penalty term expressing the complexity of the model to a negative twofold maximum log-likelihood. By clarifying the condition of the penalty term to ensure the consistency, we derive conditions for consistency of the AIC, BIC and other information criteria under the high-dimensional asymptotic framework.
linear and quadratic prediction problems infinite populations have become of great interest to Many authors recently. In the present paper, we mainly aim to extend the problem Of quadratic prediction from a general li...
详细信息
linear and quadratic prediction problems infinite populations have become of great interest to Many authors recently. In the present paper, we mainly aim to extend the problem Of quadratic prediction from a general linearmodel, of form y = X beta + e, e similar to N(0, sigma V-2), to a multivariate linear model, denoted by Y = XB + E. Vec(E) similar to N(0, Sigma circle times V) with Y = (Y-ij)(nxq) = (y(1),...., y(q)). Firstly, the optimal invariant quadratic unbiased (OIQU) predictor and the optimal invariant quadratic (potentially) biased (OIQB) predictor of Y'HY for any particular symmetric nonnegative definite matrix H satisfying HX = 0 are derived. Secondly, we consider predicting a'Y'HYb and tr(Y'HY).The corresponding restricted OIQU predictor and restricted OIQB predictor for them are given. In addition, we also offer four concluding remarks. One concerns the generalization of predicting a'Y'HYb and tr(Y'HY), and the others are concerned with three possible extensions from OlUltivariate linearmodels to growth curve models, to restricted multivariate linear models, and to matrix elliptical linearmodels. (C) 2008 Elsevier Inc. All rights reserved.
For multivariate linear model Y=XΘ+ε, ~N(0,σ2∑(×)V), this paper is concerned with the admissibility of linear estimators of estimable function SXΘ in the class of ali estimators. All admissible linear estim...
详细信息
For multivariate linear model Y=XΘ+ε, ~N(0,σ2∑(×)V), this paper is concerned with the admissibility of linear estimators of estimable function SXΘ in the class of ali estimators. All admissible linear estimators of SXΘ are given under each of four defioitions of admissibility.
We present a novel approach for modeling mandible shape variations in Temporomandibular Joint Osteoarthritis (TMJ OA) patients. We have employed weighted spherical harmonic (SPHARM) representation to parameterize and ...
详细信息
ISBN:
(纸本)9781479974184
We present a novel approach for modeling mandible shape variations in Temporomandibular Joint Osteoarthritis (TMJ OA) patients. We have employed weighted spherical harmonic (SPHARM) representation to parameterize and normalize mandible surfaces. This representation is fed to multivariate linear models which account for nuisance covariates such as age and mandible size. multivariate linear models are implemented using SurfStat package and using this implementation one can avoid the complexity of specifying design matrices. In several multivariate shape models, the Hotelling's T-square has been used as a test statistic. In Hotelling's T-square statistic, we can test the equality of vector means without considering redundant covariates. Thus we have used SurfStat package in which Hotelling's T-square framework is generalized to incorporate additional covariates. Our proposed methodology has been applied for investigating Mandibular condyle shape variations in 19 TMJ OA subjects. Promising results have been demonstrated in lesion localization which is an important step in surgical planning and treatment.
In this article, we present a framework of estimating patterned covariance of interest in the multivariate linear models. The main idea in it is to estimate a patterned covariance by minimizing a trace distance functi...
详细信息
In this article, we present a framework of estimating patterned covariance of interest in the multivariate linear models. The main idea in it is to estimate a patterned covariance by minimizing a trace distance function between outer product of residuals and its expected value. The proposed framework can provide us explicit estimators, called outer product least-squares estimators, for parameters in the patterned covariance of the multivariate linear model without or with restrictions on regression coefficients. The outer product least-squares estimators enjoy the desired properties in finite and large samples, including unbiasedness, invariance, consistency and asymptotic normality. We still apply the framework to three special situations where their patterned covariances are the uniform correlation, a generalized uniform correlation and a generalq-dependence structure, respectively. Simulation studies for three special cases illustrate that the proposed method is a competent alternative of the maximum likelihood method in finite size samples. [ABSTRACT FROM AUTHOR]
In this paper, the authors consider an adaptive recursive algorithm by selecting an adaptive sequence for computing M-estimators in multivariatelinear regression models. Its asymptotic property is investigated. The r...
详细信息
In this paper, the authors consider an adaptive recursive algorithm by selecting an adaptive sequence for computing M-estimators in multivariatelinear regression models. Its asymptotic property is investigated. The recursive algorithm given by Miao and Wu (1996) is modified accordingly. Simu- lation studies of the Mgorithm is also provided. In addition, the Newton-Raphson iterative algorithm is considered for the purpose of comparison.
暂无评论