Limitations in the use and interpretation of continuous relative phase

Continuous relative phase (CRP), a variable used to quantify intersegmental coordination, is difficult to interpret if care is not taken regarding the assumptions and limitations of the measure. Specifically, CRP is often interpreted as a higher resolution form of discrete relative phase (DRP). DRP, however, yields information regarding the relative dispersion of events in oscillatory signals while CRP describes their relationship in a higher order phase-plane domain. In this paper we address issues surrounding the calculation of CRP and suggest a new interpretation based on the aforementioned methodological issues. Through the use of test signals, with known properties, it was found that the CRP information will be arbitrary if no normalization procedures are used to account for frequency differences in the component oscillators. In addition, signals with non-sinusoidal trajectories will produce patterns in CRP that are not equivalent to discrete relative phase (DRP) measures. The implications of these issues are discussed.     

Water velocity shapes fish movement behavior

Stream and river ecosystems present fluvial fishes with a dynamic energy landscape because moving water generates heterogeneous flow fields that are rarely static in space and time. Fish movement behavior should be consistent with conserving energy in these dynamic flowing environments, but little evidence supporting this hypothesis exists. Here, we tested experimentally whether three general movement behaviors—against the current, with the current, or holding position (i.e., staying in one position and location)—were performed in a way consistent with minimizing the cost of swimming in a heterogeneous flow field. We tested the effects of water velocity on movement behavior across three age classes (0, 1, and 5 years) of two different fluvial specialist fishes, the pallid sturgeon (Scaphirhynchus albus) and shovelnose sturgeon (Scaphirhynchus platorynchus). Individuals from the three age classes were exposed to a continuous and dynamic velocity field ranging from 0.02 to 0.53 m s⁻¹, which represented natural benthic flow regimes occupied by these species in rivers. Both sturgeon species exhibited the same pattern with regard to their tendency to hold position, move upstream, or move downstream. Moving downstream was positively associated with velocity for all age groups. Moving upstream was inversely related to velocity for young fish, but as the fish aged, moving upstream was not related to water velocity. The oldest fish (age 5) moved upstream more frequently compared to the younger age classes. Holding position within a water current was the most frequent behavior and occurred with similar probability across the range of experimental velocity for youngest fish (age 0), but was inversely related to velocity in older fish. Our experiment across age classes suggests that the suite of swimming behaviors exhibited by fluvial specialists might have evolved to mitigate the energetic costs of complex energy landscapes generated by moving water to ultimately maximize net energy gain.     

Hilbert-Huang transform for analysis of heart rate variability in cardiac health

This paper introduces a modified technique based on Hilbert-Huang transform (HHT) to improve the spectrum estimates of heart rate variability (HRV). In order to make the beat-to-beat (RR) interval be a function of time and produce an evenly sampled time series, we first adopt a preprocessing method to interpolate and resample the original RR interval. Then, the HHT, which is based on the empirical mode decomposition (EMD) approach to decompose the HRV signal into several monocomponent signals that become analytic signals by means of Hilbert transform, is proposed to extract the features of preprocessed time series and to characterize the dynamic behaviors of parasympathetic and sympathetic nervous system of heart. At last, the frequency behaviors of the Hilbert spectrum and Hilbert marginal spectrum (HMS) are studied to estimate the spectral traits of HRV signals. In this paper, two kinds of experiment data are used to compare our method with the conventional power spectral density (PSD) estimation. The analysis results of the simulated HRV series show that interpolation and resampling are basic requirements for HRV data processing, and HMS is superior to PSD estimation. On the other hand, in order to further prove the superiority of our approach, real HRV signals are collected from seven young health subjects under the condition that autonomic nervous system (ANS) is blocked by certain acute selective blocking drugs: atropine and metoprolol. The high-frequency power/total power ratio and low-frequency power/high-frequency power ratio indicate that compared with the Fourier spectrum based on principal dynamic mode, our method is more sensitive and effective to identify the low-frequency and high-frequency bands of HRV.     

