C293. A note on using a hypothesis-dependent variance estimate

Detecting the number of signals and estimating the parameters of the signals is an important problem in signal processing. Quite a number of papers appeared in the last twenty years regarding the estimation of the parameters of the sinusoidal components but not that much of attention has been given in estimating the number of terms present in a sinusoidal signal. Fuchs developed a criterion based on the perturbation analysis of the data auto correlation matrix to estimate the number of sinusoids, which is in some sense a subjective-based method. Recently Reddy and Biradar proposed two criteria based on AIC and MDL and developed an analytical framework for analyzing the performance of these criteria. In this paper we develop a method using the extended order modelling and singular value decomposition technique similar to that of Reddy and Biradar. We use penalty function technique but instead of using any fixed penalty function like AIC or MDL, a class of penalty functions satisfying some special properties has been used. We prove that any penalty function from that special class will give consistent estimate under the assumptions that the error random variables are independent and identically distributed with mean zero and finite variance. We also obtain the probabilities of wrong detection for any particular penalty function under somewhat weaker assumptions than that of Reddy and Biradar of Kaveh et al. It gives some idea to choose the proper penalty function for any particular model. Simulations are performed to verify the usefulness of the analysis and to compare our methods with the existing ones.     

Analysis of Variance for Bayesian Inference

This paper develops a multiway analysis of variance for non-Gaussian multivariate distributions and provides a practical simulation algorithm to estimate the corresponding components of variance. It specifically addresses variance in Bayesian predictive distributions, showing that it may be decomposed into the sum of extrinsic variance, arising from posterior uncertainty about parameters, and intrinsic variance, which would exist even if parameters were known. Depending on the application at hand, further decomposition of extrinsic or intrinsic variance (or both) may be useful. The paper shows how to produce simulation-consistent estimates of all of these components, and the method demands little additional effort or computing time beyond that already invested in the posterior simulator. It illustrates the methods using a dynamic stochastic general equilibrium model of the US economy, both before and during the global financial crisis.     

The relative impact of interviewer effects and sample design effects on survey precision

One of the principal sources of error in data collected from structured face-to-face interviews is the interviewer. The other major component of imprecision in survey estimates is sampling variance. It is rare, however, to find studies in which the complex sampling variance and the complex interviewer variance are both computed. This paper compares the relative impact of interviewer effects and sample design effects on survey precision by making use of an interpenetrated primary sampling unit–interviewer experiment which was designed by the authors for implementation in the second wave of the British Household Panel Study as part of its scientific programme. It also illustrates the use of a multilevel (hierarchical) approach in which the interviewer and sample design effects are estimated simultaneously while being incorporated in a substantive model of interest.     

Maximum likelihood analysis of the mixed model: the balanced case

The primary purpose of this paper is to develop, analytically, the inverse of the covariance matrix for the mixed analysis-of-variance model with balanced data. The use of this matrix in the identification of minimal sufficient statistics and in developing the likelihood equations is illustrated.     

Visualization for Large‐scale Gaussian Updates

In geostatistics and also in other applications in science and engineering, it is now common to perform updates on Gaussian process models with many thousands or even millions of components. These large‐scale inferences involve modelling, representational and computational challenges. We describe a visualization tool for large‐scale Gaussian updates, the ‘medal plot’. The medal plot shows the updated uncertainty at each observation location and also summarizes the sharing of information across observations, as a proxy for the sharing of information across the state vector (or latent process). As such, it reflects characteristics of both the observations and the statistical model. We illustrate with an application to assess mass trends in the Antarctic Ice Sheet, for which there are strong constraints from the observations and the physics.     

A generalized improvement procedure for variance bounds for minimum variance unbiased estimator in inverse sampling

New bounds are obtained for the variance of the minimum variance unbiased estimator of p i n inverse sampling. A generalized procedure for further improving the bounds is also discussed.     

On the variance estimate of a wilcoxon-mann whitney statistic for group ordered data

Two consistent estimators for the non-null variance of Wil-coxon-Mann-Whitney’s statistic applied to grouped ordered data, are considered. The first is based on U-statistics and the sec-ond is obtained by the Delta method. Some examples are given to demonstrate the extent of error when using a null variance esti-mate for constructing confidence intervals. It appears that the two consistent estimates are very close, but may both be disting-uishably larger or smaller than the null variance estimate.     

Multivariate Process Variability Monitoring

Understanding multivariate variability is a difficult task because there is no single measure that can be properly used. This article presents a new measure that features good properties. If this measure is simultaneously used with generalized variance, it will give a better understanding of multivariate variability. It can also efficiently be used for large data sets with high dimensions. Furthermore, when it is used for constructing a Shewhart-type chart to monitor multivariate variability, the resulting chart has a much better out-of-control ARL than the generalized variance chart. An example illustrates its advantage.     

On the asymptotic normality of linear rank statistics

