Model Switching and Model Averaging in Time- Varying Parameter Regression Models

This paper investigates the usefulness of switching Gaussian state space models as a tool for implementing dynamic model selecting (DMS) or averaging (DMA) in time-varying parameter regression models. DMS methods allow for model switching, where a different model can be chosen at each point in time. Thus, they allow for the explanatory variables in the time-varying parameter regression model to change over time. DMA will carry out model averaging in a time-varying manner. We compare our exact approach to DMA/DMS to a popular existing procedure which relies on the use of forgetting factor approximations. In an application, we use DMS to select different predictors in an in ation forecasting application. We also compare different ways of implementing DMA/DMS and investigate whether they lead to similar results.

Investigation of field and diffusion time dependence of the diffusion-weighted signal at ultrahigh magnetic fields.

Over the last decade, there has been a significant increase in the number of high-magnetic-field MRI magnets. However, the exact effect of a high magnetic field strength (B0 ) on diffusion-weighted MR signals is not yet fully understood. The goal of this study was to investigate the influence of different high magnetic field strengths (9.4 T and 14.1 T) and diffusion times (9, 11, 13, 15, 17 and 24 ms) on the diffusion-weighted signal in rat brain white matter. At a short diffusion time (9 ms), fractional anisotropy values were found to be lower at 14.1 T than at 9.4 T, but this difference disappeared at longer diffusion times. A simple two-pool model was used to explain these findings. The model describes the white matter as a first hindered compartment (often associated with the extra-axonal space), characterized by a faster orthogonal diffusion and a lower fractional anisotropy, and a second restricted compartment (often associated with the intra-axonal space), characterized by a slower orthogonal diffusion (i.e. orthogonal to the axon direction) and a higher fractional anisotropy. Apparent T2 relaxation time measurements of the hindered and restricted pools were performed. The shortening of the pseudo-T2 value from the restricted compartment with B0 is likely to be more pronounced than the apparent T2 changes in the hindered compartment. This study suggests that the observed differences in diffusion tensor imaging parameters between the two magnetic field strengths at short diffusion time may be related to differences in the apparent T2 values between the pools. Copyright © 2013 John Wiley & Sons, Ltd.

Optimització i implementació en Simulink de l'adquisició d'un senyal GPS amb blocs FFT

Aquest projecte es tracta de la optimització i la implementació de l’etapa d’adquisició d’un receptor GPS. També inclou una revisió breu del sistema GPS i els seus principis de funcionament. El procés d’adquisició s’ha estudiat amb detall i programat en els entorns de treball Matlab i Simulink. El fet d’implementar aquesta etapa en dos entorns diferents ha estat molt útil tant de cara a l’aprenentatge com també per la comprovació dels resultats obtinguts. El principal objectiu del treball és el disseny d’un model Simulink que es capaç d’adquirir una senyal capturada amb hardware real. En realitat, s’han fet dues implementacions: una que utilitza blocs propis de Simulink i l’altra que utilitza blocs de la llibreria Xilinx. D’aquesta manera, posteriorment, es facilitaria la transició del model a la FPGA utilitzant l’entorn ISE de Xilinx. La implementació de l’etapa d’adquisició es basa en el mètode de cerca de fase de codi en paral·lel, el qual empra la operació correlació creuada mitjançant la transformada ràpida de Fourier (FFT). Per aquest procés es necessari realitzar dues transformades (per a la senyal entrant i el codi de referència) i una antitransformada de Fourier (per al resultat de la correlació). Per tal d’optimitzar el disseny s’utilitza un bloc FFT, ja que tres blocs consumeixen gran part dels recursos d’una FPGA. En lloc de replicar el bloc FFT, en el model el bloc és compartit en el temps gràcies a l’ús de buffers i commutadors, com a resultat la quantitat de recursos requerits per una implementació en una FPGA es podria reduir considerablement.

Gaussian estimates for the density of the non-linear stochastic heat equation in any space dimension

In this paper, we establish lower and upper Gaussian bounds for the probability density of the mild solution to the stochastic heat equation with multiplicative noise and in any space dimension. The driving perturbation is a Gaussian noise which is white in time with some spatially homogeneous covariance. These estimates are obtained using tools of the Malliavin calculus. The most challenging part is the lower bound, which is obtained by adapting a general method developed by Kohatsu-Higa to the underlying spatially homogeneous Gaussian setting. Both lower and upper estimates have the same form: a Gaussian density with a variance which is equal to that of the mild solution of the corresponding linear equation with additive noise.

