Sequential Monte Carlo for Bayesian sequentially designed experiments for discrete data

In this paper we present a sequential Monte Carlo algorithm for Bayesian sequential experimental design applied to generalised non-linear models for discrete data. The approach is computationally convenient in that the information of newly observed data can be incorporated through a simple re-weighting step. We also consider a flexible parametric model for the stimulus-response relationship together with a newly developed hybrid design utility that can produce more robust estimates of the target stimulus in the presence of substantial model and parameter uncertainty. The algorithm is applied to hypothetical clinical trial or bioassay scenarios. In the discussion, potential generalisations of the algorithm are suggested to possibly extend its applicability to a wide variety of scenarios

Anisotropy of articular cartilage reflects the ECM gradient architecture : Hough-Radon Transform Analysis

Tissue-specific extracellular matrix (ECM) is known to be an ideal bioscaffold to inspire the future of regenerative medicine. It holds the secret of how nature has developed such an organization of molecules into a unique functional complexity. This work exploited an innovative image processing algorithm and high resolution microscopy associated with mechanical analysis to establish a correlation between the gradient organization of cartiligous ECM and its anisotropic biomechanical response. This was hypothesized to be a reliable determinant that can elucidate how microarchitecture interrelates with biomechanical properties. Hough-Radon transform of the ECM cross-section images revealed its conformational variation from tangential interface down to subchondral region. As the orientation varied layer by layer, the anisotropic mechanical response deviated relatively. Although, results were in good agreement (Kendall's tau-b > 90%), there were evidences proposing that alignment of the fibrous network, specifically in middle zone, is not as random as it was previously thought.

Practical stability of approximating discrete-time filters with respect to model mismatch

This paper establishes practical stability results for an important range of approximate discrete-time filtering problems involving mismatch between the true system and the approximating filter model. Practical stability is established in the sense of an asymptotic bound on the amount of bias introduced by the model approximation. Our analysis applies to a wide range of estimation problems and justifies the common practice of approximating intractable infinite dimensional nonlinear filters by simpler computationally tractable filters.

An extension to the rank constraint for improving the reliability of stereo matching algorithms using the rank transform

The rank transform is one non-parametric transform which has been applied to the stereo matching problem The advantages of this transform include its invariance to radio metric distortion and its amenability to hardware implementation. This paper describes the derivation of the rank constraint for matching using the rank transform Previous work has shown that this constraint was capable of resolving ambiguous matches thereby improving match reliability A new matching algorithm incorporating this constraint was also proposed. This paper extends on this previous work by proposing a matching algorithm which uses a dimensional match surface in which the match score is computed for every possible template and match window combination. The principal advantage of this algorithm is that the use of the match surface enforces the left�right consistency and uniqueness constraints thus improving the algorithms ability to remove invalid matches Experimental results for a number of test stereo pairs show that the new algorithm is capable of identifying and removing a large number of in incorrect matches particularly in the case of occlusions

A constraint to improve the reliability of stereo matching using the rank transform

The rank transform is a non-parametric technique which has been recently proposed for the stereo matching problem. The motivation behind its application to the matching problem is its invariance to certain types of image distortion and noise, as well as its amenability to real-time implementation. This paper derives an analytic expression for the process of matching using the rank transform, and then goes on to derive one constraint which must be satisfied for a correct match. This has been dubbed the rank order constraint or simply the rank constraint. Experimental work has shown that this constraint is capable of resolving ambiguous matches, thereby improving matching reliability. This constraint was incorporated into a new algorithm for matching using the rank transform. This modified algorithm resulted in an increased proportion of correct matches, for all test imagery used.

Estimating the probability of an invalid match for matching using the rank transform

A fundamental problem faced by stereo matching algorithms is the matching or correspondence problem. A wide range of algorithms have been proposed for the correspondence problem. For all matching algorithms, it would be useful to be able to compute a measure of the probability of correctness, or reliability of a match. This paper focuses in particular on one class for matching algorithms, which are based on the rank transform. The interest in these algorithms for stereo matching stems from their invariance to radiometric distortion, and their amenability to fast hardware implementation. This work differs from previous work in that it derives, from first principles, an expression for the probability of a correct match. This method was based on an enumeration of all possible symbols for matching. The theoretical results for disparity error prediction, obtained using this method, were found to agree well with experimental results. However, disadvantages of the technique developed in this chapter are that it is not easily applicable to real images, and also that it is too computationally expensive for practical window sizes. Nevertheless, the exercise provides an interesting and novel analysis of match reliability.

Synthesis of discrete allyl-functional polyethylene glycol for efficient hydrophilization of dihydr (polydimethylsiloxane)

Amphiphilic poly(ethylene glycol)-block-pol (dimethylsiloxane)-block-poly(ethylene glycol)(PEG-block-PDMS block-PEG) triblock copolymers have been successfully prepared via hydrosilylation using discrete and polydisperse PEG of various chain lengths. Facile synthesis of discrete PEG (dPEG) is achieved via systematic tosylation and etherification of lower glycols. Amphiphilicity of the dPEG block-PDMS-block-dPEG triblock copolymer is illustrated by dynamic light scattering (DLS) and measurement of the critical micelle concentration (CMC).

Discrete loading of the foot during hurdle drills: A pilot investigation

Introduction: Lower limb function in hurdling is patently asymmetrical. The lead limb undertakes the preparatory and landing steps while the trail limb contends with the hurdle and recovery steps. Discrete loading profiles of these steps will reflect the asymmetrical function and may provide useful insight into injury mechanisms. A pilot study was undertaken to determine the loading profiles of the hurdle, landing and recovery steps of elite male hurdlers. Equivalent data for steps between hurdles, where the running action is more symmetrical, were used for the purpose of comparison, simultaneously minimising the confounding effect of speed. Methodology: In-shoe pressures were recorded (FScan, 200 Hz) for four elite male hurdlers while they completed a routine hurdle drill at a self-selected fast but sub-race speed. The drill comprised of three consecutive hurdles. Data for the hurdle, landing and recovery steps of the first and second hurdles, along with data for the running steps between hurdles 1 and 2, and 2 and 3, were used for the purpose of analysis. Peak pressures within 1cm2 masks were determined for the hallux, first, central and fifth metatarsals (T1, M1, M2–4 and M5 respectively). Peak pressure (kPa) and loading duration (ms) for the hurdle, landing and recovery steps are reported as a percentage of the respective limb-matched values for between-hurdle steps. Results/discussion: For between-hurdle steps, T1, M1 and M2–4 peak pressures were 312/357, 356/306 and 362/368 kPa, lead/trail limbs respectively. For the hurdle, landing and recovery steps, pressures at T1 and M1 increased. For T1 the increases were in the order of 17%, 36% and 8% (hurdle, landing and recovery steps, respectively) while the corresponding increases at M1 were 7%, 54% and 20%. Pressures at M2–4 were similar for all steps, while M5 loaded erratically. For the between-hurdle steps, the loading durations at T1, M1 and M2–4, were 160/162, 170/142 and 190/191 ms, respectively. For the landing step, loading duration decreased for T1, M1and M2–4 (−8%, −19% and −18%, respectively). In the hurdle step, loading duration decreased for the metatarsals but not for T1. Conclusions: The hurdling action leads to regional pressure increases that act for shorter durations in comparison to the between-hurdle running steps. These changes are most notable at the first metatarsal, a common site of foot injury.

Rapid quantification of methamphetamine: using Attenuated Total Reflectance Fourier Transform Infrared Spectroscopy (ATR-FTIR) and Chemometrics

In Australia and increasingly worldwide, methamphetamine is one of the most commonly seized drugs analysed by forensic chemists. The current well-established GC/MS methods used to identify and quantify methamphetamine are lengthy, expensive processes, but often rapid analysis is requested by undercover police leading to an interest in developing this new analytical technique. Ninety six illicit drug seizures containing methamphetamine (0.1% - 78.6%) were analysed using Fourier Transform Infrared Spectroscopy with an Attenuated Total Reflectance attachment and Chemometrics. Two Partial Least Squares models were developed, one using the principal Infrared Spectroscopy peaks of methamphetamine and the other a Hierarchical Partial Least Squares model. Both of these models were refined to choose the variables that were most closely associated with the methamphetamine % vector. Both of the models were excellent, with the principal peaks in the Partial Least Squares model having Root Mean Square Error of Prediction 3.8, R2 0.9779 and lower limit of quantification 7% methamphetamine. The Hierarchical Partial Least Squares model had lower limit of quantification 0.3% methamphetamine, Root Mean Square Error of Prediction 5.2 and R2 0.9637. Such models offer rapid and effective methods for screening illicit drug samples to determine the percentage of methamphetamine they contain.

Lattice-free models of cell invasion : discrete simulations and travelling waves

Invasion waves of cells play an important role in development, disease and repair. Standard discrete models of such processes typically involve simulating cell motility, cell proliferation and cell-to-cell crowding effects in a lattice-based framework. The continuum-limit description is often given by a reaction–diffusion equation that is related to the Fisher–Kolmogorov equation. One of the limitations of a standard lattice-based approach is that real cells move and proliferate in continuous space and are not restricted to a predefined lattice structure. We present a lattice-free model of cell motility and proliferation, with cell-to-cell crowding effects, and we use the model to replicate invasion wave-type behaviour. The continuum-limit description of the discrete model is a reaction–diffusion equation with a proliferation term that is different from lattice-based models. Comparing lattice based and lattice-free simulations indicates that both models lead to invasion fronts that are similar at the leading edge, where the cell density is low. Conversely, the two models make different predictions in the high density region of the domain, well behind the leading edge. We analyse the continuum-limit description of the lattice based and lattice-free models to show that both give rise to invasion wave type solutions that move with the same speed but have very different shapes. We explore the significance of these differences by calibrating the parameters in the standard Fisher–Kolmogorov equation using data from the lattice-free model. We conclude that estimating parameters using this kind of standard procedure can produce misleading results.

Evidence for the adsorption of molecules at special sites located at copper/zinc oxide interfaces. Part 2. -- A fourier-transform infrared spectroscopy study of methanol adsorption on reduced and oxidised Cu/ZnO/SiO2 catalysts

FTIR spectra are reported of methanol adsorbed at 295 K on ZnO/SiO 2, on reduced Cu/ZnO/SiO2 and on Cu/ZnO/SiO2 which had been preoxidised by exposure to nitrous oxide. Methanol on ZnO/SiO2 gave methoxy species on ZnO and SiO, in addition to both strongly and weakly physisorbed methanol on SiO2. The corresponding adsorption of methanol on reduced Cu/ZnO/SiO2 also gave methoxy species on Cu and a small amount of bridging formate. Reaction of methanol with a reoxidised Cu/ZnO/SiO2 catalyst resulted in an enhanced quantity of methoxy species on Cu. Heating adsorbed species on Cu/ZnO/SiO2 at 393 K led to the loss of methoxy groups on Cu and the concomitant formation of formate species on both ZnO and Cu. The comparable reaction on a reoxidised Cu/ZnO/SiO2 catalyst gave an increased amount of formate species on ZnO and this correlated with an increased quantity of methoxy groups lost from Cu. An explanation is given in terms of adsorption of formate and formaldehyde species at special sites located at the copper/zinc oxide interface.

Evidence for the adsorption of molecules at special sites located at copper/zinc oxide interfaces: Part 1. -- A Fourier-transform infrared study of formic acid and formaldehyde adsorption on reduced and oxidised Cu/ZnO/SiO2 catalysts

Fourier-transform infrared (FTIR) spectra are reported of formic acid and formaldehyde on ZnO/SiO2, reduced Cu/ZnO/SiO2 and reoxidised Cu/ZnO/SiO2 catalyst. Formic acid adsorption on ZnO/SiO2 produced mainly bidentate zinc formate species with a lesser quantity of unidentate zinc formate. Formic acid on reduced Cu/ZnO/SiO2 catalyst resulted not only in the formation of bridging copper formate structures but also in an enhanced amount of formate relative to that for ZnO/SiO2 catalyst. Formic acid on reoxidised Cu/ZnO/SiO2 gave unidentate formate species on copper in addition to zinc formate moieties. The interaction of formaldehyde with ZnO/SiO2 catalyst resulted in the formation of zinc formate species. The same reaction on reduced Cu/ZnO/SiO2 catalyst gave bridging formate on copper and a remarkable increase in the quantity of formate species associated with the zinc oxide. Adsorption of formaldehyde on a reoxidised Cu/ZnO/SiO2 catalyst produced bridging copper formate and again an apparent increase in the concentration of zinc formate species. An explanation in terms of the adsorption of molecules at special sites located at the interface between copper and zinc oxide is given.

A combined temperature-programmed reaction spectroscopy and Fourier-transform infrared spectroscopy study of CO2/H2 and CO/CO2/H2 interactions with model ZnO/SiO2, Cu/SiO2 and Cu/ZnO/SiO2 methanol-synthesis catalysts

The reaction of CO2 and H2 with ZnO/SiO2 catalyst at 295 K gave predominantly hydrogencarbonate on zinc oxide and a small quantity of formate was evolved after heating at 393 K. Elevation of the reaction temperature to 503 K enhanced the rate of formation of zinc formate species. Significantly these formate species decomposed at 573 K almost entirely to CO2 and H2. Even after exposure of CO2-H2 or CO-CO2-H2 mixtures to highly defected ZnO/SiO2 catalyst, the formate species produced still decomposed to give CO2 and H2. It was concluded that carboxylate species which were formed at oxygen anion vacancies on polar Zn planes were not significantly hydrogenated to formate. Consequently it was proposed that the non-polar planes on zinc oxide contained sites which were specific for the synthesis of methanol. The interaction of CO2 and H2 with reduced Cu/ZnO/SiO2 catalyst at 393 K gave copper formate species in addition to substantial quantities of formate created at interfacial sites between copper and zinc oxide. It was deduced that interfacial formate species were produced from the hydrogenation of interfacial bidentate carbonate structures. The relevance of interfacial formate species in the methanol synthesis reaction is discussed. Experiments concerning the reaction of CO2-H2 with physical mixtures of Cu/SiO2 and ZnO/SiO2 gave results which were simply characteristic of the individual components. By careful consideration of previous data a detailed proposal regarding the role of spillover hydrogen is outlined. Admission of CO to a gaseous CO2-H2 feedstock resulted in a considerably diminished amount of formate species on copper. This was ascribed to a combination of over-reduction of the surface and site-blockage.