Representation errors and retrievals in linear and nonlinear data assimilation

This article shows how one can formulate the representation problem starting from Bayes’ theorem. The purpose of this article is to raise awareness of the formal solutions,so that approximations can be placed in a proper context. The representation errors appear in the likelihood, and the different possibilities for the representation of reality in model and observations are discussed, including nonlinear representation probability density functions. Specifically, the assumptions needed in the usual procedure to add a representation error covariance to the error covariance of the observations are discussed,and it is shown that, when several sub-grid observations are present, their mean still has a representation error ; socalled ‘superobbing’ does not resolve the issue. Connection is made to the off-line or on-line retrieval problem, providing a new simple proof of the equivalence of assimilating linear retrievals and original observations. Furthermore, it is shown how nonlinear retrievals can be assimilated without loss of information. Finally we discuss how errors in the observation operator model can be treated consistently in the Bayesian framework, connecting to previous work in this area.

Analysis of wall plasters and natural sediments from the Neolithic town of Çatalhöyük (Turkey) by a range of analytical techniques

Wall plaster sequences from the Neolithic town of Çatalhöyük have been analysed and compared to three types of natural sediment found in the vicinity of the site, using a range of analytical techniques. Block samples containing the plaster sequences were removed from the walls of several different buildings on the East Mound. Sub-samples were examined by IR spectroscopy, X-ray diffraction and X-ray fluorescence to determine the overall mineralogical and elemental composition, whilst thin sections were studied using optical polarising microscopy, IR Microscopy and Environmental Scanning Electron Microscopy with Energy Dispersive X-ray analysis. The results of this study have shown that there are two types of wall plaster found in the sequences and that the sediments used to produce these were obtained from at least two distinct sources. In particular, the presence of clay, calcite and magnesian calcite in the foundation plasters suggested that these were prepared predominantly from a marl source. On the other hand, the finishing plasters were found to contain dolomite with a small amount of clay and no calcite, revealing that softlime was used in their preparation. Whilst marl is located directly below and around Çatalhöyük, the nearest source of softlime is 6.5 km away, an indication that the latter was important to the Neolithic people, possibly due to the whiter colour (5Y 8/1) of this sediment. Furthermore, the same two plaster types were found on each wall of Building 49, the main building studied in this research, and in all five buildings investigated, suggesting that the use of these sources was an established practice for the inhabitants of several different households across the site.

Sensory profiles and consumer acceptability of a range of sugar-reduced products on the UK market

Current UK intake of non-milk extrinsic sugars (NMES) is above recommendations. Reducing the sugar content of processed high sugar foods through reformulation is one option for reducing consumption of NMES at a population level. However, reformulation can alter the sensory attributes of food products and influence consumer liking. This study evaluated consumer acceptance of a selection of products that are commercially-available in the UK; these included regular and sugar-reduced baked beans, strawberry jam, milk chocolate, cola and cranberry & raspberry juice. Sweeteners were present in the reformulated chocolate (maltitol), cola (aspartame and acesulfame-K) and juice (sucralose) samples. Healthy, non-smoking consumers (n = 116; 55 men, 61 women, age: 33 ± 9 years; BMI: 25.7 ± 4.6 kg/m2) rated the products for overall liking and on liking of appearance, flavor and texture using a nine-point hedonic scale. There were significant differences between standard and reduced sugar products in consumers’ overall liking and on liking of each modality (appearance, flavor and texture; all P < 0.0001). For overall liking, only the regular beans and cola were significantly more liked than their reformulated counterparts (P < 0.0001). Cluster analysis identified three consumer clusters that were representative of different patterns of consumer liking. For the largest cluster (cluster 3: 45%), there was a significant difference in mean liking scores across all products, except jam. Differences in liking were predominantly driven by sweet taste in 2 out of 3 clusters. The current research has demonstrated that a high proportion of consumers prefer conventional products over sugar-reduced products across a wide range of product types (45%) or across selected products (27%), when tasted unbranded, and so there is room for further optimization of commercial reduced sugar products that were evaluated in the current study. Future work should evaluate strategies to facilitate compliance to dietary recommendations on NMES and free sugars, such as the impact of sugar-reduced food exposure on their acceptance.