Phylogenetic Code in the Cloud - Can it Meet the Expectations?

Cloud computing has recently become very popular, and several bioinformatics applications exist already in that domain. The aim of this article is to analyse a current cloud system with respect to usability, benchmark its performance and compare its user friendliness with a conventional cluster job submission system. Given the current hype on the theme, user expectations are rather high, but current results show that neither the price/performance ratio nor the usage model is very satisfactory for large-scale embarrassingly parallel applications. However, for small to medium scale applications that require CPU time at certain peak times the cloud is a suitable alternative.

Spin-labeling coronary MR angiography with steady-state free precession and radial k-space sampling: initial results in healthy volunteers.

The purpose of this study was to prospectively compare free-breathing navigator-gated cardiac-triggered three-dimensional steady-state free precession (SSFP) spin-labeling coronary magnetic resonance (MR) angiography performed by using Cartesian k-space sampling with that performed by using radial k-space sampling. A new dedicated placement of the two-dimensional selective labeling pulse and an individually adjusted labeling delay time approved by the institutional review board were used. In 14 volunteers (eight men, six women; mean age, 28.8 years) who gave informed consent, signal-to-noise ratio (SNR), contrast-to-noise ratio (CNR), vessel sharpness, vessel length, and subjective image quality were investigated. Differences between groups were analyzed with nonparametric tests (Wilcoxon, Pearson chi2). Radial imaging, as compared with Cartesian imaging, resulted in a significant reduction in the severity of motion artifacts, as well as an increase in SNR (26.9 vs 12.0, P < .05) in the coronary arteries and CNR (23.1 vs 8.8, P < .05) between the coronary arteries and the myocardium. A tendency toward improved vessel sharpness and vessel length was also found with radial imaging. Radial SSFP imaging is a promising technique for spin-labeling coronary MR angiography.

GLUT2 surface expression and intracellular transport via the constitutive pathway in pancreatic beta cells and insulinoma: evidence for a block in trans-Golgi network exit by brefeldin A.

The biosynthesis, intracellular transport, and surface expression of the beta cell glucose transporter GLUT2 was investigated in isolated islets and insulinoma cells. Using a trypsin sensitivity assay to measure cell surface expression, we determined that: (a) greater than 95% of GLUT2 was expressed on the plasma membrane; (b) GLUT2 did not recycle in intracellular vesicles; and (c) after trypsin treatment, reexpression of the intact transporter occurred with a t1/2 of approximately 7 h. Kinetics of intracellular transport of GLUT2 was investigated in pulse-labeling experiments combined with glycosidase treatment and the trypsin sensitivity assay. We determined that transport from the endoplasmic reticulum to the trans-Golgi network (TGN) occurred with a t1/2 of 15 min and that transport from the TGN to the plasma membrane required a similar half-time. When added at the start of a pulse-labeling experiment, brefeldin A prevented exit of GLUT2 from the endoplasmic reticulum. When the transporter was first accumulated in the TGN during a 15-min period of chase, but not following a low temperature (22 degrees C) incubation, addition of brefeldin A (BFA) prevented subsequent surface expression of the transporter. This indicated that brefeldin A prevented GLUT2 exit from the TGN by acting at a site proximal to the 22 degrees C block. Together, these data demonstrate that GLUT2 surface expression in beta cells is via the constitutive pathway, that transport can be blocked by BFA at two distinct steps and that once on the surface, GLUT2 does not recycle in intracellular vesicles.

Isolated congenital atrioventricular block diagnosed in utero: natural history and outcome.