Three-dimensional head angular velocity detection from otolith afferent signals

Afferent signals from the otolith organs can produce compensatory eye position and velocity signals which has been described as linear vestibulo-ocular reflex (LVOR). The afferent otolith signals carry information about head orientation and changes of head orientation relative to gravity. A head orientation (tilt) related position signal can be obtained from population vector coding of tonic otolith afferent signals during static or dynamic head tilts, which in turn could produce compensatory eye position signals in the LVOR. On the other hand, eye angular velocity signals may be extracted, as proposed in this study, from the population response of tilt-velocity sensitive otolith afferents. Such afferents are shown to encode instantaneous head orientation relative to gravity at onset of a head movement and, as the movement continues, the projection of head angular velocity onto the earth-horizontal plane, indicating the instantaneous direction of movement relative to gravity. Angular velocity components along the earth-vertical direction which are not directly encoded by otolith afferents can be detected by central signal processing. Central reconstruction of 3D head angular velocity allows to obtain information about absolute head orientation in space even in the absence of semicircular canal related information. Such information is important for generating compensatory eye movements as well as for dynamic control of posture.     

Estimation of muscle fiber conduction velocity of doublet discharges

A doublet is defined as two consecutive discharges of a motor unit occurring at short time interval between each other (e.g., <20 ms). In this paper, we propose a method for the estimation of muscle fiber conduction velocity (CV) from two partly overlapping action potentials generated by the same motor units. The method is based on the minimization of the mean square error between time-filtered versions of two surface EMG signals recorded along the direction of muscle fibers. The minimization is performed over the filter parameters that define the two propagation delays. The method was tested on simulated and experimental signals. Simulation results showed that the method is only in some cases superior to the simpler peak approach, due to limitations in the ideal model used for the algorithm development. However, application to experimental signals that mimic doublet motor unit discharges showed a substantial improvement in estimation quality of the new method with respect to the peak method.     

Spectral Analysis of Binary Time Series: Square Waves vs. Sinusoidal Functions

The analysis of signals consisting of discrete and irregular data causes methodological problems for the Fourier spectral Analysis: Since it is based on sinusoidal functions, rectangular signals with unequal periodicities cannot easily be replicated. The Walsh spectral Analysis is based on the so called "Walsh functions", a complete set of orthonormal, rectangular waves and thus seems to be the method of choice for analysing signals consisting of binary or ordinal data. The paper compares the Walsh spectral analysis and the Fourier spectral analysis on the basis of simulated and real binary data sets of various length. Simulated data were derived from signals with defined cyclic patterns that were noised by randomly generated signals of the same length. The Walsh and Fourier spectra of each set were determined and up to 25% of the periodogram coefficients were utilized as input for an inverse transform. Mean square approximation error (MSE) was calculated for each of the series in order to compare the goodness of fit between the original and the reconstructed signal. The same procedure was performed with real data derived from a behavioral observation in pigs. The comparison of the two methods revealed that, in the analysis of discrete and binary time series, Walsh spectral analysis is the more appropriate method, if the time series is rather short. If the length of the signal increases, the difference between the two methods is less substantial.     

Spectral Analysis of Binary Time Series: Square Waves vs. Sinusoidal Functions

The analysis of signals consisting of discrete and irregular data causes methodological problems for the Fourier spectral Analysis: Since it is based on sinusoidal functions, rectangular signals with unequal periodicities cannot easily be replicated. The Walsh spectral Analysis is based on the so called "Walsh functions", a complete set of orthonormal, rectangular waves and thus seems to be the method of choice for analysing signals consisting of binary or ordinal data. The paper compares the Walsh spectral analysis and the Fourier spectral analysis on the basis of simulated and real binary data sets of various length. Simulated data were derived from signals with defined cyclic patterns that were noised by randomly generated signals of the same length. The Walsh and Fourier spectra of each set were determined and up to 25% of the periodogram coefficients were utilized as input for an inverse transform. Mean square approximation error (MSE) was calculated for each of the series in order to compare the goodness of fit between the original and the reconstructed signal. The same procedure was performed with real data derived from a behavioral observation in pigs. The comparison of the two methods revealed that, in the analysis of discrete and binary time series, Walsh spectral analysis is the more appropriate method, if the time series is rather short. If the length of the signal increases, the difference between the two methods is less substantial.     

