Numerical Analysis and Scientific Computing Commons^™

Pstat: An R Package To Assess Population Differentiation In Phenotypic Traits, Stéphane Blondeau Da Silva, Anne Da Silva

The R Journal

The package Pstat calculates P_ST values to assess differentiation among populations from a set of quantitative traits and provides bootstrapped distributions and confidence intervals for P_ST. Variations of P_ST as a function of the parameter c/h² are studied as well. The package implements different transformations of the measured phenotypic traits to eliminate variation resulting from allometric growth, including calculation of residuals from linear regression, Reist standardization, and the Aitchison transformation.

Residuals And Diagnostics For Binary And Ordinal Regression Models: An Introduction To The Sure Package, Brandon M. Greenwell, Andrew J. Mccarthy, Bradley C. Boehmke, Dungang Liu

The R Journal

Residual diagnostics is an important topic in the classroom, but it is less often used in practice when the response is binary or ordinal. Part of the reason for this is that generalized models for discrete data, like cumulative link models and logistic regression, do not produce standard residuals that are easily interpreted as those in ordinary linear regression. In this paper, we introduce the R package sure, which implements a recently developed idea of SUrrogate REsiduals. We demonstrate the utility of the package in detection of cumulative link model misspecification with respect to mean structures, link functions, …

Hrm: An R Package For Analysing High-Dimensional Multi-Factor Repeated Measures Authors: Martin Happ, Solomon W. Harrar And Arne C. Bathke, Martin Happ, Solomon W. Harrar, Arne C. Bathke

The R Journal

High-dimensional longitudinal data pose a serious challenge for statistical inference as many test statistics cannot be computed for high-dimensional data, or they do not maintain the nominal type-I error rate, or have very low power. Therefore, it is necessary to derive new inference methods capable of dealing with high dimensionality, and to make them available to statistics practitioners. One such method is implemented in the package HRM described in this article. This new method uses a similar approach as the Welch-Satterthwaite t-test approximation and works very well for high-dimensional data as long as the data distribution is not too skewed …

R Foundation News, Torsten Hothorn

The R Journal

Donations and members

Donations

Supporting benefactors

Supporting institutions

Supporting members

Inventorymodel: An R Package For Centralized Inventory Problems, Alejandro Saavedra-Nieves

The R Journal

Inventory management of goods is an integral part of logistics systems; common to various economic sectors such as industry, agriculture and trade; and independent of production volume. In general, as companies seek to minimize economic losses, studies on problems of multi-agent inventory have increased in recent years. A multi-agent inventory problem is a situation in which several agents face individual inventory problems and agree to coordinate their orders with the objective of reducing their costs. The R package Inventorymodel allows the determination of both the optimal policy for some inventory situations with deterministic demands and the allocation of costs from …

Rpostgis: Linking R With A Postgis Spatial Database, David Bucklin, Mathieu Basille

The R Journal

