Cross-spectral analysis of physiological tremor and muscle activity. I. Theory and application to unsynchronized EMG

We investigate the relationship between the extensor electromyogram (EMG) and tremor time series in physiological hand tremor by cross-spectral analysis. Special attention is directed to the phase spectrum and the effects of observational noise. We calculate the theoretical phase spectrum for a second order linear stochastic process and compare the results to measured tremor data recorded from subjects who did not show a synchronized EMG activity in the corresponding extensor muscle. The results show that physiological tremor is well described by the proposed model and that the measured EMG represents a Newtonian force by which the muscle acts on the hand.Comment: 9 pages, 6 figures, to appear in Biological Cybernetic

Linear and nonlinear time series analysis of the black hole candidate Cygnus X-1

We analyze the variability in the X-ray lightcurves of the black hole candidate Cygnus X-1 by linear and nonlinear time series analysis methods. While a linear model describes the over-all second order properties of the observed data well, surrogate data analysis reveals a significant deviation from linearity. We discuss the relation between shot noise models usually applied to analyze these data and linear stochastic autoregressive models. We debate statistical and interpretational issues of surrogate data testing for the present context. Finally, we suggest a combination of tools from linear andnonlinear time series analysis methods as a procedure to test the predictions of astrophysical models on observed data.Comment: 15 pages, to appear in Phys. Rev.

Talking quiescence: a rigorous theory that supports parallel composition, action hiding and determinisation

The notion of quiescence - the absence of outputs - is vital in both behavioural modelling and testing theory. Although the need for quiescence was already recognised in the 90s, it has only been treated as a second-class citizen thus far. This paper moves quiescence into the foreground and introduces the notion of quiescent transition systems (QTSs): an extension of regular input-output transition systems (IOTSs) in which quiescence is represented explicitly, via quiescent transitions. Four carefully crafted rules on the use of quiescent transitions ensure that our QTSs naturally capture quiescent behaviour. We present the building blocks for a comprehensive theory on QTSs supporting parallel composition, action hiding and determinisation. In particular, we prove that these operations preserve all the aforementioned rules. Additionally, we provide a way to transform existing IOTSs into QTSs, allowing even IOTSs as input that already contain some quiescent transitions. As an important application, we show how our QTS framework simplifies the fundamental model-based testing theory formalised around ioco.Comment: In Proceedings MBT 2012, arXiv:1202.582

Reporting of prognostic studies of tumour markers: a review of published articles in relation to REMARK guidelines

Background: Poor reporting compromises the reliability and clinical value of prognostic tumour marker studies. We review articles to assess the reporting of patients and events using REMARK guidelines, at the time of guideline publication. Methods: We sampled 50 prognostic tumour marker studies from higher impact cancer journals between 2006 and 2007. The inclusion criteria were cancer; focus on single biological tumour marker; survival analysis; multivariable analysis; and not gene array or proteomic data. Articles were assessed for the REMARK profile and other REMARK guideline items. We propose a reporting aid, the REMARK profile, motivated by the CONSORT flowchart. Results: In 50 studies assessed for the REMARK profile, the number of eligible patients (56% of articles), excluded patients (54%) and patients in analyses (98%) was reported. Only 50% of articles reported the number of outcome events. In multivariable analyses, 54% and 30% of articles reported patient and event numbers for all variables. Of the studies, 66% used archival samples, indicating a potentially biased patient selection. Only 36% of studies reported clearly defined outcomes. Conclusions: Good reporting is critical for the interpretability and clinical applicability of prognostic studies. Current reporting of key information, such as the number of outcome events in all patients and subgroups, is poor. Use of the REMARK profile would greatly improve reporting and enhance prognostic research

Phase synchronization from noisy univariate signals

We present methods for detecting phase synchronization of two unidirectionally coupled, self-sustained noisy oscillators from a signal of the driven oscillator alone. One method detects soft, another hard phase locking. Both are applied to the problem of detecting phase synchronization in von Karman vortex flow meters.Comment: 4 pages, 4 figure

A Multi-Wavelength Study of Sgr A*: The Role of Near-IR Flares in Production of X-ray, Soft γ\gamma-ray and Sub-millimeter Emission

