Maximum Likelihood Estimation of Mixture Densities for Binned and Truncated Multivariate Data

Binning and truncation of data are common in data analysis and machine learning. This paper addresses the problem of fitting mixture densities to multivariate binned and truncated data. The EM approach proposed by McLachlan and Jones (Biometrics, 44: 2, 571-578, 1988) for the univariate case is generalized to multivariate measurements. The multivariate solution requires the evaluation of multidimensional integrals over each bin at each iteration of the EM procedure. Naive implementation of the procedure can lead to computationally inefficient results. To reduce the computational cost a number of straightforward numerical techniques are proposed. Results on simulated data indicate that the proposed methods can achieve significant computational gains with no loss in the accuracy of the final parameter estimates. Furthermore, experimental results suggest that with a sufficient number of bins and data points it is possible to estimate the true underlying density almost as well as if the data were not binned. The paper concludes with a brief description of an application of this approach to diagnosis of iron deficiency anemia, in the context of binned and truncated bivariate measurements of volume and hemoglobin concentration from an individual's red blood cells.

Robust Mixture Modelling Using the t Distribution

Normal mixture models are being increasingly used to model the distributions of a wide variety of random phenomena and to cluster sets of continuous multivariate data. However, for a set of data containing a group or groups of observations with longer than normal tails or atypical observations, the use of normal components may unduly affect the fit of the mixture model. In this paper, we consider a more robust approach by modelling the data by a mixture of t distributions. The use of the ECM algorithm to fit this t mixture model is described and examples of its use are given in the context of clustering multivariate data in the presence of atypical observations in the form of background noise.

Robust algorithms for solving stochastic partial differential equations

A robust semi-implicit central partial difference algorithm for the numerical solution of coupled stochastic parabolic partial differential equations (PDEs) is described. This can be used for calculating correlation functions of systems of interacting stochastic fields. Such field equations can arise in the description of Hamiltonian and open systems in the physics of nonlinear processes, and may include multiplicative noise sources. The algorithm can be used for studying the properties of nonlinear quantum or classical field theories. The general approach is outlined and applied to a specific example, namely the quantum statistical fluctuations of ultra-short optical pulses in chi((2)) parametric waveguides. This example uses a non-diagonal coherent state representation, and correctly predicts the sub-shot noise level spectral fluctuations observed in homodyne detection measurements. It is expected that the methods used wilt be applicable for higher-order correlation functions and other physical problems as well. A stochastic differencing technique for reducing sampling errors is also introduced. This involves solving nonlinear stochastic parabolic PDEs in combination with a reference process, which uses the Wigner representation in the example presented here. A computer implementation on MIMD parallel architectures is discussed. (C) 1997 Academic Press.

Conditional production of superpositions of coherent states with inefficient photon detection

It is shown that a linear superposition of two macroscopically distinguishable optical coherent states can be generated using a single photon source and simple all-optical operations. Weak squeezing on a single photon, beam mixing with an auxiliary coherent state, and photon detecting with imperfect threshold detectors are enough to generate a coherent state superposition in a free propagating optical field with a large coherent amplitude (alpha>2) and high fidelity (F>0.99). In contrast to all previous schemes to generate such a state, our scheme does not need photon number resolving measurements nor Kerr-type nonlinear interactions. Furthermore, it is robust to detection inefficiency and exhibits some resilience to photon production inefficiency.

Genetic variability of Tomato spotted wilt virus in Australia and validation of real time RT-PCR for its detection in single and bulked leaf samples

The potential for large-scale use of a sensitive real time reverse transcription polymerase chain reaction (RT-PCR) assay was evaluated for the detection of Tomato spotted wilt virus (TSWV) in single and bulked leaf samples by comparing its sensitivity with that of DAS-ELISA. Using total RNA extracted with RNeasy (R) or leaf soak methods, real time RT-PCR detected TSWV in all infected samples collected from 16 horticultural crop species (including flowers, herbs and vegetables), two arable crop species, and four weed species by both assays. In samples in which DAS-ELISA had previously detected TSWV, real time RT-PCR was effective at detecting it in leaf tissues of all 22 plant species tested at a wide range of concentrations. Bulk samples required more robust and extensive extraction methods with real time RT-PCR, but it generally detected one infected sample in 1000 uninfected ones. By contrast, ELISA was less sensitive when used to test bulked samples, once detecting up to I infected in 800 samples with pepper but never detecting more than I infected in 200 samples in tomato and lettuce. It was also less reliable than real time RT-PCR when used to test samples from parts of the leaf where the virus concentration was low. The genetic variability among Australian isolates of TSWV was small. Direct sequencing of a 587 bp region of the nucleoprotein gene (S RNA) of 29 isolates from diverse crops and geographical locations yielded a maximum of only 4.3% nucleotide sequence difference. Phylogenetic analysis revealed no obvious groupings of isolates according to geographic origin or host species. TSWV isolates, that break TSWV resistance genes in tomato or pepper did not differ significantly in the N gene region studied, indicating that a different region of the virus genome is responsible for this trait.