With the proliferation of sensors and the ease of data collection from online sources, large datasets have become the norm in many scientific disciplines, and efficient data storage, management, and retrival is imperative for large research projects. Relational databases provide a solution, but in order to be useful, must be able to be linked to analysis and visualization tools, such as R. Here, we present a package intended to facilitate integration of R with the open-source database software PostgreSQL, with a focus on its spatial extension, PostGIS. The package rpostgis (version 1.4.1) provides methods for spatial data handling (vector and …

R Package Imputetestbench To Compare Imputation Methods For Univariate Time Series, Marcus W. Beck, Neeraj Bokde, Gualberto Asencio-Cortés, Kishore Kulat

The R Journal

Missing observations are common in time series data and several methods are available to impute these values prior to analysis. Variation in statistical characteristics of univariate time series can have a profound effect on characteristics of missing observations and, therefore, the accuracy of different imputation methods. The imputeTestbench package can be used to compare the prediction accuracy of different methods as related to the amount and type of missing data for a user-supplied dataset. Missing data are simulated by removing observations completely at random or in blocks of different sizes depending on characteristics of the data. Several imputation algorithms are …

Advanced Bayesian Multilevel Modeling With The R Package Brms, Paul-Christian Bürkner

The R Journal

The brms package allows R users to easily specify a wide range of Bayesian single-level and multilevel models which are fit with the probabilistic programming language Stan behind the scenes. Several response distributions are supported, of which all parameters (e.g., location, scale, and shape) can be predicted. Non-linear relationships may be specified using non-linear predictor terms or semi-parametric approaches such as splines or Gaussian processes. Multivariate models can be fit as well. To make all of these modeling options possible in a multilevel framework, brms provides an intuitive and powerful formula syntax, which extends the well known formula syntax of …

Welfare, Inequality And Poverty Analysis With Rtip: An Approach Based On Stochastic Dominance, Angel Berihuete, Carmen D. Ramos, Miguel A. Sordo

The R Journal

Disparities in economic welfare, inequality and poverty across and within countries are of great interest to sociologists, economists, researchers, social organizations and political scientists. Information about these topics is commonly based on surveys. We present a package called rtip that implements techniques based on stochastic dominance to make unambiguous comparisons, in terms of welfare, poverty and inequality, among income distributions. Besides providing point estimates and confidence intervals for the most commonly used indicators of these characteristics, the package rtip estimates the usual Lorenz curve, the generalized Lorenz curve, the TIP (Three I’s of Poverty) curve and allows to test statistically …

Support Vector Machines For Survival Analysis With R, Césaire J.K. Fouodo, Inke R. König, Claus Weihs, Andreas Ziegler, Marvin N. Wright

The R Journal

This article introduces the R package survivalsvm, implementing support vector machines for survival analysis. Three approaches are available in the package: The regression approach takes censoring into account when formulating the inequality constraints of the support vector problem. In the ranking approach, the inequality constraints set the objective to maximize the concordance index for comparable pairs of observations. The hybrid approach combines the regression and ranking constraints in a single model. We describe survival support vector machines and their implementation, provide examples and compare the prediction performance with the Cox proportional hazards model, random survival forests and gradient boosting using …

Mglm: An R Package For Multivariate Categorical Data Analysis, Juhyun Kim, Yiwen Zhang, Joshua Day, Hua Zhou

The R Journal

Data with multiple responses is ubiquitous in modern applications. However, few tools are available for regression analysis of multivariate counts. The most popular multinomial-logit model has a very restrictive mean-variance structure, limiting its applicability to many data sets. This article introduces an R package MGLM, short for multivariate response generalized linear models, that expands the current tools for regression analysis of polytomous data. Distribution fitting, random number generation, regression, and sparse regression are treated in a unifying framework. The algorithm, usage, and implementation details are discussed.

Cchs: An R Package For Stratified Case-Cohort Studies, Edmund Jones

The R Journal

The cchs package contains a function, also called cchs, for analyzing data from a stratified case-cohort study, as used in epidemiology. For data from this type of study, cchs calculates Estimator III of Borgan et al. (2000), which is a score-unbiased estimator for the regression coefficients in the Cox proportional hazards model. From the user’s point of view, the function is similar to coxph (in the survival package) and other widely used model-fitting functions. Convenient software has not previously been available for Estimator III since it is complicated to calculate. SAS and S-Plus code-fragments for the calculation have been published, …

Small Area Disease Risk Estimation And Visualization Using R, Paula Moraga

The R Journal

Small area disease risk estimation is essential for disease prevention and control. In this paper, we demonstrate how R can be used to obtain disease risk estimates and quantify risk factors using areal data. We explain how to define disease risk models and how to perform Bayesian inference using the INLA package. We also show how to make interactive maps of estimates using the leaflet package to better understand the disease spatial patterns and communicate the results. We show an example of lung cancer risk in Pennsylvania, United States, in year 2002, and demonstrate that R represents an excellent tool …

Semiparametric Generalized Linear Models With The Gldrm Package, Michael J. Wurm, Paul J. Rathouz

The R Journal

This paper introduces a new algorithm to estimate and perform inferences on a recently proposed and developed semiparametric generalized linear model (glm). Rather than selecting a particular parametric exponential family model, such as the Poisson distribution, this semiparametric glm assumes that the response is drawn from the more general exponential tilt family. The regression coefficients and unspecified reference distribution are estimated by maximizing a semiparametric likelihood. The new algorithm incorporates several computational stability and efficiency improvements over the algorithm originally proposed. In particular, the new algorithm performs well for either small or large support for the nonparametric response distribution. The …

Ratingscalereduction Package: Stepwise Rating Scale Item Reduction Without Predictability Loss, Waldemar W. Koczkodaj, Feng Li, Alicja Wolny–Dominiak

The R Journal

This study presents an innovative method for reducing the number of rating scale items without predictability loss. The “area under the receiver operator curve” method (AUC ROC) is used for the stepwise method of reducing items of a rating scale. RatingScaleReduction R package contains the presented implementation. Differential evolution (a metaheuristic for optimization) was applied to one of the analyzed datasets to illustrate that the presented stepwise method can be used with other classifiers to reduce the number of rating scale items (variables). The targeted areas of application are decision making, data mining, machine learning, and psychometrics.

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.7 was released on 1 May, 2018. It is compatible with R 3.5.1 and consists of 1560 software packages, 342 experiment data packages, and 919 up-to-date annotation packages. The release announcement includes descriptions of 98 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

A System For An Accountable Data Analysis Process In R, Jonathan Gelfond, Martin Goros, Brian Hernandez, Alex Bokov

The R Journal

Efficiently producing transparent analyses may be difficult for beginners or tedious for the experienced. This implies a need for computing systems and environments that can efficiently satisfy reproducibility and accountability standards. To this end, we have developed a system, R package, and R Shiny application called adapr (Accountable Data Analysis Process in R) that is built on the principle of accountable units. An accountable unit is a data file (statistic, table or graphic) that can be associated with a provenance, meaning how it was created, when it was created and who created it, and this is similar to the ’verifiable …

Mmpf: Monte-Carlo Methods For Prediction Functions, Zachary M. Jones

The R Journal

Machine learning methods can often learn high-dimensional functions which generalize well but are not human interpretable. The mmpf package marginalizes prediction functions using Monte-Carlo methods, allowing users to investigate the behavior of these learned functions, as on a lower dimensional subset of input features: partial dependence and variations thereof. This makes machine learning methods more useful in situations where accurate prediction is not the only goal, such as in the social sciences where linear models are commonly used because of their interpretability.

Many methods for estimating prediction functions produce estimated functions which are not directly human-interpretable because of their complexity: …

Generalized Additive Model Multiple Imputation By Chained Equations With Package Imputerobust, Daniel Salfran, Martin Spiess

The R Journal

Data analysis, common to all empirical sciences, often requires complete data sets. Unfortunately, real world data collection will usually result in data values not being observed. We present a package for robust multiple imputation (the ImputeRobust package) that allows the use of generalized additive models for location, scale, and shape in the context of chained equations. The paper describes the basics of the imputation technique which builds on a semi-parametric regression model (GAMLSS) and the algorithms and functions provided with the corresponding package. Furthermore, some illustrative examples are provided.

Fhdi: An R Package For Fractional Hot Deck Imputation, Jongho Im, In Ho Cho, Jae Kwang Kim

The R Journal

Fractional hot deck imputation (FHDI), proposed by Kalton and Kish (1984) and investigated by Kim and Fuller (2004), is a tool for handling item nonresponse in survey sampling. In FHDI, each missing item is filled with multiple observed values yielding a single completed data set for subsequent analyses. An R package FHDI is developed to perform FHDI and also the fully efficient fractional imputation (FEFI) method of (Fuller and Kim, 2005) to impute multivariate missing data with arbitrary missing patterns. FHDI substitutes missing items with a few observed values jointly obtained from a set of donors whereas the FEFI uses …

The R Journal (July 2018) 10(1): Complete Issue, The R Foundation

The R Journal

Editorial, John Verzani

Contributed Research Articles

A System for an Accountable Data Analysis Process in R, Jonathan Gelfond, Martin Goros, Brian Hernandez and Alex Bokov

RealVAMS: An R Package for Fitting a Multivariate Value-added Model (VAM), Jennifer Broatch, Jennifer Green, and Andrew Karl

InfoTrad: An R Package for Estimating the Probability of Informed Trading, Duygu Çelik and Murat Tiniç

RatingScaleReduction Package: Stepwise Rating Scale Item Reduction without Predictability Loss, Waldemar W. Koczkodaj, Feng Li, and Alicja Wolny-Dominiak

mmpf: Monte-Carlo Methods for Prediction Functions, Zachary M. Jones

Generalized Additive Model Multiple Imputation by Chained Equations with Package ImputeRobust, Daniel Salfran and …

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 …