Method for exploratory cluster analysis and visualisation of single-trial ERP ensembles

Background: The validity of ensemble averaging on event-related potential (ERP) data has been questioned, due to its assumption that the ERP is identical across trials. Thus, there is a need for preliminary testing for cluster structure in the data. New method: We propose a complete pipeline for the cluster analysis of ERP data. To increase the signalto-noise (SNR) ratio of the raw single-trials, we used a denoising method based on Empirical Mode Decomposition (EMD). Next, we used a bootstrap-based method to determine the number of clusters, through a measure called the Stability Index (SI). We then used a clustering algorithm based on a Genetic Algorithm (GA)to define initial cluster centroids for subsequent k-means clustering. Finally, we visualised the clustering results through a scheme based on Principal Component Analysis (PCA). Results: After validating the pipeline on simulated data, we tested it on data from two experiments – a P300 speller paradigm on a single subject and a language processing study on 25 subjects. Results revealed evidence for the existence of 6 clusters in one experimental condition from the language processing study. Further, a two-way chi-square test revealed an influence of subject on cluster membership.

Extracting planetary waves from geomagnetic time series using Empirical Mode Decomposition

Empirical Mode Decomposition is presented as an alternative to traditional analysis methods to decompose geomagnetic time series into spectral components. Important comments on the algorithm and its variations will be given. Using this technique, planetary wave modes of 5-, 10-, and 16-day mean periods can be extracted from magnetic field components of three different stations in Germany. In a second step, the amplitude modulation functions of these wave modes can be shown to contain significant contribution from solar cycle variation through correlation with smoothed sunspot numbers. Additionally, the data indicate connections with geomagnetic jerk occurrences, supported by a second set of data providing reconstructed near-Earth magnetic field for 150 years. Usually attributed to internal dynamo processes within the Earth's outer core, the question of who is impacting whom will be briefly discussed here.

Human versus animal: contrasting decomposition dynamics of mammalian analogues in experimental taphonomy

Taphonomic studies regularly employ animal analogues for human decomposition due to ethical restrictions relating to the use of human tissue. However, the validity of using animal analogues in soil decomposition studies is still questioned. This study compared the decomposition of skeletal muscle tissues (SMTs) from human (Homo sapiens), pork (Sus scrofa), beef (Bos taurus), and lamb (Ovis aries) interred in soil microcosms. Fixed interval samples were collected from the SMT for microbial activity and mass tissue loss determination; samples were also taken from the underlying soil for pH, electrical conductivity, and nutrient (potassium, phosphate, ammonium, and nitrate) analysis. The overall patterns of nutrient fluxes and chemical changes in nonhuman SMT and the underlying soil followed that of human SMT. Ovine tissue was the most similar to human tissue in many of the measured parameters. Although no single analogue was a precise predictor of human decomposition in soil, all models offered close approximations in decomposition dynamics.

Freezing skeletal muscle tissue does not affect its decomposition in soil: Evidence from temporal changes in tissue mass, microbial activity and soil chemistry based on excised samples

The study of decaying organisms and death assemblages is referred to as forensic taphonomy, or more simply the study of graves. This field is dominated by the fields of entomology, anthropology and archaeology. Forensic taphonomy also includes the study of the ecology and chemistry of the burial environment. Studies in forensic taphonomy often require the use of analogues for human cadavers or their component parts. These might include animal cadavers or skeletal muscle tissue. However, sufficient supplies of cadavers or analogues may require periodic freezing of test material prior to experimental inhumation in the soil. This study was carried out to ascertain the effect of freezing on skeletal muscle tissue prior to inhumation and decomposition in a soil environment under controlled laboratory conditions. Changes in soil chemistry were also measured. In order to test the impact of freezing, skeletal muscle tissue (Sus scrofa) was frozen (−20 °C) or refrigerated (4 °C). Portions of skeletal muscle tissue (∼1.5 g) were interred in microcosms (72 mm diameter × 120 mm height) containing sieved (2 mm) soil (sand) adjusted to 50% water holding capacity. The experiment had three treatments: control with no skeletal muscle tissue, microcosms containing frozen skeletal muscle tissue and those containing refrigerated tissue. The microcosms were destructively harvested at sequential periods of 2, 4, 6, 8, 12, 16, 23, 30 and 37 days after interment of skeletal muscle tissue. These harvests were replicated 6 times for each treatment. Microbial activity (carbon dioxide respiration) was monitored throughout the experiment. At harvest the skeletal muscle tissue was removed and the detritosphere soil was sampled for chemical analysis. Freezing was found to have no significant impact on decomposition or soil chemistry compared to unfrozen samples in the current study using skeletal muscle tissue. However, the interment of skeletal muscle tissue had a significant impact on the microbial activity (carbon dioxide respiration) and chemistry of the surrounding soil including: pH, electroconductivity, ammonium, nitrate, phosphate and potassium. This is the first laboratory controlled study to measure changes in inorganic chemistry in soil associated with the decomposition of skeletal muscle tissue in combination with microbial activity.

