3-d residual eddy current field characterisation: applied to diffusion weighted magnetic resonance imaging.

Clinical use of the Stejskal-Tanner diffusion weighted images is hampered by the geometric distortions that result from the large residual 3-D eddy current field induced. In this work, we aimed to predict, using linear response theory, the residual 3-D eddy current field required for geometric distortion correction based on phantom eddy current field measurements. The predicted 3-D eddy current field induced by the diffusion-weighting gradients was able to reduce the root mean square error of the residual eddy current field to ~1 Hz. The model's performance was tested on diffusion weighted images of four normal volunteers, following distortion correction, the quality of the Stejskal-Tanner diffusion-weighted images was found to have comparable quality to image registration based corrections (FSL) at low b-values. Unlike registration techniques the correction was not hindered by low SNR at high b-values, and results in improved image quality relative to FSL. Characterization of the 3-D eddy current field with linear response theory enables the prediction of the 3-D eddy current field required to correct eddy current induced geometric distortions for a wide range of clinical and high b-value protocols.

Adaptively weighted maximum likelihood estimation of discrete distributions

SummaryDiscrete data arise in various research fields, typically when the observations are count data.I propose a robust and efficient parametric procedure for estimation of discrete distributions. The estimation is done in two phases. First, a very robust, but possibly inefficient, estimate of the model parameters is computed and used to indentify outliers. Then the outliers are either removed from the sample or given low weights, and a weighted maximum likelihood estimate (WML) is computed.The weights are determined via an adaptive process such that if the data follow the model, then asymptotically no observation is downweighted.I prove that the final estimator inherits the breakdown point of the initial one, and that its influence function at the model is the same as the influence function of the maximum likelihood estimator, which strongly suggests that it is asymptotically fully efficient.The initial estimator is a minimum disparity estimator (MDE). MDEs can be shown to have full asymptotic efficiency, and some MDEs have very high breakdown points and very low bias under contamination. Several initial estimators are considered, and the performances of the WMLs based on each of them are studied.It results that in a great variety of situations the WML substantially improves the initial estimator, both in terms of finite sample mean square error and in terms of bias under contamination. Besides, the performances of the WML are rather stable under a change of the MDE even if the MDEs have very different behaviors.Two examples of application of the WML to real data are considered. In both of them, the necessity for a robust estimator is clear: the maximum likelihood estimator is badly corrupted by the presence of a few outliers.This procedure is particularly natural in the discrete distribution setting, but could be extended to the continuous case, for which a possible procedure is sketched.RésuméLes données discrètes sont présentes dans différents domaines de recherche, en particulier lorsque les observations sont des comptages.Je propose une méthode paramétrique robuste et efficace pour l'estimation de distributions discrètes. L'estimation est faite en deux phases. Tout d'abord, un estimateur très robuste des paramètres du modèle est calculé, et utilisé pour la détection des données aberrantes (outliers). Cet estimateur n'est pas nécessairement efficace. Ensuite, soit les outliers sont retirés de l'échantillon, soit des faibles poids leur sont attribués, et un estimateur du maximum de vraisemblance pondéré (WML) est calculé.Les poids sont déterminés via un processus adaptif, tel qu'asymptotiquement, si les données suivent le modèle, aucune observation n'est dépondérée.Je prouve que le point de rupture de l'estimateur final est au moins aussi élevé que celui de l'estimateur initial, et que sa fonction d'influence au modèle est la même que celle du maximum de vraisemblance, ce qui suggère que cet estimateur est pleinement efficace asymptotiquement.L'estimateur initial est un estimateur de disparité minimale (MDE). Les MDE sont asymptotiquement pleinement efficaces, et certains d'entre eux ont un point de rupture très élevé et un très faible biais sous contamination. J'étudie les performances du WML basé sur différents MDEs.Le résultat est que dans une grande variété de situations le WML améliore largement les performances de l'estimateur initial, autant en terme du carré moyen de l'erreur que du biais sous contamination. De plus, les performances du WML restent assez stables lorsqu'on change l'estimateur initial, même si les différents MDEs ont des comportements très différents.Je considère deux exemples d'application du WML à des données réelles, où la nécessité d'un estimateur robuste est manifeste : l'estimateur du maximum de vraisemblance est fortement corrompu par la présence de quelques outliers.La méthode proposée est particulièrement naturelle dans le cadre des distributions discrètes, mais pourrait être étendue au cas continu.

