Numerical Analysis and Scientific Computing Commons^™

Automated Species Classification Methods For Passive Acoustic Monitoring Of Beaked Whales, John Lebien

LSU New Orleans Theses and Dissertations

The Littoral Acoustic Demonstration Center has collected passive acoustic monitoring data in the northern Gulf of Mexico since 2001. Recordings were made in 2007 near the Deepwater Horizon oil spill that provide a baseline for an extensive study of regional marine mammal populations in response to the disaster. Animal density estimates can be derived from detections of echolocation signals in the acoustic data. Beaked whales are of particular interest as they remain one of the least understood groups of marine mammals, and relatively few abundance estimates exist. Efficient methods for classifying detected echolocation transients are essential for mining long-term passive …

Underwater Acoustic Signal Analysis Toolkit, Kirk Bienvenu Jr

LSU New Orleans Theses and Dissertations

This project started early in the summer of 2016 when it became evident there was a need for an effective and efficient signal analysis toolkit for the Littoral Acoustic Demonstration Center Gulf Ecological Monitoring and Modeling (LADC-GEMM) Research Consortium. LADC-GEMM collected underwater acoustic data in the northern Gulf of Mexico during the summer of 2015 using Environmental Acoustic Recording Systems (EARS) buoys. Much of the visualization of data was handled through short scripts and executed through terminal commands, each time requiring the data to be loaded into memory and parameters to be fed through arguments. The vision was to develop …

News From The Bioconductor Project, Bioconductor Core Team

The R Journal

The Bioconductor project provides tools for the analysis and comprehension of high throughput genomic data. Bioconductor 3.6 was released on 31 October, 2017. It is compatible with R 3.4.3 and consists of 1473 software packages, 326 experiment data packages, and 911 up-to-date annotation packages. The release announcement includes descriptions of 100 new software packages, and updated NEWS files for many additional packages. Start using Bioconductor by installing the most recent version of R and evaluating the commands

Changes On Cran, Kurt Hornik, Uwe Ligges, Achim Zeileis

The R Journal

In the past 6 months,1244 new packages were added to the CRAN package repository. 19 packages were unarchived, 55 archived and 3 removed. The following shows the growth of the number of active packages in the CRAN package repository

R Teaching Column, Matthias Gehrke, Reed Davis, Norman Matloff, Paul Thompson, Tiffany Chen, Emily Watkins, Laurel Beckett

The R Journal

The revisit package, developed as a collaborative tool for scientists, also serves as a tool for teaching statistics, in a manner that can be highly motivating for students. Using either the included case studies or datasets/code provided by the instructor, students can explore several alternate paths of analysis, such as the effects of including/excluding certain variables, employing different types of statistical methodology and so on. The package includes features that help students follow modern statistical standards and avoid various statistical errors, such as “p-hacking” and lack of attention to outlier data.

Forwards Column, Stella Bollmann, Dianne Cook, Jasmine Dumas, John Fox, Julie Josse, Oliver Keyes, Carolin Strobl, Heather Turner, Rudolf Debelak

The R Journal

Forwards is a task force that was set up by the R Foundation in 2015 to address the under representation of women that has since widened its scope to encompass other under represented groups. The task force is organised as a core team comprising leaders from a number of sub-teams that focus on particular aspects:

An Introduction To Rocker: Docker Containers For R, Carl Boettiger, Dirk Eddelbuettel

The R Journal

We describe the Rocker project, which provides a widely-used suite of Docker images with customized R environments for particular tasks. We discuss how this suite is organized, and how these tools can increase portability, scaling, reproducibility, and convenience of R users and developers.

Openebgm: An R Implementation Of The Gamma-Poisson Shrinker Data Mining Model, Travis Canida, John Ihrie

The R Journal

We introduce the R package openEBGM, an implementation of the Gamma-Poisson Shrinker (GPS) model for identifying unexpected counts in large contingency tables using an empirical Bayes approach. The Empirical Bayes Geometric Mean (EBGM) and quantile scores are obtained from the GPS model estimates. openEBGM provides for the evaluation of counts using a number of different methods, including the model-based disproportionality scores, the relative reporting ratio (RR), and the proportional reporting ratio (PRR). Data squashing for computational efficiency and stratification for confounding variable adjustment are included. Application to adverse event detection is discussed.