Heterogeneous tensor decomposition for clustering via manifold optimization

Tensor clustering is an important tool that exploits intrinsically rich structures in real-world multiarray or Tensor datasets. Often in dealing with those datasets, standard practice is to use subspace clustering that is based on vectorizing multiarray data. However, vectorization of tensorial data does not exploit complete structure information. In this paper, we propose a subspace clustering algorithm without adopting any vectorization process. Our approach is based on a novel heterogeneous Tucker decomposition model taking into account cluster membership information. We propose a new clustering algorithm that alternates between different modes of the proposed heterogeneous tensor model. All but the last mode have closed-form updates. Updating the last mode reduces to optimizing over the multinomial manifold for which we investigate second order Riemannian geometry and propose a trust-region algorithm. Numerical experiments show that our proposed algorithm compete effectively with state-of-the-art clustering algorithms that are based on tensor factorization.

Black-box decomposition approach for computational hemodynamics: One-dimensional models

Increasing efforts exist in integrating different levels of detail in models of the cardiovascular system. For instance, one-dimensional representations are employed to model the systemic circulation. In this context, effective and black-box-type decomposition strategies for one-dimensional networks are needed, so as to: (i) employ domain decomposition strategies for large systemic models (1D-1D coupling) and (ii) provide the conceptual basis for dimensionally-heterogeneous representations (1D-3D coupling, among various possibilities). The strategy proposed in this article works for both of these two scenarios, though the several applications shown to illustrate its performance focus on the 1D-1D coupling case. A one-dimensional network is decomposed in such a way that each coupling point connects two (and not more) of the sub-networks. At each of the M connection points two unknowns are defined: the flow rate and pressure. These 2M unknowns are determined by 2M equations, since each sub-network provides one (non-linear) equation per coupling point. It is shown how to build the 2M x 2M non-linear system with arbitrary and independent choice of boundary conditions for each of the sub-networks. The idea is then to solve this non-linear system until convergence, which guarantees strong coupling of the complete network. In other words, if the non-linear solver converges at each time step, the solution coincides with what would be obtained by monolithically modeling the whole network. The decomposition thus imposes no stability restriction on the choice of the time step size. Effective iterative strategies for the non-linear system that preserve the black-box character of the decomposition are then explored. Several variants of matrix-free Broyden`s and Newton-GMRES algorithms are assessed as numerical solvers by comparing their performance on sub-critical wave propagation problems which range from academic test cases to realistic cardiovascular applications. A specific variant of Broyden`s algorithm is identified and recommended on the basis of its computer cost and reliability. (C) 2010 Elsevier B.V. All rights reserved.

Fiedler trees for multiscale surface analysis