Ecology models and newton vector fields

We generate phase plane portraits of autonomous differential equations by constructing integrals of motion—even in nonconservative systems, where no integrals produced by quadrature exist. These phase plane portraits are produced with no propagated error. Phase portraits do not usually display the parametrization along orbits. The use of high resolution color graphics allows us to display this parametrization conveniently. We examine three ecology models—the logistic equation in the complex domain, and the predator-prey and competing species equations in the real domain.     

Hilbert-Huang transform in MHD plasma diagnostics

A new method for processing experimental data from MHD diagnostics is discussed that provides a more detailed study of the dynamics of large-scale MHD instabilities. The method is based on the Hilbert-Huang transform method and includes an empirical mode decomposition algorithm, which is used to decompose the experimental MHD diagnostic signals into a set of frequency-and amplitude-modulated harmonics in order to construct the time evolutions of the amplitudes and frequencies of these harmonics with the help of the Hilbert transform. The method can also be applied to analyze data from other diagnostics that measure unsteady oscillating signals.     

Comparison of the SCoRE and HA methods for locating in vivo the glenohumeral joint centre

A recent paper has described a new functional method, the symmetrical centre of rotation (SCoRE), for locating joint centre position [Ehrig, R.M., Taylor, W.R., Duda, G.N., Heller, M.O., 2006. A survey of formal methods for determining the centre of rotation of ball joints. Journal of Biomechanics 39 (15), 2798-2809]. For in vitro analyses, the SCoRE method showed better precision than helical axis (HA) or sphere fitting methods. Despites HA determination is very sensitive to small angular velocity, the International Society of Biomechanics has recommended to use HA for locating the glenohumeral joint centre. This paper aims at comparing the SCoRE method with the HA method for locating in vivo the glenohumeral joint centre according to the movement characteristics. Nine subjects performed 10 cycles of three different movements at two different velocities. For each test (combination of movements) the location of the centre of rotation was estimated with both methods (SCoRE and HA). Analyses focused on the 3D location of the glenohumeral joint centre and on the repeatability of location (standard deviation). This study showed that SCoRE and HA methods yielded the same GH location. Nevertheless, with SCoRE method, the location of the glenohumeral joint centre was different according to the test. This study evidenced that the SCoRE method was more precise than HA method (error of 3 mm versus 4.6 mm) and that the GH location with the SCoRE method was not affected by movements with slow velocities.     

Discrimination of Velocity and Spectral Composition of a Sound Stimulus during Its Movement in the Vertical Plane

Differential sensitivity to the velocity of a sound source image in the vertical plane was studied by using two signals with different spectral bandwidths: 0.25–4 kHz (signal 1) and 4–12.5 kHz (signal 2). Five subjects were tested. Sequential switching of loudspeakers with similar frequency characteristics simulated movement of a sound source. Differential velocity thresholds were determined for two reference velocities: 58 and 115°/s. Significant differences in the absolute values of these thresholds were found for signals with different spectral compositions. The threshold for signal 1 (0.25–4.0 kHz) exceeded that for signal 2 (4.0–12.5 kHz) twofold at 58°/s and 1.6-fold at 115°/s.     

Compensation for changing motor uncertainty

When movement outcome differs consistently from the intended movement, errors are used to correct subsequent movements (e.g., adaptation to displacing prisms or force fields) by updating an internal model of motor and/or sensory systems. Here, we examine changes to an internal model of the motor system under changes in the variance structure of movement errors lacking an overall bias. We introduced a horizontal visuomotor perturbation to change the statistical distribution of movement errors anisotropically, while monetary gains/losses were awarded based on movement outcomes. We derive predictions for simulated movement planners, each differing in its internal model of the motor system. We find that humans optimally respond to the overall change in error magnitude, but ignore the anisotropy of the error distribution. Through comparison with simulated movement planners, we found that aimpoints corresponded quantitatively to an ideal movement planner that updates a strictly isotropic (circular) internal model of the error distribution. Aimpoints were planned in a manner that ignored the direction-dependence of error magnitudes, despite the continuous availability of unambiguous information regarding the anisotropic distribution of actual motor errors.     

