Seasonal temperature variability of the Neoglacial (3300–2500 BP) and Roman Warm Period (2500–1600 BP) reconstructed from oxygen isotope ratios of limpet shells (Patella vulgata), Northwest Scotland

Seasonal sea-surface temperaturevariability for the Neoglacial (3300–2500 BP) and Roman WarmPeriod (RWP; 2500–1600 BP), which correspond to the Bronze and Iron Ages, respectively, was estimated using oxygen isotope ratios obtained from high-resolution samples micromilled from radiocarbon-dated, archaeological limpet (Patella vulgata) shells. The coldest winter months recorded in Neoglacial shells averaged 6.6 ± 0.3 °C, and the warmest summer months averaged 14.7 ± 0.4 °C. One Neoglacial shell captured a year without a summer, which may have resulted from a dust veil from a volcanic eruption in the Katla volcanic system in Iceland. RWP shells record average winter and summer monthly temperatures of 6.3 ± 0.1 °C and 13.3 ± 0.3 °C, respectively. These results capture a cooling transition from the Neoglacial to RWP, which is further supported by earlier studies of pine history in Scotland, pollen type analyses in northeast Scotland, and European glacial events. The cooling transition observed at the boundary between the Neoglacial and RWP in our study also agrees with the abrupt climate deterioration at 2800–2700 BP (also referred to as the Subboreal/Subatlantic transition) and therefore may have been driven by decreased solar radiation and weakened North Atlantic Oscillation conditions.

Oxygen isotope ratios in the shell of Mytilus edulis: archives of glacier meltwater in Greenland?

Melting of the Greenland Ice Sheet (GrIS) is accelerating and will contribute significantly to global sea level rise during the 21st century. Instrumental data on GrIS melting only cover the last few decades, and proxy data extending our knowledge into the past are vital for validating models predicting the influence of ongoing climate change. We investigated a potential meltwater proxy in Godthåbsfjord (West Greenland), where glacier meltwater causes seasonal excursions with lower oxygen isotope water (δ18Ow) values and salinity. The blue mussel (Mytilus edulis) potentially records these variations, because it precipitates its shell calcite in oxygen isotopic equilibrium with ambient seawater. As M. edulis shells are known to occur in raised shorelines and archaeological shell middens from previous Holocene warm periods, this species may be ideal in reconstructing past meltwater dynamics. We investigate its potential as a palaeo-meltwater proxy. First, we confirmed that M. edulis shell calcite oxygen isotope (δ18Oc) values are in equilibrium with ambient water and generally reflect meltwater conditions. Subsequently we investigated if this species recorded the full range of δ18Ow values occurring during the years 2007 to 2010. Results show that δ18Ow values were not recorded at very low salinities (< ~ 19), because the mussels appear to cease growing. This implies that Mytilus edulis δ18Oc values are suitable in reconstructing past meltwater amounts in most cases, but care has to be taken that shells are collected not too close to a glacier, but rather in the mid-region or mouth of the fjord. The focus of future research will expand on the geographical and temporal range of the shell measurements by sampling mussels in other fjords in Greenland along a south–north gradient, and by sampling shells from raised shorelines and archaeological shell middens from prehistoric settlements in Greenland.

Error, reproducibility and sensitivity: a pipeline for data processing of Agilent oligonucleotide expression arrays

Background: Expression microarrays are increasingly used to obtain large scale transcriptomic information on a wide range of biological samples. Nevertheless, there is still much debate on the best ways to process data, to design experiments and analyse the output. Furthermore, many of the more sophisticated mathematical approaches to data analysis in the literature remain inaccessible to much of the biological research community. In this study we examine ways of extracting and analysing a large data set obtained using the Agilent long oligonucleotide transcriptomics platform, applied to a set of human macrophage and dendritic cell samples. Results: We describe and validate a series of data extraction, transformation and normalisation steps which are implemented via a new R function. Analysis of replicate normalised reference data demonstrate that intrarray variability is small (only around 2 of the mean log signal), while interarray variability from replicate array measurements has a standard deviation (SD) of around 0.5 log(2) units (6 of mean). The common practise of working with ratios of Cy5/Cy3 signal offers little further improvement in terms of reducing error. Comparison to expression data obtained using Arabidopsis samples demonstrates that the large number of genes in each sample showing a low level of transcription reflect the real complexity of the cellular transcriptome. Multidimensional scaling is used to show that the processed data identifies an underlying structure which reflect some of the key biological variables which define the data set. This structure is robust, allowing reliable comparison of samples collected over a number of years and collected by a variety of operators. Conclusions: This study outlines a robust and easily implemented pipeline for extracting, transforming normalising and visualising transcriptomic array data from Agilent expression platform. The analysis is used to obtain quantitative estimates of the SD arising from experimental (non biological) intra- and interarray variability, and for a lower threshold for determining whether an individual gene is expressed. The study provides a reliable basis for further more extensive studies of the systems biology of eukaryotic cells.

