Characterization of microcirculation in multiple sclerosis lesions by dynamic texture parameter analysis (DTPA)

OBJECTIVE Texture analysis is an alternative method to quantitatively assess MR-images. In this study, we introduce dynamic texture parameter analysis (DTPA), a novel technique to investigate the temporal evolution of texture parameters using dynamic susceptibility contrast enhanced (DSCE) imaging. Here, we aim to introduce the method and its application on enhancing lesions (EL), non-enhancing lesions (NEL) and normal appearing white matter (NAWM) in multiple sclerosis (MS). METHODS We investigated 18 patients with MS and clinical isolated syndrome (CIS), according to the 2010 McDonald's criteria using DSCE imaging at different field strengths (1.5 and 3 Tesla). Tissues of interest (TOIs) were defined within 27 EL, 29 NEL and 37 NAWM areas after normalization and eight histogram-based texture parameter maps (TPMs) were computed. TPMs quantify the heterogeneity of the TOI. For every TOI, the average, variance, skewness, kurtosis and variance-of-the-variance statistical parameters were calculated. These TOI parameters were further analyzed using one-way ANOVA followed by multiple Wilcoxon sum rank testing corrected for multiple comparisons. RESULTS Tissue- and time-dependent differences were observed in the dynamics of computed texture parameters. Sixteen parameters discriminated between EL, NEL and NAWM (pAVG = 0.0005). Significant differences in the DTPA texture maps were found during inflow (52 parameters), outflow (40 parameters) and reperfusion (62 parameters). The strongest discriminators among the TPMs were observed in the variance-related parameters, while skewness and kurtosis TPMs were in general less sensitive to detect differences between the tissues. CONCLUSION DTPA of DSCE image time series revealed characteristic time responses for ELs, NELs and NAWM. This may be further used for a refined quantitative grading of MS lesions during their evolution from acute to chronic state. DTPA discriminates lesions beyond features of enhancement or T2-hypersignal, on a numeric scale allowing for a more subtle grading of MS-lesions.

Effect of nonnormal random components on level and power in group-randomized trials

The use of group-randomized trials is particularly widespread in the evaluation of health care, educational, and screening strategies. Group-randomized trials represent a subset of a larger class of designs often labeled nested, hierarchical, or multilevel and are characterized by the randomization of intact social units or groups, rather than individuals. The application of random effects models to group-randomized trials requires the specification of fixed and random components of the model. The underlying assumption is usually that these random components are normally distributed. This research is intended to determine if the Type I error rate and power are affected when the assumption of normality for the random component representing the group effect is violated. ^ In this study, simulated data are used to examine the Type I error rate, power, bias and mean squared error of the estimates of the fixed effect and the observed intraclass correlation coefficient (ICC) when the random component representing the group effect possess distributions with non-normal characteristics, such as heavy tails or severe skewness. The simulated data are generated with various characteristics (e.g. number of schools per condition, number of students per school, and several within school ICCs) observed in most small, school-based, group-randomized trials. The analysis is carried out using SAS PROC MIXED, Version 6.12, with random effects specified in a random statement and restricted maximum likelihood (REML) estimation specified. The results from the non-normally distributed data are compared to the results obtained from the analysis of data with similar design characteristics but normally distributed random effects. ^ The results suggest that the violation of the normality assumption for the group component by a skewed or heavy-tailed distribution does not appear to influence the estimation of the fixed effect, Type I error, and power. Negative biases were detected when estimating the sample ICC and dramatically increased in magnitude as the true ICC increased. These biases were not as pronounced when the true ICC was within the range observed in most group-randomized trials (i.e. 0.00 to 0.05). The normally distributed group effect also resulted in bias ICC estimates when the true ICC was greater than 0.05. However, this may be a result of higher correlation within the data. ^

Can we infer plant facilitation from remote sensing? a test across global drylands