In this work we introduce a new hierarchical surface decomposition method for multiscale analysis of surface meshes. In contrast to other multiresolution methods, our approach relies on spectral properties of the surface to build a binary hierarchical decomposition. Namely, we utilize the first nontrivial eigenfunction of the Laplace-Beltrami operator to recursively decompose the surface. For this reason we coin our surface decomposition the Fiedler tree. Using the Fiedler tree ensures a number of attractive properties, including: mesh-independent decomposition, well-formed and nearly equi-areal surface patches, and noise robustness. We show how the evenly distributed patches can be exploited for generating multiresolution high quality uniform meshes. Additionally, our decomposition permits a natural means for carrying out wavelet methods, resulting in an intuitive method for producing feature-sensitive meshes at multiple scales. Published by Elsevier Ltd.

The role of citrate precursors on the morphology of lanthanide oxides obtained by thermal decomposition

Two series of lanthanide oxides with different morphologies were synthesized through calcinations of two types of citrate polymeric precursors. These oxides were characterized by XRD patterns, SEM electronic microscopy, and N(2) adsorption isotherms. SEM microscopy analysis showed that the calcination of crystalline fibrous precursors [Ln(2)(LH)(3)center dot 2H(2)O] (L = citrate) originated fibrous shaped particles. On the other hand, the calcination of irregular shaped particles of precursors [LnL center dot xH(2)O] originated irregular shaped particles of oxide, pointing out a morphological template effect of precursors on the formation of the respective oxides.

Many-body energy decomposition of hydrogen-bonded glycine clusters in gas-phase

A detailed analysis of the many-body contribution to the interaction energies of the gas-phase hydrogen-bonded glycine clusters, (Gly)(N), N = 1-4 is presented. The energetics of the hydrogen-bonded dimer, trimer and tetramer complexes have been analyzed using density-functional theory. The magnitude of the two-through four-body energy terms have been calculated and compared. The relaxation energy and the two-body energy terms are the principal contributors to the total binding energy. Four-body contribution is negligible. However, the three-body contribution is found to be sizable and the formation of the cyclic glycine trimer presents geometric strains that make it less favorable. (C) 2010 Elsevier B.V. All rights reserved.

Leverage analysis for linear mixed models

We consider a generalized leverage matrix useful for the identification of influential units and observations in linear mixed models and show how a decomposition of this matrix may be employed to identify high leverage points for both the marginal fitted values and the random effect component of the conditional fitted values. We illustrate the different uses of the two components of the decomposition with a simulated example as well as with a real data set.

Towards a maximal extension of Pelczynski`s decomposition method in Banach spaces

l Suppose that X, Y. A and B are Banach spaces such that X is isomorphic to Y E) A and Y is isomorphic to X circle plus B. Are X and Y necessarily isomorphic? In this generality. the answer is no, as proved by W.T. Cowers in 1996. In the present paper, we provide a very simple necessary and sufficient condition on the 10-tuples (k, l, m, n. p, q, r, s, u, v) in N with p+q+u >= 3, r+s+v >= 3, uv >= 1, (p,q)$(0,0), (r,s)not equal(0,0) and u=1 or v=1 or (p. q) = (1, 0) or (r, s) = (0, 1), which guarantees that X is isomorphic to Y whenever these Banach spaces satisfy X(u) similar to X(p)circle plus Y(q), Y(u) similar to X(r)circle plus Y(s), and A(k) circle plus B(l) similar to A(m) circle plus B(n). Namely, delta = +/- 1 or lozenge not equal 0, gcd(lozenge, delta (p + q - u)) divides p + q - u and gcd(lozenge, delta(r + s - v)) divides r + s - v, where 3 = k - I - in + n is the characteristic number of the 4-tuple (k, l, m, n) and lozenge = (p - u)(s - v) - rq is the discriminant of the 6-tuple (p, q, r, s, U, v). We conjecture that this result is in some sense a maximal extension of the classical Pelczynski`s decomposition method in Banach spaces: the case (1, 0. 1, 0, 2. 0, 0, 2. 1. 1). (C) 2009 Elsevier Inc. All rights reserved.

Flow injection analysis of ethyl xanthate by gas diffusion and UV detection as CS(2) for process monitoring of sulfide ore flotation

