Orthogonal sets of balanced incomplete block designs with nested rows and column

Known series of balanced incomplete block designs with nested rows and columns are used to find orthogonal sets of these designs, producing main effects plans in nested rows and columns. Two infinite series are so constructed and shown to be universally optimum for the analysis with recovery of row and column information, a benefit produced by the additional higher strata orthogonality they enjoy. One of these series achieves orthogonality with just v − 1 replicates of v treatments, fewer than required by Latin squares.     

Block Design with Nested Rows and Columns for Research on Food Acceptability Limitation

Block designs with nested rows and columns were introduced in Srivastava (1978 Srivastava , J. N. ( 1978 ). Statistical design of agricultural experiments . J. Ind. Soc. Agric. Statist. 30 : 1 – 10 . [Google Scholar]) and more specifically defined in Singh and Dey (1979 Singh , M. , Dey , A. ( 1979 ). Block designs with nested rows and columns . Biometrika 66 ( 2 ): 321 – 326 .[Crossref], [Web of Science ®] , [Google Scholar]). We present properties of some partially balanced block designs with nested rows and columns. The designs are considered for near-factorial experiments, when there are a levels of experimental factor A and b levels of experimental factor B and there is one control treatment added. In this article we extend the theory of block designs with nested rows and columns described by \L acka and Koz?owska (2009 ?acka , A. , Koz?owska , M. ( 2009 ). Planning of factorial experiments in a block design with nested rows and columns for environmental research . Environmetrics 20 ( 6 ): 730 – 742 .[Crossref], [Web of Science ®] , [Google Scholar]) and Koz?owska et al. (2011 Koz?owska , M. , ?acka , A. , Krawczyk , R. , Koz?owski , R. J. ( 2011 ). Some block designs with nested rows and columns for research on pesticide dose limitation . Environmetrics 22 ( 6 ): 781 – 788 .[Crossref], [Web of Science ®] , [Google Scholar]). We carry out our consideration by the derived mixed linear model resulting from randomization of blocks, rows and columns. For this model of observations, some properties of estimation of treatment contrasts are discussed. We calculate the efficiency factors of estimation of treatment contrasts. We formulate theorems for some partially balanced block designs with nested rows and columns. Plant protection experiment on limitation of food acceptability for T. urticae is given to show how the obtained results can be applied.     

Equineighboured balanced nested row-column designs

This paper presents equineighboured balanced nested row-column designs for v treatments arranged in b blocks each comprising pq units further grouped into p rows and q columns. These are two-dimensional designs with the property that all pairs of treatments are neighbours equally frequently at all levels in both rows and columns. Methods of construction are given both for designs based on Latin squares and those where pq≤v. Cyclic equineighboured designs are defined and tabulated for 5≤v≤10, p≤3, q≤5, where r=bpq/v is the number of replications of each treatment.     

Bayesian Hierarchical Kernelized Probabilistic Matrix Factorization

We present a Bayesian analysis framework for matrix-variate normal data with dependency structures induced by rows and columns. This framework of matrix normal models includes prior specifications, posterior computation using Markov chain Monte Carlo methods, evaluation of prediction uncertainty, model structure search, and extensions to multidimensional arrays. Compared with Bayesian probabilistic matrix factorization, which integrates a Gaussian prior for single row of the data matrix, our proposed model, namely Bayesian hierarchical kernelized probabilistic matrix factorization, imposes Gaussian Process priors over multiple rows of the matrix. Hence, the learned model explicitly captures the underlying correlation among the rows and the columns. In addition, our method requires no specific assumptions like independence of latent factors for rows and columns, which obtains more flexibility for modeling real data compared to existing works. Finally, the proposed framework can be adapted to a wide range of applications, including multivariate analysis, times series, and spatial modeling. Experiments highlight the superiority of the proposed model in handling model uncertainty and model optimization.     

A survey of nested designs

A nest with parameters (r,k,λ)→(r′,k′,λ′) is a BIBD on (b,v,r,k,λ) where each block has a distinguished sublock of cardinality k, the sublocks forming a (b,v,r,k,λ)-design.These designs are ‘nested’ in the sense of W.T. Federer (1972), who recommended the use of these designs for the sequential addition of periods in marketing experiments in order to retain Youden design properties as rows are added. Note that for a Youden design, the b columns and v treatments are in an SBIBD arrangement with parameters v=b, k=r, and λ.     

Optimal row–column design for three treatments

The A-optimality problem is solved for three treatments in a row–column layout. Depending on the numbers of rows and columns, the requirements for optimality can be decidedly counterintuitive: replication numbers need not be as equal as possible, and trace of the information matrix need not be maximal. General rules for comparing 3×3$3\times 3$ information matrices for their A-behavior are also developed, and the A-optimality problem is also solved for three treatments in simple block designs.     