BACKGROUND: Isolated congenital atrioventricular block (CAVB) diagnosed in utero is associated with a high morbidity and mortality. Prognosis is especially poor when heart rate drops below 55 beats per minute (bpm) and when fetal hydrops develops. We describe the natural history and outcome of 24 infants with isolated CAVB diagnosed in utero, review the literature, and assess the risk factors that could predict outcome. METHODS: This was a retrospective multicenter study of 24 patients with isolated CAVB diagnosed in utero. RESULTS: CAVB was detected at a mean gestational age (GA) of 24.7 +/- 5.1 weeks. Ten fetuses initially presented with complete heart block. Low heart rate or incomplete heart block was the first documentation of bradyarrhythmia in the other 14 fetuses. In 11 of them, CAVB developed during pregnancy after a median time of 3 (range 1-16) weeks. Fetal hydrops developed in 10 of 24 (42%) fetuses at a mean GA of 27.6 +/- 5.1 weeks. Hydropic fetuses showed lower heart rates during pregnancy (47 +/- 10 bpm) than non-hydropic fetuses (57 +/- 10 bpm). There were three intrauterine deaths; all were hydropic and female. Nine viable females and 12 males were born at a mean GA of 37.1 +/- 6.1 weeks with an average birth weight of 3097 +/- 852 g. Fifteen CAVB patients required pacemaker (PM) intervention, 10 of them immediately after birth. Dilated cardiomyopathy (DCM) developed in three infants of whom two died of congestive heart failure, shortly after the diagnosis was made; one is still alive. Mortality before or after birth was 21%, and was associated with heart rates below 50 bpm and development of fetal hydrops. Poor outcome, defined as death, PM implantation, or development of DCM, occurred in 83% of cases and was associated with heart rates below 60 bpm during pregnancy. CONCLUSIONS: Isolated CAVB diagnosed in utero is associated with high morbidity and mortality. Patients who develop fetal hydrops show lower heart rates during pregnancy than patients who do not. A fetal heart rate below 50 bpm and development of fetal hydrops is associated with increased mortality. Rates below 60 bpm are associated with PM requirement and/or DCM.

Documenting a time-bound, circular view of hierarchies: a microanalysis of parent-infant dyadic interaction.

This article presents a new theory that separates the levels of communication and relates them circularly, namely, by separating time from space/meaning variables. Documenting this proposition requires sequential microdescriptions--a far-out project in the field of family therapy. In an extensive study of clinical and nonclinical families, starting with available microanalytic data on nonverbal parent-infant dialogue, distinct time organizations have been found to modify the degree of circularity between the levels of interaction according to the observed types of engagement, that is, consensual, conflictual, and paradoxical. The double description of the dyad as a totality versus the dyad as a framing/developing organization imparts crucial information on how development proceeds in dyadic, co-evolutive systems, and presumably in larger ones too. In this perspective, a model is elaborated and then applied to a case description in our therapeutic consultation.

Single-injection or continuous femoral nerve block for total knee arthroplasty?

BACKGROUND: The ideal local anesthetic regime for femoral nerve block that balances analgesia with mobility after total knee arthroplasty (TKA) remains undefined. QUESTIONS/PURPOSES: We compared two volumes and concentrations of a fixed dose of ropivacaine for continuous femoral nerve block after TKA to a single injection femoral nerve block with ropivacaine to determine (1) time to discharge readiness; (2) early pain scores and analgesic consumption; and (3) functional outcomes, including range of motion and WOMAC scores at the time of recovery. METHODS: Ninety-nine patients were allocated to one of three continuous femoral nerve block groups for this randomized, placebo-controlled, double-blind trial: a high concentration group (ropivacaine 0.2% infusion), a low concentration group (ropivacaine 0.1% infusion), or a placebo infusion group (saline 0.9% infusion). Infusions were discontinued on postoperative Day (POD) 2. The primary outcome was time to discharge readiness. Secondary outcomes included opioid consumption, pain, and functional outcomes. Ninety-three patients completed the study protocol; the study was halted early because of unanticipated changes to pain protocols at the host institution, by which time only 61% of the required number of patients had been enrolled. RESULTS: With the numbers available, the mean time to discharge readiness was not different between groups (high concentration group, 62 hours [95% confidence interval [CI], 51-72 hours]; low concentration group, 73 hours [95% CI, 63-83 hours]; placebo infusion group 65 hours [95% CI, 56-75 hours]; p = 0.27). Patients in the low concentration group consumed significantly less morphine during the period of infusion (POD 1, high concentration group, 56 mg [95% CI, 42-70 mg]; low concentration group, 35 mg [95% CI, 27-43 mg]; placebo infusion group, 48 mg [95% CI, 38-59 mg], p = 0.02; POD 2, high concentration group, 50 mg [95% CI, 41-60 mg]; low concentration group, 33 mg [95% CI, 24-42 mg]; placebo infusion group, 39 mg [95% CI, 30-48 mg], p = 0.04); however, there were no important differences in pain scores or opioid-related side effects with the numbers available. Likewise, there were no important differences in functional outcomes between groups. CONCLUSIONS: Based on this study, which was terminated prematurely before the desired sample size could be achieved, we were unable to demonstrate that varying the concentration and volume of a fixed-dose ropivacaine infusion for continuous femoral nerve block influences time to discharge readiness when compared with a conventional single-injection femoral nerve block after TKA. A low concentration of ropivacaine infusion can reduce postoperative opioid consumption but without any important differences in pain scores, side effects, or functional outcomes. These pilot data may be used to inform the statistical power of future randomized trials. LEVEL OF EVIDENCE: Level II, therapeutic study. See Guidelines for Authors for a complete description of levels of evidence.

