Spatial modelling and prediction on river networks: up model, down model or hybrid?

Preservation of rivers and water resources is crucial in most environmental policies and many efforts are made to assess water quality. Environmental monitoring of large river networks are based on measurement stations. Compared to the total length of river networks, their number is often limited and there is a need to extend environmental variables that are measured locally to the whole river network. The objective of this paper is to propose several relevant geostatistical models for river modeling. These models use river distance and are based on two contrasting assumptions about dependency along a river network. Inference using maximum likelihood, model selection criterion and prediction by kriging are then developed. We illustrate our approach on two variables that differ by their distributional and spatial characteristics: summer water temperature and nitrate concentration. The data come from 141 to 187 monitoring stations in a network on a large river located in the Northeast of France that is more than 5000 km long and includes Meuse and Moselle basins. We first evaluated different spatial models and then gave prediction maps and error variance maps for the whole stream network.

Spatial statistical models that use flow and stream distance

We develop spatial statistical models for stream networks that can estimate relationships between a response variable and other covariates, make predictions at unsampled locations, and predict an average or total for a stream or a stream segment. There have been very few attempts to develop valid spatial covariance models that incorporate flow, stream distance, or both. The application of typical spatial autocovariance functions based on Euclidean distance, such as the spherical covariance model, are not valid when using stream distance. In this paper we develop a large class of valid models that incorporate flow and stream distance by using spatial moving averages. These methods integrate a moving average function, or kernel, against a white noise process. By running the moving average function upstream from a location, we develop models that use flow, and by construction they are valid models based on stream distance. We show that with proper weighting, many of the usual spatial models based on Euclidean distance have a counterpart for stream networks. Using sulfate concentrations from an example data set, the Maryland Biological Stream Survey (MBSS), we show that models using flow may be more appropriate than models that only use stream distance. For the MBSS data set, we use restricted maximum likelihood to fit a valid covariance matrix that uses flow and stream distance, and then we use this covariance matrix to estimate fixed effects and make kriging and block kriging predictions.

Applications of Linear Programming to Coding Theory

Maximum-likelihood decoding is often the optimal decoding rule one can use, but it is very costly to implement in a general setting. Much effort has therefore been dedicated to find efficient decoding algorithms that either achieve or approximate the error-correcting performance of the maximum-likelihood decoder. This dissertation examines two approaches to this problem. In 2003 Feldman and his collaborators defined the linear programming decoder, which operates by solving a linear programming relaxation of the maximum-likelihood decoding problem. As with many modern decoding algorithms, is possible for the linear programming decoder to output vectors that do not correspond to codewords; such vectors are known as pseudocodewords. In this work, we completely classify the set of linear programming pseudocodewords for the family of cycle codes. For the case of the binary symmetric channel, another approximation of maximum-likelihood decoding was introduced by Omura in 1972. This decoder employs an iterative algorithm whose behavior closely mimics that of the simplex algorithm. We generalize Omura's decoder to operate on any binary-input memoryless channel, thus obtaining a soft-decision decoding algorithm. Further, we prove that the probability of the generalized algorithm returning the maximum-likelihood codeword approaches 1 as the number of iterations goes to infinity.

Análise multivariada de características reprodutivas e produtivas de vacas da raça Holandesa

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Schizophrenia, brain disease and meta-analyses: Integrating the pieces and testing Fusar-Poli`s hypothesis

This paper aims to discuss and test the hypothesis raised by Fusar-Poli [Fusar-Poli P. Can neuroimaging prove that schizophrenia is a brain disease? A radical hypothesis. Medical Hypotheses in press, corrected proof] that ""on the basis of the available imaging literature there is no consistent evidence to reject the radical and provocative hypothesis that schizophrenia is not a brain disease"". To achieve this goal, all meta-analyses on `fMRI and schizophrenia` published during the current decade and indexed in Pubmed were summarized, as much as some other useful information, e.g., meta-analyses on genetic risk factors. Our main conclusion is that the literature fully supports the hypothesis that schizophrenia is a syndrome (not a disease) associated with brain abnormalities, despite the fact that there is no singular and reductionist pathway from the nosographic entity (schizophrenia) to its causes. This irreducibility is due to the fact that the syndrome has more than one dimension (e.g., cognitive, psychotic and negative) and each of them is related to abnormalities in specific neuronal networks. A psychiatric diagnosis is a statistical procedure; these dimensions are not identically represented in each diagnosticated case and this explains the existence of more than one pattern of brain abnormalities related to schizophrenia. For example, chronification is associated with negativism while the first psychotic episode is not; in that sense, the same person living with schizophrenia may reveal different symptoms and fMRI patterns along the course of his life, and this is precisely what defines schizophrenia since the time when it was called Dementia Praecox (first by pick then by Kraepelin). It is notable that 100% of the collected meta-analyses on `fMRI and schizophrenia` reveal positive findings. Moreover, all meta-analyses that found positive associations between schizophrenia and genetic risk factors have to do with genes (SNPs) especially activated in neuronal tissue of the central nervous system (CNS), suggesting that, to the extent these polymorphisms are related to schizophrenia`s etiology, they are also related to abnormal brain activity. (C) 2009 Elsevier Ltd. All rights reserved.