Influence of hip-flexion angle on hamstrings isokinetic activity in sprinters.

CONTEXT: Hamstrings strains are common and debilitating injuries in many sports. Most hamstrings exercises are performed at an inadequately low hip-flexion angle because this angle surpasses 70° at the end of the sprinting leg's swing phase, when most injuries occur. OBJECTIVE: To evaluate the influence of various hip-flexion angles on peak torques of knee flexors in isometric, concentric, and eccentric contractions and on the hamstrings-to-quadriceps ratio. DESIGN: Descriptive laboratory study. SETTING: Research laboratory. Patients and Other Participants: Ten national-level sprinters (5 men, 5 women; age = 21.2 ± 3.6 years, height = 175 ± 6 cm, mass = 63.8 ± 9.9 kg). Intervention(s): For each hip position (0°, 30°, 60°, and 90° of flexion), participants used the right leg to perform (1) 5 seconds of maximal isometric hamstrings contraction at 45° of knee flexion, (2) 5 maximal concentric knee flexion-extensions at 60° per second, (3) 5 maximal eccentric knee flexion-extensions at 60° per second, and (4) 5 maximal eccentric knee flexionextensions at 150° per second. Main Outcome Measure(s): Hamstrings and quadriceps peak torque, hamstrings-to-quadriceps ratio, lateral and medial hamstrings root mean square. RESULTS: We found no difference in quadriceps peak torque for any condition across all hip-flexion angles, whereas hamstrings peak torque was lower at 0° of hip flexion than at any other angle (P < .001) and greater at 90° of hip flexion than at 30° and 60° (P < .05), especially in eccentric conditions. As hip flexion increased, the hamstrings-to-quadriceps ratio increased. No difference in lateral or medial hamstrings root mean square was found for any condition across all hip-flexion angles (P > .05). CONCLUSIONS: Hip-flexion angle influenced hamstrings peak torque in all muscular contraction types; as hip flexion increased, hamstrings peak torque increased. Researchers should investigate further whether an eccentric resistance training program at sprint-specific hip-flexion angles (70° to 80°) could help prevent hamstrings injuries in sprinters. Moreover, hamstrings-to-quadriceps ratio assessment should be standardized at 80° of hip flexion.

Optimum power allocation for parallel Gaussian channels with arbitrary input distributions

The mutual information of independent parallel Gaussian-noise channels is maximized, under an average power constraint, by independent Gaussian inputs whose power is allocated according to the waterfilling policy. In practice, discrete signalling constellations with limited peak-to-average ratios (m-PSK, m-QAM, etc) are used in lieu of the ideal Gaussian signals. This paper gives the power allocation policy that maximizes the mutual information over parallel channels with arbitrary input distributions. Such policy admits a graphical interpretation, referred to as mercury/waterfilling, which generalizes the waterfilling solution and allows retaining some of its intuition. The relationship between mutual information of Gaussian channels and nonlinear minimum mean-square error proves key to solving the power allocation problem.

Using composite estimators to improve both domain and total area estimation

In this article we propose using small area estimators to improve the estimatesof both the small and large area parameters. When the objective is to estimateparameters at both levels accurately, optimality is achieved by a mixed sampledesign of fixed and proportional allocations. In the mixed sample design, oncea sample size has been determined, one fraction of it is distributedproportionally among the different small areas while the rest is evenlydistributed among them. We use Monte Carlo simulations to assess theperformance of the direct estimator and two composite covariant-freesmall area estimators, for different sample sizes and different sampledistributions. Performance is measured in terms of Mean Squared Errors(MSE) of both small and large area parameters. It is found that the adoptionof small area composite estimators open the possibility of 1) reducingsample size when precision is given, or 2) improving precision for a givensample size.

