On discrimination algorithms for ill-posed problems with an application to magnetic tomography

In this paper we explore classification techniques for ill-posed problems. Two classes are linearly separable in some Hilbert space X if they can be separated by a hyperplane. We investigate stable separability, i.e. the case where we have a positive distance between two separating hyperplanes. When the data in the space Y is generated by a compact operator A applied to the system states ∈ X, we will show that in general we do not obtain stable separability in Y even if the problem in X is stably separable. In particular, we show this for the case where a nonlinear classification is generated from a non-convergent family of linear classes in X. We apply our results to the problem of quality control of fuel cells where we classify fuel cells according to their efficiency. We can potentially classify a fuel cell using either some external measured magnetic field or some internal current. However we cannot measure the current directly since we cannot access the fuel cell in operation. The first possibility is to apply discrimination techniques directly to the measured magnetic fields. The second approach first reconstructs currents and then carries out the classification on the current distributions. We show that both approaches need regularization and that the regularized classifications are not equivalent in general. Finally, we investigate a widely used linear classification algorithm Fisher's linear discriminant with respect to its ill-posedness when applied to data generated via a compact integral operator. We show that the method cannot stay stable when the number of measurement points becomes large.

Problems with binary pattern measures for flood model evaluation

As the calibration and evaluation of flood inundation models are a prerequisite for their successful application, there is a clear need to ensure that the performance measures that quantify how well models match the available observations are fit for purpose. This paper evaluates the binary pattern performance measures that are frequently used to compare flood inundation models with observations of flood extent. This evaluation considers whether these measures are able to calibrate and evaluate model predictions in a credible and consistent way, i.e. identifying the underlying model behaviour for a number of different purposes such as comparing models of floods of different magnitudes or on different catchments. Through theoretical examples, it is shown that the binary pattern measures are not consistent for floods of different sizes, such that for the same vertical error in water level, a model of a flood of large magnitude appears to perform better than a model of a smaller magnitude flood. Further, the commonly used Critical Success Index (usually referred to as F<2 >) is biased in favour of overprediction of the flood extent, and is also biased towards correctly predicting areas of the domain with smaller topographic gradients. Consequently, it is recommended that future studies consider carefully the implications of reporting conclusions using these performance measures. Additionally, future research should consider whether a more robust and consistent analysis could be achieved by using elevation comparison methods instead.

Application of complex extreme learning machine to multiclass classification problems with high dimensionality: a THz spectra classification problem

We extend extreme learning machine (ELM) classifiers to complex Reproducing Kernel Hilbert Spaces (RKHS) where the input/output variables as well as the optimization variables are complex-valued. A new family of classifiers, called complex-valued ELM (CELM) suitable for complex-valued multiple-input–multiple-output processing is introduced. In the proposed method, the associated Lagrangian is computed using induced RKHS kernels, adopting a Wirtinger calculus approach formulated as a constrained optimization problem similarly to the conventional ELM classifier formulation. When training the CELM, the Karush–Khun–Tuker (KKT) theorem is used to solve the dual optimization problem that consists of satisfying simultaneously smallest training error as well as smallest norm of output weights criteria. The proposed formulation also addresses aspects of quaternary classification within a Clifford algebra context. For 2D complex-valued inputs, user-defined complex-coupled hyper-planes divide the classifier input space into four partitions. For 3D complex-valued inputs, the formulation generates three pairs of complex-coupled hyper-planes through orthogonal projections. The six hyper-planes then divide the 3D space into eight partitions. It is shown that the CELM problem formulation is equivalent to solving six real-valued ELM tasks, which are induced by projecting the chosen complex kernel across the different user-defined coordinate planes. A classification example of powdered samples on the basis of their terahertz spectral signatures is used to demonstrate the advantages of the CELM classifiers compared to their SVM counterparts. The proposed classifiers retain the advantages of their ELM counterparts, in that they can perform multiclass classification with lower computational complexity than SVM classifiers. Furthermore, because of their ability to perform classification tasks fast, the proposed formulations are of interest to real-time applications.

