Testing equality of two normal covariance matrices with monotone missing data

Abstract

The problem of testing equality of two multivariate normal covariance matrices is considered. Assuming that the incomplete data are of monotone pattern, a quantity similar to the Likelihood Ratio Test Statistic is proposed. A satisfactory approximation to the distribution of the quantity is derived. Hypothesis testing based on the approximate distribution is outlined. The merits of the test are investigated using Monte Carlo simulation. Monte Carlo studies indicate that the test is very satisfactory even for moderately small samples. The proposed methods are illustrated using an example.     

A Bayesian Analysis of Abundance,Trend, and Population Viability for Harbor Seals in Iliamna Lake,Alaska

Harbor seals in Iliamna Lake, Alaska, are a small, isolated population, and one of only two freshwater populations of harbor seals in the world, yet little is known about their abundance or risk for extinction. Bayesian hierarchical models were used to estimate abundance and trend of this population. Observational models were developed from aerial survey and harvest data, and they included effects for time of year and time of day on survey counts. Underlying models of abundance and trend were based on a Leslie matrix model that used prior information on vital rates from the literature. We developed three scenarios for variability in the priors and used them as part of a sensitivity analysis. The models were fitted using Markov chain Monte Carlo methods. The population production rate implied by the vital rate estimates was about 5% per year, very similar to the average annual harvest rate. After a period of growth in the 1980s, the population appears to be relatively stable at around 400 individuals. A population viability analysis assessing the risk of quasi‐extinction, defined as any reduction to 50 animals or below in the next 100 years, ranged from 1% to 3%, depending on the prior scenario. Although this is moderately low risk, it does not include genetic or catastrophic environmental events, which may have occurred to the population in the past, so our results should be applied cautiously.     

Heteroscedastic Global Tests that the Regression Parameters for Two or More Independent Groups are Identical

Simulation results are reported on methods that allow both within group and between group heteroscedasticity when testing the hypothesis that independent groups have identical regression parameters. The methods are based on a combination of extant techniques, but their finite-sample properties have not been studied. Included are results on the impact of removing all leverage points or just bad leverage points. The method used to identify leverage points can be important and can improve control over the Type I error probability. Results are illustrated using data from the Well Elderly II study.     

Restricted minimax estimation in semiparametric linear models

Abstract

In this article we develop the minimax estimation approach of general linear models to the semiparametric linear models when the parameters are simultaneously constrained by an ellipsoid and linear restrictions. Combining sample information and prior constraints the minimax estimator is obtained and compared with partially least square estimator by theoretical and simulation methods.     

The residual effect of a growth-decay mechanism and the distributions of covariance structures

The product of two independent or dependent scalar normal variables, sums of products, sample covariances, and general bilinear forms are considered. Their distributions are shown to belong to a class called generalized Laplacian. A growth-decay mechanism is also shown to produce such a generalized Laplacian. Sets of necessary and sufficient conditions are derived for bilinear forms to belong to this class. As a generalization, the distributions of rectangular matrices associated with multivariate normal random vectors are also discussed.     

复亚半正定矩阵

给出了复亚半正定矩阵的概念，研究了它的基本性质及行列式理论，将Hermite阵的Schur定理，华罗庚定理，Minkowski不等式，凸性不等式，Ostrowski-Taussky不等式推广到了较广泛的复矩阵类，扩大了Minkowski不等式的指数范围，削弱了华罗庚不等式的条件。     

Computing maximum likelihood estimates from type II doubly censored exponential data

It is well-known that, under Type II double censoring, the maximum likelihood (ML) estimators of the location and scale parameters, θ and δ, of a twoparameter exponential distribution are linear functions of the order statistics. In contrast, when θ is known, theML estimator of δ does not admit a closed form expression. It is shown, however, that theML estimator of the scale parameter exists and is unique. Moreover, it has good large-sample properties. In addition, sharp lower and upper bounds for this estimator are provided, which can serve as starting points for iterative interpolation methods such as regula falsi. Explicit expressions for the expected Fisher information and Cramér-Rao lower bound are also derived. In the Bayesian context, assuming an inverted gamma prior on δ, the uniqueness, boundedness and asymptotics of the highest posterior density estimator of δ can be deduced in a similar way. Finally, an illustrative example is included.     

关于第二类级联构造Bent序列的讨论(英文)

考虑用 4个长为 2 n - 2 的Bent序列或线性序列级联构造Bent序列的问题 ,给出了第二类线性基Bent序列的类型以及级联序列为第二类Bent基Bent序列的一个充分条件 .     

增长曲线模型非齐次线性估计的可容许性

讨论增长曲线模型Y =X1BX2 +ε中回归矩阵B的函数C1BC2 的估计L1YL2 +A ,在矩阵损失 (LT2 L1)Y +A - (ST2 XT2 S1X1)B (LT2 L1)Y +A - (ST2 XT2 S1X1)B T 下 ,我们得到了非齐次线性估计L1YL2 +A在非齐次线性估计类Г ={L1YL2 +A|L1:t×p ,L2 ;n×n ,A :t×s均为已知实阵 }中可容许的充要条件 :L1YL2在Г0 ={L1YL2 |L1:t×p ,L2 :n×s均为已知实阵 }中容许且当LT2 XT2 L1X1=ST2 XT2 S1X1时有A =0。     

Examples of the use of an equivalence theorem in constructing optimum experimental designs for random-effects nonlinear regression models

The problem considered is that of finding an optimum measurement schedule to estimate population parameters in a nonlinear model when the patient effects are random. The paper presents examples of the use of sensitivity functions, derived from the General Equivalence Theorem for D-optimality, in the construction of optimum population designs for such schedules. With independent observations, the theorem applies to the potential inclusion of a single observation. However, in population designs the observations are correlated and the theorem applies to the inclusion of an additional measurement schedule. In one example, three groups of patients of differing size are subject to distinct schedules. Numerical, as opposed to analytical, calculation of the sensitivity function is advocated. The required covariances of the observations are found by simulation.