Real-Time implementation of a blind authentication method using self-synchronous speech watermarking

A blind speech watermarking scheme that meets hard real-time deadlines is presented and implemented. In addition, one of the key issues in these block-oriented watermarking techniques is to preserve the synchronization. Namely, to recover the exact position of each block in the mark extract process. In fact, the presented scheme can be split up into two distinguished parts, the synchronization and the information mark methods. The former is embedded into the time domain and it is fast enough to be run meeting real-time requirements. The latter contains the authentication information and it is embedded into the wavelet domain. The synchronization and information mark techniques are both tunable in order to allow a con gurable method. Thus, capacity, transparency and robustness can be con gured depending on the needs. It makes the scheme useful for professional applications, such telephony authentication or even sending information throw radio applications.

Discovering similarities in the time use patterns of the Spanish Autonomous Communities by Fuzzy techniques

In 2000 the European Statistical Office published the guidelines for developing theHarmonized European Time Use Surveys system. Under such a unified framework,the first Time Use Survey of national scope was conducted in Spain during 2002–03. The aim of these surveys is to understand human behavior and the lifestyle ofpeople. Time allocation data are of compositional nature in origin, that is, they aresubject to non-negativity and constant-sum constraints. Thus, standard multivariatetechniques cannot be directly applied to analyze them. The goal of this work is toidentify homogeneous Spanish Autonomous Communities with regard to the typicalactivity pattern of their respective populations. To this end, fuzzy clustering approachis followed. Rather than the hard partitioning of classical clustering, where objects areallocated to only a single group, fuzzy method identify overlapping groups of objectsby allowing them to belong to more than one group. Concretely, the probabilistic fuzzyc-means algorithm is conveniently adapted to deal with the Spanish Time Use Surveymicrodata. As a result, a map distinguishing Autonomous Communities with similaractivity pattern is drawn.Key words: Time use data, Fuzzy clustering; FCM; simplex space; Aitchison distance

Cell size distribution in a random tessellation of space governed by the Kolmogorov-Johnson-Mehl-Avrami model: Grain size distribution in crystallization

The space subdivision in cells resulting from a process of random nucleation and growth is a subject of interest in many scientific fields. In this paper, we deduce the expected value and variance of these distributions while assuming that the space subdivision process is in accordance with the premises of the Kolmogorov-Johnson-Mehl-Avrami model. We have not imposed restrictions on the time dependency of nucleation and growth rates. We have also developed an approximate analytical cell size probability density function. Finally, we have applied our approach to the distributions resulting from solid phase crystallization under isochronal heating conditions

A new efficient approach to the direct restricted active space self-consistent field method

An implicitly parallel method for integral-block driven restricted active space self-consistent field (RASSCF) algorithms is presented. The approach is based on a model space representation of the RAS active orbitals with an efficient expansion of the model subspaces. The applicability of the method is demonstrated with a RASSCF investigation of the first two excited states of indole

What can spike train distances tell us about the neural code?

Time scale parametric spike train distances like the Victor and the van Rossum distancesare often applied to study the neural code based on neural stimuli discrimination.Different neural coding hypotheses, such as rate or coincidence coding,can be assessed by combining a time scale parametric spike train distance with aclassifier in order to obtain the optimal discrimination performance. The time scalefor which the responses to different stimuli are distinguished best is assumed to bethe discriminative precision of the neural code. The relevance of temporal codingis evaluated by comparing the optimal discrimination performance with the oneachieved when assuming a rate code.We here characterize the measures quantifying the discrimination performance,the discriminative precision, and the relevance of temporal coding. Furthermore,we evaluate the information these quantities provide about the neural code. Weshow that the discriminative precision is too unspecific to be interpreted in termsof the time scales relevant for encoding. Accordingly, the time scale parametricnature of the distances is mainly an advantage because it allows maximizing thediscrimination performance across a whole set of measures with different sensitivitiesdetermined by the time scale parameter, but not due to the possibility toexamine the temporal properties of the neural code.