Empirical likelihood based longitudinal studies

In longitudinal data analysis, our primary interest is in the regression parameters for the marginal expectations of the longitudinal responses; the longitudinal correlation parameters are of secondary interest. The joint likelihood function for longitudinal data is challenging, particularly for correlated discrete outcome data. Marginal modeling approaches such as generalized estimating equations (GEEs) have received much attention in the context of longitudinal regression. These methods are based on the estimates of the first two moments of the data and the working correlation structure. The confidence regions and hypothesis tests are based on the asymptotic normality. The methods are sensitive to misspecification of the variance function and the working correlation structure. Because of such misspecifications, the estimates can be inefficient and inconsistent, and inference may give incorrect results. To overcome this problem, we propose an empirical likelihood (EL) procedure based on a set of estimating equations for the parameter of interest and discuss its characteristics and asymptotic properties. We also provide an algorithm based on EL principles for the estimation of the regression parameters and the construction of a confidence region for the parameter of interest. We extend our approach to variable selection for highdimensional longitudinal data with many covariates. In this situation it is necessary to identify a submodel that adequately represents the data. Including redundant variables may impact the model’s accuracy and efficiency for inference. We propose a penalized empirical likelihood (PEL) variable selection based on GEEs; the variable selection and the estimation of the coefficients are carried out simultaneously. We discuss its characteristics and asymptotic properties, and present an algorithm for optimizing PEL. Simulation studies show that when the model assumptions are correct, our method performs as well as existing methods, and when the model is misspecified, it has clear advantages. We have applied the method to two case examples.

Territorial strategies and regions

The notion of territorial strategy emerged in the 1990s and has become more and more popular since. It refers to that combination of factors purposely assembled by governments, private and public companies, universities, and industrial associations to exploit a specific geographic competitive advantage in order to boost economic growth through the development of entrepreneurial activity and innovation. Three factors are generally considered to be the building blocks of a territorial strategy: natural resources, human capital, and industrial capabilities. Natural resources derive from environ­ mental conditions and represent raw materials or land available in a region. The presence of natural resources characterizes the typology of an industry (related to tourism, oil, wood, fish, and so forth) that exists or could exist in a certain area. Human capital refers to the stock of competences available in a certain region resulting from education and work experience. Industrial capabilities relate to complex constructs of specialized expertise, the confidence to apply knowledge and skills in various contexts and under changing conditions, and an ability repeatedly to improve methods and processes in a specific industry.

Integrating sorghum whole genome sequence information with a compendium of sorghum QTL studies reveals uneven distribution of QTL and of gene-rich regions with significant implications for crop improvement.

A comprehensive analysis was conducted using 48 sorghum QTL studies published from 1995 to 2010 to make information from historical sorghum QTL experiments available in a form that could be more readily used by sorghum researchers and plant breeders. In total, 771 QTL relating to 161 unique traits from 44 studies were projected onto a sorghum consensus map. Confidence intervals (CI) of QTL were estimated so that valid comparisons could be made between studies. The method accounted for the number of lines used and the phenotypic variation explained by individual QTL from each study. In addition, estimated centimorgan (cM) locations were calculated for the predicted sorghum gene models identified in Phytozome (JGI GeneModels SBI v1.4) and compared with QTL distribution genome-wide, both on genetic linkage (cM) and physical (base-pair/bp) map scales. QTL and genes were distributed unevenly across the genome. Heterochromatic enrichment for QTL was observed, with approximately 22% of QTL either entirely or partially located in the heterochromatic regions. Heterochromatic gene enrichment was also observed based on their predicted cM locations on the sorghum consensus map, due to suppressed recombination in heterochromatic regions, in contrast to the euchromatic gene enrichment observed on the physical, sequence-based map. The finding of high gene density in recombination-poor regions, coupled with the association with increased QTL density, has implications for the development of more efficient breeding systems in sorghum to better exploit heterosis. The projected QTL information described, combined with the physical locations of sorghum sequence-based markers and predicted gene models, provides sorghum researchers with a useful resource for more detailed analysis of traits and development of efficient marker-assisted breeding strategies.

Patterns of nucleotide substitution in angiosperm cpDNA trnL (UAA)-trnF (GAA) regions