Riskregression: Predicting The Risk Of An Event Using Cox Regression Models, Brice Ozenne, Anne Lyngholm Sørensen, Thomas Scheike, Christian Torp-Pedersen, Thomas Alexander Gerds

The R Journal

In the presence of competing risks a prediction of the time-dynamic absolute risk of an event can be based on cause-specific Cox regression models for the event and the competing risks (Benichou and Gail, 1990). We present computationally fast and memory optimized C++functions with an R inter face for predicting the covariate specific absolute risks, their confidence intervals, and their confidence bands based on right censored time to event data. We provide explicit formulas for our implementation of the estimator of the (stratified) baseline hazard function in the presence of tied event times. As a by-product we obtain fast access …

Partial Rank Data With The Hyper2 Package: Likelihood Functions For Generalized Bradley-Terry Models, Robin K. S Hankin

The R Journal

Here I present the hyper2 package for generalized Bradley-Terry models and give examples from two competitive situations: single scull rowing, and the competitive cooking game show Master Chef Australia. A number of natural statistical hypotheses may be tested straightforwardly using the software.

The R Journal (December 2017) 9(2): Complete Issue, The R Foundation

The R Journal

Editorial, Roger Bivand

Contributed Research Articles

anchoredDistr: A Package for the Bayesian Inversion of Geostatistical Parameters with Multi-type and Multi-scale Data, Heather Savoy, Falk Heße, and Yoram Rubin

dGAselID: An R Package for Selecting a Variable Number of Features in High Dimensional Data, Nicolae Teodor Melita and Stefan Holban

Allele Imputation and Haplotype Determination from Databases Composed of Nuclear Families, Nathan Medina-Rodríguez and Ángelo Santana

Visualization of Regression Models Using visreg, Patrick Breheny and Woodrow Burchett

fourierin: An R package to compute Fourier integrals, Guillermo Basulto-Elias, Alicia Carriquiry, Kris De Brabanter, and Daniel J. Nordman

Discrete Time Markov Chains with …

Anomalydetection: Implementation Of Augmented Network Log Anomaly Detection Procedures, Robert J. Gutierrez, Bradley C. Boehmke, Air Force Institute Of Technology, Cade M. Saie, Trevor J. Bihl

The R Journal

As the number of cyber-attacks continues to grow on a daily basis, so does the delay in threat detection. For instance, in 2015, the Office of Personnel Management discovered that approximately 21.5 million individual records of Federal employees and contractors had been stolen. On average, the time between an attack and its discovery is more than 200 days. In the case of the OPM breach, the attack had been going on for almost a year. Currently, cyber analysts inspect numerous potential incidents on a daily basis, but have neither the time nor the resources available to perform such a task. …

The Welchadf Package For Robust Hypothesis Testing In Unbalanced Multivariate Mixed Models With Heteroscedastic And Non-Normal Data, Pablo J. Villacorta

The R Journal

A new R package is presented for dealing with non-normality and variance heterogeneity of sample data when conducting hypothesis tests of main effects and interactions in mixed models. The proposal departs from an existing SAS program which implements Johansen’s general formulation of Welch-James’s statistic with approximate degrees of freedom, which makes it suitable for testing any linear hypothesis concerning cell means in univariate and multivariate mixed model designs when the data pose non-normality and non-homogeneous variance. Improved type I error rate control is obtained using bootstrapping for calculating an empirical critical value, whereas robustness against non-normality is achieved through trimmed …

Mle.Tools: An R Package For Maximum Likelihood Bias Correction, Josmar Mazucheli, André Felipe B. Menezes, Saralees Nadarajah

The R Journal

Recently, Mazucheli (2017) uploaded the package mle.tools to CRAN. It can be used for bias corrections of maximum likelihood estimates through the methodology proposed by Cox and Snell (1968). The main function of the package, coxsnell.bc(), computes the bias corrected maximum likelihood estimates. Although in general, the bias corrected estimators may be expected to have better sampling properties than the uncorrected estimators, analytical expressions from the formula proposed by Cox and Snell (1968) are either tedious or impossible to obtain. The purpose of this paper is twofolded: to introduce the mle.tools package, especially the coxsnell.bc() function; secondly, to compare, for …

Liureg: A Comprehensive R Package For The Liu Estimation Of Linear Regression Model With Collinear Regressors, Muhammad Imdadullah, Muhammad Aslam, Saima Altaf