An implementation of nonlinear multiuser detection in Rayleigh fading channel

A blind nonlinear interference cancellation receiver for code-division multiple-access- (CDMA-) based communication systems operating over Rayleigh flat-fading channels is proposed. The receiver which assumes knowledge of the signature waveforms of all the users is implemented in an asynchronous CDMA environment. Unlike the conventional MMSE receiver, the proposed blind ICA multiuser detector is shown to be robust without training sequences and with only knowledge of the signature waveforms. It has achieved nearly the same performance of the conventional training-based MMSE receiver. Several comparisons and experiments are performed based on examining BER performance in AWGN and Rayleigh fading in order to verify the validity of the proposed blind ICA multiuser detector.

A method to lower the detection limit for optical beat signals from a frequency-dithered stabilized laser

A narrow absorption feature in an atomic or molecular gas (such as iodine or methane) is used as the frequency reference in many stabilized lasers. As part of the stabilization scheme an optical frequency dither is applied to the laser. In optical heterodyne experiments, this dither is transferred to the RF beat signal, reducing the spectral power density and hence the signal to noise ratio over that in the absence of dither. We removed the dither by mixing the raw beat signal with a dithered local oscillator signal. When the dither waveform is matched to that of the reference laser the output signal from the mixer is rendered dither free. Application of this method to a Winters iodine-stabilized helium-neon laser reduced the bandwidth of the beat signal from 6 MHz to 390 kHz, thereby lowering the detection threshold from 5 pW of laser power to 3 pW. In addition, a simple signal detection model is developed which predicts similar threshold reductions.

Development and evaluation of neural network freeway incident detection models using field data

This paper discusses a multi-layer feedforward (MLF) neural network incident detection model that was developed and evaluated using field data. In contrast to published neural network incident detection models which relied on simulated or limited field data for model development and testing, the model described in this paper was trained and tested on a real-world data set of 100 incidents. The model uses speed, flow and occupancy data measured at dual stations, averaged across all lanes and only from time interval t. The off-line performance of the model is reported under both incident and non-incident conditions. The incident detection performance of the model is reported based on a validation-test data set of 40 incidents that were independent of the 60 incidents used for training. The false alarm rates of the model are evaluated based on non-incident data that were collected from a freeway section which was video-taped for a period of 33 days. A comparative evaluation between the neural network model and the incident detection model in operation on Melbourne's freeways is also presented. The results of the comparative performance evaluation clearly demonstrate the substantial improvement in incident detection performance obtained by the neural network model. The paper also presents additional results that demonstrate how improvements in model performance can be achieved using variable decision thresholds. Finally, the model's fault-tolerance under conditions of corrupt or missing data is investigated and the impact of loop detector failure/malfunction on the performance of the trained model is evaluated and discussed. The results presented in this paper provide a comprehensive evaluation of the developed model and confirm that neural network models can provide fast and reliable incident detection on freeways. (C) 1997 Elsevier Science Ltd. All rights reserved.

Specific PCR based detection of Phytophthora medicaginis using the intergenic spacer region of the ribosomal DNA

A technique based on the polymerase chain reaction (PCR) for the specific detection of Phytophthora medicaginis was developed using nucleotide sequence information of the ribosomal DNA (rDNA) regions. The complete IGS 2 region between the 5 S gene of one rDNA repeat and the small subunit of the adjacent repeat was sequenced for P. medicaginis and related species. The entire nucleotide sequence length of the IGS 2 of P. medicaginis was 3566 bp. A pair of oligonucleotide primers (PPED04 and PPED05), which allowed amplification of a specific fragment (364 bp) within the IGS 2 of P. medicaginis using the PCR, was designed. Specific amplification of this fragment from P. medicaginis was highly sensitive, detecting template DNA as low as 4 ng and in a host-pathogen DNA ratio of 1000000:1. Specific PCR amplification using PPED04 and PPED05 was successful in detecting P. medicaginis in lucerne stems infected under glasshouse conditions and field infected lucerne roots. The procedures developed in this work have application to improved identification and detection of a wide range of Phytophthora spp. in plants and soil.

Improved one-step solid-phase extraction method for morphine, morphine-3-glucuronide, and morphine-6-glucuronide from plasma and quantitation using high-performance liquid chromatography with electrochemical detection

This communication describes an improved one-step solid-phase extraction method for the recovery of morphine (M), morphine-3-glucuronide (M3G), and morphine-6-glucuronide (M6G) from human plasma with reduced coextraction of endogenous plasma constituents, compared to that of the authors' previously reported method. The magnitude of the peak caused by endogenous plasma components in the chromatogram that eluted immediately before the retention time of M3G has been reduced (similar to 80%) significantly (p < 0.01) while achieving high extraction efficiencies for the compounds of interest, viz morphine, M6G, and M3G (93.8 +/- 2.5, 91.7 +/- 1.7, and 93.1 +/- 2.2%, respectively). Furthermore, when the improved solid-phase extraction method was used, the extraction cartridge-derived late-eluting peak (retention time 90 to 100 minutes) reported in our previous method, was no longer present in the plasma extracts. Therefore the combined effect of reducing the recovery of the endogenous components of plasma that chromatographed just before the retention time of M3G and the removal of the late-eluting, extraction cartridge-derived peak has resulted in a decrease in the chromatographic run-time to 20 minutes, thereby increasing the sample throughput by up to 100%.