(abridged) We describe highlights of the results of two observing campaigns in 2004 to investigate the correlation of flare activity in Sgr A* in different wavelength regimes, using a total of nine ground and space-based telescopes. We report the detection of several new near-IR flares during the campaign based on {\it HST} observations. The level of near-IR flare activity can be as low as $\sim0.15$ mJy at 1.6 $\mu$m and continuous up to about 40% of the total observing time. Using the NICMOS instrument on the {\it HST}, the {\it XMM-Newton} and CSO observatories, we also detect simultaneous bright X-ray and near-IR flare in which we observe for the first time correlated substructures as well as simultaneous submillimeter and near-IR flaring. X-ray emission is arising from the population of near-IR-synchrotron-emitting relativistic particles which scatter submillimeter seed photons within the inner 10 Schwarzschild radii of Sgr A* up to X-ray energies. In addition, using the inverse Compton scattering picture, we explain the high energy 20-120 keV emission from the direction toward Sgr A*, and the lack of one-to-one X-ray counterparts to near-IR flares, by the variation of the magnetic field and the spectral index distributions of this population of nonthermal particles. In this picture, the evidence for the variability of submillimeter emission during a near-IR flare is produced by the low-energy component of the population of particles emitting synchrotron near-IR emission. Based on the measurements of the duration of flares in near-IR and submillimeter wavelengths, we argue that the cooling could be due to adiabatic expansion with the implication that flare activity may drive an outflow.Comment: 48 pages, 12 figures, ApJ (in press

Characterization of the Crab Pulsar's Timing Noise

We present a power spectral analysis of the Crab pulsar's timing noise, mainly using radio measurements from Jodrell Bank taken over the period 1982-1989. The power spectral analysis is complicated by nonuniform data sampling and the presence of a steep red power spectrum that can distort power spectra measurement by causing severe power ``leakage''. We develop a simple windowing method for computing red noise power spectra of uniformly sampled data sets and test it on Monte Carlo generated sample realizations of red power-law noise. We generalize time-domain methods of generating power-law red noise with even integer spectral indices to the case of noninteger spectral indices. The Jodrell Bank pulse phase residuals are dense and smooth enough that an interpolation onto a uniform time series is possible. A windowed power spectrum is computed revealing a periodic or nearly periodic component with a period of about 568 days and a 1/f^3 power-law noise component with a noise strength of 1.24 +/- 0.067 10^{-16} cycles^2/sec^2 over the analysis frequency range 0.003 - 0.1 cycles/day. This result deviates from past analyses which characterized the pulse phase timing residuals as either 1/f^4 power-law noise or a quasiperiodic process. The analysis was checked using the Deeter polynomial method of power spectrum estimation that was developed for the case of nonuniform sampling, but has lower spectral resolution. The timing noise is consistent with a torque noise spectrum rising with analysis frequency as f implying blue torque noise, a result not predicted by current models of pulsar timing noise. If the periodic or nearly periodic component is due to a binary companion, we find a companion mass > 3.2 Earth masses.Comment: 53 pages, 9 figures, submitted to MNRAS, abstract condense

Neuron dynamics in the presence of 1/f noise

Interest in understanding the interplay between noise and the response of a non-linear device cuts across disciplinary boundaries. It is as relevant for unmasking the dynamics of neurons in noisy environments as it is for designing reliable nanoscale logic circuit elements and sensors. Most studies of noise in non-linear devices are limited to either time-correlated noise with a Lorentzian spectrum (of which the white noise is a limiting case) or just white noise. We use analytical theory and numerical simulations to study the impact of the more ubiquitous "natural" noise with a 1/f frequency spectrum. Specifically, we study the impact of the 1/f noise on a leaky integrate and fire model of a neuron. The impact of noise is considered on two quantities of interest to neuron function: The spike count Fano factor and the speed of neuron response to a small step-like stimulus. For the perfect (non-leaky) integrate and fire model, we show that the Fano factor can be expressed as an integral over noise spectrum weighted by a (low pass) filter function. This result elucidates the connection between low frequency noise and disorder in neuron dynamics. We compare our results to experimental data of single neurons in vivo, and show how the 1/f noise model provides much better agreement than the usual approximations based on Lorentzian noise. The low frequency noise, however, complicates the case for information coding scheme based on interspike intervals by introducing variability in the neuron response time. On a positive note, the neuron response time to a step stimulus is, remarkably, nearly optimal in the presence of 1/f noise. An explanation of this effect elucidates how the brain can take advantage of noise to prime a subset of the neurons to respond almost instantly to sudden stimuli.Comment: Phys. Rev. E in pres

Inzicht en toezicht: controle in de kennissamenleving

Nieuwe technieken maken het opslaan en verwerken van informatie eenvoudiger. Sterker\ud nog: ze maken de weg vrij om gegevens vast te leggen voordat duidelijk is welk inzicht\ud we eigenlijk nastreven. Toezicht wint zo terrein ten opzichte van de kennisverwerving.\ud \ud Voorbeelden van toenemend toezicht zijn elektronische dossiers over burgers,\ud informatievergaring door de politie en in private initiatieven zoals Google Earth. Deze\ud systemen maken het mogelijk om op grote schaal verbanden te ontdekken en afwijkingen op\ud te sporen, zonder dat er duidelijke kennisvragen aan vooraf gaan.\ud \ud De relatie tussen inzicht en toezicht is bepalend voor hoe we met de verzamelde kennis\ud omgaan. Ligt de nadruk daarbij op inzicht, overzicht, toezicht, of nog iets anders? Dit\ud Jaarboek Kennissamenleving 2010 stelt kritische vragen bij de alomtegenwoordige blik van\ud toezicht. Het draagt ideeën aan om deze in de toekomst anders te richten, zodat inzicht en\ud toezicht in balans kunnen blijven