A decomposition formula for option prices in the Heston model and applications to option pricing approximation

By means of classical Itô's calculus we decompose option prices asthe sum of the classical Black-Scholes formula with volatility parameterequal to the root-mean-square future average volatility plus a term dueby correlation and a term due to the volatility of the volatility. Thisdecomposition allows us to develop first and second-order approximationformulas for option prices and implied volatilities in the Heston volatilityframework, as well as to study their accuracy. Numerical examples aregiven.

Estimadores compuestos en estadística regional: aplicación para la tasa de variación de la ocupación en la industria

This work is part of a project studying the performance of model basedestimators in a small area context. We have chosen a simple statisticalapplication in which we estimate the growth rate of accupation for severalregions of Spain. We compare three estimators: the direct one based onstraightforward results from the survey (which is unbiassed), and a thirdone which is based in a statistical model and that minimizes the mean squareerror.

An empirical evaluation of small area estimators

This paper investigates the comparative performance of five small areaestimators. We use Monte Carlo simulation in the context of boththeoretical and empirical populations. In addition to the direct andindirect estimators, we consider the optimal composite estimator withpopulation weights, and two composite estimators with estimatedweights: one that assumes homogeneity of within area variance andsquare bias, and another one that uses area specific estimates ofvariance and square bias. It is found that among the feasibleestimators, the best choice is the one that uses area specificestimates of variance and square bias.

Portfolio choice beyond the traditional approach

This paper surveys asset allocation methods that extend the traditional approach. An important feature of the the traditional approach is that measures the risk and return tradeoff in terms of mean and variance of final wealth. However, there are also other important features that are not always made explicit in terms of investor s wealth, information, and horizon: The investor makes a single portfolio choice based only on the mean and variance of her final financial wealth and she knows the relevant parameters in that computation. First, the paper describes traditional portfolio choice based on four basic assumptions, while the rest of the sections extend those assumptions. Each section will describe the corresponding equilibrium implications in terms of portfolio advice and asset pricing.

Combining heat stress and moderate hypoxia reduces cycling time to exhaustion without modifying neuromuscular fatigue characteristics.

PURPOSE: This study investigated the isolated and combined effects of heat [temperate (22 °C/30 % rH) vs. hot (35 °C/40 % rH)] and hypoxia [sea level (FiO2 0.21) vs. moderate altitude (FiO2 0.15)] on exercise capacity and neuromuscular fatigue characteristics. METHODS: Eleven physically active subjects cycled to exhaustion at constant workload (66 % of the power output associated with their maximal oxygen uptake in temperate conditions) in four different environmental conditions [temperate/sea level (control), hot/sea level (hot), temperate/moderate altitude (hypoxia) and hot/moderate altitude (hot + hypoxia)]. Torque and electromyography (EMG) responses following electrical stimulation of the tibial nerve (plantar-flexion; soleus) were recorded before and 5 min after exercise. RESULTS: Time to exhaustion was reduced (P < 0.05) in hot (-35 ± 15 %) or hypoxia (-36 ± 14 %) compared to control (61 ± 28 min), while hot + hypoxia (-51 ± 20 %) further compromised exercise capacity (P < 0.05). However, the effect of temperature or altitude on end-exercise core temperature (P = 0.089 and P = 0.070, respectively) and rating of perceived exertion (P > 0.05) did not reach significance. Maximal voluntary contraction torque, voluntary activation (twitch interpolation) and peak twitch torque decreased from pre- to post-exercise (-9 ± 1, -4 ± 1 and -6 ± 1 % all trials compounded, respectively; P < 0.05), with no effect of the temperature or altitude. M-wave amplitude and root mean square activity were reduced (P < 0.05) in hot compared to temperate conditions, while normalized maximal EMG activity did not change. Altitude had no effect on any measured parameters. CONCLUSION: Moderate hypoxia in combination with heat stress reduces cycling time to exhaustion without modifying neuromuscular fatigue characteristics. Impaired oxygen delivery or increased cardiovascular strain, increasing relative exercise intensity, may have also contributed to earlier exercise cessation.