Facilitation is a major force shaping the structure and diversity of plant communities in terrestrial ecosystems. Detecting positive plant–plant interactions relies on the combination of field experimentation and the demonstration of spatial association between neighboring plants. This has often restricted the study of facilitation to particular sites, limiting the development of systematic assessments of facilitation over regional and global scales. Here we explore whether the frequency of plant spatial associations detected from high-resolution remotely sensed images can be used to infer plant facilitation at the community level in drylands around the globe. We correlated the information from remotely sensed images freely available through Google Earth with detailed field assessments, and used a simple individual-based model to generate patch-size distributions using different assumptions about the type and strength of plant–plant interactions. Most of the patterns found from the remotely sensed images were more right skewed than the patterns from the null model simulating a random distribution. This suggests that the plants in the studied drylands show stronger spatial clustering than expected by chance. We found that positive plant co-occurrence, as measured in the field, was significantly related to the skewness of vegetation patch-size distribution measured using Google Earth images. Our findings suggest that the relative frequency of facilitation may be inferred from spatial pattern signals measured from remotely sensed images, since facilitation often determines positive co-occurrence among neighboring plants. They pave the road for a systematic global assessment of the role of facilitation in terrestrial ecosystems. Read More: http://www.esajournals.org/doi/10.1890/14-2358.1

Modeling absolute differences in life expectancy with a censored skew-normal regression approach.

Parameter estimates from commonly used multivariable parametric survival regression models do not directly quantify differences in years of life expectancy. Gaussian linear regression models give results in terms of absolute mean differences, but are not appropriate in modeling life expectancy, because in many situations time to death has a negative skewed distribution. A regression approach using a skew-normal distribution would be an alternative to parametric survival models in the modeling of life expectancy, because parameter estimates can be interpreted in terms of survival time differences while allowing for skewness of the distribution. In this paper we show how to use the skew-normal regression so that censored and left-truncated observations are accounted for. With this we model differences in life expectancy using data from the Swiss National Cohort Study and from official life expectancy estimates and compare the results with those derived from commonly used survival regression models. We conclude that a censored skew-normal survival regression approach for left-truncated observations can be used to model differences in life expectancy across covariates of interest.

Automatic quality control in clinical (1) H MRSI of brain cancer.

MRSI grids frequently show spectra with poor quality, mainly because of the high sensitivity of MRS to field inhomogeneities. These poor quality spectra are prone to quantification and/or interpretation errors that can have a significant impact on the clinical use of spectroscopic data. Therefore, quality control of the spectra should always precede their clinical use. When performed manually, quality assessment of MRSI spectra is not only a tedious and time-consuming task, but is also affected by human subjectivity. Consequently, automatic, fast and reliable methods for spectral quality assessment are of utmost interest. In this article, we present a new random forest-based method for automatic quality assessment of (1) H MRSI brain spectra, which uses a new set of MRS signal features. The random forest classifier was trained on spectra from 40 MRSI grids that were classified as acceptable or non-acceptable by two expert spectroscopists. To account for the effects of intra-rater reliability, each spectrum was rated for quality three times by each rater. The automatic method classified these spectra with an area under the curve (AUC) of 0.976. Furthermore, in the subset of spectra containing only the cases that were classified every time in the same way by the spectroscopists, an AUC of 0.998 was obtained. Feature importance for the classification was also evaluated. Frequency domain skewness and kurtosis, as well as time domain signal-to-noise ratios (SNRs) in the ranges 50-75 ms and 75-100 ms, were the most important features. Given that the method is able to assess a whole MRSI grid faster than a spectroscopist (approximately 3 s versus approximately 3 min), and without loss of accuracy (agreement between classifier trained with just one session and any of the other labelling sessions, 89.88%; agreement between any two labelling sessions, 89.03%), the authors suggest its implementation in the clinical routine. The method presented in this article was implemented in jMRUI's SpectrIm plugin. Copyright © 2016 John Wiley & Sons, Ltd.

Reliability studies of the skew normal distribution

It has been observed in various practical applications that data do not conform to the normal distribution, which is symmetric with no skewness. The skew normal distribution proposed by Azzalini(1985) is appropriate for the analysis of data which is unimodal but exhibits some skewness. The skew normal distribution includes the normal distribution as a special case where the skewness parameter is zero. In this thesis, we study the structural properties of the skew normal distribution, with an emphasis on the reliability properties of the model. More specifically, we obtain the failure rate, the mean residual life function, and the reliability function of a skew normal random variable. We also compare it with the normal distribution with respect to certain stochastic orderings. Appropriate machinery is developed to obtain the reliability of a component when the strength and stress follow the skew normal distribution. Finally, IQ score data from Roberts (1988) is analyzed to illustrate the procedure.

