Heat kernel smoothing via Laplace-Beltrami eigenfunctions and its application to subcortical structure modeling

We present a new subcortical structure shape modeling framework using heat kernel smoothing constructed with the Laplace-Beltrami eigenfunctions. The cotan discretization is used to numerically obtain the eigenfunctions of the Laplace-Beltrami operator along the surface of subcortical structures of the brain. The eigenfunctions are then used to construct the heat kernel and used in smoothing out measurements noise along the surface. The proposed framework is applied in investigating the influence of age (38-79 years) and gender on amygdala and hippocampus shape. We detected a significant age effect on hippocampus in accordance with the previous studies. In addition, we also detected a significant gender effect on amygdala. Since we did not find any such differences in the traditional volumetric methods, our results demonstrate the benefit of the current framework over traditional volumetric methods.

Gender differences in the relationship between depression and suicidal ideation in young adolescents

Objective: This study examined the risk relationship between depressive symptomatology and suicidal ideation for young adolescent males and females. Method: A large cohort of students in their first year of high school completed the Center for Epidemiological Studies Depression Scale (CES-D) and the Adolescent Suicide Questionnaire. The risk relationship between depressive symptomatology and suicidal ideation was modelled using non-parametric kernel-smoothing techniques. Results: Suicidal ideation was more frequently reported by females compared with males which was partly explained by females having higher mean depression scores. At moderate levels of depression females also had a significantly higher risk of suicidal ideation compared with males and this increased risk contributed to the overall higher levels of female ideation. Conclusions: The risk relationship between depressive symptomatology and suicidal ideation is different for young adolescent males and females. The results indicate that moderate levels of depressive symptomatology can be associated with suicidal ideation (especially among young females) and that for these young people a suicide risk assessment is required.

Estimation of dynamic latent variable models using simulated nonparametric moments

Given a model that can be simulated, conditional moments at a trial parameter value can be calculated with high accuracy by applying kernel smoothing methods to a long simulation. With such conditional moments in hand, standard method of moments techniques can be used to estimate the parameter. Since conditional moments are calculated using kernel smoothing rather than simple averaging, it is not necessary that the model be simulable subject to the conditioning information that is used to define the moment conditions. For this reason, the proposed estimator is applicable to general dynamic latent variable models. Monte Carlo results show that the estimator performs well in comparison to other estimators that have been proposed for estimation of general DLV models.

Estimation of dynamic latent variable models using simulated nonparametric moments

Abstract. Given a model that can be simulated, conditional moments at a trial parameter value can be calculated with high accuracy by applying kernel smoothing methods to a long simulation. With such conditional moments in hand, standard method of moments techniques can be used to estimate the parameter. Because conditional moments are calculated using kernel smoothing rather than simple averaging, it is not necessary that the model be simulable subject to the conditioning information that is used to define the moment conditions. For this reason, the proposed estimator is applicable to general dynamic latent variable models. It is shown that as the number of simulations diverges, the estimator is consistent and a higher-order expansion reveals the stochastic difference between the infeasible GMM estimator based on the same moment conditions and the simulated version. In particular, we show how to adjust standard errors to account for the simulations. Monte Carlo results show how the estimator may be applied to a range of dynamic latent variable (DLV) models, and that it performs well in comparison to several other estimators that have been proposed for DLV models.

Spatial distribution of Schistosoma mansoni infection before and after chemotherapy with two praziquantel doses in a community of Pernambuco, Brazil

Praziquantel chemotherapy has been the focus of the Schistosomiasis Control Program in Brazil for the past two decades. Nevertheless, information on the impact of selective chemotherapy against Schistosoma mansoni infection under the conditions confronted by the health teams in endemic municipalities remains scarce. This paper compares the spatial pattern of infection before and after treatment with either a 40 mg/kg or 60 mg/kg dose of praziquantel by determining the intensity of spatial cluster among patients at 180 and 360 days after treatment. The spatial-temporal distribution of egg-positive patients was analysed in a Geographic Information System using the kernel smoothing technique. While all patients became egg-negative after 21 days, 17.9% and 30.9% reverted to an egg-positive condition after 180 and 360 days, respectively. Both the prevalence and intensity of infection after treatment were significantly lower in the 60 mg/kg than in the 40 mg/kg treatment group. The higher intensity of the kernel in the 40 mg/kg group compared to the 60 mg/kg group, at both 180 and 360 days, reflects the higher number of reverted cases in the lower dose group. Auxiliary, preventive measures to control transmission should be integrated with chemotherapy to achieve a more enduring impact.

A matter of scale: apparent niche differentiation of diploid and tetraploid plants may depend on extent and grain of analysis