Continuity and Analysis of Sequences of Principal Components

Double arrays of n rows and p columns can be regarded as n drawings from some p-dimensional population. A sequence of such arrays is considered. Principal component analysis for each array forms sequences of sample principal components and eigenvalues. The continuity of these sequences, in the sense of convergence with probability one and convergence in probability, is investigated, that appears to be informative for pattern study and prediction of principal components. Various features of paths of sequences of population principal components are highlighted through an example.     

Block clustering with collapsed latent block models

We introduce a Bayesian extension of the latent block model for model-based block clustering of data matrices. Our approach considers a block model where block parameters may be integrated out. The result is a posterior defined over the number of clusters in rows and columns and cluster memberships. The number of row and column clusters need not be known in advance as these are sampled along with cluster memberhips using Markov chain Monte Carlo. This differs from existing work on latent block models, where the number of clusters is assumed known or is chosen using some information criteria. We analyze both simulated and real data to validate the technique.     

MBMUDs: a combinatorial extension of BIBDs showing good optimality behaviour

The construction of a balanced incomplete block design (BIBD) is formulated in terms of combinatorial optimization by defining a cost function that reaches its lower bound on all and only those configurations corresponding to a BIBD. This cost function is a linear combination of distribution measures for each of the properties of a block design (number of plots, uniformity of rows, uniformity of columns, and balance). The approach generalizes naturally to a super-class BIBDs, which we call maximally balanced maximally uniform designs (MBMUDs), that allow two consecutive values for their design parameters [r,r+1;k,k+1;λ,λ+1]. In terms of combinatorial balance, MBMUDs are the closest possible approximation to BIBDs for all experimental settings where no set of admissible parameters exists. Thus, other design classes previously proposed with the same approximation aim—such as RDGs, SRDGs and NBIBDs of type I—can be viewed as particular cases of MBMUDs. Interestingly, experimental results show that the proposed combinatorial cost function has a monotonic relation with A- and D-statistical optimality in the space of designs with uniform rows and columns, while its computational cost is much lower.     

LIMIT THEOREMS FOR STANDARDIZED PARTIAL SUMS OF WEAKLY EXCHANGEABLE ARRAYS

A symmetric array of random variables is weakly exchangeable if, when the same arbitrary permutation is applied to rows and columns, the joint distribution remains the same. We consider the asymptotic distribution of the standardized sums of the elements of the upper left hand corner of the partitioned array, generalizing results on U-statistics. In general, the asymptotic distribution is normal, but if the array is first standardized by subtracting row and column means, then it is a linear form in a normal variable and independent squares of normal variables with coefficients depending on the limits of the eigenvalues of the array.     

Two graphical displays for the detection of potentially influential subsets in regression

In the context of the general linear model Y=Xβ+ε, the matrix P_z =Z(Z^TZ)^?1 Z^T , where Z=(X: Y), plays an important role in determining least squares results. In this article we propose two graphical displays for the off-diagonal as well as the diagonal elements of P_Z . The two graphs are based on simple ideas and are useful in the detection of potentially influential subsets of observations in regression. Since P_Z is invariant with respect to permutations of the columns of Z, an added advantage of these graphs is that they can be used to detect outliers in multivariate data where the rows of Z are usually regarded as a random sample from a multivariate population. We also suggest two calibration points, one for the diagonal elements of P_Z and the other for the off-diagonal elements. The advantage of these calibration points is that they take into consideration the variability of the off-diagonal as well as the diagonal elements of P_Z . They also do not suffer from masking.     

Prediction intervals for growth curves

The growth curve model Y^n×p = A^{n×p ξ mtimes;k}B^k×p+ E^nxp, where Y is an observation matrix, &sigma is a matrix of unknown parameters, A is a known matrix of rank m, B is a known matrix of rank k with 1'= (1, …, 1) as its first row, and the rows of E are independent each distributed as N^p(0,Σ,) is considered. The problem of constructing the prediction intervals for future observations using the above model is considered and approximate intervals assuming different structures on σ are derived. The results are illustrated with several data sets.     

Order-preserving Estimators and an Inequality on the Integration of Zonal Polynomial

ABSTRACT

