Some properties of BLUE in a linear model and canonical correlations associated with linear transformations

Let (x, Xβ, V) be a linear model and let A′ = (A′₁, A′₂) be a p × p nonsingular matrix such that A₂X = 0, Rank A₂ = p − Rank X. We represent the BLUE and its covariance matrix in alternative forms under the conditions that the number of unit canonical correlations between y₁ ( = A₁x) and y₂ ( = A₂x) is zero. For the second problem, let x′ = (x′₁, x′₂) and let a g-inverse V⁻ of V be written as (V⁻)′ = (A′₁, A′₂). We investigate the reations (if any) between the nonzero canonical correlations {1 ₁ … ₁ > 0} due to y₁ ( = A₁x) and y₂ ( = A₂x), and the nonzero canonical correlations {1 λ₁ … λ_v+r > 0} due to x₁ and x₂. We answer some of the questions raised by Latour et al. (1987, in Proceedings, 2nd Int. Tampere Conf. Statist. (T. Pukkila and S. Puntanen, Eds.), Univ. of Tampere, Finland) in the case of the Moore-Penrose inverse V⁺ = (A′₁, A′₂) of V.     

一般Gauss-Markov模型的线性充分性

This paper provides further contributions to the theory of linear sufficiency in the general Gauss-Markov model E(y)= Xt3, Var (y)= V. The notion of linear sufficiency introduced by Baksalary and Kala(1981) and Drygas(1983) is extended for any specific estimable function c‘β. Some general results with respect to the extended concept are obtained. An essential result concerning the former notion is a direct consequence of this paper.     

Likelihood-Based Local Polynomial Fitting for Single-Index Models

The parametric generalized linear model assumes that the conditional distribution of a response Y given a d-dimensional covariate X belongs to an exponential family and that a known transformation of the regression function is linear in X. In this paper we relax the latter assumption by considering a nonparametric function of the linear combination β^TX, say η₀(β^TX). To estimate the coefficient vector β and the nonparametric component η₀ we consider local polynomial fits based on kernel weighted conditional likelihoods. We then obtain an estimator of the regression function by simply replacing β and η₀ in η₀(β^TX) by these estimators. We derive the asymptotic distributions of these estimators and give the results of some numerical experiments.     

The blue and minque in Gauss-Markoff model with linear transformation of the observable variables

For a singular linear model A = (y, Xβ, σ2 V) and its transformed model AF = (Fy, FXβ, σ2FVF'), where V is nonnegative definite and X can be rank-deficient,the expressions for the differences of the estimates for the vector of FXβ and the variance factor σ2 are given. Moreover, the necessary and sufficient conditions for the equalities of the estimates for the vector of FXβ and the variance factor σ2 are also established. In the meantime, works in Baksalary and Kala (1981) are strengthened and consequences in Puntanen and Nurhonen (1992), and Puntanen (1996) are extended.     

Common Principal Components for Dependent Random Vectors

Let the kp-variate random vector X be partitioned into k subvectors X_i of dimension p each, and let the covariance matrix Ψ of X be partitioned analogously into submatrices Ψ_ij. The common principal component (CPC) model for dependent random vectors assumes the existence of an orthogonal p by p matrix β such that β^tΨ_ijβ is diagonal for all (i, j). After a formal definition of the model, normal theory maximum likelihood estimators are obtained. The asymptotic theory for the estimated orthogonal matrix is derived by a new technique of choosing proper subsets of functionally independent parameters.     

More on Stochastic Comparisons and Dependence among Concomitants of Order Statistics

For a sample of iid observations {(X_i, Y_i)} from an absolutely continuous distribution, the multivariate dependence of concomitants Y_[]=(Y_[1], Y_[2], …, Y_[n]) and the stochastic order of subsets of Y_[] are studied. If (X, Y) is totally positive dependent of order 2, Y_[] is multivariate totally positive dependent of order 2. If the conditional hazard rate function of Y given X, h_YX(yx), is decreasing in x for every y, Y_[] is multivariate right corner set increasing. And if Y is stochastically increasing in X, the concomitants are increasing in multivariate stochastic order.     

On stochastic integral representation of stable processes with sample paths in Banach spaces

Certain path properties of a symmetric α-stable process X(t) = ∫_Sh(t, s) dM(s), t T, are studied in terms of the kernel h. The existence of an appropriate modification of the kernel h enables one to use results from stable measures on Banach spaces in studying X. Bounds for the moments of the norm of sample paths of X are obtained. This yields definite bounds for the moments of a double α-stable integral. Also, necessary and sufficient conditions for the absolute continuity of sample paths of X are given. Along with the above stochastic integral representation of stable processes, the representation of stable random vectors due to[13], Ann. Probab.9, 624–632) is extensively used and the relationship between these two representations is discussed.     

