Likelihood-based estimation of the regression model with scrambled responses

A significant problem in the collection of responses to potentially sensitive questions, such as relating to illegal, immoral or embarrassing activities, is non-sampling error due to refusal to respond or false responses. Eichhorn & Hayre (1983) suggested the use of scrambled responses to reduce this form of bias. This paper considers a linear regression model in which the dependent variable is unobserved but for which the sum or product with a scrambling random variable of known distribution, is known. The performance of two likelihood-based estimators is investigated, namely of a Bayesian estimator achieved through a Markov chain Monte Carlo (MCMC) sampling scheme, and a classical maximum-likelihood estimator. These two estimators and an estimator suggested by Singh, Joarder & King (1996) are compared. Monte Carlo results show that the Bayesian estimator outperforms the classical estimators in almost all cases, and the relative performance of the Bayesian estimator improves as the responses become more scrambled.

Calculating current densities and fields due to shielded bi-planar radio-frequency coils

A method for the accurate computation of the current densities produced in a wide-runged bi-planar radio-frequency coil is presented. The device has applications in magnetic resonance imaging. There is a set of opposing primary rungs, symmetrically placed on parallel planes and a similar arrangement of rungs on two parallel planes surrounding the primary serves as a shield. Current densities induced in these primary and shielding rungs are calculated to a high degree of accuracy using an integral-equation approach, combined with the inverse finite Hilbert transform. Once these densities are known, accurate electrical and magnetic fields are then computed without difficulty. Some test results are shown. The method is so rapid that it can be incorporated into optimization software. Some preliminary fields produced from optimized coils are presented.

A numerical study of large sparse matrix exponentials arising in Markov chains

Krylov subspace techniques have been shown to yield robust methods for the numerical computation of large sparse matrix exponentials and especially the transient solutions of Markov Chains. The attractiveness of these methods results from the fact that they allow us to compute the action of a matrix exponential operator on an operand vector without having to compute, explicitly, the matrix exponential in isolation. In this paper we compare a Krylov-based method with some of the current approaches used for computing transient solutions of Markov chains. After a brief synthesis of the features of the methods used, wide-ranging numerical comparisons are performed on a power challenge array supercomputer on three different models. (C) 1999 Elsevier Science B.V. All rights reserved.AMS Classification: 65F99; 65L05; 65U05.

Atoms in squeezed light fields

Squeezed light is of interest as an example of a non-classical state of the electromagnetic field and because of its applications both in technology and in fundamental quantum physics. This review concentrates on one aspect of squeezed light, namely its application in atomic spectroscopy. The general properties, detection and application of squeezed light are first reviewed. The basic features of the main theoretical methods (master equations, quantum Langevin equations, coupled systems) used to treat squeezed light spectroscopy are then outlined. The physics of squeezed light interactions with atomic systems is dealt with first for the simpler case of two-level atoms and then for the more complex situation of multi-level atoms and multi-atom systems. Finally the specific applications of squeezed light spectroscopy are reviewed.

Polymerization of diglycidyl ether of bisphenol A/diaminodiphenyl sulfone epoxy resins using radio frequency fields

Curing of diglycidyl ether of bisphenol A/diaminodiphenyl sulfone (DGEBA/DDS) epoxy resin has been effected by heating with radio frequency (RF) radiation at frequencies of 30-99 MHz. The epoxy resins can be cured rapidly at low RF power levels. Comparison of the kinetics of the RF curing with thermal curing while maintaining the same curing temperature revealed no differences. Previous differences in rates of thermal and microwave curing are believed to be due to lack of temperature control during microwave curing. For RF curing,the rate of cure, at constant power level, increases at lower RF frequency, thus emphasizing one of the principal advantages of RF curing over microwave curing. (C) 1999 John Wiley & Sons, Inc.

Periodic orbit resonances in layered metals in tilted magnetic fields

The frequency dependence of the interlayer conductivity of a layered Fermi liquid in a magnetic field that is tilted away from the normal to the layers is considered. For both quasi-one- and quasi-two-dimensional systems resonances occur when the frequency is a harmonic of the frequency at which the magnetic field causes the electrons to oscillate on the Fermi surface within the layers. The intensity of the different harmonic resonances varies significantly with the direction of the field. The resonances occur for both coherent and weakly incoherent interlayer transport and so their observation does not imply the existence of a three-dimensional Fermi surface. [S0163-1829(99)51240-X].

Sausage: protein threading with flexible force fields

