Random orthogonal matrix simulation

This paper introduces a method for simulating multivariate samples that have exact means, covariances, skewness and kurtosis. We introduce a new class of rectangular orthogonal matrix which is fundamental to the methodology and we call these matrices L matrices. They may be deterministic, parametric or data specific in nature. The target moments determine the L matrix then infinitely many random samples with the same exact moments may be generated by multiplying the L matrix by arbitrary random orthogonal matrices. This methodology is thus termed “ROM simulation”. Considering certain elementary types of random orthogonal matrices we demonstrate that they generate samples with different characteristics. ROM simulation has applications to many problems that are resolved using standard Monte Carlo methods. But no parametric assumptions are required (unless parametric L matrices are used) so there is no sampling error caused by the discrete approximation of a continuous distribution, which is a major source of error in standard Monte Carlo simulations. For illustration, we apply ROM simulation to determine the value-at-risk of a stock portfolio.

A sequential two meal challenge reveals abnormalities in postprandial TAG but not glucose in men with increasing numbers of metabolic syndrome components

Objective To examine the impact of increasing numbers of metabolic syndrome (MetS) components on postprandial lipaemia. Methods Healthy men (n = 112) underwent a sequential meal postprandial investigation, in which blood samples were taken at regular intervals after a test breakfast (0 min) and lunch (330 min). Lipids and glucose were measured in the fasting sample, with triacylglycerol (TAG), non-esterified fatty acids and glucose analysed in the postprandial samples. Results Subjects were grouped according to the number of MetS components regardless of the combinations of components (0/1, 2, 3 and 4/5). As expected, there was a trend for an increase in body mass index, blood pressure, fasting TAG, glucose and insulin, and a decrease in fasting high-density lipoprotein cholesterol with increasing numbers of MetS components (P≤0.0004). A similar trend was observed for the summary measures of the postprandial TAG and glucose responses. For TAG, the area under the curve (AUC) and maximum concentration (maxC) were significantly greater in men with ≥ 3 than < 3 components (P < 0.001), whereas incremental AUC was greater in those with 3 than 0/1 and 2, and 4/5 compared with 2 components (P < 0.04). For glucose, maxC after the test breakfast (0-330 min) and total AUC (0-480 min) were higher in men with ≥ 3 than < 3 components (P≤0.001). Conclusions Our data analysis has revealed a linear trend between increasing numbers of MetS components and magnitude (AUC) of the postprandial TAG and glucose responses. Furthermore, the two meal challenge discriminated a worsening of postprandial lipaemic control in subjects with ≥ 3 MetS components.

Crystallization of random aromatic copolyesters containing flexible spacer chains and side-groups

The crystallization behaviour of a series of random copolymers of varying chemical composition is reported. For polymers containing a high proportion of alternating rigid aromatic units and flexible spacers, conventional liquid crystalline and crystalline phase behaviour is observed. The introduction of a substantial fraction of a second shorter rigid unit containing side-chains leads to a broad endotherm in the d.s.c. scan covering some 150°C. Subsequent isothermal crystallization at any point within the broad endotherm leads to the generation of sharp endotherms at temperatures just above the recrystallization temperature. We attribute this behaviour to the crystallization of clusters of molecules containing similar random sequences. Such crystals are non-periodic along the chain direction.

Biaxial optical properties of thermotropic random copolyesters

The optical microstructures of thin sections of two liquid crystalline polymers are examined in the polarizing microscope. The polymers are random copolyesters based on hydroxybenzoic and hydroxynaphthoic acids (B-N), and hydroxybenzoic acid and ethylene terephthalate (B-ET). Sections cut from oriented samples, so as to include the extrusion direction, show microstructures in which there is no apparent preferred orientation of the axes describing the local optical anisotropy. The absence of preferred orientation in the microstructure, despite marked axial alignment of molecular chain segments as demonstrated by X-Ray diffraction, is interpreted in terms of the polymer having biaxial optical properties. The implication of optical biaxiality is that, although the mesophases are nematic, the orientation of the molecules is correlated about three (orthogonal) axes over distances greater than a micron. The structure is classified as a multiaxial nematic.

Three-dimensional random Voronoi tessellations: From cubic crystal lattices to Poisson point processes