Bayesian model comparison with un-normalised likelihoods

Models for which the likelihood function can be evaluated only up to a parameter-dependent unknown normalizing constant, such as Markov random field models, are used widely in computer science, statistical physics, spatial statistics, and network analysis. However, Bayesian analysis of these models using standard Monte Carlo methods is not possible due to the intractability of their likelihood functions. Several methods that permit exact, or close to exact, simulation from the posterior distribution have recently been developed. However, estimating the evidence and Bayes’ factors for these models remains challenging in general. This paper describes new random weight importance sampling and sequential Monte Carlo methods for estimating BFs that use simulation to circumvent the evaluation of the intractable likelihood, and compares them to existing methods. In some cases we observe an advantage in the use of biased weight estimates. An initial investigation into the theoretical and empirical properties of this class of methods is presented. Some support for the use of biased estimates is presented, but we advocate caution in the use of such estimates.

Objective Assessment of Internal Nasal Dimensions and Speech Resonance in Individuals With Repaired Unilateral Cleft Lip and Palate After Rhinoseptoplasty

The objective of the current study was to analyze the effects of rhinoseptoplasty on internal nasal dimensions and speech resonance of individuals with unilateral cleft lip and palate, estimated by acoustic rhinometry and nasometry, respectively. Twenty-one individuals (aged 15-46 years) with previously repaired unilateral cleft lip and palate were analyzed before (PRE), and 6 to 9 (POST1) and 12 to 18 months (POST2) after surgery. Acoustic rhinometry was used to measure the cross-sectional areas (CSAs) of segments corresponding to the nasal valve (CSA1), anterior portion (CSA2), and posterior portion (CSA3) of the lower turbinate, and the volumes at the nasal valve (V1) and turbinate (V2) regions at cleft and noncleft sides, before and after nasal decongestion with a topical vasoconstrictor. Nasometry was used to evaluate speech nasalance during the reading of a set of sentences containing nasal sounds and other devoid of nasal sounds. At the cleft side, before nasal decongestion, there was a significant increase (P < 0.05) in mean CSA1 and V1 values at POST1 and POST2 compared with PRE. After decongestion, increased values were also observed for CSA2 and V2 at POST2. No significant changes were observed at the noncleft side. Mean nasalance values at PRE, POST1, an POST2 were not different from each other in both oral and nasal sentences. The measurement of CSAs and volumes by acoustic rhinometry revealed that rhinoseptoplasty provided, in most cases analyzed, a significant increase in nasal patency, without concomitant changes in speech resonance, as estimated by nasalance assessment.

Genetic algorithm inversion of the average 1D crustal structure using local and regional earthquakes

Knowing the best 1D model of the crustal and upper mantle structure is useful not only for routine hypocenter determination, but also for linearized joint inversions of hypocenters and 3D crustal structure, where a good choice of the initial model can be very important. Here, we tested the combination of a simple GA inversion with the widely used HYPO71 program to find the best three-layer model (upper crust, lower crust, and upper mantle) by minimizing the overall P- and S-arrival residuals, using local and regional earthquakes in two areas of the Brazilian shield. Results from the Tocantins Province (Central Brazil) and the southern border of the Sao Francisco craton (SE Brazil) indicated an average crustal thickness of 38 and 43 km, respectively, consistent with previous estimates from receiver functions and seismic refraction lines. The GA + HYPO71 inversion produced correct Vp/Vs ratios (1.73 and 1.71, respectively), as expected from Wadati diagrams. Tests with synthetic data showed that the method is robust for the crustal thickness, Pn velocity, and Vp/Vs ratio when using events with distance up to about 400 km, despite the small number of events available (7 and 22, respectively). The velocities of the upper and lower crusts, however, are less well constrained. Interestingly, in the Tocantins Province, the GA + HYPO71 inversion showed a secondary solution (local minimum) for the average crustal thickness, besides the global minimum solution, which was caused by the existence of two distinct domains in the Central Brazil with very different crustal thicknesses. (C) 2010 Elsevier Ltd. All rights reserved.