Spectral analysis of the elliptic sine-Gordon equation in the quarter plane

We study the elliptic sine-Gordon equation in the quarter plane using a spectral transform approach. We determine the Riemann-Hilbert problem associated with well-posed boundary value problems in this domain and use it to derive a formal representation of the solution. Our analysis is based on a generalization of the usual inverse scattering transform recently introduced by Fokas for studying linear elliptic problems.

Spectral analysis of diffusions with jump boundary

In this paper we consider one-dimensional diffusions with constant coefficients in a finite interval with jump boundary and a certain deterministic jump distribution. We use coupling methods in order to identify the spectral gap in the case of a large drift and prove that there is a threshold drift above which the bottom of the spectrum no longer depends on the drift. As a corollary to our result we are able to answer two questions concerning elliptic eigenvalue problems with non-local boundary conditions formulated previously by Iddo Ben-Ari and Ross Pinsky.

A comparison of duality and energy a posteriori estimates for $\mathrm {L}_{\infty }(0,T;\mathrm {L}_2(\varOmega ))$ in parabolic problems

We use the elliptic reconstruction technique in combination with a duality approach to prove a posteriori error estimates for fully discrete backward Euler scheme for linear parabolic equations. As an application, we combine our result with the residual based estimators from the a posteriori estimation for elliptic problems to derive space-error indicators and thus a fully practical version of the estimators bounding the error in the $ \mathrm {L}_{\infty }(0,T;\mathrm {L}_2(\varOmega ))$ norm. These estimators, which are of optimal order, extend those introduced by Eriksson and Johnson in 1991 by taking into account the error induced by the mesh changes and allowing for a more flexible use of the elliptic estimators. For comparison with previous results we derive also an energy-based a posteriori estimate for the $ \mathrm {L}_{\infty }(0,T;\mathrm {L}_2(\varOmega ))$-error which simplifies a previous one given by Lakkis and Makridakis in 2006. We then compare both estimators (duality vs. energy) in practical situations and draw conclusions.

"I don't know what you are trying to establish": The socio-environmental 'problems' of visual impaired children

This paper uses the social model of disability to examine visually impaired children's experiences of their housing and neighbourhoods and finds that they did not experience any significant problems with the design of them. The source of their problems was within these environments, and was caused by factors such as the intensity of movement, for example, from flows of traffic. We conclude by discussing the social policy implications of these findings.

Emergence of variants with altered survival properties in stationary phase cultures of Campylobacter jejuni

During the stationary phase of Campylobacter jejuni NCTC 11351 viable numbers fluctuate in a characteristic fashion. After reaching the maximum cell count (ca. 2 X 10(9) CFU/ml) in early stationary phase (denoted phase 1), viable numbers subsequently decrease to about 10(6) CFU/ml after 48 h and then increase again to about 10(8) CFU/ml (denoted phase 2) before decreasing once more to a value intermediate between the previous maximum and minimum values. To investigate whether the increase in viable numbers following the initial decline was due to the emergence of a new strain with a growth advantage in stationary phase analogous to the 'GASP' phenotype described in Escherichia coli [Science 259 (1993) 1757], we conducted mixed culture experiments with cells from the original culture and antibiotic-resistant marked organisms isolated from the re-growth phase. In many experiments of this type, strains isolated from phase 2 failed to out-compete the original strain and we have thus been unable to demonstrate a convincing GASP phenotype. However, strains isolated from phase 2 showed a much lower rate of viability loss in early stationary phase and a small increase in resistance to aeration, peroxide challenge and heat, indicating that the emergent strain was different from the parent. These results support the view that dynamic population changes occur during the stationary phase of C jejuni that may play a role in the survival of this organism. (C) 2003 Published by Elsevier B.V.

Upper limb robot mediated stroke therapy - GENTLE/s approach

