Bayesian analysis of an inverse Gaussian correlated frailty model

In survival analysis frailty is often used to model heterogeneity between individuals or correlation within clusters. Typically frailty is taken to be a continuous random effect, yielding a continuous mixture distribution for survival times. A Bayesian analysis of a correlated frailty model is discussed in the context of inverse Gaussian frailty. An MCMC approach is adopted and the deviance information criterion is used to compare models. As an illustration of the approach a bivariate data set of corneal graft survival times is analysed. (C) 2006 Elsevier B.V. All rights reserved.

Sample-induced beam distortions in terahertz time domain spectroscopy and imaging systems

It is demonstrated that distortion of the terahertz beam profile and generation of a cross-polarised component occur when the beam in terahertz time domain spectroscopy and imaging systems interacts with the sample under test. These distortions modify the detected signal, leading to spectral and image artefacts. The degree of distortion depends on the optical design of the system as well as the properties of the sample.

Upper limb tremor suppression in ADL via an orthosis incorporating a controllable double viscous beam actuator

Neuromuscular disorders affect millions of people world-wide. Upper limb tremor is a common symptom, and due to its complex aetiology it is difficult to compensate for except, in particular cases by surgical intervention or drug therapy. Wearable devices that mechanically compensate for limb tremor could benefit a considerable number of patients, but the technology to assist suffers in this way is under-developed. In this paper we propose an innovative orthosis that can dynamically suppress pathological tremor, by applying viscous damping to the affected limb in a controlled manner. The orthosis design utilises a new actuator design based on Magneto-Rheological Fluids that efficiently deliver damping action in response to the instantaneous tremor frequency and amplitude.

Application of Gaussian multi-scale representation to feature tracking in meteorological satellite imagery

Gaussian multi-scale representation is a mathematical framework that allows to analyse images at different scales in a consistent manner, and to handle derivatives in a way deeply connected to scale. This paper uses Gaussian multi-scale representation to investigate several aspects of the derivation of atmospheric motion vectors (AMVs) from water vapour imagery. The contribution of different spatial frequencies to the tracking is studied, for a range of tracer sizes, and a number of tracer selection methods are presented and compared, using WV 6.2 images from the geostationary satellite MSG-2.

Comparative study of Barium Fluoride thin films produced by ion-beam sputtering and thermal evaporation

Coatings and filters for spaceflight far-infrared components require a robust, non-absorptive low-index thin film material to contrast with the typically higher refractive index infrared materials. Barium fluoride is one such material for the 10 to 20µm wavelength infrared region, however its optical and mechanical properties vary depending on the process used to deposit it in thin film form. Thin films of dielectric produced by thermal evaporation are well documented as having a lower packing density and refractive index than bulk material. The porous and columnar micro structure of these films causes possible deterioration of their performance in varied environmental conditions, primarily because of the moisture absorption. Dielectric thin films produced by the more novel technique of ion-beam sputtering are denser with no columnar micro structure and have a packing density and refractive index similar to the bulk material. A comparative study of Barium Fluoride (BaF2) thin films made by conventional thermal evaporation and ion-beam sputtering is reported. Films of similar thicknesses are deposited on Cadmium Telluride and Germanium substrates. The optical and mechanical properties of these films are then examined. The refractive index n of the films is obtained from applying the modified Manifacier's evvelope method to the spectral measurements made on a Perkin Elmer 580 spectrophotometer. An estimate is also made of the value of the extinction coefficient k in the infrared wavelength transparent region of the thin film. In order to study the mechanical properties of the BaF2 films, and evaluate their usefulness in spaceflight infrared filters and coatings, the thin film samples are subjected to MIL-F-48616 environmental tests. Comparisons are made of their performance under these tests.

Robust initialisation of Gaussian radial basis function networks using partitioned k-means clustering

Radial basis function networks can be trained quickly using linear optimisation once centres and other associated parameters have been initialised. The authors propose a small adjustment to a well accepted initialisation algorithm which improves the network accuracy over a range of problems. The algorithm is described and results are presented.

3D scanning by multiple fan beam X-ray sources and sensors

*** Purpose – Computer tomography (CT) for 3D reconstruction entails a huge number of coplanar fan-beam projections for each of a large number of 2D slice images, and excessive radiation intensities and dosages. For some applications its rate of throughput is also inadequate. A technique for overcoming these limitations is outlined. *** Design methodology/approach – A novel method to reconstruct 3D surface models of objects is presented, using, typically, ten, 2D projective images. These images are generated by relative motion between this set of objects and a set of ten fanbeam X-ray sources and sensors, with their viewing axes suitably distributed in 2D angular space. *** Findings – The method entails a radiation dosage several orders of magnitude lower than CT, and requires far less computational power. Experimental results are given to illustrate the capability of the technique *** Practical implications – The substantially lower cost of the method and, more particularly, its dramatically lower irradiation make it relevant to many applications precluded by current techniques *** Originality/value – The method can be used in many applications such as aircraft hold-luggage screening, 3D industrial modelling and measurement, and it should also have important applications to medical diagnosis and surgery.