Patterns of substitution in chloroplast encoded trnL_F regions were compared between species of Actaea (Ranunculales), Digitalis (Scrophulariales), Drosera (Caryophyllales), Panicoideae (Poales), the small chromosome species clade of Pelargonium (Geraniales), each representing a different order of flowering plants, and Huperzia (Lycopodiales). In total, the study included 265 taxa, each with > 900-bp sequences, totaling 0.24 Mb. Both pairwise and phylogeny-based comparisons were used to assess nucleotide substitution patterns. In all six groups, we found that transition/transversion ratios, as estimated by maximum likelihood on most-parsimonious trees, ranged between 0.8 and 1.0 for ingroups. These values occurred both at low sequence divergences, where substitutional saturation, i.e., multiple substitutions having occurred at the same (homologous) nucleotide position, was not expected, and at higher levels of divergence. This suggests that the angiosperm trnL-F regions evolve in a pattern different from that generally observed for nuclear and animal mtDNA (transitional/transversion ratio > or = 2). Transition/transversion ratios in the intron and the spacer region differed in all alignments compared, yet base compositions between the regions were highly similar in all six groups. A>-C transversions were significantly less frequent than the other four substitution types. This correlates with results from studies on fidelity mechanisms in DNA replication that predict A<->T and G<->C transversions to be least likely to occur. It therefore strengthens confidence in the link between mutation bias at the polymerase level and the actual fixation of substitutions as recorded on evolutionary trees, and concomitantly, in the neutrality of nucleotide substitutions as phylogenetic markers.

Mortality rate in children caused by traffic accidents according to geographical regions : Brazil, 1997 - 2005

Lesões fatais em crianças causadas por acidentes de trânsito representam um problema em muitos países. Este estudo analisou a taxa de mortalidade em crianças passageiras de automóveis menores de 10 anos de idade no Brasil, entre 1997 e 2005. Para isso, o número de mortes foi obtido diretamente no banco de dados do Sistema de Informação sobre Mortalidade (SIM) e os dados da população são projeções intercensitárias a partir censo demográfico do Instituto Brasileiro de Geografia e Estatística (IBGE) disponíveis pelo site do DATASUS. Foram calculadas, para os triênios compreendidos no período em estudo, as taxas de mortalidade por acidente de trânsito entre crianças passageiras de automóveis segundo faixa etária (menor que 1 ano, 1 a 4 e 5 a 9) e região geográfica. Os resultados mostraram taxas de mortalidade de 5,68, 7,32 e 6,78 (por 1.000.000), respectivamente, para os períodos 1997-1999, 2000-2002 e 2003-2005 para todo o Brasil. Crianças menores de 1 ano de idade apresentam taxa de mortalidade de 10,18 (por 1,000,000), maior que as observadas para as outras faixas etárias. Para o período 1997-2005, as maiores taxas foram observadas nas regiões Centro-Oeste e Sul, representando, respectivamente, 13,88 e 11,47 (por 1.000.000). Tais resultados mostram a situação de risco da criança em relação a acidentes de trânsito como passageiras de automóveis e contribuem para a elaboração de campanhas educativas de prevenção de lesões.

POOR PERFORMANCE OF BOOTSTRAP CONFIDENCE INTERVALS FOR THE LOCATION OF A QUANTITATIVE TRAIT LOUCS

The aim of many genetic studies is to locate the genomic regions (called quantitative trait loci, QTLs) that contribute to variation in a quantitative trait (such as body weight). Confidence intervals for the locations of QTLs are particularly important for the design of further experiments to identify the gene or genes responsible for the effect. Likelihood support intervals are the most widely used method to obtain confidence intervals for QTL location, but the non-parametric bootstrap has also been recommended. Through extensive computer simulation, we show that bootstrap confidence intervals are poorly behaved and so should not be used in this context. The profile likelihood (or LOD curve) for QTL location has a tendency to peak at genetic markers, and so the distribution of the maximum likelihood estimate (MLE) of QTL location has the unusual feature of point masses at genetic markers; this contributes to the poor behavior of the bootstrap. Likelihood support intervals and approximate Bayes credible intervals, on the other hand, are shown to behave appropriately.

Descriptive epidemiology of cancer in three regions of Kazakhstan