Stroke is a leading cause of disability in particular affecting older people. Although the causes of stroke are well known and it is possible to reduce these risks, there is still a need to improve rehabilitation techniques. Early studies in the literature suggest that early intensive therapies can enhance a patient's recovery. According to physiotherapy literature, attention and motivation are key factors for motor relearning following stroke. Machine mediated therapy offers the potential to improve the outcome of stroke patients engaged on rehabilitation for upper limb motor impairment. Haptic interfaces are a particular group of robots that are attractive due to their ability to safely interact with humans. They can enhance traditional therapy tools, provide therapy "on demand" and can present accurate objective measurements of a patient's progression. Our recent studies suggest the use of tele-presence and VR-based systems can potentially motivate patients to exercise for longer periods of time. The creation of human-like trajectories is essential for retraining upper limb movements of people that have lost manipulation functions following stroke. By coupling models for human arm movement with haptic interfaces and VR technology it is possible to create a new class of robot mediated neuro rehabilitation tools. This paper provides an overview on different approaches to robot mediated therapy and describes a system based on haptics and virtual reality visualisation techniques, where particular emphasis is given to different control strategies for interaction derived from minimum jerk theory and the aid of virtual and mixed reality based exercises.

Comments on “An ensemble cumulus convection parameterization with explicit cloud treatment”

Wagner and Graf (2010) derive a population evolution equation for an ensemble of convective plumes, an analogue with the Lotka–Volterra equation, from the energy equations for convective plumes provided by Arakawa and Schubert (1974). Although their proposal is interesting, as the present note shows, there are some problems with their derivation.

The consistency of private and public real estate within mixed-asset portfolios

This study considers the consistency of the role of both the private and public real estate markets within a mixed-asset context. While a vast literature has developed that has examined the potential role of both the private and public real estate markets, most studies have largely relied on both single time horizons and single sample periods. This paper builds upon the analysis of Lee and Stevenson (2005) who examined the consistency of REITs in a US capital market portfolio. The current paper extends that by also analyzing the role of the private market. To address the question, the allocation of both the private and traded markets is evaluated over different holding periods varying from 5- to 20-years. In general the results show that optimum mixed-asset portfolios already containing private real estate have little place for public real estate securities, especially in low risk portfolios and for longer investment horizons. Additionally, mixed-asset portfolios with public real estate either see the allocations to REITs diminished or eliminated if private real estate is also considered. The results demonstrate that there is a still a strong case for private real estate in the mixed-asset portfolio on the basis of an increase in risk-adjusted performance, even if the investor is already holding REITs, but that the reverse is not always the case.

Disorder-specific functional abnormalities during sustained attention in youth with Attention Deficit Hyperactivity Disorder (ADHD) and with Autism

Attention Deficit Hyperactivity Disorder (ADHD) and Autism Spectrum Disorder (ASD) are often comorbid and share behavioural-cognitive abnormalities in sustained attention. A key question is whether this shared cognitive phenotype is based on common or different underlying pathophysiologies. To elucidate this question, we compared 20 boys with ADHD to 20 age and IQ matched ASD and 20 healthy boys using functional magnetic resonance imaging (fMRI) during a parametrically modulated vigilance task with a progressively increasing load of sustained attention. ADHD and ASD boys had significantly reduced activation relative to controls in bilateral striato–thalamic regions, left dorsolateral prefrontal cortex (DLPFC) and superior parietal cortex. Both groups also displayed significantly increased precuneus activation relative to controls. Precuneus was negatively correlated with the DLPFC activation, and progressively more deactivated with increasing attention load in controls, but not patients, suggesting problems with deactivation of a task-related default mode network in both disorders. However, left DLPFC underactivation was significantly more pronounced in ADHD relative to ASD boys, which furthermore was associated with sustained performance measures that were only impaired in ADHD patients. ASD boys, on the other hand, had disorder-specific enhanced cerebellar activation relative to both ADHD and control boys, presumably reflecting compensation. The findings show that ADHD and ASD boys have both shared and disorder-specific abnormalities in brain function during sustained attention. Shared deficits were in fronto–striato–parietal activation and default mode suppression. Differences were a more severe DLPFC dysfunction in ADHD and a disorder-specific fronto–striato–cerebellar dysregulation in ASD.