Sensitive detection of nitric oxide using seeded parametric four-wave mixing

A sensitive near-resonant four-wave mixing technique based on two-photon parametric four-wave mixing has been developed. Seeded parametric four-wave mixing requires only a single laser as an additional phase matched seeder field is generated via parametric four-wave mixing of the pump beam in a high gain cell. The seeder field travels collinearly with the pump beam providing efficient nondegenerate four-wave mixing in a second medium. This simple arrangement facilitates the detection of complex molecular spectra by simply scanning the pump laser. Seeded parametric four-wave mixing is demonstrated in both a low pressure cell and an air/acetylene flame with detection of the two-photon C (2) Pi(upsilon'=0)<--X (2) Pi(upsilon =0) spectrum of nitric oxide. From the cell data a detection limit of 10(12) molecules/cm(3) is established. A theoretical model of seeded parametric four-wave mixing is developed from existing parametric four-wave mixing theory. The addition of the seeder field significantly modifies the parametric four-wave mixing behaviour such that in the small signal regime, the signal intensity can readily be made to scale as the cube of the laser pump power while the density dependence follows a more familiar square law dependence, In general, we find excellent agreement between theory and experiment. Limitations to the process result from an ac Stark shift of the two-photon resonance in the high pressure seeder cell caused by the generation of a strong seeder field, as well as a reduction in phase matching efficiency due to the presence of certain buffer species. Various optimizations are suggested which should overcome these limitations, providing even greater detection sensitivity. (C) 1998 American Institute of Physics, [S0021-9606(98)01014-9].

Sensitive detection of sodium in a flame using parametric four-wave mixing and seeded parametric four-wave mixing

Two-photon resonant parametric four-wave mixing and a newly developed variant called seeded parametric four-wave mixing are used to detect trace quantities of sodium in a flame. Both techniques are simple, requiring only a single laser to generate a signal beam at a different wavelength which propagates collinearly with the pump beam, allowing efficient signal recovery. A comparison of the two techniques reveals that seeded parametric four-wave mixing is more than two orders of magnitude more sensitive than parametric four-wave mixing, with an estimated detection sensitivity of 5 x 10(9) atoms/cm(3). Seeded parametric four-wave mixing is achieved by cascading two parametric four-wave mixing media such that one of the parametric fields generated in the first high-density medium is then used to seed the same four-wave mixing process in a second medium in order to increase the four-wave mixing gain. The behavior of this seeded parametric four-wave mixing is described using semiclassical perturbation theory. A simplified small-signal theory is found to model most of the data satisfactorily. However, an anomalous saturationlike behavior is observed in the large signal regime. The full perturbation treatment, which includes the competition between two different four-wave mixing processes coupled via the signal field, accounts for this apparently anomalous behavior.

Solid-phase extraction method with high-performance liquid chromatography and electrochemical detection for the quantitative analysis of oxycodone in human plasma

A sensitive and reproducible solid-phase extraction (SPE) method for the quantification of oxycodone in human plasma was developed. Varian Certify SPE cartridges containing both C-8 and benzoic acid functional groups were the most suitable for the extraction of oxycodone and codeine (internal standard), with consistently high (greater than or equal to 80%) and reproducible recoveries. The elution mobile phase consisted of 1.2 ml of butyl chloride-isopropanol (80:20, v/v) containing 2% ammonia. The quantification limit for oxycodone was 5.3 pmol on-column. Within-day and inter-day coefficients of variation were 1.2% and 6.8% respectively for 284 nM oxycodone and 9.5% and 6.2% respectively for 28.4 nM oxycodone using 0.5-ml plasma aliquots. (C) 1998 Elsevier Science BN. All rights reserved.

Sensitive and specific detection of Clavibacter xyli subsp. xyli, causal agent of ratoon stunting disease of sugarcane, with a polymerase chain reaction-based assay

A sensitive, specific polymerase chain reaction-based assay was developed for the detection of the causal agent of ratoon stunting disease of sugarcane, Clavibacter xyli subsp. xyli. This assay uses oligonucleotide primers derived from the internal transcribed spacer region between the 16S and 23S rRNA genes of the bacterial rRNA operon. The assay is specific for C. xyli subsp. xyli and does not produce an amplification product from the template of the closely related bacterium C. xyli subsp. cynodontis, nor from other bacterial species. The assay was successfully applied to the detection of C. xyli subsp. xyli in fibrovascular fluid extracted from sugarcane and was sensitive to approximately 22 cells per PCR assay. A multiplex PCR test was also developed which identified and differentiated C. xyli subsp. xyli and C. xyli subsp. cynodontis in a single PCR assay.