Variable selection in monotone single‐index models via the adaptive LASSO

We consider the problem of variable selection for monotone single‐index models. A single‐index model assumes that the expectation of the outcome is an unknown function of a linear combination of covariates. Assuming monotonicity of the unknown function is often reasonable and allows for more straightforward inference. We present an adaptive LASSO penalized least squares approach to estimating the index parameter and the unknown function in these models for continuous outcome. Monotone function estimates are achieved using the pooled adjacent violators algorithm, followed by kernel regression. In the iterative estimation process, a linear approximation to the unknown function is used, therefore reducing the situation to that of linear regression and allowing for the use of standard LASSO algorithms, such as coordinate descent. Results of a simulation study indicate that the proposed methods perform well under a variety of circumstances and that an assumption of monotonicity, when appropriate, noticeably improves performance. The proposed methods are applied to data from a randomized clinical trial for the treatment of a critical illness in the intensive care unit. Copyright © 2013 John Wiley & Sons, Ltd.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/100172/1/sim5834.pd

Assessment of DNA extracted from FTA® cards for use on the Illumina iSelect BeadChip

<p>Abstract</p> <p>Background</p> <p>As FTA^®cards provide an ideal medium for the field collection of DNA we sought to assess the quality of genomic DNA extracted from this source for use on the Illumina BovineSNP50 iSelect BeadChip which requires unbound, relatively intact (fragment sizes ≥ 2 kb), and high-quality DNA. Bovine blood and nasal swab samples collected on FTA cards were extracted using the commercially available GenSolve kit with a minor modification. The call rate and concordance of genotypes from each sample were compared to those obtained from whole blood samples extracted by standard PCI extraction.</p> <p>Findings</p> <p>An ANOVA analysis indicated no significant difference (P > 0.72) in BovineSNP50 genotype call rate between DNA extracted from FTA cards by the GenSolve kit or extracted from whole blood by PCI. Two sample t-tests demonstrated that the DNA extracted from the FTA cards produced genotype call and concordance rates that were not different to those produced by assaying DNA samples extracted by PCI from whole blood.</p> <p>Conclusion</p> <p>We conclude that DNA extracted from FTA cards by the GenSolve kit is of sufficiently high quality to produce results comparable to those obtained from DNA extracted from whole blood when assayed by the Illumina iSelect technology. Additionally, we validate the use of nasal swabs as an alternative to venous blood or buccal samples from animal subjects for reliably producing high quality genotypes on this platform.</p

Analysis of time-to-event for observational studies: Guidance to the use of intensity models

This paper provides guidance for researchers with some mathematical background on the conduct of time-to-event analysis in observational studies based on intensity (hazard) models. Discussions of basic concepts like time axis, event definition and censoring are given. Hazard models are introduced, with special emphasis on the Cox proportional hazards regression model. We provide check lists that may be useful both when fitting the model and assessing its goodness of fit and when interpreting the results. Special attention is paid to how to avoid problems with immortal time bias by introducing time-dependent covariates. We discuss prediction based on hazard models and difficulties when attempting to draw proper causal conclusions from such models. Finally, we present a series of examples where the methods and check lists are exemplified. Computational details and implementation using the freely available R software are documented in Supplementary Material. The paper was prepared as part of the STRATOS initiative.Comment: 28 pages, 12 figures. For associated Supplementary material, see http://publicifsv.sund.ku.dk/~pka/STRATOSTG8

On the survivability and detectability of terrestrial meteorites on the moon

Materials blasted into space from the surface of early Earth may preserve a unique record of our planet's early surface environment. Armstrong et al. (2002) pointed out that such materials, in the form of terrestrial meteorites, may exist on the Moon and be of considerable astrobiological interest if biomarkers from early Earth are preserved within them. Here, we report results obtained via the AUTODYN hydrocode to calculate the peak pressures within terrestrial meteorites on the lunar surface to assess their likelihood of surviving the impact. Our results confirm the order-of-magnitude estimates of Armstrong et al. (2002) that substantial survivability is to be expected, especially in the case of relatively low velocity (ca. 2.5 km/s) or oblique (≤45°) impacts, or both. We outline possible mechanisms for locating such materials on the Moon and conclude that searching for them would be a scientifically valuable activity for future lunar exploration

What's in your next-generation sequence data? An exploration of unmapped DNA and RNA sequence reads from the bovine reference individual.