Suppression of directional interference for STBC MIMO system based on beam-forming

In civil applications, many researches on MIMO technique have achieved great progress. However, we consider military applications here. Differing from civil applications, military MIMO system may face many kinds of interferences, and the interference source may even not be equipped with multiple antennas. So the military MIMO system may receive some kind of strong interference coming from certain direction. Therefore, the military MIMO system must have capability to suppress directional interference. This paper presents a scheme to suppress directional interference for STBC MIMO system based on beam-forming. Simulation result shows that the scheme is valid to suppress directional strong interference for STBC MIMO system although with some performance loss compared with the ideal case of non-interference.

Gaussian quadrature 4D-VAR

A new incremental four-dimensional variational (4D-Var) data assimilation algorithm is introduced. The algorithm does not require the computationally expensive integrations with the nonlinear model in the outer loops. Nonlinearity is accounted for by modifying the linearization trajectory of the observation operator based on integrations with the tangent linear (TL) model. This allows us to update the linearization trajectory of the observation operator in the inner loops at negligible computational cost. As a result the distinction between inner and outer loops is no longer necessary. The key idea on which the proposed 4D-Var method is based is that by using Gaussian quadrature it is possible to get an exact correspondence between the nonlinear time evolution of perturbations and the time evolution in the TL model. It is shown that J-point Gaussian quadrature can be used to derive the exact adjoint-based observation impact equations and furthermore that it is straightforward to account for the effect of multiple outer loops in these equations if the proposed 4D-Var method is used. The method is illustrated using a three-level quasi-geostrophic model and the Lorenz (1996) model.

Measures of observation impact in non-Gaussian data assimilation

ABSTRACT Non-Gaussian/non-linear data assimilation is becoming an increasingly important area of research in the Geosciences as the resolution and non-linearity of models are increased and more and more non-linear observation operators are being used. In this study, we look at the effect of relaxing the assumption of a Gaussian prior on the impact of observations within the data assimilation system. Three different measures of observation impact are studied: the sensitivity of the posterior mean to the observations, mutual information and relative entropy. The sensitivity of the posterior mean is derived analytically when the prior is modelled by a simplified Gaussian mixture and the observation errors are Gaussian. It is found that the sensitivity is a strong function of the value of the observation and proportional to the posterior variance. Similarly, relative entropy is found to be a strong function of the value of the observation. However, the errors in estimating these two measures using a Gaussian approximation to the prior can differ significantly. This hampers conclusions about the effect of the non-Gaussian prior on observation impact. Mutual information does not depend on the value of the observation and is seen to be close to its Gaussian approximation. These findings are illustrated with the particle filter applied to the Lorenz ’63 system. This article is concluded with a discussion of the appropriateness of these measures of observation impact for different situations.

Extending Murty interferometry to the Terahertz part of the spectrum

We discuss some novel technologies that enable the implementation of shearing interferometry at the terahertz part of the spectrum. Possible applications include the direct measurement of lens parameters, the measurement of refractive index of materials that are transparent to terahertz frequencies, determination of homogeneity of samples, measurement of optical distortions and the non-contact evaluation of thermal expansion coefficient of materials buried inside media that are opaque to optical or infrared frequencies but transparent to THz frequencies. The introduction of a shear to a Gaussian free-space propagating terahertz beam in a controlled manner also makes possible a range of new encoding and optical signal processing modalities.

An assessment of a three-beam Doppler lidar wind profiling method for use in urban areas

Currently there are few observations of the urban wind field at heights other than rooftop level. Remote sensing instruments such as Doppler lidars provide wind speed data at many heights, which would be useful in determining wind loadings of tall buildings, and predicting local air quality. Studies comparing remote sensing with traditional anemometers carried out in flat, homogeneous terrain often use scan patterns which take several minutes. In an urban context the flow changes quickly in space and time, so faster scans are required to ensure little change in the flow over the scan period. We compare 3993 h of wind speed data collected using a three-beam Doppler lidar wind profiling method with data from a sonic anemometer (190 m). Both instruments are located in central London, UK; a highly built-up area. Based on wind profile measurements every 2 min, the uncertainty in the hourly mean wind speed due to the sampling frequency is 0.05–0.11 m s−1. The lidar tended to overestimate the wind speed by ≈0.5 m s−1 for wind speeds below 20 m s−1. Accuracy may be improved by increasing the scanning frequency of the lidar. This method is considered suitable for use in urban areas.