We perturb the SC, BCC, and FCC crystal structures with a spatial Gaussian noise whose adimensional strength is controlled by the parameter a, and analyze the topological and metrical properties of the resulting Voronoi Tessellations (VT). The topological properties of the VT of the SC and FCC crystals are unstable with respect to the introduction of noise, because the corresponding polyhedra are geometrically degenerate, whereas the tessellation of the BCC crystal is topologically stable even against noise of small but finite intensity. For weak noise, the mean area of the perturbed BCC and FCC crystals VT increases quadratically with a. In the case of perturbed SCC crystals, there is an optimal amount of noise that minimizes the mean area of the cells. Already for a moderate noise (a>0.5), the properties of the three perturbed VT are indistinguishable, and for intense noise (a>2), results converge to the Poisson-VT limit. Notably, 2-parameter gamma distributions are an excellent model for the empirical of of all considered properties. The VT of the perturbed BCC and FCC structures are local maxima for the isoperimetric quotient, which measures the degre of sphericity of the cells, among space filling VT. In the BCC case, this suggests a weaker form of the recentluy disproved Kelvin conjecture. Due to the fluctuations of the shape of the cells, anomalous scalings with exponents >3/2 is observed between the area and the volumes of the cells, and, except for the FCC case, also for a->0. In the Poisson-VT limit, the exponent is about 1.67. As the number of faces is positively correlated with the sphericity of the cells, the anomalous scaling is heavily reduced when we perform powerlaw fits separately on cells with a specific number of faces.

Spectrum of a Feinberg-Zee random hopping matrix

This paper provides a new proof of a theorem of Chandler-Wilde, Chonchaiya, and Lindner that the spectra of a certain class of infinite, random, tridiagonal matrices contain the unit disc almost surely. It also obtains an analogous result for a more general class of random matrices whose spectra contain a hole around the origin. The presence of the hole forces substantial changes to the analysis.

Random-walk-based characterisation of multiple access interference in a DS/SSMA system

The problem of calculating the probability of error in a DS/SSMA system has been extensively studied for more than two decades. When random sequences are employed some conditioning must be done before the application of the central limit theorem is attempted, leading to a Gaussian distribution. The authors seek to characterise the multiple access interference as a random-walk with a random number of steps, for random and deterministic sequences. Using results from random-walk theory, they model the interference as a K-distributed random variable and use it to calculate the probability of error in the form of a series, for a DS/SSMA system with a coherent correlation receiver and BPSK modulation under Gaussian noise. The asymptotic properties of the proposed distribution agree with other analyses. This is, to the best of the authors' knowledge, the first attempt to propose a non-Gaussian distribution for the interference. The modelling can be extended to consider multipath fading and general modulation

Bayesian parameter estimation for latent Markov random fields and social networks

Undirected graphical models are widely used in statistics, physics and machine vision. However Bayesian parameter estimation for undirected models is extremely challenging, since evaluation of the posterior typically involves the calculation of an intractable normalising constant. This problem has received much attention, but very little of this has focussed on the important practical case where the data consists of noisy or incomplete observations of the underlying hidden structure. This paper specifically addresses this problem, comparing two alternative methodologies. In the first of these approaches particle Markov chain Monte Carlo (Andrieu et al., 2010) is used to efficiently explore the parameter space, combined with the exchange algorithm (Murray et al., 2006) for avoiding the calculation of the intractable normalising constant (a proof showing that this combination targets the correct distribution in found in a supplementary appendix online). This approach is compared with approximate Bayesian computation (Pritchard et al., 1999). Applications to estimating the parameters of Ising models and exponential random graphs from noisy data are presented. Each algorithm used in the paper targets an approximation to the true posterior due to the use of MCMC to simulate from the latent graphical model, in lieu of being able to do this exactly in general. The supplementary appendix also describes the nature of the resulting approximation.

Sources of variation in the ampicillin-resistant Escherichia coli concentration in the feces of organic broiler chickens

Currently, there are limited published data for the population dynamics of antimicrobial-resistant commensal bacteria. This study was designed to evaluate both the proportions of the Escherichia coli populations that are resistant to ampicillin at the level of the individual chicken on commercial broiler farms and the feasibility of obtaining repeated measures of fecal E. coli concentrations. Short-term temporal variation in the concentration of fecal E. coli was investigated, and a preliminary assessment was made of potential factors involved in the shedding of high numbers of ampicillin-resistant E. coli by growing birds in the absence of the use of antimicrobial drugs. Multilevel linear regression modeling revealed that the largest component of random variation in log-transformed fecal E. coli concentrations was seen between sampling occasions for individual birds. The incorporation of fixed effects into the model demonstrated that the older, heavier birds in the study were significantly more likely (P = 0.0003) to shed higher numbers of ampicillin-resistant E. coli. This association between increasing weight and high shedding was not seen for the total fecal E. coli population (P = 0.71). This implies that, in the absence of the administration of antimicrobial drugs, the proportion of fecal E. coli that was resistant to ampicillin increased as the birds grew. This study has shown that it is possible to collect quantitative microbiological data on broiler farms and that such data could make valuable contributions to risk assessments concerning the transfer of resistant bacteria between animal and human populations.