Model selection for convolutive ICA with an application to spatiotemporal analysis of EEG

We present a new algorithm for maximum likelihood convolutive independent component analysis (ICA) in which components are unmixed using stable autoregressive filters determined implicitly by estimating a convolutive model of the mixing process. By introducing a convolutive mixing model for the components, we show how the order of the filters in the model can be correctly detected using Bayesian model selection. We demonstrate a framework for deconvolving a subspace of independent components in electroencephalography (EEG). Initial results suggest that in some cases, convolutive mixing may be a more realistic model for EEG signals than the instantaneous ICA model.     

tweezercalib 2.0: Faster version of MatLab package for precise calibration of optical tweezers

We present a vectorized version of the MatLab (MathWorks Inc.) package tweezercalib for calibration of optical tweezers with precision. The calibration is based on the power spectrum of the Brownian motion of a dielectric bead trapped in the tweezers. Precision is achieved by accounting for a number of factors that affect this power spectrum, as described in vs. 1 of the package [I.M. Toli?-Nørrelykke, K. Berg-Sørensen, H. Flyvbjerg, Matlab program for precision calibration of optical tweezers, Comput. Phys. Comm. 159 (2004) 225-240]. The graphical user interface allows the user to include or leave out each of these factors. Several “health tests” are applied to the experimental data during calibration, and test results are displayed graphically. Thus, the user can easily see whether the data comply with the theory used for their interpretation. Final calibration results are given with statistical errors and covariance matrix.

New version program summary

Title of program: tweezercalibCatalogue identifier: ADTV_v2_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/ADTV_v2_0Program obtainable from: CPC Program Library, Queen's University of Belfast, N. IrelandReference in CPC to previous version: I.M. Toli?-Nørrelykke, K. Berg-Sørensen, H. Flyvbjerg, Comput. Phys. Comm. 159 (2004) 225Catalogue identifier of previous version: ADTVDoes the new version supersede the original program: YesComputer for which the program is designed and others on which it has been tested: General computer running MatLab (Mathworks Inc.)Operating systems under with the program has been tested: Windows2000, Windows-XP, LinuxProgramming language used: MatLab (Mathworks Inc.), standard licenseMemory required to execute with typical data: Of order four times the size of the data fileHigh speed storage required: noneNo. of lines in distributed program, including test data, etc.: 135 989No. of bytes in distributed program, including test data, etc.: 1 527 611Distribution format: tar. gzNature of physical problem: Calibrate optical tweezers with precision by fitting theory to experimental power spectrum of position of bead doing Brownian motion in incompressible fluid, possibly near microscope cover slip, while trapped in optical tweezers. Thereby determine spring constant of optical trap and conversion factor for arbitrary-units-to-nanometers for detection system.Method of solution: Elimination of cross-talk between quadrant photo-diode's output channels for positions (optional). Check that distribution of recorded positions agrees with Boltzmann distribution of bead in harmonic trap. Data compression and noise reduction by blocking method applied to power spectrum. Full accounting for hydrodynamic effects: Frequency-dependent drag force and interaction with nearby cover slip (optional). Full accounting for electronic filters (optional), for “virtual filtering” caused by detection system (optional). Full accounting for aliasing caused by finite sampling rate (optional). Standard non-linear least-squares fitting. Statistical support for fit is given, with several plots facilitating inspection of consistency and quality of data and fit.Summary of revisions: A faster fitting routine, adapted from [J. Nocedal, Y.x. Yuan, Combining trust region and line search techniques, Technical Report OTC 98/04, Optimization Technology Center, 1998; W.H. Press, B.P. Flannery, S.A. Teukolsky, W.T. Vetterling, Numerical Recipes. The Art of Scientific Computing, Cambridge University Press, Cambridge, 1986], is applied. It uses fewer function evaluations, and the remaining function evaluations have been vectorized. Calls to routines in Toolboxes not included with a standard MatLab license have been replaced by calls to routines that are included in the present package. Fitting parameters are rescaled to ensure that they are all of roughly the same size (of order 1) while being fitted. Generally, the program package has been updated to comply with MatLab, vs. 7.0, and optimized for speed.Restrictions on the complexity of the problem: Data should be positions of bead doing Brownian motion while held by optical tweezers. For high precision in final results, data should be time series measured over a long time, with sufficiently high experimental sampling rate: The sampling rate should be well above the characteristic frequency of the trap, the so-called corner frequency. Thus, the sampling frequency should typically be larger than 10 kHz. The Fast Fourier Transform used works optimally when the time series contain n² data points, and long measurement time is obtained with n>12-15. Finally, the optics should be set to ensure a harmonic trapping potential in the range of positions visited by the bead. The fitting procedure checks for harmonic potential.Typical running time: SecondsUnusual features of the program: NoneReferences: The theoretical underpinnings for the procedure are found in [K. Berg-Sørensen, H. Flyvbjerg, Power spectrum analysis for optical tweezers, Rev. Sci. Ins. 75 (2004) 594-612].     