Identification of the structure of large-scale MHD perturbations in a tokamak from Mirnov signals

An algorithm is proposed that allows one to identify the MHD mode structure in toroidal plasmas by processing signals from Mirnov probes measuring plasma MHD activity. The algorithm differs fundamentally from the diagnostic methods presently used in tokamaks, being simpler and more efficient. The algorithm is based on constructing an analytic signal using the Hilbert transformation of the Mirnov signals at a given instant. The phase and amplitude dependences obtained take into account the toroidal effects and allow one to determine the number and amplitude of the excited MHD mode. The algorithm was approbated with both test signals and actual signals from MHD diagnostics in the T-10 tokamak. It is demonstrated that the algorithm can be used to analyze single-mode MHD instabilities in toroidal plasmas.     

基于EMD的左右手运动想象脑电信号分析研究

基于经验模态分解（EMD）理论,提出一种左右手运动想象脑电信号分析方法。首先利用时间窗对脑电信号数据进行划分,对每段数据通过经验模态分解法将其分解为一组固有模态函数IMF,提取主要信号所在的IMF层去除信号中的噪声。对含有主要信号的几层IMF进行Hilbert变换,得到瞬时频率与对应的瞬时幅值。再提取左右手想象的特定频段mu节律和beta节律的能量信号作为特征,分别利用支持向量机（SVM）和Fisher进行了分类比较。对EMD和小波包在去噪和特征提取进行了比较。结果表明,EMD是一种很有效的去噪方法,经过EMD分解后提取的能量信号在区分左右手想象上更具有优势,识别率高。     

Evoked potentials of the cat inferior colliculus to acoustic stimuli simulating sound source movement with different velocities in opposite directions

Amplitude changes of inferior colliculus evoked potentials (EPs) in anaesthetized adult cats were studied under presentation of acoustic stimuli simulating both azimuth-moving and stationary sound source. The movement was simulated with gradual changes of interaural time delay between binaurally presented click trains. It was shown that the amplitude of EPs elicited by "moving" signals depended on the velocity of movement. Amplitude differences between EPs to "moving" and stationary stimuli were observed under motion velocities up to 320 deg./s. The greatest response amplitudes in different experiments took place under velocities within the range of 67-320 deg./s with most of them recorded under velocities of 170 and 125 deg./s. Amplitude of the responses to lateral-medial movement with any velocity were always greater than those to opposite direction of movement with the same velocity.     

Control of a one-link arm by burst signal generators

 The focus of this paper is the study of stability and point-to-point movement of a one-link arm. The sagittal arm has two musculotendon actuators, two neural oscillators that generate burst signals as motoneuron inputs, and spindles and Golgi tendon organs for extrinsic reflex feedbacks. It is shown that coactivation leads to intrinsic position and velocity feedback, and that the tendons introduce intrinsic force and rate of force feedback. In addition, the integrating effects of the tendons are studied when the actuator is constructed from a large number of identical fibers that are excited by alpha signals whose arrival times at the fiber are randomly distributed. Each of the musculotendon actuators receives two input signals – a burst signal analogous to alpha inputs and a conventional analogue signal that represents fusimotor input to the spindles. The process of combining burst signals and conventional analogue signals is studied. Simulation results show that the movement of the system with burst signals as inputs has overshoot and speed similar to the system with analogue signals. Received: 30 May 1994/Accepted in revised form: 13 January 1995     

A hypothetical explanation of congenital nystagmus