Bayesian analysis and constraints on kinematic models from union SNIa

The kinematic expansion history of the universe is investigated by using the 307 supernovae type Ia from the Union Compilation set. Three simple model parameterizations for the deceleration parameter ( constant, linear and abrupt transition) and two different models that are explicitly parametrized by the cosmic jerk parameter ( constant and variable) are considered. Likelihood and Bayesian analyses are employed to find best fit parameters and compare models among themselves and with the flat Lambda CDM model. Analytical expressions and estimates for the deceleration and cosmic jerk parameters today (q(0) and j(0)) and for the transition redshift (z(t)) between a past phase of cosmic deceleration to a current phase of acceleration are given. All models characterize an accelerated expansion for the universe today and largely indicate that it was decelerating in the past, having a transition redshift around 0.5. The cosmic jerk is not strongly constrained by the present supernovae data. For the most realistic kinematic models the 1 sigma confidence limits imply the following ranges of values: q(0) is an element of [-0.96, -0.46], j(0) is an element of [-3.2,-0.3] and z(t) is an element of [0.36, 0.84], which are compatible with the Lambda CDM predictions, q(0) = -0.57 +/- 0.04, j(0) = -1 and z(t) = 0.71 +/- 0.08. We find that even very simple kinematic models are equally good to describe the data compared to the concordance Lambda CDM model, and that the current observations are not powerful enough to discriminate among all of them.

Mitochondrial sequence data and Dipsacales phylogeny: Mixed models, partitioned Bayesian analyses, and model selection

Phylogenetic analyses of chloroplast DNA sequences, morphology, and combined data have provided consistent support for many of the major branches within the angiosperm, clade Dipsacales. Here we use sequences from three mitochondrial loci to test the existing broad scale phylogeny and in an attempt to resolve several relationships that have remained uncertain. Parsimony, maximum likelihood, and Bayesian analyses of a combined mitochondrial data set recover trees broadly consistent with previous studies, although resolution and support are lower than in the largest chloroplast analyses. Combining chloroplast and mitochondrial data results in a generally well-resolved and very strongly supported topology but the previously recognized problem areas remain. To investigate why these relationships have been difficult to resolve we conducted a series of experiments using different data partitions and heterogeneous substitution models. Usually more complex modeling schemes are favored regardless of the partitions recognized but model choice had little effect on topology or support values. In contrast there are consistent but weakly supported differences in the topologies recovered from coding and non-coding matrices. These conflicts directly correspond to relationships that were poorly resolved in analyses of the full combined chloroplast-mitochondrial data set. We suggest incongruent signal has contributed to our inability to confidently resolve these problem areas. (c) 2007 Elsevier Inc. All rights reserved.

Efficiency of wear and decalcification technique for estimating the age of estuarine dolphin Sotalia guianensis

Most techniques used for estimating the age of Sotalia guianensis (van B,n,den, 1864) (Cetacea; Delphinidae) are very expensive, and require sophisticated equipment for preparing histological sections of teeth. The objective of this study was to test a more affordable and much simpler method, involving of the manual wear of teeth followed by decalcification and observation under a stereomicroscope. This technique has been employed successfully with larger species of Odontoceti. Twenty-six specimens were selected, and one tooth of each specimen was worn and demineralized for growth layers reading. Growth layers were evidenced in all specimens; however, in 4 of the 26 teeth, not all the layers could be clearly observed. In these teeth, there was a significant decrease of growth layer group thickness, thus hindering the layers count. The juxtaposition of layers hindered the reading of larger numbers of layers by the wear and decalcification technique. Analysis of more than 17 layers in a single tooth proved inconclusive. The method applied here proved to be efficient in estimating the age of Sotalia guianensis individuals younger than 18 years. This method could simplify the study of the age structure of the overall population, and allows the use of the more expensive methodologies to be confined to more specific studies of older specimens. It also enables the classification of the calf, young and adult classes, which is important for general population studies.