BackgroundNext-generation sequencing projects commonly commence by aligning reads to a reference genome assembly. While improvements in alignment algorithms and computational hardware have greatly enhanced the efficiency and accuracy of alignments, a significant percentage of reads often remain unmapped.ResultsWe generated de novo assemblies of unmapped reads from the DNA and RNA sequencing of the Bos taurus reference individual and identified the closest matching sequence to each contig by alignment to the NCBI non-redundant nucleotide database using BLAST. As expected, many of these contigs represent vertebrate sequence that is absent, incomplete, or misassembled in the UMD3.1 reference assembly. However, numerous additional contigs represent invertebrate species. Most prominent were several species of Spirurid nematodes and a blood-borne parasite, Babesia bigemina. These species are either not present in the US or are not known to infect taurine cattle and the reference animal appears to have been host to unsequenced sister species.ConclusionsWe demonstrate the importance of exploring unmapped reads to ascertain sequences that are either absent or misassembled in the reference assembly and for detecting sequences indicative of parasitic or commensal organisms

Moving forward in circles: challenges and opportunities in modelling population cycles

Population cycling is a widespread phenomenon, observed across a multitude of taxa in both laboratory and natural conditions. Historically, the theory associated with population cycles was tightly linked to pairwise consumer–resource interactions and studied via deterministic models, but current empirical and theoretical research reveals a much richer basis for ecological cycles. Stochasticity and seasonality can modulate or create cyclic behaviour in non-intuitive ways, the high-dimensionality in ecological systems can profoundly influence cycling, and so can demographic structure and eco-evolutionary dynamics. An inclusive theory for population cycles, ranging from ecosystem-level to demographic modelling, grounded in observational or experimental data, is therefore necessary to better understand observed cyclical patterns. In turn, by gaining better insight into the drivers of population cycles, we can begin to understand the causes of cycle gain and loss, how biodiversity interacts with population cycling, and how to effectively manage wildly fluctuating populations, all of which are growing domains of ecological research

Tissue Tropism in Host Transcriptional Response to Members of the Bovine Respiratory Disease Complex.

Bovine respiratory disease (BRD) is the most common infectious disease of beef and dairy cattle and is characterized by a complex infectious etiology that includes a variety of viral and bacterial pathogens. We examined the global changes in mRNA abundance in healthy lung and lung lesions and in the lymphoid tissues bronchial lymph node, retropharyngeal lymph node, nasopharyngeal lymph node and pharyngeal tonsil collected at the peak of clinical disease from beef cattle experimentally challenged with either bovine respiratory syncytial virus, infectious bovine rhinotracheitis, bovine viral diarrhea virus, Mannheimia haemolytica or Mycoplasma bovis. We identified signatures of tissue-specific transcriptional responses indicative of tropism in the coordination of host's immune tissue responses to infection by viral or bacterial infections. Furthermore, our study shows that this tissue tropism in host transcriptional response to BRD pathogens results in the activation of different networks of response genes. The differential crosstalk among genes expressed in lymphoid tissues was predicted to be orchestrated by specific immune genes that act as 'key players' within expression networks. The results of this study serve as a basis for the development of innovative therapeutic strategies and for the selection of cattle with enhanced resistance to BRD

Impacts of lower speed limits in South Australia

ustralia has recently undergone a change of urban speed limits in most of its jurisdictions. The political and social shift in attitudes required for this change is significant in a country that has a strong reliance on the use of the private motor vehicle. At present five states and one territory have lowered the speed limit in urban areas from 60km/h to 50km/h. Of the remaining two states and territories, one has already implemented lower speed limits (40km/h) in small areas within cities. In the majority of cases, the only criterion for judging the success of such schemes is an observed reduction in both speeds and crash numbers. This paper reports on a more holistic assessment of such schemes taking into account factors in addition to speed and crashes including: traffic volume displacement, physical road network characteristics, environmental factors, community ownership and acceptance, enforcement effort and impact on travel times. The research work has included the analysis of extensive traffic data, community surveys and focus groups, the collection of environmental and travel time data from an instrumented probe vehicle and the computer modelling of road networks. The work reported is based on over 10 years of working with a 40km/h lower urban speed limit area in South Australia. The paper expands the notion of using speed and crash outcomes as the only criteria for measuring the success of lower speed limit schemes. Document type: Articl