Congenital nystagmus (CN) is a conjugate, rhythmic, eye movement disorder characterized by a wide variety of waveforms ranging from jerk to pendular types. No detailed mechanisms have been proposed to explain the generation of the CN wave-form This paper proposes a hypothetical mechanism for CN, and shows with computer simulations that a model based on this hypothesis can account for a variety of disparate waveforms. The basis of this model is a gaze-holding network, or neural integrator, that has both position and velocity feedback loops. The signals carried in these loops could arise from either afference or efference. In normal subjects, the position feedback would be positive and the velocity feedback would be negative. Both would help to increase the time constant of an imperfect neural integrator in the brain stem. We propose that in patients with CN the sign of the velocity pathway is reversed, making the neural integrator unstable. This instability could manifest as many different CN waveforms, depending on the direction and velocity of post-saccadic ocular drift and actions of nonlinearities within the position and velocity feedback loops. Thus a single underlying abnormality may be responsible for a variety of CN waveforms.     

Pedigree simulations reveal that maternity assignment is reliable in populations with conspecific brood parasitism,incomplete parental sampling and kin structure

Modern genetic parentage methods reveal that alternative reproductive strategies are common in both males and females. Under ideal conditions, genetic methods accurately connect the parents to offspring produced by extra-pair matings or conspecific brood parasitism. However, some breeding systems and sampling scenarios present significant complications for accurate parentage assignment. We used simulated genetic pedigrees to assess the reliability of parentage assignment for a series of challenging sampling regimes that reflect realistic conditions for many brood-parasitic birds: absence of genetic samples from sires, absence of samples from brood parasites and female kin-structured populations. Using 18 microsatellite markers and empirical allele frequencies from two populations of a conspecific brood parasite, the wood duck (Aix sponsa), we simulated brood parasitism and determined maternity using two widely used programs, cervus and colony . Errors in assignment were generally modest for most sampling scenarios but differed by program: cervus suffered from false assignment of parasitic offspring, whereas colony sometimes failed to assign offspring to their known mothers. Notably, colony was able to accurately infer unsampled parents. Reducing the number of markers (nine loci rather than 18) caused the assignment error to slightly worsen with colony but balloon with cervus . One potential error with important biological implications was rare in all cases—few nesting females were incorrectly excluded as the mother of their own offspring, an error that could falsely indicate brood parasitism. We consider the implications of our findings for both a retrospective assessment of previous studies and suggestions for best practices for future studies.     

Rapid determination of particle velocity from space-time images using the Radon transform

Laser-scanning methods are a means to observe streaming particles, such as the flow of red blood cells in a blood vessel. Typically, particle velocity is extracted from images formed from cyclically repeated line-scan data that is obtained along the center-line of the vessel; motion leads to streaks whose angle is a function of the velocity. Past methods made use of shearing or rotation of the images and a Singular Value Decomposition (SVD) to automatically estimate the average velocity in a temporal window of data. Here we present an alternative method that makes use of the Radon transform to calculate the velocity of streaming particles. We show that this method is over an order of magnitude faster than the SVD-based algorithm and is more robust to noise.     

Comparison between simulated annealing algorithms and rapid chain delineation in the construction of genetic maps

The efficiency of simulated annealing algorithms and rapid chain delineation in establishing the best linkage order, when constructing genetic maps, was evaluated. Linkage refers to the phenomenon by which two or more genes, or even more molecular markers, can be present in the same chromosome or linkage group. In order to evaluate the capacity of algorithms, four F(2) co-dominant populations, 50, 100, 200 and 1000 in size, were simulated. For each population, a genome with four linkage groups (100 cM) was generated. The linkage groups possessed 51, 21, 11 and 6 marks, respectively, and a corresponding distance of 2, 5, 10 and 20 cM between adjacent marks, thereby causing various degrees of saturation. For very saturated groups, with an adjacent distance between marks of 2 cM and in greater number, i.e., 51, the method based upon stochastic simulation by simulated annealing presented orders with distances equivalent to or lower than rapid chain delineation. Otherwise, the two methods were commensurate through presenting the same SARF distance.