Grain size analysis of sediment core CRP-3 from the Ross Sea, Antarctica

Grain-size analyses by sieve and Sedigraph are presented for 115 samples of core from CRP-3, 12 km off the coast of south Victoria Land. The data provide a useful check on visual core descriptions. The geographic setting for the strata sampled, some 790 m of early Oligocene nearshore marine sediments with a persistent glacial influence, is reviewed, and sediment textures interpreted in that context. Sand textures from the CRP-3 samples in the lower part of the core suggest that deposition was initially primarily wave-dominated, but that at times the influence of the waves was over-ridden by episodes of rapid sedimentation. Sedimentary cycles, recognised in the visual description of the core above 485 mbsf, show an increasing proportion of mudstone in the middle of each cycle above 330 mbsf that is interpreted to record periodic sedimentation in deeper water. Sandstone textures in the lower and upper parts of each cycle are interpreted to record departure from and return to shoreface deposition with changes in sea level. Mudstone textures above 176 mbsf indicate sedimentation below wave base. Many of the textures in both sand and mud samples show the coarse 'tail' characteristic of ice-rafted debris, but others do not, indicating ice-free periods. Many sandstones below c. 200 mbsf have virtually no silt, but significant amounts of clay (6 to 17%) that is thought to be of post-depositional origin.

Grain size analysis of sediment coring site CRP-2 from the Ross Sea, Antarctica

The purpose of this note is to present results of grain size analyses from 118 samples of the CRP-2/2A core using sieve and Sedigraph techniques. The samples were selected to represent the range of facies encountered, and tend to become more widely spaced with depth. Fifteen came from the upper 27 m of Quaternary and Pliocene sediments, 62 from the early Miocene-late Oligocene strata (27 to 307 mbsf), and 41 from the early Oligocene strata beneath (307 to 624 mbsf). The results are intended to provide reference data for lithological descriptions in the core logs (Cape Roberts Science Team, 1999), and to help with facies interpretation. The analytical technique used for determining size frequency of the sand fraction in our samples (sieving) is simple, physical and widely practised for over a century. Thus it provides a useful reference point for analyses produced by other faster and more sophisticated techniques, such as the Malvern laser particle size analysis system (Woolfe et al., 2000), and estimates derived from measurements taken with down-hole logging tools (Bücker, pers. com., 1999).

Description and chemical composition of ferromanganese deposits from Trout Lake, northern Wisconsin

Sediments were collected with Eckman and Petersen dredges from the bottom of Trout Lake, northern Wisconsin, at 221 stations. Sampling was done with a spud sampler at 32 stations, and core samples were obtained with a Jenkins and Mortimer and a Twenhofel sampler at 17 stations. The shore and offshore deposits of the shores of Trout Lake and the shores of the islands are described. Megascopic descriptions are given of the samples collected with the Eckman and Petersen dredges. Sediments on bottoms of about 10 meters or deeper are mainly gyttja, or crusts composed of mixtures of organic matter, ferric hydroxide, and some form of manganese oxide. The latter deposits are extensive. Detailed descriptions of some of the samples of sands are given, and generalizations respecting size and distribution are made. Tables showing quartiles, medians, and coefficients of sorting and skewness of the coarse sediments collected from the bottom are given in tables. Mechanical analyses of all fine sediments, mainly gyttja, were not made, as previous experience seems to have demonstrated that results have no sedimentational value. Organic matter of the gyttja was determined and also the percentages of lignin in the organic matter. Core samples are composed almost entirely of fine materials, mainly gyttja, and determinations were made on these samples in the same way as on the samples obtained with the Eckman and Petersen dredges. Studies of the core samples show that the fine sediments usually contain in excess of 90 per cent moisture and there is very little change in the moisture content from top to bottom of cores. A map shows the distribution of the iron and manganese deposits. These deposits were found to contain 10 to 20 per cent of organic matter, 11 to 16 per cent of metallic iron, and 12 to 30 per cent of metallic manganese. No stratification of any kind was found in any of the deep-water sediments of Trout Lake except in the iron and manganese crusts. Absence of stratification is considered to be due to the slow rate of deposition and the mixing of sediments by organisms which dwell in them. The data indicate that the rate of deposition in the deep waters of Trout Lake is of the order of 1 foot in 15,000 years.