Effects of prior repeated-sprints with different recovery duration on oxygen uptake kinetics during subsequent heavy-exercise

Introduction: Prior repeated-sprints (6) has become an interesting method to resolve the debate surrounding the principal factors that limits the oxygen uptake (V'O2) kinetics at the onset of exercise [i.e., muscle O2 delivery (5) or metabolic inertia (3)]. The aim of this study was to compare the effects of two repeated-sprints sets of 6x6s separated by different recovery duration between the sprints on V'O2 and muscular de-oxygenation [HHb] kinetics during a subsequent heavy-intensity exercise. Methods: 10 male subjects performed a 6-min constant-load cycling test (T50) at intensity corresponding to half of the difference between V'O2max and the ventilatory threshold. Then, they performed two repeated-sprints sets of 6x6s all-out separated by different recovery duration between the sprints (S1:30s and S2:3min) followed, after 7-min-recovery, by the T50 (S1T50 and S2T50, respectively). V'O2, [HHb] of the vastus lateralis (VL) and surface electromyography activity [i.e., root-mean-square (RMS) and the median frequency of the power density spectrum (MDF)] from VL and vastus medialis (VM) were recorded throughout T50. Models using a bi-exponential function for the overall T50 and a mono-exponential for the first 90s of T50 were used to define V'O2 and [HHb] kinetics respectively. Results: V'O2 mean value was higher in S1 (2.9±0.3l.min-1) than in S2 (1.2±0.3l.min-1); (p<0.001). The peripheral blood flow was increased after sprints as attested by a higher basal heart rate (HRbaseline) (S1T50: +22%; S2T50: +17%; p≤0.008). Time delay [HHb] was shorter for S1T50 and S2T50 than for T50 (-22% for both; p≤0.007) whereas the mean response time of V'O2 was accelerated only after S1 (S1T50: 32.3±2.5s; S2T50: 34.4±2.6s; T50: 35.7±5.4s; p=0.031). There were no significant differences in RMS between the three conditions (p>0.05). MDF of VM was higher during the first 3-min in S1T50 than in T50 (+6%; p≤0.05). Conclusion: The study show that V'O2 kinetics was speeded by prior repeated-sprints with a short (30s) but not a long (3min) inter-sprints-recovery even though the [HHb] kinetics was accelerated and the peripheral blood flow was enhanced after both sprints. S1, inducing a greater PCr depletion (1) and change in the pattern of the fibres recruitment (increase in MDF) compared with S2, may decrease metabolic inertia (2), stimulate the oxidative phosphorylation activation (4) and accelerate V'O2 kinetics at the beginning of the subsequent high-intensity exercise.

Estimation of surface roughness in a semiarid region from C-band ERS-1 synthetic aperture radar data

In this study, we investigated the feasibility of using the C-band European Remote Sensing Satellite (ERS-1) synthetic aperture radar (SAR) data to estimate surface soil roughness in a semiarid rangeland. Radar backscattering coefficients were extracted from a dry and a wet season SAR image and were compared with 47 in situ soil roughness measurements obtained in the rocky soils of the Walnut Gulch Experimental Watershed, southeastern Arizona, USA. Both the dry and the wet season SAR data showed exponential relationships with root mean square (RMS) height measurements. The dry C-band ERS-1 SAR data were strongly correlated (R² = 0.80), while the wet season SAR data have somewhat higher secondary variation (R² = 0.59). This lower correlation was probably provoked by the stronger influence of soil moisture, which may not be negligible in the wet season SAR data. We concluded that the single configuration C-band SAR data is useful to estimate surface roughness of rocky soils in a semiarid rangeland.

Effects of a 5-h hilly running on ankle plantar and dorsal flexor force and fatigability

This study aimed to examine the effects of a 5-h hilly run on ankle plantar (PF) and dorsal flexor (DF) force and fatigability. It was hypothesised that DF fatigue/fatigability would be greater than PF fatigue/fatigability. Eight male trail long distance runners (42.5 ± 5.9 years) were tested for ankle PF and DF maximal voluntary isokinetic contraction strength and fatigue resistance tests (percent decrement score), maximal voluntary and electrically evoked isometric contraction strength before and after the run. Maximal EMG root mean square (RMS(max)) and mean power frequency (MPF) values of the tibialis anterior (TA), gastrocnemius lateralis (GL) and soleus (SOL) EMG activity were calculated. The peak torque of the potentiated high- and low-frequency doublets and the ratio of paired stimulation peak torques at 10 Hz over 100 Hz (Db10:100) were analysed for PF. Maximal voluntary isometric contraction strength of PF decreased from pre- to post-run (-17.0 ± 6.2%; P < 0.05), but no significant decrease was evident for DF (-7.9 ± 6.2%). Maximal voluntary isokinetic contraction strength and fatigue resistance remained unchanged for both PF and DF. RMS(max) SOL during maximal voluntary isometric contraction and RMS(max) TA during maximal voluntary isokinetic contraction were decreased (P < 0.05) after the run. For MPF, a significant decrease for TA (P < 0.05) was found and the ratio Db10:100 decreased for PF (-6.5 ± 6.0%; P < 0.05). In conclusion, significant isometric strength loss was only detected for PF after a 5-h hilly run and was partly due to low-frequency fatigue. This study contradicted the hypothesis that neuromuscular alterations due to prolonged hilly running are predominant for DF.

Breakpoints in ventilation, cerebral and muscle oxygenation, and muscle activity during an incremental cycling exercise.

The aim of this study was to locate the breakpoints of cerebral and muscle oxygenation and muscle electrical activity during a ramp exercise in reference to the first and second ventilatory thresholds. Twenty-five cyclists completed a maximal ramp test on an electromagnetically braked cycle-ergometer with a rate of increment of 25 W/min. Expired gazes (breath-by-breath), prefrontal cortex and vastus lateralis (VL) oxygenation [Near-infrared spectroscopy (NIRS)] together with electromyographic (EMG) Root Mean Square (RMS) activity for the VL, rectus femoris (RF), and biceps femoris (BF) muscles were continuously assessed. There was a non-linear increase in both cerebral deoxyhemoglobin (at 56 ± 13% of the exercise) and oxyhemoglobin (56 ± 8% of exercise) concomitantly to the first ventilatory threshold (57 ± 6% of exercise, p > 0.86, Cohen's d < 0.1). Cerebral deoxyhemoglobin further increased (87 ± 10% of exercise) while oxyhemoglobin reached a plateau/decreased (86 ± 8% of exercise) after the second ventilatory threshold (81 ± 6% of exercise, p < 0.05, d > 0.8). We identified one threshold only for muscle parameters with a non-linear decrease in muscle oxyhemoglobin (78 ± 9% of exercise), attenuation in muscle deoxyhemoglobin (80 ± 8% of exercise), and increase in EMG activity of VL (89 ± 5% of exercise), RF (82 ± 14% of exercise), and BF (85 ± 9% of exercise). The thresholds in BF and VL EMG activity occurred after the second ventilatory threshold (p < 0.05, d > 0.6). Our results suggest that the metabolic and ventilatory events characterizing this latter cardiopulmonary threshold may affect both cerebral and muscle oxygenation levels, and in turn, muscle recruitment responses.