Suppose X , p × p p.d. random matrix, has the distribution which depends on a p × p p.d. parameter matrix Σ and this distribution is orthogonally invariant. The orthogonally invariant estimator of Σ which has the eigenvalues of the same order as the eigenvalues of X is called order-preserving. We conjecture that a non-order-preserving estimator is dominated by modified order-preserving estimators with respect to the entropy (Stein's) loss function. We show that an inequality on the integration of zonal polynomial is sufficient for this conjecture. We also prove this inequality for the case p = 2.     

Interblock information in multidimensional block designs

In many experimental situations, d-way heterogeneity among experimental units may be controlled through use of multiple blocking criteria. In some cases it is reasonable to regard some or all of the block effects as random. Then the model is mixed and observations within blocks are correlated. Very general estimators of treatment effects and their dispersion matrix with recovery of interblock information are provided. They apply to designs with d > 1 blocking criteria that may be crossed, nested, or a combination thereof. These general results may be specialized to provide analyses of new classes of MBD's or used directly for numerical analyses of designs in the general class, perhaps through use as the basis for very general computer programs. Estimation of variance components is discussed, and an example is provided to illustrate adaptation of the general results.     

TWO-WAY CHANGE-POINT MODEL AND ITS APPLICATION

Exact null and alternative distributions of the two-way maximally selected x² for interaction between the ordered rows and columns are derived for each of the normal and Poisson models, respectively. The method is one of the multiple comparison procedures for ordered parameters and is useful for defining a block interaction or a two-way change-point model as a simple alternative to the two-way additive model. The construction of a confidence region for the two-way change-point is then described. An important application is found in a dose-response clinical trial with ordered categorical responses, where detecting the dose level which gives significantly higher responses than the lower doses can be formulated as a problem of detecting a change in the interaction effects.     

On the use of diagnostic biplots in model screening for genotype by environment tables

Popular rank-2 and rank-3 models for two-way tables have geometrical properties which can be used as diagnostic keys in screening for an appropriate model. Row and column levels of two-way tables are represented by points in two or three dimensional space, whereupon collinearity and coplanarity of row and column points provide diagnostic keys for informal model choice. Coordinates are obtained from a factorization of the two-way table Y in the matrix product UV ^T. The rows of U then contain row-point coordinates and the rows of V column-point coordinates. Illustrations of applications of diagnostic biplots in the literature were restricted to data from chemistry and physics with little or no noise. In plant breeding, two-way tables containing substantial amounts of noise regularly arise in the form of genotype by environment tables. To investigate the usefulness of diagnostic biplots for model screening for genotype by environment tables, data tables were generated from a range of two-way models under the addition of various amounts of noise. Chances for correct diagnosis of the generating model depended on the type of model. Diagnostic biplots on their own do not seem to provide a sufficient means for model selection for genotype by environment tables, but in combination with other methods they certainly can provide extra insight into the structure of the data.     

Resolvable row–column designs

This paper discusses resolvable row-column designs for v treatments arranged in b blocks each comprising pq units further grouped into p rows and q columns. A resolvable row-column design has v =pqs treatments set out in r groups of s blocks. Each rectangular block has p rows and q columns.     

Goodness-of-fit tests for the Gompertz distribution

Abstract

While the Gompertz distribution is often fitted to lifespan data, testing whether the fit satisfies theoretical criteria is being neglected. Here four goodness-of-fit measures – the Anderson–Darling statistic, the correlation coefficient test, a statistic using moments, and a nested test against the generalized extreme value distributions – are discussed. Along with an application to laboratory rat data, critical values calculated by the empirical distribution of the test statistics are also presented.     

On the Distribution of Matrix Quadratic Forms

A characterization of the distribution of the multivariate quadratic form given by X A X′, where X is a p × n normally distributed matrix and A is an n × n symmetric real matrix, is presented. We show that the distribution of the quadratic form is the same as the distribution of a weighted sum of non central Wishart distributed matrices. This is applied to derive the distribution of the sample covariance between the rows of X when the expectation is the same for every column and is estimated with the regular mean.     

Crisis cycles in the slightly unstable block random model

ABSTRACT

The paper proposes a new approach for studying the time to time appearing breakdowns in economy. Block random model can describe stability of large complicated systems with variable number of participants. Theoretical background of the model is given by a theorem about the eigenvalues of block random matrices [Juhász F. On the characteristic values of non-symmetric block random matrices. J Theoret Probab. 1990;67:199–205; On the structural eigenvalues of block random matrices. Linear Algebra Appl. 1996;246:225–231]. The model takes into account not only effects of participants but of groups formed from them as well. Slight instability means group level stability and participant level instability [Juhász F. On the turbulence of slightly unstable block random systems. In: Taylor C, et al., editors. Numerical methods for laminar and turbulent flow. Atlanta; 1995. p. 113–121]. Lability index of block random systems is introduced for measuring instability. It is showed that lability index of a slightly unstable block random model is growing while number of participants increases. Alteration in the number of participants makes it possible to describe crisis cycles.