The Kumaraswamy Gumbel distribution

The Gumbel distribution is perhaps the most widely applied statistical distribution for problems in engineering. We propose a generalization-referred to as the Kumaraswamy Gumbel distribution-and provide a comprehensive treatment of its structural properties. We obtain the analytical shapes of the density and hazard rate functions. We calculate explicit expressions for the moments and generating function. The variation of the skewness and kurtosis measures is examined and the asymptotic distribution of the extreme values is investigated. Explicit expressions are also derived for the moments of order statistics. The methods of maximum likelihood and parametric bootstrap and a Bayesian procedure are proposed for estimating the model parameters. We obtain the expected information matrix. An application of the new model to a real dataset illustrates the potentiality of the proposed model. Two bivariate generalizations of the model are proposed.

Evolution of extrafloral nectaries: adaptive process and selective regime changes from forest to savanna

Much effort has been devoted to understanding the function of extrafloral nectaries (EFNs) for antplantherbivore interactions. However, the pattern of evolution of such structures throughout the history of plant lineages remains unexplored. In this study, we used empirical knowledge on plant defences mediated by ants as a theoretical framework to test specific hypotheses about the adaptive role of EFNs during plant evolution. Emphasis was given to different processes (neutral or adaptive) and factors (habitat change and trade-offs with new trichomes) that may have affected the evolution of antplant associations. We measured seven EFN quantitative traits in all 105 species included in a well-supported phylogeny of the tribe Bignonieae (Bignoniaceae) and collected field data on antEFN interactions in 32 species. We identified a positive association between ant visitation (a surrogate of ant guarding) and the abundance of EFNs in vegetative plant parts and rejected the hypothesis of phylogenetic conservatism of EFNs, with most traits presenting K-values < 1. Modelling the evolution of EFN traits using maximum likelihood approaches further suggested adaptive evolution, with static-optimum models showing a better fit than purely drift models. In addition, the abundance of EFNs was associated with habitat shifts (with a decrease in the abundance of EFNs from forest to savannas), and a potential trade-off was detected between the abundance of EFNs and estipitate glandular trichomes (i.e. trichomes with sticky secretion). These evolutionary associations suggest divergent selection between species as well as explains K-values < 1. Experimental studies with multiple lineages of forest and savanna taxa may improve our understanding of the role of nectaries in plants. Overall, our results suggest that the evolution of EFNs was likely associated with the adaptive process which probably played an important role in the diversification of this plant group.

Heritability Estimate and Genetic Correlations of Reproductive Features in Nellore Bulls, Offspring of Super Precocious, Precocious and Normal Cows Under Extensive Farming Conditions

The present work aimed to estimate heritability and genetic correlations of reproductive features of Nellore bulls, offspring of mothers classified as superprecocious (M1), precocious (M2) and normal (M3). Twenty one thousand hundred and eighty-six animals with average age of 21.29 months were used, evaluated through the breeding soundness evaluation from 1999 to 2008. The breeding soundness features included physical semen evaluation (progressive sperm motility and sperm vigour), semen morphology (major, minor and total sperm defects), scrotal circumference (SC), testicular volume (TV) and SC at 18 months of age (SC18). The components of variance, heritability and genetic correlations for and between the features were estimated simultaneously by restricted maximum likelihood, with the use of the vce software system vs 6. The heritability estimates were high for SC18, SC and TV (0.43, 0.63 and 0.54; 0.45, 0.45 and 0.44; 0.42, 0.45 and 0.41, respectively for the categories of mothers M1, M2 and M3) and low for physical and morphological semen aspects. The genetic correlations between SC18 and SC were high, as well as between these variables with TV. High and positive genetic correlations were recorded among SC18, SC and TV with the physical aspects of the semen, although no favourable association was verified with the morphological aspects, for the three categories of mothers. It can be concluded that the mothers sexual precocity did not affect the heritability of their offspring reproduction features.

A novel linkage map of sugarcane with evidence for clustering of retrotransposon-based markers