The R Journal

The Liu regression estimator is now a commonly used alternative to the conventional ordinary least squares estimator that avoids the adverse effects in the situations when there exists a considerable degree of multicollinearity among the regressors. There are only a few software packages available for estimation of the Liu regression coefficients, though with limited methods to estimate the Liu biasing parameter without addressing testing procedures. Our liureg package can be used to estimate the Liu regression coefficients utilizing a range of different existing biasing parameters, to test these coefficients with more than 15 Liu related statistics, and to present different …

Ider: Intrinsic Dimension Estimation With R, Hideitsu Hino

The R Journal

In many data analyses, the dimensionality of the observed data is high while its intrinsic dimension remains quite low. Estimating the intrinsic dimension of an observed dataset is an essential preliminary step for dimensionality reduction, manifold learning, and visualization. This paper introduces an R package, named ider, that implements eight intrinsic dimension estimation methods, including a recently proposed method based on a second-order expansion of a probability mass function and a generalized linear model. The usage of each function in the package is explained with datasets generated using a function that is also included in the package

Carx: An R Package To Estimate Censored Autoregressive Time Series With Exogenous Covariates, Chao Wang, Kung-Sik Chan

The R Journal

We implement in the R package carx a novel and computationally efficient quasi-likelihood method for estimating a censored autoregressive model with exogenous covariates. The proposed quasi-likelihood method reduces to maximum likelihood estimation in absence of censoring. The carx package contains many useful functions for practical data analysis with censored stochastic regression, including functions for outlier detection, model diagnostics, and prediction with censored time series data. We illustrate the capabilities of the carx package with simulations and an elaborate real data analysis.

Furniture For Quantitative Scientists, Tyson S. Barrett, Emily Brignone

The R Journal

A basic understanding of the distributions of study variables and the relationships among them is essential to inform statistical modeling. This understanding is achieved through the computation of summary statistics and exploratory data analysis. Unfortunately, this step tends to be under-emphasized in the research process, in part because of the often tedious nature of thorough exploratory data analysis. The table1() function in the furniture package streamlines much of the exploratory data analysis process, making the computation and communication of summary statistics simple and beautiful while offering significant time-savings to the researcher

Ctmcd: An R Package For Estimating The Parameters Of A Continuous-Time Markov Chain From Discrete-Time Data, Marius Pfeuffer

The R Journal

This article introduces the R package ctmcd, which provides an implementation of methods for the estimation of the parameters of a continuous-time Markov chain given that data are only available on a discrete-time basis. This data consists of partial observations of the state of the chain, which are made without error at discrete times, an issue also known as the embedding problem for Markovchains. The functions provided comprise matrix logarithm based approximations as described in Israel et al. (2001), as well as Kreinin and Sidelnikova (2001), an expectation-maximization algorithm and a Gibbs sampling approach, both introduced by Bladt and …

Discrete Time Markov Chains With R, Giorgio Alfredo Spedicato

The R Journal

The markovchain package aims to provide S4 classes and methods to easily handle Discrete Time Markov Chains (DTMCs), filling the gap with what is currently available in the CRAN repository. In this work, I provide an exhaustive description of the main functions included in the package, as well as hands-on examples.

Fourierin: An R Package To Compute Fourier Integrals, Guillermo Basulto-Elias, Alicia Carriquiry, Kris De Brabanter, Daniel J. Nordman

The R Journal

We present the R package fourierin (Basulto-Elias, 2017) for evaluating functions defined as Fourier-type integrals over a collection of argument values. The integrals are finitely supported with integrands involving continuous functions of one or two variables. As an important application, such Fourier integrals arise in so-called “inversion formulas”, where one seeks to evaluate a probability density at a series of points from a given characteristic function (or vice versa) through Fourier transforms. This paper intends to fill a gap in current R software, where tools for repeated evaluation of functions as Fourier integrals are not directly available. We implement two …

Allele Imputation And Haplotype Determination From Databases Composed Of Nuclear Families, Nathan Medina-Rodríguez, Ángelo Santana

The R Journal