Evolutionary history of Ramphastos toucans: Molecular phylogenetics, temporal diversification, and biogeography

The toucan genus Ramphastos (Piciformes: Ramphastidae) has been a model in the formulation of Neotropical paleobiogeographic hypotheses. Weckstein (2005) reported on the phylogenetic history of this genus based on three mitochondrial genes, but some relationships were weakly supported and one of the subspecies of R. vitellinus (citreolaemus) was unsampled. This study expands on Weckstein (2005) by adding more DNA sequence data (including a nuclear marker) and more samples, including R v. citreolaemus. Maximum parsimony, maximum likelihood, and Bayesian methods recovered similar trees, with nodes showing high support. A monophyletic R. vitellinus complex was strongly supported as the sister-group to R. brevis. The results also confirmed that the southeastern and northern populations of R. vitellinus ariel are paraphyletic. X v. citreolaemus is sister to the Amazonian subspecies of the vitellinus complex. Using three protein-coding genes (COI, cytochrome-b and ND2) and interval-calibrated nodes under a Bayesian relaxed-clock framework, we infer that ramphastid genera originated in the middle Miocene to early Pliocene, Ramphastos species originated between late Miocene and early Pleistocene, and intra-specific divergences took place throughout the Pleistocene. Parsimony-based reconstruction of ancestral areas indicated that evolution of the four trans-Andean Ramphastos taxa (R. v. citreolaemus, R. a. swainsonii, R. brevis and R. sulfuratus) was associated with four independent dispersals from the cis-Andean region. The last pulse of Andean uplift may have been important for the evolution of R. sulfuratus, whereas the origin of the other trans-Andean Ramphastos taxa is consistent with vicariance due to drying events in the lowland forests north of the Andes. Estimated rates of molecular evolution were higher than the ""standard"" bird rate of 2% substitutions/site/million years for two of the three genes analyzed (cytochrome-b and ND2). (C) 2009 Elsevier Inc. All rights reserved.

Samba, a Xenopus hnRNP Expressed in Neural and Neural Crest Tissues

RNA binding proteins regulate gene expression at the posttranscriptional level and play important roles in embryonic development. Here, we report the cloning and expression of Samba, a Xenopus hnRNP that is maternally expressed and persists at least until tail bud stages. During gastrula stages, Samba is enriched in the dorsal regions. Subsequently, its expression is elevated only in neural and neural crest tissues. In the latter, Samba expression overlaps with that of Slug in migratory neural crest cells. Thereafter, Samba is maintained in the neural crest derivatives, as well as other neural tissues, including the anterior and posterior neural tube and the eyes. Overexpression of Samba in the animal pole leads to defects in neural crest migration and cranial cartilage development. Thus, Samba encodes a Xenopus hnRNP that is expressed early in neural and neural crest derivatives and may regulate crest cells migratory behavior. Developmental Dynamics 238:204-209, 2009. (C) 2008 Wiley-Liss, Inc.

A Phylogenetic Lineage of Closely Related Trypanosomes (Trypanosomatidae, Kinetoplastida) of Anurans and Sand Flies (Psychodidae, Diptera) Sharing the Same Ecotopes in Brazilian Amazonia

Analysis of the phylogenetic relationships among trypanosomes from vertebrates and invertebrates disclosed a new lineage of trypanosomes circulating among anurans and sand flies that share the same ecotopes in Brazilian Amazonia. This assemblage of closely related trypanosomes was determined by comparing whole SSU rDNA sequences of anuran trypanosomes from the Brazilian biomes of Amazonia, the Pantanal, and the Atlantic Forest and from Europe, North America, and Africa, and from trypanosomes of sand flies from Amazonia. Phylogenetic trees based on maximum likelihood and parsimony corroborated the positioning of all new anuran trypanosomes in the aquatic clade but did not support the monophyly of anuran trypanosomes. However, all analyses always supported four major clades (An01-04) of anuran trypanosomes. Clade An04 is composed of trypanosomes from exotic anurans. Isolates in clades An01 and An02 were from Brazilian frogs and toads captured in the three biomes studied, Amazonia, the Pantanal and the Atlantic Forest. Clade An01 contains mostly isolates from Hylidae whereas clade An02 comprises mostly isolates from Bufonidae; and clade An03 contains trypanosomes from sand flies and anurans of Bufonidae, Leptodactylidae, and Leiuperidae exclusively from Amazonia. To our knowledge, this is the first study describing morphological and growth features, and molecular phylogenetic affiliation of trypanosomes from anurans and phlebotomines, incriminating these flies as invertebrate hosts and probably also as important vectors of Amazonian terrestrial anuran trypanosomes.