Linear vs. log-linear unit-root specification: an application of mis-specification encompassing

The objective of this paper is to apply the mis-specification (M-S) encompassing perspective to the problem of choosing between linear and log-linear unit-root models. A simple M-S encompassing test, based on an auxiliary regression stemming from the conditional second moment, is proposed and its empirical size and power are investigated using Monte Carlo simulations. It is shown that by focusing on the conditional process the sampling distributions of the relevant statistics are well behaved under both the null and alternative hypotheses. The proposed M-S encompassing test is illustrated using US total disposable income quarterly data.

1H-13C HSQC NMR spectroscopy for estimating procyanidin/prodelphinidin and cis/trans-flavan-3-ol ratios of condensed tannin samples: correlation with thiolysis

Studies with a diverse array of 22 purified condensed tannin (CT) samples from nine plant species demonstrated that procyanidin/prodelphinidin (PC/PD) and cis/trans-flavan-3-ol ratios can be appraised by 1H-13C HSQC NMR spectroscopy. The method was developed from samples containing 44 to ~100% CT, PC/PD ratios ranging from 0/100 to 99/1, and cis/trans ratios from 58/42 to 95/5 as determined by thiolysis with benzyl mercaptan. Integration of cross-peak contours of H/C-6' signals from PC and of H/C-2',6' signals from PD yielded nuclei adjusted estimates that were highly correlated with PC/PD ratios obtained by thiolysis (R2 = 0.99). Cis/trans-flavan-3-ol ratios, obtained by integration of the respective H/C-4 cross-peak contours, were also related to determinations made by thiolysis (R2 = 0.89). Overall, 1H-13C HSQC NMR spectroscopy appears to be a viable alternative to thiolysis for estimating PC/PD and cis/trans ratios of CT, if precautions are taken to avoid integration of cross-peak contours of contaminants.

Generalised golden ratios over integer alphabets

It is a well known result that for β ∈ (1,1+√52) and x ∈ (0,1β−1) there exists uncountably many (ǫi)∞i=1 ∈ {0,1}N such that x = P∞i=1ǫiβ−i. When β ∈ (1+√52,2] there exists x ∈ (0,1β−1) for which there exists a unique (ǫi)∞i=1 ∈ {0,1}N such that x=P∞i=1ǫiβ−i. In this paper we consider the more general case when our sequences are elements of {0, . . . , m}N. We show that an analogue of the golden ratio exists and give an explicit formula for it.

Ratios of UV, PAR and NIR components to global solar radiation measured at Botucatu site in Brazil

The relationships between the four radiant fluxes are analyzed based on a 4 year data archive of hourly and daily global ultraviolet (I(UV)), photosynthetically active-PAR (I(PAR)), near infrared (I(NIR)) and broadband global solar radiation (I(G)) collected at Botucatu, Brazil. These data are used to establish both the fractions of spectral components to global solar radiation and the proposed linear regression models. Verification results indicated that the proposed regression models predict accurately the spectral radiant fluxes at least for the Brazilian environment. Finally, results obtained in this analysis agreed well with most published results in the literature. (c) 2010 Elsevier Ltd. All rights reserved.

Isometric immersions into a homogeneous Lorentzian Heisenberg group and rigidity

In this paper we prove an existence result for local and global isometric immersions of semi-Riemannian surfaces into the three dimensional Heisenberg group endowed with a homogeneous left-invariant Lorentzian metric. As a corollary, we prove a rigidity result for such immersions.

On estimation and influence diagnostics for log-Birnbaum-Saunders Student-t regression models: Full Bayesian analysis

The purpose of this paper is to develop a Bayesian approach for log-Birnbaum-Saunders Student-t regression models under right-censored survival data. Markov chain Monte Carlo (MCMC) methods are used to develop a Bayesian procedure for the considered model. In order to attenuate the influence of the outlying observations on the parameter estimates, we present in this paper Birnbaum-Saunders models in which a Student-t distribution is assumed to explain the cumulative damage. Also, some discussions on the model selection to compare the fitted models are given and case deletion influence diagnostics are developed for the joint posterior distribution based on the Kullback-Leibler divergence. The developed procedures are illustrated with a real data set. (C) 2010 Elsevier B.V. All rights reserved.

The log-exponentiated Weibull regression model for interval-censored data