(Table 1) Deuterium excess record during different time periods in Lomonosovfonna ice core

In this paper we present a deuterium excess (d) record from an ice core drilled on a small ice cap in Svalbard in 1997. The core site is located at Lomonosovfonna at 1255 m asl, and the analyzed time series spans the period 1400-1990 A.D. The record shows pronounced multidecadal to centennial-scale variations coherent with sea surface temperature changes registered in the subtropical to southern middle-latitude North Atlantic during the instrumental period. We interpret the negative trend in the deuterium excess during the 1400s and 1500s as an indication of cooling in the North Atlantic associated with the onset of the Little Ice Age. Consistently positive anomalies of d after 1900, peaking at about 1950, correspond with well-documented contemporary warming. Yet the maximum values of deuterium excess during 1900-1990 are not as high as in the early part of the record (pre-1550). This suggests that the sea surface temperatures during this earlier period of time in the North Atlantic to the south of approximately 45°N were at least comparable with those registered in the 20th century before the end of the 1980s. We examine the potential for a cold bias to exist in the deuterium excess record due to increased evaporation from the local colder sources of moisture having isotopically cold signature. It is argued that despite a recent oceanic warming, the contribution from this local moisture to the Lomonosovfonna precipitation budget is still insufficient to interfere with the isotopic signal from the primary moisture region in the midlatitude North Atlantic.

Sedimentology and mineralogy of surface sediments from the Chilean continental margin

The regional patterns of texture and composition of modern continental slope and pelagic sediments off Chile between 25°S and 43°S reflect the latitudinal segmentation of geological, morphological, and climatic features of the continental hinterland. Grain-size characteristics are controlled by the grain-size of source rocks, the weathering regime, and mode of sediment input (eolian off northern Chile vs fluvial further south). Bulk-mineral assemblages reveal a low grade of maturity. Regional variations are governed by the source-rock composition of the different geological terranes and the relative source-rock contribution of the Coastal Range and Andes, as controlled by the continental hydrology. The relative abundance of clay minerals is also predominantly influenced by the source-rock composition and partly by continental smectite neoformation. Latitudinal variations of illite crystallinities along the Chilean continental slope (and west of the Peru-Chile trench) clearly reflect modifications of the weathering regime which correspond to the strong climatic zonation of Chile.

Micro-granulometry of sedimentary dykes, Kathmandu basin, Nepal

Soft-sediment deformation structures have been analyzed at six sites of the Kathmandu valley. Microgranulometric study (this Supplement and Fig. 3B of Mugnier et al., Tectonophysics, 2011) reveals that silty levels (60 to 80% silt) favor the development of soft-sediment deformation structures, while sandy levels (60 to 80% sand) are passively deformed. Nonetheless well sorted sand levels (more than 80% sand) generate over-fluid pressure during compaction if located beneath a silty cap, leading to fluidization and dike development. 3-D geometry of seismites indicates a very strong horizontal shearing during their development. Using a physical approach based on soil liquefaction during horizontal acceleration, we show that the fluidization zone progressively grows down-section during the shaking, but does not exactly begin at the surface. The comparison of bed-thickness and strength/depth evolution indicates three cases: i) no soft-sediment deformation occurs for thin (few centimeters) silty beds; ii) the thickness of soft-sediment deformation above sandy beds is controlled by the lithological contrast; iii) the thickness of soft-sediment deformation depends on the shaking intensity for very thick silty beds. These 3 cases are evidenced in the Kathmandu basin. We use the 30 cm-thick soft-sediment deformation level formed during the 1833 earthquake as a reference: the 1833 earthquake rupture zone extended very close to Kathmandu, inducing there MMI IX-X damages. A 90 cm-thick sediment deformation has therefore to be induced by an event greater than MMI X. From a compilation of paleo and historic seismology studies, it is found that the great (M ~ 8.1) historical earthquakes are not characteristic of the greatest earthquakes of Himalaya; hence earthquakes greater than M ~ 8.6 occurred. Kathmandu is located above one of the asperities that laterally limits the extent of mega-earthquake ruptures and two successive catastrophic events already affected Kathmandu, in 1255 located to the west of this asperity and in ~ 1100 to the east.