A sensitive and robust analytical method for spectrophotometric determination of ethyl xanthate, CH(3)CH(2)OCS(2)(-) at trace concentrations in pulp solutions from froth flotation process is proposed. The analytical method is based on the decomposition of ethyl xanthate. EtX(-), with 2.0 mol L(-1) HCl generating ethanol and carbon disulfide. CS(2). A gas diffusion cell assures that only the volatile compounds diffuse through a PTFE membrane towards an acceptor stream of deionized water, thus avoiding the interferences of non-volatile compounds and suspended particles. The CS(2) is selectively detected by UV absorbance at 206 nm (epsilon = 65,000 L mol(-1) cm(-1)). The measured absorbance is directly proportional to EtX(-) concentration present in the sample solutions. The Beer`s law is obeyed in a 1 x 10(-6) to 2 x 10(-4) mol L(-1) concentration range of ethyl xanthate in the pulp with an excellent correlation coefficient (r = 0.999) and a detection limit of 3.1 x 10(-7) mol L(-1), corresponding to 38 mu g L. At flow rates of 200 mu L min(-1) of the donor stream and 100 mu L min(-1) of the acceptor channel a sampling rate of 15 injections per hour could be achieved with RSD < 2.3% (n = 10, 300 mu L injections of 1 x 10(-5) mol L(-1) EtX(-)). Two practical applications demonstrate the versatility of the FIA method: (i) evaluation the free EtX(-) concentration during a laboratory study of the EtX(-) adsorption capacity on pulverized sulfide ore (pyrite) and (ii) monitoring of EtX(-) at different stages (from starting load to washing effluents) of a flotation pilot plant processing a Cu-Zn sulfide ore. (C) 2010 Elsevier By. All rights reserved.

Thermal analysis and compatibility studies of prednicarbate with excipients used in semi solid pharmaceutical form

Differential Scanning Calorimetry (DSC), thermogravimetry/derivative thermogravimetry (TG/DTG) and infrared spectroscopy (IR) techniques were used to investigate the compatibility between prednicarbate and several excipients commonly used in semi solid pharmaceutical form. The thermoanalytical studies of 1:1 (m/m) drug/excipient physical mixtures showed that the beginning of the first thermal decomposition stage of the prednicarbate (T (onset) value) was decreased in the presence of stearyl alcohol and glyceryl stearate compared to the drug alone. For the binary mixture of drug/sodium pirrolidone carboxilate the first thermal decomposition stage was not changed, however the DTG peak temperature (T (peak DTG)) decreased. The comparison of the IR spectra of the drug, the physical mixtures and of the thermally treated samples confirmed the thermal decomposition of prednicarbate. By the comparison of the thermal profiles of 1:1 prednicarbate:excipients mixtures (methylparaben, propylparaben, carbomer 940, acrylate crosspolymer, lactic acid, light liquid paraffin, isopropyl palmitate, myristyl lactate and cetyl alcohol) no interaction was observed.

THERMAL BEHAVIOR STUDY AND DECOMPOSITION KINETICS OF SALBUTAMOL UNDER ISOTHERMAL AND NON-ISOTHERMAL CONDITIONS

The thermal decomposition of salbutamol (beta(2) - selective adrenoreceptor) was studied using differential scanning calorimetry (DSC) and thermogravimetry/derivative thermogravimetry (TG/DTG). It was observed that the commercial sample showed a different thermal profile than the standard sample caused by the presence of excipients. These compounds increase the thermal stability of the drug. Moreover, higher activation energy was calculated for the pharmaceutical sample, which was estimated by isothermal and non-isothermal methods for the first stage of the thermal decomposition process. For isothermal experiments the average values were E(act) = 130 kJ mol(-1) (for standard sample) and E(act) = 252 kJ mol(-1) (for pharmaceutical sample) in a dynamic nitrogen atmosphere (50 mL min(-1)). For non-isothermal method, activation energy was obtained from the plot of log heating rates vs. 1/T in dynamic air atmosphere (50 mL min(-1)). The calculated values were E(act) = 134 kJ mol(-1) (for standard sample) and E(act) (=) 139 kJ mol(-1) (for pharmaceutical sample).