The alleHap package is designed for imputing genetic missing data and reconstruct non recombinant haplotypes from pedigree databases in a deterministic way. When genotypes of related individuals are available in a number of linked genetic markers, the program starts by identifying haplotypes compatible with the observed genotypes in those markers without missing values. If haplotypes are identified in parents or offspring, missing alleles can be imputed in subjects containing missing values. Several scenarios are analyzed: family completely genotyped, children partially genotyped and parents completely genotyped, children fully genotyped and parents containing entirely or partially missing genotypes, and founders and their …

Anchoreddistr: A Package For The Bayesian Inversion Of Geostatistical Parameters With Multi-Type And Multi-Scale Data, Heather Savoy, Falk Heße, Yoram Rubin

The R Journal

The Method of Anchored Distributions (MAD) is a method for Bayesian inversion designed for inferring both local (e.g. point values) and global properties (e.g. mean and variogram parameters) of spatially heterogenous fields using multi-type and multi-scale data. Software implementations of MAD exist in C++ and C# to import data, execute an ensemble of forward model simulations, and perform basic post-processing of calculating likelihood and posterior distributions for a given application. This article describes the R package anchoredDistr that has been built to provide an R based environment for this method. In particular, anchoredDistr provides a range of post-processing capabilities for …

Editorial, Roger Bivand

The R Journal

In my editorial for the 2017–1 issue, I concentrated on tabulating the status of this journal with respect to its authors and reviewers (updated tables below). This time, I was prompted by an interesting blog posting by Jan Wijffels of BNOSAC, describing the use of the udpipe package to apply natural language processing (NLP) to the CRAN package database available from tools::CRAN_package_db() since the release of R 3.4. The interactive NLP searcher is a dashboard permitting exploration of annotated CRAN package title and description NLP data.

Visualization Of Regression Models Using Visreg, Patrick Breheny, Woodrow Burchett

The R Journal

Regression models allow one to isolate the relationship between the outcome and an explanatory variable while the other variables are held constant. Here, we introduce an R package, visreg, for the convenient visualization of this relationship via short, simple function calls. In addition to estimates of this relationship, the package also provides pointwise confidence bands and partial residuals to allow assessment of variability as well as outliers and other deviations from modeling assumptions. The package provides several options for visualizing models with interactions, including lattice plots, contour plots, and both static and interactive perspective plots. The implementation of the …

Rqgis: Integrating R With Qgis For Statistical Geocomputing, Jannes Muenchow, Patrick Schratz, Alexander Brenning

The R Journal

Integrating R with Geographic Information Systems (GIS) extends R’s statistical capabilities with numerous geoprocessing and data handling tools available in a GIS. QGIS is one of the most popular open-source GIS, and it furthermore integrates other GIS programs such as the System for Automated Geoscientific Analyses (SAGA) GIS and the Geographic Resources Analysis Support System (GRASS) GIS within a single software environment. This and its QGIS Python API makes it a perfect candidate for console-based geoprocessing. By establishing an interface, the R package RQGIS makes it possible to use QGIS as a geoprocessing workhorse from within R. Compared to other …

A Tidy Data Model For Natural Language Processing Using Cleannlp, Taylor Arnold

The R Journal

Recent advances in natural language processing have produced libraries that extract low level features from a collection of raw texts. These features, known as annotations, are usually stored internally in hierarchical, tree-based data structures. This paper proposes a data model to represent annotations as a collection of normalized relational data tables optimized for exploratory data analysis and predictive modeling. The R package cleanNLP, which calls one of two state of the art NLP libraries (CoreNLP or spaCy), is presented as an implementation of this data model. It takes raw text as an input and returns a list of normalized …

R Foundation News, Torsten Hothorn

The R Journal

Donations and members

Donations

Supporting benefactors

Supporting members

Splitting It Up: The Spduration Split-Population Duration Regression Package For Time-Varying Covariates, Andreas Beger, Daniel W. Hill Jr, Nils W. Metternich, Shahryar Minhas, Michael D. Ward

The R Journal

We present an implementation of split-population duration regression in the spduration (Beger et al., 2017) package for R that allows for time-varying covariates. The statistical model accounts for units that are immune to a certain outcome and are not part of the duration process the researcher is primarily interested in. We provide insights for when immune units exist, that can significantly increase the predictive performance compared to standard duration models. The package includes estimation and several post-estimation methods for split-population Weibull and log-logistic models. Weprovide an empirical application to data on military coups.

Changes In R, R Core Team

The R Journal

CHANGES IN R 3.4.3