A compound class of Weibull and power series distributions

In this paper we introduce the Weibull power series (WPS) class of distributions which is obtained by compounding Weibull and power series distributions where the compounding procedure follows same way that was previously carried out by Adamidis and Loukas (1998) This new class of distributions has as a particular case the two-parameter exponential power series (EPS) class of distributions (Chahkandi and Gawk 2009) which contains several lifetime models such as exponential geometric (Adamidis and Loukas 1998) exponential Poisson (Kus 2007) and exponential logarithmic (Tahmasbi and Rezaei 2008) distributions The hazard function of our class can be increasing decreasing and upside down bathtub shaped among others while the hazard function of an EPS distribution is only decreasing We obtain several properties of the WPS distributions such as moments order statistics estimation by maximum likelihood and inference for a large sample Furthermore the EM algorithm is also used to determine the maximum likelihood estimates of the parameters and we discuss maximum entropy characterizations under suitable constraints Special distributions are studied in some detail Applications to two real data sets are given to show the flexibility and potentiality of the new class of distributions (C) 2010 Elsevier B V All rights reserved

Relationship between structure and photoactivity of porphyrins derived from protoporphyrin IX

Protoporphyrin (Pp IX) derivatives were prepared to study the relationship between photosensitizer structure and photoactivity, with an emphasis on understanding the role of membrane interactions in the efficiency of photosensitizers used in photodynamic therapy (PDT). The synthetic strategies described here aimed at changing protoporphyrin periferic groups, varying overall charge and oil/water partition, while maintaining their photochemical properties. Three synthetic routes were used: (1) modification of Pp IX at positions 3(1) and 8(1) by addition of alkyl amine groups of different lengths (compounds 2-5), (2) change of Pp IX at positions 13(3) and 17(3), generating alkyl amines (compounds 6 and 7), a phosphate amine (compound 8), and quarternary ammonium compounds (compounds 9 and 10), and (3) amine-alkylation of Hematoporphyrin IX (Hp IX) at positions 3(1), 8(1), 13(3) and 17(3) (compound 12). Strategy 1 leads to hydrophobic compounds with low photocytotoxicity. Strategy 2 leads to compounds 6-10 that have high levels of binding/incorporation in vesicles, mitochondria and cells, which are indicative of high bioavailability. Addition of the phosphate group (compound 8), generates an anionic compound that has low liposome and cell incorporation, plus low photocytotoxicity. Compound 12 has intermediate incorporation and photocytotoxic properties. Compound modification is also associated with changes in their sub-cellular localization: 30% of 8 (anionic) is found in mitochondria as compared to 95% of compound 10 (cationic). Photocytotoxicity was shown to be highly correlated with membrane affinity, which depends on the asymmetrical and amphiphilic characters of sens, as well as with sub-cellular localization.

A Study on China's Income Inequality and the Relationship with Economic Growth

The purpose of this paper is to study China’s income inequality under rapid economic growth.Does the relationship between economic growth and income inequality in China follow theKuznets hypothesis? What is the main cause and trend of China’s income inequality? We usedata which covers the period 1980-2005 to analyze the overall inequality, and data coveringthe period 1980-2002 to analyze the inequality inside rural and urban areas. The derivedresults doubt the validity of Kuznets hypothesis on explaining the relationship betweeneconomic growth and income inequality in China. Also we derive the trend of China’sincreased income inequality and find that the urban-rural income disparity is the main causeof China’s income inequality.