分块Z-循环矩阵的逆与Moore-Penrose逆的一个算法

本文给出了分块Z-循环矩阵的逆与Moorc-Pcnrosc逆的一个算法。     

Adjusting overall survival for treatment switches: commonly used methods and practical application

In parallel group trials, long‐term efficacy endpoints may be affected if some patients switch or cross over to the alternative treatment arm prior to the event. In oncology trials, switch to the experimental treatment can occur in the control arm following disease progression and potentially impact overall survival. It may be a clinically relevant question to estimate the efficacy that would have been observed if no patients had switched, for example, to estimate ‘real‐life’ clinical effectiveness for a health technology assessment. Several commonly used statistical methods are available that try to adjust time‐to‐event data to account for treatment switching, ranging from naive exclusion and censoring approaches to more complex inverse probability of censoring weighting and rank‐preserving structural failure time models. These are described, along with their key assumptions, strengths, and limitations. Best practice guidance is provided for both trial design and analysis when switching is anticipated. Available statistical software is summarized, and examples are provided of the application of these methods in health technology assessments of oncology trials. Key considerations include having a clearly articulated rationale and research question and a well‐designed trial with sufficient good quality data collection to enable robust statistical analysis. No analysis method is universally suitable in all situations, and each makes strong untestable assumptions. There is a need for further research into new or improved techniques. This information should aid statisticians and their colleagues to improve the design and analysis of clinical trials where treatment switch is anticipated. Copyright © 2013 John Wiley & Sons, Ltd.     

On the matrix-variate generalized hyperbolic distribution and its Bayesian applications

In the first part of the paper, we introduce the matrix-variate generalized hyperbolic distribution by mixing the matrix normal distribution with the matrix generalized inverse Gaussian density. The p-dimensional generalized hyperbolic distribution of [Barndorff-Nielsen, O. (1978). Hyperbolic distributions and distributions on hyperbolae. Scand. J. Stat., 5, 151–157], the matrix-T distribution and many well-known distributions are shown to be special cases of the new distribution. Some properties of the distribution are also studied. The second part of the paper deals with the application of the distribution in the Bayesian analysis of the normal multivariate linear model.     

Correcting the bias due to dependent censoring of the survival estimator by conditioning

We consider the conditional estimation of the survival function of the time T₂ to a second event as a function of the time T₁ to a first event when there is a censoring mechanism acting on their sum T₁+T₂. The problem has been motivated by a treatment interruption study aimed at improving the quality of life of HIV-infected patients. We base the analysis on the survival function of T₂ given that T₁∈I, where I represents a period of scientific interest (1 trimester, 1 year, 2 years, etc.) and propose a non-parametric estimator for the survival function of T₂ given that T₁∈I, which takes into account both the selection bias and the heterogeneity due to the dependent censoring. The proposed estimator for the survival function uses the risk group of T₂ conditioned on the categories of T₁ and corrects for the dependent censoring using weights defined by the observed values of T₁. The estimator, properly normalized, converges weakly to a zero-mean Gaussian process. We estimate the variance of the limiting process via a bootstrap methodology. Properties of the proposed estimator are illustrated by an extensive simulation study. The motivating data set is analysed by means of this new methodology.     

On asymptotic approximation of inverse moments for a class of nonnegative random variables

Sung [On inverse moments for a class of nonnegative random variables. J Inequal Appl. 2010;2010:1–13. Article ID 823767, doi:10.1155/2010/823767] obtained the asymptotic approximation of inverse moments for a class of nonnegative random variables with finite second moments and satisfying a Rosenthal-type inequality. In the paper, we further study the asymptotic approximation of inverse moments for a class of nonnegative random variables with finite first moments, which generalizes and improves the corresponding ones of Wu et al. [Asymptotic approximation of inverse moments of nonnegative random variables. Statist Probab Lett. 2009;79:1366–1371], Wang et al. [Exponential inequalities and inverse moment for NOD sequence. Statist Probab Lett. 2010;80:452–461; On complete convergence for weighted sums of ? mixing random variables. J Inequal Appl. 2010;2010:1–13, Article ID 372390, doi:10.1155/2010/372390], Sung (2010) and Hu et al. [A note on the inverse moment for the nonnegative random variables. Commun Statist Theory Methods. 2012. Article ID 673677, doi:10.1080/03610926.2012.673677].     

Bayesian analysis in multivariate regression models with conjugate priors

In this paper, we consider the full rank multivariate regression model with matrix elliptically contoured distributed errors. We formulate a conjugate prior distribution for matrix elliptical models and derive the posterior distributions of mean and scale matrices. In the sequel, some characteristics of regression matrix parameters are also proposed.     

A Short-Cut Derivation for the Solution of Autoregressive Models from Sharp Algebraic Arguments

This article uses algebraic arguments to cast light on the solution of vector autoregressive models in the presence of unit roots. First, the linear case and then the multi-lag specification are investigated. Clear-cut representations of the model solutions are obtained, closed-form expressions of the coefficient matrices are provided, and integration features and cointegration mechanisms for stationarity recovery are elucidated.     

Gamma-Generalized Inverse Gaussian Class of Distributions with Applications

In this article, a new family of probability distributions with domain in ?⁺ is introduced. This class can be considered as a natural extension of the exponential-inverse Gaussian distribution in Bhattacharya and Kumar (1986 Bhattacharya , S. K. , Kumar , S. ( 1986 ). E-IG model in life testing . Calcutta Statist. Assoc. Bull. 35 : 85 – 90 . [Google Scholar]) and Frangos and Karlis (2004 Frangos , N. , Karlis , D. ( 2004 ). Modelling losses using an exponential-inverse Gaussian distribution . Insur. Math. Econo. 35 : 53 – 67 .[Crossref], [Web of Science ®] , [Google Scholar]). This new family is obtained through the mixture of gamma distribution with generalized inverse Gaussian distribution. We also show some important features such as expressions of probability density function, moments, etc. Special attention is paid to the mixture with the inverse Gaussian distribution, as a particular case of the generalized inverse Gaussian distribution. From the exponential-inverse Gaussian distribution two one-parameter family of distributions are obtained to derive risk measures and credibility expressions. The versatility of this family has been proven in numerical examples.     

On bivariate transformation of scale distributions

ABSTRACT

Elsewhere, I have promoted (univariate continuous) “transformation of scale” (ToS) distributions having densities of the form 2g(Π^?1(x)) where g is a symmetric distribution and Π is a transformation function with a special property. Here, I develop bivariate (readily multivariate) ToS distributions. Univariate ToS distributions have a transformation of random variable relationship with Azzalini-type skew-symmetric distributions; the bivariate ToS distribution here arises from marginal variable transformation of a particular form of bivariate skew-symmetric distribution. Examples are given, as are basic properties—unimodality, a covariance property, random variate generation—and connections with a bivariate inverse Gaussian distribution are pointed out.