A method for proving the asymptotic normality of linear rank statistics is presented. Use is made of the fact that the variances of many rank statistics met in practice admit a factorization expressed in terms of the sample size. The method of proof consists of two steps corresponding to the limiting behavior of the modulus and argument of the characteristic function of the rank statistic.     

几种方差概念的比较

当分布数列中各标志值代表的含义不完全一样时，方差仅仅描述总体的离散程度，而半方差概念则能更准确地描述所关注的标志值的离散程度；相对方差反映的是分布数列中各标志值的相对离散程度，与标准差系数一样，相对方差可用于不同性质、不同单位或具有不同水平的数列的离散程度的比较，而方差却不适合用于这种比较；半相对方差是综合了半方差和相对方差优点的一个新指标。上述几种方差概念，从不同视角描述了数列的离散程度，各有优劣点，而且相互间有着紧密的联系。     

On the coefficient of determination for mixed regression models

For mixed regression models, we define a variance decomposition including three terms, explained individual variance, unexplained individual variance and noise variance. In contrast to traditional variance decomposition, it is thus the unexplained  , not the explained, variance that is split. It gives rise to a coefficient of individual determination (CID) defined as the estimated fraction of explained individual variance. We argue that in many applications CID is a valuable complement to R²$R^2$, since it excludes noise variance (which can never be explained) and thus has one as a natural upper bound.     

Umvu estimators for the population mean and variance based on random effects models for lognormal data

Cross-classified data are often obtained in controlled experimental situations and in epidemiologic studies. As an example of the latter, occupational health studies sometimes require personal exposure measurements on a random sample of workers from one or more job groups, in one or more plant locations, on several different sampling dates. Because the marginal distributions of exposure data from such studies are generally right-skewed and well-approximated as lognormal, researchers in this area often consider the use of ANOVA models after a logarithmic transformation. While it is then of interest to estimate original-scale population parameters (e.g., the overall mean and variance), standard candidates such as maximum likelihood estimators (MLEs) can be unstable and highly biased. Uniformly minimum variance unbiased (UMVU) cstiniators offer a viable alternative, and are adaptable to sampling schemes that are typiral of experimental or epidemiologic studies. In this paper, we provide UMVU estimators for the mean and variance under two random effects ANOVA models for logtransformed data. We illustrate substantial mean squared error gains relative to the MLE when estimating the mean under a one-way classification. We illustrate that the results can readily be extended to encompass a useful class of purely random effects models, provided that the study data are balanced.     

Counterexamples to a result of broemeling on simultaneous inferences for variance ratios of some mixed linear models

In a recent paper5 Broemeling (1978) extended his earlier work on one-sided confidence regions for the variance ratios of balanced random-effects models to the two-sided case. The extension depends on a probability Inequality which was claimed to be tru We show here that it is false, hence the proof of the main result given in Ms parer is in error W also show Lhat the ntatement of his result remains true in certain special cases.     

Efficient variance reduction approach based on the variation of the input importance

ABSTRACT

To reduce the output variance, the variance-based importance analysis can provide an efficient way by reducing the variance of the ‘important’ inputs. But with the reduction of the variance of those ‘important’ inputs, the input importance will change and it is no longer the most efficient way to reduce the variance of those ‘important’ inputs alone. Thus, analyst needs to consider reducing the variance of other inputs to obtain a more efficient way. This work provides a graphical solution for analyst to decide how to reduce the input variance to achieve the targeted reduction of the output variance efficiently. Furthermore, by the importance sampling-based approach, the graphical solution can be obtained with only a single group of samples, which can decrease the computational cost greatly.     

An approximation to the distribution of the sample variance

A simple three-moment approximation is introduced for the distribution of the sample variance. Comparisons are given with other approximations discussed by Tan and Wong (1977) and with an approximation developed very recently by Mudholkar and Trivedi (1981).     

Model selection for infinite variance time series

In this we consider the problem of model selection for infinite variance time series. We introduce a group of model selection critera based on a general loss function Ψ. This family includes various generalizations of predictive least square and AIC Parameter estimation is carried out using Ψ. We use two loss functions commonly used in robust estimation and show that certain criteria out perform the conventional approach based on least squares or Yule-Walker estima­tion for heavy tailed innovations. Our conclusions are based on a comprehensive study of the performance of competing criteria for a wide selection of AR(2) models. We also consider the performance of these techniques when the ‘true’ model is not contained in the family of candidate models.     

Optimal allocation with fixed power in a completely randomized design with levels of subsampling

In many completely randomized design experiments, levels of subsampling may be performed on each experimental unit. In such cases the expected mean square error E(MSE) for testing among treatment groups is comprised of variance components analogour to those associated with the primary sampling unit is nested sampling Marcuse (1949) gives a procedure to minimize the cost of obtaining the samples if a desired degree of precision in the E(MSE) is fixed. However, her method gives no consideration to the resulting power of the test for differences among the treatment groups. Our method stipulates that the power, rather than the precision, is fixed at a critical level and the total cost is minimized subject to this constraint.