Molecular dynamics simulations of oscillatory flows in microfluidic channels

In this paper we apply the direct non-equilibrium molecular dynamics technique to oscillatory flows of fluids in microscopic channels. Initially, we show that the microscopic simulations resemble the macroscopic predictions based on the Navier–Stokes equation very well for large channel width, high density and low temperature. Further simulations for high temperature and low density show that the non-slip boundary condition traditionally used in the macroscopic equation is greatly compromised when the fluid–wall interactions are the same as the fluid–fluid interactions. Simulations of a system with very narrow channel width confirm earlier findings of Poiseuille flow, namely, that the velocity profiles are modulated. We find that these modulations cannot be explained by the local area density model.

Interactive display of isosurfaces with global illumination

In many applications, volumetric data sets are examined by displaying isosurfaces, surfaces where the data, or some function of the data, takes on a given value. Interactive applications typically use local lighting models to render such surfaces. This work introduces a method to precompute or lazily compute global illumination to improve interactive isosurface renderings. The precompiled illumination resides in a separate volume and includes direct light, shadows, and intersections. Using this volume, interactive globally illuminated renderings of isosurfaces become feasible while still allowing dynamic manipulation of lighting, viewpoint and isovalue.     

Heroin related overdose problems

During the first ten months of 1992 the Oslo Ambulance Department registered 716 incidences of assumed drug-related intoxications. 80% happened in down town Oslo. 19 cases of asystoly were recorded, 13 of the patients recovered after treatment, without sequelae. Five of these patients left the location after emergency help and they refused hospitalization. 432 of the patients were unconscious when the ambulance personal arrived, 472 were treated with naloxone both by the intramuscular and the intravenous route. Most of the persons refused further observation. A team of specially trained out reach workers offers help after acute medical treatment by means of "streetwork". The intervention is directed at addicts who have experienced an overdose.     

Low temperature hopping conduction in neutron transmutation doped isotopically enriched70Ge:Ga single crystals

The temperature dependence of variable range hopping resistivity in neutron transmutation doped (NTD) isotopically enriched⁷⁰Ge:Ga samples is reported. Five samples with compensation ratios K less than 0.001 and Ga concentrations between 3×10¹⁶ and 1.77×10¹⁷ cm^–3 were studied. All samples investigated show the ln T^–1/2 dependence in the temperature range below 1.5K. As thermistor materials NTD⁷⁰Ge:Ga samples are found to have more than factor of two higher sensitivity than commonly used natural NTD Ge in the temperature range between 0.2K and 1K. Our results are compared with theoretical predictions for variable range hopping conduction.     

Bicriterion cluster analysis

Cluster analysis is concerned with the problem of partitioning a given set of entities into homogeneous and well-separated subsets called clusters. The concepts of homogeneity and of separation can be made precise when a measure of dissimilarity between the entities is given. Let us define the diameter of a partition of the given set of entities into clusters as the maximum dissimilarity between any pair of entities in the same cluster and the split of a partition as the minimum dissimilarity between entities in different clusters. The problems of determining a partition into a given number of clusters with minimum diameter (i.e., a partition of maximum homogeneity) or with maximum split (i.e., a partition of maximum separation) are first considered. It is shown that the latter problem can be solved by the classical single-link clustering algorithm, while the former can be solved by a graph-theoretic algorithm involving the optimal coloration of a sequence of partial graphs, described in more detail in a previous paper. A partition into a given number of clusters will be called efficient if and only if there exists no partition into at most the same number of clusters with smaller diameter and not smaller split or with larger split and not larger diameter. Two efficient partitions are called equivalent if and only if they have the same values for the split and for the diameter.     