Geochemistry, clay mineralogy and grain size distribution in various sediment cores drilled during ODP Leg 119 on the Antarctic continental shelf off Prydz Bay

Drilling was undertaken at five sites (739-743) on ODP Leg 119 on a transect across the continental shelf of Prydz Bay, East Antarctica, to elucidate the long-term glacial history of the area and to examine the importance of the area with respect to the development of the East Antarctic ice sheet as a whole. In addition to providing a record of glaciation spanning 36 m.y. or more, Leg 119 has provided information concerning the development of a continental margin under the prolonged influence of a major ice sheet. This has allowed the development of a sedimentary model that may be applicable not only to other parts of the Antarctic continental margin, but also to northern high-latitude continental shelves. The cored glacial sedimentary record in Prydz Bay consists of three major sequences, dominated by diamictite: 1. An upper flat-lying sequence that ranges in thickness from a few meters in inner and western Prydz Bay to nearly 250 m in the outer or eastern parts of the bay. The uppermost few meters consist of Holocene diatom ooze and diatomaceous mud with a minor ice-rafted component overlying diamicton and diamictite of late Miocene to Quaternary age. The diamictite is mainly massive, but stratified varieties and minor mudstone and diatomite also occur. 2. An upper prograding sequence cored at Sites 739 and 743, unconformly below the flat-lying sequence. This consists of a relatively steep (4° inclination) prograding wedge with a number of discrete sedimentary packages. At Sites 739 and 743 the sequence is dominated by massive and stratified diamictite, some of which shows evidence of slumping and minor debris flowage. 3. A lower, more gently inclined, prograding sequence lies unconformably below the flat-lying sequence at Site 742 and the upper prograding sequence at Site 739. This extends to the base of both sites, to 316 and 487 mbsf, respectively. It is dominated by massive, relatively clast-poor diamictite which is kaolinite-rich, light in color, and contains sporadic carbonate-cemented layers. The lower part of Site 742 includes well-stratified diamictites and very poorly sorted mudstones. The base of this site has indications of large-scale soft-sediment deformation and probably represents proximity to the base of the glacial sequence. Facies analysis of the Prydz Bay glacial sequence indicates a range of depositional environments. Massive diamictite is interpreted largely as waterlain till, deposited close to the grounding line of a floating glacier margin, although basal till and debris flow facies are also present. Weakly stratified diamictite is interpreted as having formed close to or under the floating ice margin and influenced by the input of marine diatomaceous sediment (proximal glaciomarine setting). Well-stratified diamictite has a stronger marine input, being more diatom-rich, and probably represents a proximal-distal glaciomarine sediment with the glaciogenic component being supplied by icebergs. Other facies include a variety of mudstones and diatom-rich sediments of marine origin, in which an ice-rafted component is still significant. None of the recovered sediments are devoid of a glacial influence. The overall depositional setting of the prograding sequence is one in which the grounded ice margin is situated close to the shelf edge. Progradation was achieved primarily by deposition of waterlain till. The flat-lying sequence illustrates a complex sequence of advances and retreats across the outer part of the shelf, with intermittent phases of ice loading and erosion. The glacial chronology is based largely on diatom stratigraphy, which has limited resolution. It appears that ice reached the paleoshelf break by earliest Oligocene, suggesting full-scale development of the East Antarctic ice sheet by that time. The ice sheet probably dominated the continental margin for much of Oligocene to middle Miocene time. Retreat, but not total withdrawal of the ice sheet, took place in late Miocene to mid-Pliocene time. The late Pliocene to Pleistocene was characterized by further advances across, and progradation of, the continental shelf. Holocene time has been characterized by reduced glacial conditions and a limited influence of glacial processes on sedimentation.