Efficient arrangements of two-level orthogonal arrays in two and four blocks

When orthogonal arrays are used in practical applications, it is often difficult to perform all the designed runs of the experiment under homogeneous conditions. The arrangement of factorial runs into blocks is usually an action taken to overcome such obstacles. However, an arbitrary configuration might lead to spurious analysis results. In this work, the nice properties of two-level orthogonal arrays are taken into consideration and an effective method for arranging experimental runs into two and four blocks of the same size is proposed. This method is based on the so-called J-characteristics of the corresponding array. General theoretical results are given for studying up to four experimental factors in two blocks, as well as for studying up to three experimental factors in four blocks. Finally, we provide best blocking arrangements when the number of the factors of interest is larger, by exploiting the known lists of non-isomorphic orthogonal arrays with two levels and various run sizes.     

Seemingly unrelated systems of econometric equations

A generalization of Zellner's SUR model is derived for sets of seemingly unrelated systems of econometric equations. The resulting structural form – worked out for a set of Cowles Commission-type simultaneous equations systems – is general enough to include any SUR-type or panel-type specification of systems of econometric equations with contemporaneously correlated errors. Maximum estimation efficiency is obtained by treating all the individual subsystems at once rather than in a subsystem-by-subsystem fashion.     

Spatio‐temporal Ornstein–Uhlenbeck Processes: Theory,Simulation and Statistical Inference

Spatio‐temporal modelling is an increasingly popular topic in Statistics. Our paper contributes to this line of research by developing the theory, simulation and inference for a spatio‐temporal Ornstein–Uhlenbeck process. We conduct detailed simulation studies and demonstrate the practical relevance of these processes in an empirical study of radiation anomaly data. Finally, we describe how predictions can be carried out in the Gaussian setting.     

Marginal models for the association structure of hierarchical binary responses

Clustered binary responses are often found in ecological studies. Data analysis may include modeling the marginal probability response. However, when the association is the main scientific focus, modeling the correlation structure between pairs of responses is the key part of the analysis. Second-order generalized estimating equations (GEE) are established in the literature. Some of them are more efficient in computational terms, especially facing large clusters. Alternating logistic regression (ALR) and orthogonalized residual (ORTH) GEE methods are presented and compared in this paper. Simulation results show a slightly superiority of ALR over ORTH. Marginal probabilities and odds ratios are also estimated and compared in a real ecological study involving a three-level hierarchical clustering. ALR and ORTH models are useful for modeling complex association structure with large cluster sizes.     

Generalized correlation and kernel causality with applications in development economics

New generalized correlation measures of 2012, GMC(Y|X), use Kernel regressions to overcome the linearity of Pearson's correlation coefficients. A new matrix of generalized correlation coefficients is such that when |r*_ij| > |r*_ji|, it is more likely that the column variable X_j is what Granger called the “instantaneous cause” or what we call “kernel cause” of the row variable X_i. New partial correlations ameliorate confounding. Various examples and simulations support robustness of new causality. We include bootstrap inference, robustness checks based on the dependence between regressor and error, and on the out-of-sample forecasts. Data for 198 countries on nine development variables support growth policy over redistribution and Deaton's criticism of foreign aid. Potential applications include Big Data, since our R code is available in the online supplementary material.     

Statistical inference for the Gompertz distribution based on Type-II progressively hybrid censored data

In this article, the statistical inference for the Gompertz distribution based on Type-II progressively hybrid censored data is discussed. The estimation of the parameters for Gompertz distribution is obtained using maximum likelihood method (MLE) and Bayesian method under three different loss functions. We also proved the existence and uniqueness of the MLE. The one-sample Bayesian prediction intervals are obtained. The work is done for different values of the parameters. We apply the Monto Carlo simulation to compare the proposed methods, also an example is discussed to construct the Prediction intervals.     

A general long-term aging model with different underlying activation mechanisms: Modeling,Bayesian estimation,and case influence diagnostics

In this paper we propose a general cure rate aging model. Our approach enables different underlying activation mechanisms which lead to the event of interest. The number of competing causes of the event of interest is assumed to follow a logarithmic distribution. The model is parameterized in terms of the cured fraction which is then linked to covariates. We explore the use of Markov chain Monte Carlo methods to develop a Bayesian analysis for the proposed model. Moreover, some discussions on the model selection to compare the fitted models are given, as well as case deletion influence diagnostics are developed for the joint posterior distribution based on the ψ-divergence, which has several divergence measures as particular cases, such as the Kullback–Leibler (K-L), J-distance, L₁ norm, and χ²-square divergence measures. Simulation studies are performed and experimental results are illustrated based on a real malignant melanoma data.     

Slashed generalized Rayleigh distribution

We introduce a new family of distributions suitable for fitting positive data sets with high kurtosis which is called the Slashed Generalized Rayleigh Distribution. This distribution arises as the quotient of two independent random variables, one being a generalized Rayleigh distribution in the numerator and the other a power of the uniform distribution in the denominator. We present properties and carry out estimation of the model parameters by moment and maximum likelihood (ML) methods. Finally, we conduct a small simulation study to evaluate the performance of ML estimators and analyze real data sets to illustrate the usefulness of the new model.