Background: The development of sugarcane as a sustainable crop has unlimited applications. The crop is one of the most economically viable for renewable energy production, and CO2 balance. Linkage maps are valuable tools for understanding genetic and genomic organization, particularly in sugarcane due to its complex polyploid genome of multispecific origins. The overall objective of our study was to construct a novel sugarcane linkage map, compiling AFLP and EST-SSR markers, and to generate data on the distribution of markers anchored to sequences of scIvana_1, a complete sugarcane transposable element, and member of the Copia superfamily. Results: The mapping population parents ('IAC66-6' and 'TUC71-7') contributed equally to polymorphisms, independent of marker type, and generated markers that were distributed into nearly the same number of co-segregation groups (or CGs). Bi-parentally inherited alleles provided the integration of 19 CGs. The marker number per CG ranged from two to 39. The total map length was 4,843.19 cM, with a marker density of 8.87 cM. Markers were assembled into 92 CGs that ranged in length from 1.14 to 404.72 cM, with an estimated average length of 52.64 cM. The greatest distance between two adjacent markers was 48.25 cM. The scIvana_1-based markers (56) were positioned on 21 CGs, but were not regularly distributed. Interestingly, the distance between adjacent scIvana_1-based markers was less than 5 cM, and was observed on five CGs, suggesting a clustered organization. Conclusions: Results indicated the use of a NBS-profiling technique was efficient to develop retrotransposon-based markers in sugarcane. The simultaneous maximum-likelihood estimates of linkage and linkage phase based strategies confirmed the suitability of its approach to estimate linkage, and construct the linkage map. Interestingly, using our genetic data it was possible to calculate the number of retrotransposonscIvana_1 (similar to 60) copies in the sugarcane genome, confirming previously reported molecular results. In addition, this research possibly will have indirect implications in crop economics e. g., productivity enhancement via QTL studies, as the mapping population parents differ in response to an important fungal disease.

Computational tools for comparing asymmetric GARCH models via Bayes factors

In this paper we use Markov chain Monte Carlo (MCMC) methods in order to estimate and compare GARCH models from a Bayesian perspective. We allow for possibly heavy tailed and asymmetric distributions in the error term. We use a general method proposed in the literature to introduce skewness into a continuous unimodal and symmetric distribution. For each model we compute an approximation to the marginal likelihood, based on the MCMC output. From these approximations we compute Bayes factors and posterior model probabilities. (C) 2012 IMACS. Published by Elsevier B.V. All rights reserved.

THE LOG-BURR XII REGRESSION MODEL FOR GROUPED SURVIVAL DATA

The log-Burr XII regression model for grouped survival data is evaluated in the presence of many ties. The methodology for grouped survival data is based on life tables, where the times are grouped in k intervals, and we fit discrete lifetime regression models to the data. The model parameters are estimated by maximum likelihood and jackknife methods. To detect influential observations in the proposed model, diagnostic measures based on case deletion, so-called global influence, and influence measures based on small perturbations in the data or in the model, referred to as local influence, are used. In addition to these measures, the total local influence and influential estimates are also used. We conduct Monte Carlo simulation studies to assess the finite sample behavior of the maximum likelihood estimators of the proposed model for grouped survival. A real data set is analyzed using a regression model for grouped data.

First evidence of the emergence of novel putative infectious bronchitis virus genotypes in Cuba

The emergence of new infectious bronchitis virus (IBV) genotypes or serotypes along with the poor cross-protection observed among IBV serotypes have complicated the avian infectious bronchitis (IB) control programs in different geographic regions. In Cuba, the lack of genetic information regarding IBV and the increasing epidemiological importance of this virus in Cuban chicken flocks demand further characterization of IBV isolates. In the present work, studies of genetic diversity and phylogenetic relationships among recent IBV isolates from Cuban chicken flocks showing respiratory disorders were performed. Two putative genotypes genetically different to the Massachusetts genotype H120 strain used in the Cuban vaccination program were found in the flocks assessed. In addition, a potential nephropathogenic IBV isolate was found by first time in Cuba. (C) 2012 Elsevier Ltd. All rights reserved.

The 3D structure and function of digestive cathepsin L-like proteinases of Tenebrio molitor larval midgut