The current study investigated data quality and estimated cancer incidence and mortality rates using data provided by Pavlodar, Semipalatinsk and Ust-Kamenogorsk Regional Cancer Registries of Kazakhstan during the period of 1996–1998. Assessment of data quality was performed using standard quality indicators including internal database checks, proportion of cases verified from death certificates only, mortality:incidence ratio, data patterns, proportion of cases with unknown primary site, proportion of cases with unknown age. Crude and age-adjusted incidence and mortality rates and 95% confidence intervals were calculated, by gender, for all cancers combined and for 28 specific cancer sites for each year of the study period. The five most frequent cancers were identified and described for every population. The results of the study provide the first simultaneous assessment of data quality and standardized incidence and mortality rates for Kazakh cancer registries. ^

Statistical prediction of single-stranded regions in RNA secondary structure and application to predicting effective antisense target sites and beyond

Single-stranded regions in RNA secondary structure are important for RNA–RNA and RNA–protein interactions. We present a probability profile approach for the prediction of these regions based on a statistical algorithm for sampling RNA secondary structures. For the prediction of phylogenetically-determined single-stranded regions in secondary structures of representative RNA sequences, the probability profile offers substantial improvement over the minimum free energy structure. In designing antisense oligonucleotides, a practical problem is how to select a secondary structure for the target mRNA from the optimal structure(s) and many suboptimal structures with similar free energies. By summarizing the information from a statistical sample of probable secondary structures in a single plot, the probability profile not only presents a solution to this dilemma, but also reveals ‘well-determined’ single-stranded regions through the assignment of probabilities as measures of confidence in predictions. In antisense application to the rabbit β-globin mRNA, a significant correlation between hybridization potential predicted by the probability profile and the degree of inhibition of in vitro translation suggests that the probability profile approach is valuable for the identification of effective antisense target sites. Coupling computational design with DNA–RNA array technique provides a rational, efficient framework for antisense oligonucleotide screening. This framework has the potential for high-throughput applications to functional genomics and drug target validation.

Seed dispersal characteristics to polar regions by means of visitors

Invasive alien species are among the primary causes of biodiversity change globally, with the risks thereof broadly understood for most regions of the world. They are similarly thought to be among the most significant conservation threats to Antarctica, especially as climate change proceeds in the region. However, no comprehensive, continent-wide evaluation of the risks to Antarctica posed by such species has been undertaken. Here we do so by sampling, identifying, and mapping the vascular plant propagules carried by all categories of visitors to Antarctica during the International Polar Year's first season (2007-2008) and assessing propagule establishment likelihood based on their identity and origins and on spatial variation in Antarctica's climate. For an evaluation of the situation in 2100, we use modeled climates based on the Intergovernmental Panel on Climate Change's Special Report on Emissions Scenarios Scenario A1B [Nakicenovic N, Swart R, eds (2000) Special Report on Emissions Scenarios: A Special Report of Working Group III of the Intergovernmental Panel on Climate Change (Cambridge University Press, Cambridge, UK)]. Visitors carrying seeds average 9.5 seeds per person, although as vectors, scientists carry greater propagule loads than tourists. Annual tourist numbers (~33,054) are higher than those of scientists (~7,085), thus tempering these differences in propagule load. Alien species establishment is currently most likely for the Western Antarctic Peninsula. Recent founder populations of several alien species in this area corroborate these findings. With climate change, risks will grow in the Antarctic Peninsula, Ross Sea, and East Antarctic coastal regions. Our evidence-based assessment demonstrates which parts of Antarctica are at growing risk from alien species that may become invasive and provides the means to mitigate this threat now and into the future as the continent's climate changes.

Flow field analysis of some mixing and conveying screw element regions, within a closely intermeshing, co-rotating twin-screw extruder.