Third-order asymptomic properties of a class of test statistics under a local alternative

Suppose that {X_i; I = 1, 2, …,} is a sequence of p-dimensional random vectors forming a stochastic process. Let p_n, θ(X_n), X_n ^np, be the probability density function of X_n = (X₁, …, X_n) depending on θ Θ, where Θ is an open set of ¹. We consider to test a simple hypothesis H : θ = θ₀ against the alternative A : θ ≠ θ₀. For this testing problem we introduce a class of tests , which contains the likelihood ratio, Wald, modified Wald, and Rao tests as special cases. Then we derive the third-order asymptotic expansion of the distribution of T under a sequence of local alternatives. Using this result we elucidate various third-order asymptotic properties of T (e.g., Bartlett's adjustments, third-order asymptotically most powerful properties). Our results are very general, and can be applied to the i.i.d. case, multivariate analysis, and time series analysis. Two concrete examples will be given. One is a Gaussian ARMA process (dependent case), and the other is a nonlinear regression model (non-identically distributed case).     

Vertex Set Partitions Preserving Conservativeness

Let G be an undirected graph and ={X₁, …, X_n} be a partition of V(G). Denote by G/ the graph which has vertex set {X₁, …, X_n}, edge set E, and is obtained from G by identifying vertices in each class X_i of the partition . Given a conservative graph (G, w), we study vertex set partitions preserving conservativeness, i.e., those for which (G/ , w) is also a conservative graph. We characterize the conservative graphs (G/ , w), where is a terminal partition of V(G) (a partition preserving conservativeness which is not a refinement of any other partition of this kind). We prove that many conservative graphs admit terminal partitions with some additional properties. The results obtained are then used in new unified short proofs for a co-NP characterization of Seymour graphs by A. A. Ageev, A. V. Kostochka, and Z. Szigeti (1997, J. Graph Theory34, 357–364), a theorem of E. Korach and M. Penn (1992, Math. Programming55, 183–191), a theorem of E. Korach (1994, J. Combin. Theory Ser. B62, 1–10), and a theorem of A. V. Kostochka (1994, in “Discrete Analysis and Operations Research. Mathematics and its Applications (A. D. Korshunov, Ed.), Vol. 355, pp. 109–123, Kluwer Academic, Dordrecht).     

Normal Approximation Rate of the Kernel Smoothing Estimator in a Partial Linear Model

By establishing the asymptotic normality for the kernel smoothing estimatorβ_nof the parametric componentsβin the partial linear modelY=X′β+g(T)+, P. Speckman (1988,J. Roy. Statist. Soc. Ser. B50, 413–456) proved that the usual parametric raten^−1/2is attainable under the usual “optimal” bandwidth choice which permits the achievement of the optimal nonparametric rate for the estimation of the nonparametric componentg. In this paper we investigate the accuracy of the normal approximation forβ_nand find that, contrary to what we might expect, the optimal Berry–Esseen raten^−1/2is not attainable unlessgis undersmoothed, that is, the bandwidth is chosen with faster rate of tending to zero than the “optimal” bandwidth choice.