Cathepsin L-like proteinases (CAL) are major digestive proteinases in the beetle Tenebrio molitor. Procathepsin Ls 2 (pCAL2) and 3 (pCAL3) were expressed as recombinant proteins in Escherichia coil, purified and activated under acidic conditions. Immunoblot analyses of different T. molitor larval tissues demonstrated that a polyclonal antibody to pCAL3 recognized pCAL3 and cathepsin L 3 (CAD) only in the anterior two-thirds of midgut tissue and midgut luminal contents of T. molitor larvae. Furthermore, immunocytolocalization data indicated that pCAL3 occurs in secretory vesicles and microvilli in anterior midgut Therefore CAL3, like cathepsin L 2 (CAL2), is a digestive enzyme secreted by T. molitor anterior midgut CAD hydrolyses Z-FR-MCA and Z-RR-MCA (typical cathepsin substrates), whereas CAL2 hydrolyses only Z-FR-MCA. Active site mutants (pCAL2C25S and pCAL3C265) were constructed by replacing the catalytic cysteine with serine to prevent autocatalytic processing. Recombinant pCAL2 and pCAL3 mutants (pCAL2C25S and pCAL3C26S) were prepared, crystallized and their 3D structures determined at 1.85 and 2.1 angstrom, respectively. While the overall structure of these enzymes is similar to other members of the papain superfamily, structural differences in the S2 subsite explain their substrate specificities. The data also supported models for CAL trafficking to lysosomes and to secretory vesicles to be discharged into midgut contents. (C) 2012 Elsevier Ltd. All rights reserved.

Comparative Genome Analysis of Trichophyton rubrum and Related Dermatophytes Reveals Candidate Genes Involved in Infection

The major cause of athlete's foot is Trichophyton rubrum, a dermatophyte or fungal pathogen of human skin. To facilitate molecular analyses of the dermatophytes, we sequenced T. rubrum and four related species, Trichophyton tonsurans, Trichophyton equinum, Microsporum canis, and Microsporum gypseum. These species differ in host range, mating, and disease progression. The dermatophyte genomes are highly colinear yet contain gene family expansions not found in other human-associated fungi. Dermatophyte genomes are enriched for gene families containing the LysM domain, which binds chitin and potentially related carbohydrates. These LysM domains differ in sequence from those in other species in regions of the peptide that could affect substrate binding. The dermatophytes also encode novel sets of fungus-specific kinases with unknown specificity, including nonfunctional pseudokinases, which may inhibit phosphorylation by competing for kinase sites within substrates, acting as allosteric effectors, or acting as scaffolds for signaling. The dermatophytes are also enriched for a large number of enzymes that synthesize secondary metabolites, including dermatophyte-specific genes that could synthesize novel compounds. Finally, dermatophytes are enriched in several classes of proteases that are necessary for fungal growth and nutrient acquisition on keratinized tissues. Despite differences in mating ability, genes involved in mating and meiosis are conserved across species, suggesting the possibility of cryptic mating in species where it has not been previously detected. These genome analyses identify gene families that are important to our understanding of how dermatophytes cause chronic infections, how they interact with epithelial cells, and how they respond to the host immune response. IMPORTANCE Athlete's foot, jock itch, ringworm, and nail infections are common fungal infections, all caused by fungi known as dermatophytes (fungi that infect skin). This report presents the genome sequences of Trichophyton rubrum, the most frequent cause of athlete's foot, as well as four other common dermatophytes. Dermatophyte genomes are enriched for four gene classes that may contribute to the ability of these fungi to cause disease. These include (i) proteases secreted to degrade skin; (ii) kinases, including pseudokinases, that are involved in signaling necessary for adapting to skin; (iii) secondary metabolites, compounds that act as toxins or signals in the interactions between fungus and host; and (iv) a class of proteins (LysM) that appear to bind and mask cell wall components and carbohydrates, thus avoiding the host's immune response to the fungi. These genome sequences provide a strong foundation for future work in understanding how dermatophytes cause disease.

Genetic correlations and heritability estimates of reproductive traits of Nellore bulls

Estimates of phenotypic, genetics and residual variances for reproductive traits in 5903 Nellore bulls were obtained. The experimental model used was multiple trait derivative-free restricted maximum likelihood. The values obtained for heritability were 0.24 +/- 0.05 for scrotal circumference at 450 days of age and 0.37 +/- 0.05 at 21 months for age at the time of the breeding soundness evaluation; 0.24 +/- 0.05 and 0.26 +/- 0.05 for left and right testicle length; 0.29 +/- 0.05 and 0.31 +/- 0.05 for left and right testicle width; 0.12 +/- 0.04 for testicle format; 0.33 +/- 0.06 for testicle volume; 0.11 +/- 0.03 for gross motility; 0.08 +/- 0.03 for individual motility and 0.05 +/- 0.02 for spermatic vigor; 0.20 +/- 0.04, 0.03 +/- 0.02 and 0.19 +/- 0.04 for larger defects, smaller defects and total defects, respectively. The values for heritability for testicular biometric characteristics were moderate to high while the seminal characteristics, presented low values. Genetic correlations between scrotal circumference with all the reproductive traits were favorable, suggesting the scrotal circumference as a feature of choice in the selection of bulls.