In interval-censored survival data, the event of interest is not observed exactly but is only known to occur within some time interval. Such data appear very frequently. In this paper, we are concerned only with parametric forms, and so a location-scale regression model based on the exponentiated Weibull distribution is proposed for modeling interval-censored data. We show that the proposed log-exponentiated Weibull regression model for interval-censored data represents a parametric family of models that include other regression models that are broadly used in lifetime data analysis. Assuming the use of interval-censored data, we employ a frequentist analysis, a jackknife estimator, a parametric bootstrap and a Bayesian analysis for the parameters of the proposed model. We derive the appropriate matrices for assessing local influences on the parameter estimates under different perturbation schemes and present some ways to assess global influences. Furthermore, for different parameter settings, sample sizes and censoring percentages, various simulations are performed; in addition, the empirical distribution of some modified residuals are displayed and compared with the standard normal distribution. These studies suggest that the residual analysis usually performed in normal linear regression models can be straightforwardly extended to a modified deviance residual in log-exponentiated Weibull regression models for interval-censored data. (C) 2009 Elsevier B.V. All rights reserved.

Generalized log-gamma regression models with cure fraction

In this paper, the generalized log-gamma regression model is modified to allow the possibility that long-term survivors may be present in the data. This modification leads to a generalized log-gamma regression model with a cure rate, encompassing, as special cases, the log-exponential, log-Weibull and log-normal regression models with a cure rate typically used to model such data. The models attempt to simultaneously estimate the effects of explanatory variables on the timing acceleration/deceleration of a given event and the surviving fraction, that is, the proportion of the population for which the event never occurs. The normal curvatures of local influence are derived under some usual perturbation schemes and two martingale-type residuals are proposed to assess departures from the generalized log-gamma error assumption as well as to detect outlying observations. Finally, a data set from the medical area is analyzed.

Log-Burr XII regression models with censored data

In survival analysis applications, the failure rate function may frequently present a unimodal shape. In such case, the log-normal or log-logistic distributions are used. In this paper, we shall be concerned only with parametric forms, so a location-scale regression model based on the Burr XII distribution is proposed for modeling data with a unimodal failure rate function as an alternative to the log-logistic regression model. Assuming censored data, we consider a classic analysis, a Bayesian analysis and a jackknife estimator for the parameters of the proposed model. For different parameter settings, sample sizes and censoring percentages, various simulation studies are performed and compared to the performance of the log-logistic and log-Burr XII regression models. Besides, we use sensitivity analysis to detect influential or outlying observations, and residual analysis is used to check the assumptions in the model. Finally, we analyze a real data set under log-Buff XII regression models. (C) 2008 Published by Elsevier B.V.

First observation of the Cabibbo-suppressed decays Xi(+)(c)->Sigma(+)pi(-)pi(+) and Xi(+)(c)->Sigma(-)pi(+)pi(+) and measurement of their branching ratios

We report the first observation of two Cabibbo-suppressed ecay modes, Xi(+)(c) -> Sigma(+)pi(-)pi(+) and Xi c+ -> Sigma(-)pi(+)pi(+). We observe 59 +/- 14 over a background of 87, and 22 +/- 8 over a background of 13 events, respectively, for the signals. The data were accumulated using the SELEX spectrometer during the 1996-1997 fixed target run at Fermilab, chiefly from a 600 Gev/c Sigma(-) beam. The branching ratios of the decays relative to the Cabibbo-favored Xi c+ -> Xi(-)pi(+)pi(+) are measured to be B(Xi(+)(c) -> Sigma(+)pi(-)pi(+))/B(Xi(+)(c) -> Xi(-)pi(+)pi(+)) = 0.48 +/- 0.20, and B(Xi(+)(c) -> Sigma(-)pi(+)pi(+))/B(Xi(+)(c) -> Sigma(-)pi(+)pi(+)) = 0.18 +/- 0.09, respectively. We also report branching ratios for the same decay modes of the Delta(+)(c) relative to Delta(+)(c) -> pK(-)pi(+.) (C) 2008 Elsevier B.V. All rights reserved.

Metal location and thickness in a multilayered sheet by measuring K alpha/K beta, L alpha/L beta and L alpha/L gamma X-ray ratios

When a multilayered material is analyzed by means of energy-dispersive X-ray fluorescence analysis, then the X-ray ratios of K alpha/K beta, or L alpha/L beta and L alpha/L gamma, for an element in the multilayered material, depend on the composition and thickness of the layer in which the element is situated, and on the composition and thickness of the superimposed layer (or layers). Multilayered samples are common in archaeometry, for example, in the case of pigment layers in paintings, or in the case of gilded or silvered alloys. The latter situation is examined in detail in the present paper, with a specific reference to pre-Columbian alloys from various museums in the north of Peru. (C) 2009 Elsevier B.V. All rights reserved.