Sausage is a protein sequence threading program, but with remarkable run-time flexibility. Using different scripts, it can calculate protein sequence-structure alignments, search structure libraries, swap force fields, create models form alignments, convert file formats and analyse results. There are several different force fields which might be classed as knowledge-based, although they do not rely on Boltzmann statistics. Different force fields are used for alignment calculations and subsequent ranking of calculated models.

Quasistationarity of continuous-time Markov chains with positive drift

We shall study continuous-time Markov chains on the nonnegative integers which are both irreducible and transient, and which exhibit discernible stationarity before drift to infinity sets in. We will show how this 'quasi' stationary behaviour can be modelled using a limiting conditional distribution: specifically, the limiting state probabilities conditional on not having left 0 for the last time. By way of a dual chain, obtained by killing the original process on last exit from 0, we invoke the theory of quasistationarity for absorbing Markov chains. We prove that the conditioned state probabilities of the original chain are equal to the state probabilities of its dual conditioned on non-absorption, thus allowing us to establish the simultaneous existence and then equivalence, of their limiting conditional distributions. Although a limiting conditional distribution for the dual chain is always a quasistationary distribution in the usual sense, a similar statement is not possible for the original chain.

Further results on the relationship between mu-invariant measures and quasi-stationary distributions for absorbing continuous-time Markov chains

This note considers continuous-time Markov chains whose state space consists of an irreducible class, C, and an absorbing state which is accessible from C. The purpose is to provide results on mu-invariant and mu-subinvariant measures where absorption occurs with probability less than one. In particular, the well-known premise that the mu-invariant measure, m, for the transition rates be finite is replaced by the more natural premise that m be finite with respect to the absorption probabilities. The relationship between mu-invariant measures and quasi-stationary distributions is discussed. (C) 2000 Elsevier Science Ltd. All rights reserved.

New methods for determining quasi-stationary distributions for Markov chains

We shall be concerned with the problem of determining quasi-stationary distributions for Markovian models directly from their transition rates Q. We shall present simple conditions for a mu-invariant measure m for Q to be mu-invariant for the transition function, so that if m is finite, it can be normalized to produce a quasi-stationary distribution. (C) 2000 Elsevier Science Ltd. All rights reserved.

Signatures of the pair-coherent state

We explore in detail the possibility of generating a pair-coherent state in the nondegenerate parametric oscillator when decoherence is included. Such states are predicted in the transient regime in parametric oscillation where the pump mode is adiabatically eliminated. Two specific signatures are examined to indicate whether the state of interest has been generated, the Schrodinger cat state-like signatures, and the fidelity. Solutions in a transient regime reveal interference fringes which are indicative of the formation of a Schrodinger cat state. The fidelity indicates the purity of our prepared state compared with the ideal pair-coherent state.

Heat flow for Yang-Mills-Higgs fields, part I

The Yang-Mills-Higgs field generalizes the Yang-Mills field. The authors establish the local existence and uniqueness of the weak solution to the heat flow for the Yang-Mills-Higgs field in a vector bundle over a compact Riemannian 4-manifold, and show that the weak solution is gauge-equivalent to a smooth solution and there are at most finite singularities at the maximum existing time.

Saturation of a two-level atom in polychromatic fields

The response of a two-level atom in a strong polychromatic field composed of a large number of equidistant frequency components is investigated. We calculate numerically, as well as analytically,:the stationary population inversion and show that the saturation of the atomic transition strongly depends on whether or not there is a central (resonant) frequency component in the driving field. We find that, in the presence of the central component, the atom can remain in the ground state even for a strong Rabi frequency of the driving field. In addition, we find that the inversion is sensitive to the relative phase between the frequency components. When the central component is suppressed, the atomic transition saturates with the Rabi frequency independent of the relative phase.

Content and context of monocular regions determine perceived depth in random dot, unpaired background and phantom stereograms

Perceived depth was measured for three-types of stereograms with the colour/texture of half-occluded (monocular) regions either similar to or dissimilar to that of binocular regions or background. In a two-panel random dot stereogram the monocular region was filled with texture either similar or different to the far panel or left blank. In unpaired background stereograms the monocular region either matched the background or was different in colour or texture and in phantom stereograms the monocular region matched the partially occluded object or was a different colour or texture. In all three cases depth was considerably impaired when the monocular texture did not match either the background or the more distant surface. The content and context of monocular regions as well as their position are important in determining their role as occlusion cues and thus in three-dimensional layout. We compare coincidence and accidental view accounts of these effects. (C) 2002 Elsevier Science Ltd. All rights reserved.