Aim: Emerging polyploids may depend on environmental niche shifts for successful establishment. Using the alpine plant Ranunculus kuepferi as a model system, we explore the niche shift hypothesis at different spatial resolutions and in contrasting parts of the species range. Location: European Alps. Methods: We sampled 12 individuals from each of 102 populations of R. kuepferi across the Alps, determined their ploidy levels, derived coarse-grain (100x100m) environmental descriptors for all sampling sites by downscaling WorldClim maps, and calculated fine-scale environmental descriptors (2x2m) from indicator values of the vegetation accompanying the sampled individuals. Both coarse and fine-scale variables were further computed for 8239 vegetation plots from across the Alps. Subsequently, we compared niche optima and breadths of diploid and tetraploid cytotypes by combining principal components analysis and kernel smoothing procedures. Comparisons were done separately for coarse and fine-grain data sets and for sympatric, allopatric and the total set of populations. Results: All comparisons indicate that the niches of the two cytotypes differ in optima and/or breadths, but results vary in important details. The whole-range analysis suggests differentiation along the temperature gradient to be most important. However, sympatric comparisons indicate that this climatic shift was not a direct response to competition with diploid ancestors. Moreover, fine-grained analyses demonstrate niche contraction of tetraploids, especially in the sympatric range, that goes undetected with coarse-grained data. Main conclusions: Although the niche optima of the two cytotypes differ, separation along ecological gradients was probably less decisive for polyploid establishment than a shift towards facultative apomixis, a particularly effective strategy to avoid minority cytotype exclusion. In addition, our results suggest that coarse-grained analyses overestimate niche breadths of widely distributed taxa. Niche comparison analyses should hence be conducted at environmental data resolutions appropriate for the organism and question under study.

Nonparametric Instrumental Regression

The focus of the paper is the nonparametric estimation of an instrumental regression function P defined by conditional moment restrictions stemming from a structural econometric model : E[Y-P(Z)|W]=0 and involving endogenous variables Y and Z and instruments W. The function P is the solution of an ill-posed inverse problem and we propose an estimation procedure based on Tikhonov regularization. The paper analyses identification and overidentification of this model and presents asymptotic properties of the estimated nonparametric instrumental regression function.

Climate variability and socio-environmental changes in the northern Aegean (NE Mediterranean) during the last 1500 years

We provide new evidence on sea surface temperature (SST) variations and paleoceanographic/paleoenvironmental changes over the past 1500 years for the north Aegean Sea (NE Mediterranean). The reconstructions are based on multiproxy analyses, obtained from the high resolution (decadal to multidecadal) marine record M2 retrieved from the Athos basin. Reconstructed SSTs show an increase from ca. 850 to 950 AD and from ca. 1100 to 1300 AD. A cooling phase of almost 1.5 �C is observed from ca. 1600 AD to 1700 AD. This seems to have been the starting point of a continuous SST warming trend until the end of the reconstructed period, interrupted by two prominent cooling events at 1832 ± 15 AD and 1995 ± 1 AD. Application of an adaptive Kernel smoothing suggests that the current warming in the reconstructed SSTs of the north Aegean might be unprecedented in the context of the past 1500 years. Internal variability in atmospheric/oceanic circulations systems as well as external forcing as solar radiation and volcanic activity could have affected temperature variations in the north Aegean Sea over the past 1500 years. The marked temperature drop of approximately ~2 �C at 1832 ± 15 yr AD could be related to the 1809 АD ‘unknown’ and the 1815 AD Tambora volcanic eruptions. Paleoenvironmental proxy-indices of the M2 record show enhanced riverine/continental inputs in the northern Aegean after ca. 1450 AD. The paleoclimatic evidence derived from the M2 record is combined with a socio-environmental study of the history of the north Aegean region. We show that the cultivation of temperature-sensitive crops, i.e. walnut, vine and olive, co-occurred with stable and warmer temperatures, while its end coincided with a significant episode of cooler temperatures. Periods of agricultural growth in Macedonia coincide with periods of warmer and more stable SSTs, but further exploration is required in order to identify the causal links behind the observed phenomena. The Black Death likely caused major changes in agricultural activity in the north Aegean region, as reflected in the pollen data from land sites of Macedonia and the M2 proxy-reconstructions. Finally, we conclude that the early modern peaks in mountain vegetation in the Rhodope and Macedonia highlands, visible also in the M2 record, were very likely climate-driven.

How to not inflate population estimates? Spatial density distribution of white-lipped peccaries in a continuous Atlantic forest

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Evaluating different methods of microarray data normalization