Grafting of antioxidants and other modifiers onto polymers by reactive extrusion, has been performed successfully by the Polymer Processing and Performance Group at Aston University. Traditionally the optimum conditions for the grafting process have been established within a Brabender internal mixer. Transfer of this batch process to a continuous processor, such as an extruder, has, typically, been empirical. To have more confidence in the success of direct transfer of the process requires knowledge of, and comparison between, residence times, mixing intensities, shear rates and flow regimes in the internal mixer and in the continuous processor.The continuous processor chosen for the current work in the closely intermeshing, co-rotating twin-screw extruder (CICo-TSE). CICo-TSEs contain screw elements that convey material with a self-wiping action and are widely used for polymer compounding and blending. Of the different mixing modules contained within the CICo-TSE, the trilobal elements, which impose intensive mixing, and the mixing discs, which impose extensive mixing, are of importance when establishing the intensity of mixing. In this thesis, the flow patterns within the various regions of the single-flighted conveying screw elements and within both the trilobal element and mixing disc zones of a Betol BTS40 CICo-TSE, have been modelled using the computational fluid dynamics package Polyflow. A major obstacle encountered when solving the flow problem within all of these sets of elements, arises from both the complex geometry and the time-dependent flow boundaries as the elements rotate about their fixed axes. Simulation of the time dependent boundaries was overcome by selecting a number of sequential 2D and 3D geometries, used to represent partial mixing cycles. The flow fields were simulated using the ideal rheological properties of polypropylene and characterised in terms of velocity vectors, shear stresses generated and a parameter known as the mixing efficiency. The majority of the large 3D simulations were performed on the Cray J90 supercomputer situated at the Rutherford-Appleton laboratories, with pre- and postprocessing operations achieved via a Silicon Graphics Indy workstation. A mechanical model was constructed consisting of various CICo-TSE elements rotating within a transparent outer barrel. A technique has been developed using coloured viscous clays whereby the flow patterns and mixing characteristics within the CICo-TSE may be visualised. In order to test and verify the simulated predictions, the patterns observed within the mechanical model were compared with the flow patterns predicted by the computational model. The flow patterns within the single-flighted conveying screw elements in particular, showed good agreement between the experimental and simulated results.

Defense Expenditure: Andean and Amazon Regions

The Andean and the amazon, comprised of Venezuela, Colombia, Bolivia, Peru and Ecuador, have recently undertaken significant modernization efforts ranging from equipment, logistics, doctrine, training, deployment and the re-definition of the roles and missions of their forces. In most cases, motivations to modernize have been internal, such as continuing operations against armed groups as in the case of Colombia and Peru, enhance border control and sovereignty enforcement, as in the case of Ecuador and Brazil or regime control in Venezuela. However, they are complemented by perceptions of external threats, including traditional intra-state conventional wars. The increased tensions between Colombia and Venezuela and Ecuador as well as the historic Peru-Chile tensions are the most salient examples. Although diplomacy –especially defence diplomacy- has worked to a good degree in creating and strengthening confidence building measures, the potential for inter-state conflict is higher in this region of the Americas. This region has seen the recent emergence of long-term modernization plans, initially in Colombia followed by Venezuela and Ecuador and probably best embodied in scope and scale by the Brazilian National Defence Plan (for its long term vision). Although it has been speculated that high allocation of funds to retirement pension systems has had an impact on delaying modernization plans, this comparative study on the allocation of pension and social funds in these particular countries concludes that there is no direct linkage between the poor funding of military modernization plans and the diversion of funds to military pension systems.

Effectiveness of Depression-Suicidal Behaviour Gatekeeper training among police officers in three European regions: outcomes of the Optimising Suicide Prevention Programmes and their Implementation in Europe (OSPI-Europe) study

Background: Gatekeeper training for community facilitators, to identify and respond to those at risk of suicide, forms an important part of multi-level community-based suicide prevention programmes. Aims: This study examined the effects of gatekeeper training on attitudes, knowledge and confidence of police officers in dealing with persons at risk of suicide. Methods: A total of 828 police officers across three European regions participated in a 4-hour training programme which addressed the epidemiology of depression and suicidal behaviour, symptoms of depression, warning signs and risk factors associated with suicidal behaviour, motivating help-seeking behaviour, dealing with acute suicidal crisis and informing bereaved relatives. Participants completed internationally validated questionnaires assessing stigmatising attitudes, knowledge about depression and confidence in dealing with suicidal persons pre- and post-training. Results: There were significant differences among countries in terms of previous exposure to suicidal persons and extent of previous training. Post-training evaluation demonstrated significant improvements in stigmatising attitudes, knowledge and confidence in all three countries. Conclusion: The consistently positive effects of gatekeeper training of police officers across different regions support inclusion of this type of training as a fundamental part of multi-level community-based suicide prevention programmes and roll-out, nationally and internationally.