A VNS heuristic for escaping local extrema entrapment in normalized cut clustering

Normalized cut is one of the most popular graph clustering criteria. The main approaches proposed for its resolution are spectral clustering methods and a multilevel approach of Dhillon et al. (TPAMI 29:1944–1957, 2007), called graclus. Their aim is to obtain good solutions in a small amount of time for large instances. Metaheuristics are general frameworks for stochastic searches often employed in global optimization to improve the solutions obtained by other heuristics. Variable neighborhood search (VNS) is a metaheuristic which exploits systematically the idea of neighborhood change during the search. In this paper, we propose a VNS heuristic for normalized cut segmentation. Computational experiments show that in most cases this VNS heuristic improves significantly, and in moderate time, the solutions obtained by the current state-of-the-art algorithms, i.e., graclus and a spectral method proposed by Yu and Shi (ICCV, 2003).     

AMPHIBIAN OCCUPANCY IN FLOOD‐CREATED AND EXISTING WETLANDS OF THE LOWER MISSOURI RIVER ALLUVIAL VALLEY

During the previous century, the wetland area in the lower Missouri River alluvial valley was reduced by 39% because of river channelization and bank stabilization projects. The Great Flood of 1993 reversed the trend of wetland loss by creating 466 new wetlands in the alluvial valley between Kansas City and St. Louis, Missouri. We estimated amphibian occupancy, detection probability and number of species exhibiting evidence of reproduction in eight flood‐created and 16 pre‐flood existing wetlands from 1996 through 1998. We also evaluated whether hydroperiod (the number of days any water was present in a wetland from 20 February through 31 August) and distance to river predicted those values. Detection probabilities for adult amphibian species were relatively constant across years and ranged from 0.013 [Great Plains toad (Anaxyrus cognatus)] to 0.280 [Woodhouse's toad (Anaxyrus woodhousii woodhousii)]. Occupancy of adult amphibians differed across years and was not correlated with habitat features. Estimated occupancy probabilities for amphibian species ranged from 0.126 [Plains spadefoot (Spea bombifrons)] to 0.896 [boreal chorus frog (Pseudacris maculata)]. Almost double the number of amphibian species showed evidence of reproduction in existing wetlands (wetlands created before the Great Flood of 1993) when compared with that in flood‐created wetlands. Similarly, temporary wetlands had nearly double the number of amphibian species showing evidence of reproduction when compared with permanent wetlands. Finally, the highest number of species showed evidence of reproduction in wetlands with spring–summer hydroperiods between 135 and 140 days. All these relationships suggest that the invasion and persistence of predators in wetlands negatively influence amphibian reproduction. If the Missouri River is allowed to reconnect with the alluvial valley, more predators may be introduced into wetlands, leading to reduced amphibian occupancy and reproduction. However, this connection will not likely occur over the entire alluvial valley and, therefore, should not adversely impact amphibians that find refuge in higher‐elevation, non‐connected regions of the alluvial valley. Copyright © 2011 John Wiley & Sons, Ltd.     

Implementation of physical coordination training and cognitive behavioural training interventions at cleaning workplaces--secondary analyses of a randomised controlled trial

This study evaluates the implementation of physical coordination training (PCT) and cognitive behavioural training (CBTr) interventions in a randomised controlled trial at nine cleaners' workplaces. Female cleaners (n = 294) were randomised into a PCT, a CBTr or a reference (REF) group. Both 12-week interventions were performed in groups guided by an instructor. Records were kept on intervention dose (adherence) unanticipated events at the work place (context) and quality of intervention delivery (fidelity). Participant adherence was 37% in the PCT and 49% in the CBTr interventions. Optimal implementation was reached by only 6% in PCT and 42% in the CBTr. Analysis of the barriers to successful implementation indicated that the intervention process is sensitive to unanticipated events. In order to succeed in improving the health of high-risk populations such as cleaners and to correctly interpret intervention effects, more research on implementation is needed. Trial registration: ISRCTN96241850. PRACTITIONER SUMMARY: Both physical coordination training and cognitive behavioural training are potential effective workplace interventions among low educated job groups with high physical work demands. However, thorough consideration should be given to feasibility in the design of interventions. The optimal intervention should be tailored to closely match the implementation context and be robust and flexible to minimise susceptibility to changes in work organisation.