Numerical Methods for Stiff Two-Point Boundary Value Problems

We consider the two-point boundary value problem for stiff systems of ordinary differential equations. For systems that can be transformed to essentially diagonally dominant form with appropriate smoothness conditions, a priori estimates are obtained. Problems with turning points can be treated with this theory, and we discuss this in detail. We give robust difference approximations and present error estimates for these schemes. In particular we give a detailed description of how to transform a general system to essentially diagonally dominant form and then stretch the independent variable so that the system will satisfy the correct smoothness conditions. Numerical examples are presented for both linear and nonlinear problems.

Resolution of sharp fronts in the presence of model error in variational data assimilation

We show that the four-dimensional variational data assimilation method (4DVar) can be interpreted as a form of Tikhonov regularization, a very familiar method for solving ill-posed inverse problems. It is known from image restoration problems that L1-norm penalty regularization recovers sharp edges in the image more accurately than Tikhonov, or L2-norm, penalty regularization. We apply this idea from stationary inverse problems to 4DVar, a dynamical inverse problem, and give examples for an L1-norm penalty approach and a mixed total variation (TV) L1–L2-norm penalty approach. For problems with model error where sharp fronts are present and the background and observation error covariances are known, the mixed TV L1–L2-norm penalty performs better than either the L1-norm method or the strong constraint 4DVar (L2-norm)method. A strength of the mixed TV L1–L2-norm regularization is that in the case where a simplified form of the background error covariance matrix is used it produces a much more accurate analysis than 4DVar. The method thus has the potential in numerical weather prediction to overcome operational problems with poorly tuned background error covariance matrices.

Mixed copper, silver, and gold cyanides, (MxM′1–x)CN: tailoring chain structures to influence physical properties

Binary mixed-metal variants of the one-dimensional MCN compounds (M = Cu, Ag, and Au) have been prepared and characterized using powder X-ray diffraction, vibrational spectroscopy, and total neutron diffraction. A solid solution with the AgCN structure exists in the (CuxAg1–x)CN system over the range (0 ≤ x ≤ 1). Line phases with compositions (Cu1/2Au1/2)CN, (Cu7/12Au5/12)CN, (Cu2/3Au1/3)CN, and (Ag1/2Au1/2)CN, all of which have the AuCN structure, are found in the gold-containing systems. Infrared and Raman spectroscopies show that complete ordering of the type [M–C≡N–M′–N≡C−]n occurs only in (Cu1/2Au1/2)CN and (Ag1/2Au1/2)CN. The sense of the cyanide bonding was determined by total neutron diffraction to be [Ag–NC–Au–CN−]n in (Ag1/2Au1/2)CN and [Cu–NC–Au–CN−]n in (Cu1/2Au1/2)CN. In contrast, in (Cu0.50Ag0.50)CN, metal ordering is incomplete, and strict alternation of metals does not occur. However, there is a distinct preference (85%) for the N end of the cyanide ligand to be bonded to copper and for Ag–CN–Cu links to predominate. Contrary to expectation, aurophilic bonding does not appear to be the controlling factor which leads to (Cu1/2Au1/2)CN and (Ag1/2Au1/2)CN adopting the AuCN structure. The diffuse reflectance, photoluminescence, and 1-D negative thermal expansion (NTE) behaviors of all three systems are reported and compared with those of the parent cyanide compounds. The photophysical properties are strongly influenced both by the composition of the individual chains and by how such chains pack together. The NTE behavior is also controlled by structure type: the gold-containing mixed-metal cyanides with the AuCN structure show the smallest contraction along the chain length on heating.