Abstract Background With the development of DNA hybridization microarray technologies, nowadays it is possible to simultaneously assess the expression levels of thousands to tens of thousands of genes. Quantitative comparison of microarrays uncovers distinct patterns of gene expression, which define different cellular phenotypes or cellular responses to drugs. Due to technical biases, normalization of the intensity levels is a pre-requisite to performing further statistical analyses. Therefore, choosing a suitable approach for normalization can be critical, deserving judicious consideration. Results Here, we considered three commonly used normalization approaches, namely: Loess, Splines and Wavelets, and two non-parametric regression methods, which have yet to be used for normalization, namely, the Kernel smoothing and Support Vector Regression. The results obtained were compared using artificial microarray data and benchmark studies. The results indicate that the Support Vector Regression is the most robust to outliers and that Kernel is the worst normalization technique, while no practical differences were observed between Loess, Splines and Wavelets. Conclusion In face of our results, the Support Vector Regression is favored for microarray normalization due to its superiority when compared to the other methods for its robustness in estimating the normalization curve.

Stereological Modelling of Random Particles

Stochastic models for three-dimensional particles have many applications in applied sciences. Lévy–based particle models are a flexible approach to particle modelling. The structure of the random particles is given by a kernel smoothing of a Lévy basis. The models are easy to simulate but statistical inference procedures have not yet received much attention in the literature. The kernel is not always identifiable and we suggest one approach to remedy this problem. We propose a method to draw inference about the kernel from data often used in local stereology and study the performance of our approach in a simulation study.

A Bootstrap-Regression Procedure to Capture Unit Specific Effects in Data Envelopment Analysis

The Data Envelopment Analysis (DEA) efficiency score obtained for an individual firm is a point estimate without any confidence interval around it. In recent years, researchers have resorted to bootstrapping in order to generate empirical distributions of efficiency scores. This procedure assumes that all firms have the same probability of getting an efficiency score from any specified interval within the [0,1] range. We propose a bootstrap procedure that empirically generates the conditional distribution of efficiency for each individual firm given systematic factors that influence its efficiency. Instead of resampling directly from the pooled DEA scores, we first regress these scores on a set of explanatory variables not included at the DEA stage and bootstrap the residuals from this regression. These pseudo-efficiency scores incorporate the systematic effects of unit-specific factors along with the contribution of the randomly drawn residual. Data from the U.S. airline industry are utilized in an empirical application.

Planktic and benthic 14C reservoir ages, calibrated by a suite of 14C plateaus in the glacial-to-deglacial Suigetsu atmospheric 14C record

We here present a compilation of planktic and benthic 14C reservoir ages for the Last Glacial Maximum (LGM) and early deglacial from 11 key sites of global ocean circulation in the Atlantic and Indo-Pacific Ocean. The ages were obtained by 14C plateau tuning, a robust technique to derive both an absolute chronology for marine sediment records and a high-resolution record of changing reservoir/ventilation ages (Delta14C values) for surface and deep waters by comparing the suite of planktic 14C plateaus of a sediment record with that of the atmospheric 14C record (Sarnthein et al., 2007, doi:10.1029/173GM13). Results published thus far used as atmospheric 14C reference U/Th-dated corals, the Cariaco planktic record, and speleothems (Fairbanks et al., 2005, doi:10.1016/j.quascirev.2005.04.007; Hughen et al., 2006, doi:10.1016/j.quascirev.2006.03.014; Beck et al., 2001, doi:10.1023/A:1008175728826). We have now used the varve-counted atmospheric 14C record of Lake Suigetsu terrestrial macrofossils (Ramsey et al., 2012, doi:10.1126/science.1226660) to recalibrate the boundary ages and reservoir ages of the seven published records directly to an atmospheric 14C record. In addition, the results for four new cores and further planktic results for four published records are given. Main conclusions from the new compilation are: (1) The Suigetsu atmospheric 14C record on its varve counted time scale reflects all 14C plateaus, their internal structures and relative length previously identified, but implies a rise in the average 14C plateau age by 200-700 14C yr during LGM and early deglacial times. (2) Based on different 14C ages of coeval atmospheric and planktic 14C plateaus, marine surface water Delta14C may have temporarily dropped to an equivalent of ~0 yr in low-latitude lagoon waters, but reached >2500 14C yr both in stratified subpolar waters and in upwelled waters such as in the South China Sea. These values differ significantly from a widely assumed constant global planktic Delta14C value of 400 yr. (3) Suites of deglacial planktic Delta14C values are closely reproducible in 14C records measured at neighboring core sites. (4) Apparent deep-water 14C ventilation ages (equivalents of benthic Delta14C), deduced from the sum of planktic Delta14C and coeval benthic-planktic 14C differences, vary from 500 up to >5000 yr in LGM and deglacial ocean basins.