Numerical Analysis and Scientific Computing Commons^™

An Introduction To Principal Surrogate Evaluation With The Pseval Package, Michael C. Sachs, Erin E. Gabriel

The R Journal

We describe a new package called pseval that implements the core methods for the evaluation of principal surrogates in a single clinical trial. It provides a flexible interface for defining models for the risk given treatment and the surrogate, the models for integration over the missing counterfactual surrogate responses, and the estimation methods. Estimated maximum likelihood and pseudo-score can be used for estimation, and the bootstrap for inference. A variety of post-estimation methods are provided, including print, summary, plot, and testing. We summarize the main statistical methods that are implemented in the package and illustrate its use from the perspective …

Quantreg.Nonpar: An R Package For Performing Nonparametric Series Quantile Regression, Michael Lipsitz, Alexandre Belloni, Victor Chernozhukov, Iván Fernández-Val

The R Journal

The R package quantreg.nonpar implements nonparametric quantile regression methods to estimate and make inference on partially linear quantile models. quantreg.nonpar obtains point estimates of the conditional quantile function and its derivatives based on series approximations to the nonparametric part of the model. It also provides point-wise and uniform confidence intervals over a region of covariate values and/or quantile indices for the same functions using analytical and resampling methods. This paper serves as an introduction to the package and displays basic functionality of the functions contained within.

Two-Tier Latent Class Irt Models In R, Silvia Bacci, Francesco Bartolucci

The R Journal

In analyzing data deriving from the administration of a questionnaire to a group of individuals, Item Response Theory (IRT) models provide a flexible framework to account for several aspects involved in the response process, such as the existence of multiple latent traits. In this paper, we focus on a class of semi-parametric multidimensional IRT models, in which these traits are represented through one or more discrete latent variables; these models allow us to cluster individuals into homogeneous latent classes and, at the same time, to properly study item characteristics. In particular, we follow a within-item multidimensional formulation similar to that …

Qtools: A Collection Of Models And Tools For Quantile Inference, Marco Geraci

The R Journal

Quantiles play a fundamental role in statistics. The quantile function defines the distribution of a random variable and, thus, provides a way to describe the data that is specular but equivalent to that given by the corresponding cumulative distribution function. There are many advantages in working with quantiles, starting from their properties. The renewed interest in their usage seen in the last years is due to the theoretical, methodological, and software contributions that have broadened their applicability. This paper presents the R package Qtools, a collection of utilities for unconditional and conditional quantiles.

R Foundation News, Torsten Hothorn

The R Journal

New benefectors

Donations

New supporting institutions

New supporting members

Ggfortify: Unified Interface To Visualize Statistical Results Of Popular R Packages, Yuan Yuan, Masaaki Horikoshi, Wenxuan Li

The R Journal

The ggfortify package provides a unified interface that enables users to use one line of code to visualize statistical results of many R packages using ggplot2 idioms. With the help of ggfortify, statisticians, data scientists, and researchers can avoid the sometimes repetitive work of using the ggplot2 syntax to achieve what they need.

Ake: An R Package For Discrete And Continuous Associated Kernel Estimations, Wanbitching E. Wansouwé, Sobom M. Somé, Célestin C. Kokonendji

The R Journal

Kernel estimation is an important technique in exploratory data analysis. Its utility relies on its ease of interpretation, especially based on graphical means. The Ake package is introduced for univariate density or probability mass function estimation and also for continuous and discrete regression functions using associated kernel estimators. These associated kernels have been proposed due to their specific features of variables of interest. The package focuses on associated kernel methods appropriate for continuous (bounded, positive) or discrete (count, categorical) data often found in applied settings. Furthermore, optimal bandwidths are selected by cross-validation for any associated kernel and by Bayesian methods …

Eicompare: Comparing Ecological Inference Estimates Across Ei And Ei:Rc, Loren Collingwood, Kassra Oskooii, Sergio Garcia-Rios, Matt Barreto

The R Journal

Social scientists and statisticians often use aggregate data to predict individual-level behavior because the latter are not always available. Various statistical techniques have been developed to make inferences from one level (e.g., precinct) to another level (e.g., individual voter) that minimize errors associated with ecological inference. While ecological inference has been shown to be highly problematic in a wide array of scientific fields, many political scientists and analysis employ the techniques when studying voting patterns. Indeed, federal voting rights lawsuits now require such an analysis, yet expert reports are not consistent in which type of ecological inference is used. This …

Calculating Biological Module Enrichment Or Depletion And Visualizing Data On Large-Scale Molecular Maps With Acsnminer And Rnavicell Packages, Paul Deveau, Emmanuel Barillot, Valentina Boeva, Andrei Zinovyev, Eric Bonnet

The R Journal

Biological pathways or modules represent sets of interactions or functional relationships occurring at the molecular level in living cells. A large body of knowledge on pathways is organized in public databases such as the KEGG, Reactome, or in more specialized repositories, the Atlas of Cancer Signaling Network (ACSN) being an example. All these open biological databases facilitate analyses, improving our understanding of cellular systems. We hereby describe ACSNMineR for calculation of enrichment or depletion of lists of genes of interest in biological pathways. ACSNMineR integrates ACSNmolecular pathways gene sets, but can use any gene set encoded as a GMT file, …

Escape From Boxland, Barret Schloerke, Hadley Wickham, Dianne Cook, Heike Hofmann

The R Journal

A library of common geometric shapes can be used to train our brains for understanding data structure in high-dimensional Euclidean space. This article describes the methods for producing cubes, spheres, simplexes, and tori in multiple dimensions. It also describes new ways to define and generate high-dimensional tori. The algorithms are described, critical code chunks are given, and a large collection of generated data are provided. These are available in the R package geozoo, and selected movies and images, are available on the GeoZoo web site (http://schloerke.github.io/geozoo/)

Changes On Cran, Kurt Hornik, Achim Zeileis

The R Journal

CRAN growth

New CRAN task views

New packages in CRAN task views

Easyroc: An Interactive Web-Tool For Roc Curve Analysis Using R Language Environment, Dincer Goksuluk, Selcuk Korkmaz, Gokmen Zararsiz, A Ergun Karaagaoglu

The R Journal

ROC curve analysis is a fundamental tool for evaluating the performance of a marker in a number of research areas, e.g., biomedicine, bioinformatics, engineering etc., and is frequently used for discriminating cases from controls. There are a number of analysis tools which are used to guide researchers through their analysis. Some of these tools are commercial and provide basic methods for ROC curve analysis while others offer advanced analysis techniques and a command-based user interface, such as the R environment. The R environment includes comprehensive tools for ROC curve analysis; however, using a command-based interface might be challenging and time …

Hdm: High-Dimensional Metrics, Victor Chernozhukov, Chris Hansen, Martin Spindler

The R Journal

In this article the package High-dimensional Metrics hdm is introduced. It is a collection of statistical methods for estimation and quantification of uncertainty in high-dimensional approximately sparse models. It focuses on providing confidence intervals and significance testing for (possibly many) low-dimensional subcomponents of the high-dimensional parameter vector. Efficient estimators and uniformly valid confidence intervals for regression coefficients on target variables (e.g., treatment or policy variable) in a high-dimensional approximately sparse regression model, for average treatment effect (ATE) and average treatment effect for the treated (ATET), as well for extensions of these param eters to the endogenous setting are provided. Theory …

Distance Measures For Time Series In R: The Tsdist Package, Usue Mori, Alexander Mendiburu, Jose A. Lozano

The R Journal

The definition of a distance measure between time series is crucial for many time series data mining tasks, such as clustering and classification. For this reason, a vast portfolio of time series distance measures has been published in the past few years. In this paper, the TSdist package is presented, a complete tool which provides a unified framework to calculate the largest variety of time series dissimilarity measures available in R at the moment, to the best of our knowledge. The package implements some popular distance measures which were not previously available in R, and moreover, it also provides wrappers …

Condsurv: An R Package For The Estimation Of The Conditional Survival Function For Ordered Multivariate Failure Time Data, Luis Meira-Machado, Meira-Machado Sestelo

The R Journal

One major goal in clinical applications of time-to-event data is the estimation of survival with censored data. The usual nonparametric estimator of the survival function is the time-honored Kaplan-Meier product-limit estimator. Though this estimator has been implemented in several R packages, the development of the condSURV R package has been motivated by recent contributions that allow the estimation of the survival function for ordered multivariate failure time data. The condSURV package provides three different approaches all based on the Kaplan-Meier estimator. In one of these approaches these quantities are estimated conditionally on current or past covariate measures. Illustration of the …

Measurement Units In R, Edzer Pebesma, Thomas Mailund, James Hiebert

The R Journal

We briefly review SI units, and discuss R packages that deal with measurement units, their compatibility and conversion. Built upon udunits2 and the UNIDATA udunits library, we introduce the package units that provides a class for maintaining unit metadata. When used in expression, it automatically converts units, and simplifies units of results when possible; in case of incompatible units, errors are raised. The class flexibly allows expansion beyond predefined units. Using units may eliminate a whole class of potential scientific programming mistakes. We discuss the potential and limitations of computing with explicit units.

Tigris: An R Package To Access And Work With Geographic Data From The Us Census Bureau, Kyle Walker

The R Journal

TIGER/Line shapefiles from the United States Census Bureau are commonly used for the mapping and analysis of US demographic trends. The tigris package provides a uniform interface for R users to download and work with these shapefiles. Functions in tigris allow R users to request Census geographic datasets using familiar geographic identifiers and return those datasets as objects of class "Spatial*DataFrame". In turn, tigris ensures consistent and high-quality spatial data for R users’ cartographic and spatial analysis projects that involve US Census data. This article provides an overview of the functionality of the tigris package, and concludes with an applied …

Mixtox: An R Package For Mixture Toxicity Assessment, Xiang-Wei Zhu, Jian-Yi Chen

The R Journal

Mixture toxicity assessment is indeed necessary for humans and ecosystems that are continually exposed to a variety of chemical mixtures. This paper describes an R package, called mixtox, which offers a general framework of curve fitting, mixture experimental design, and mixture toxicity prediction for practitioners in toxicology. The unique features of mixtox include: (1) constructing a uniform table for mixture experimental design; and (2) predicting toxicity of a mixture with multiple components based on reference models such as concentration addition, independent action, and generalized concentration addition. We describe the various functions of the package and provide examples to illustrate their …

Water: Tools And Functions To Estimate Actual Evapotranspiration Using Land Surface Energy Balance Models In R, Guillermo Federico Olmedo, Samuel Ortega-Farías, Daniel De La Fuente-Sáiz, David Fonseca- Luego, Fernando Fuentes-Peñailillo

The R Journal

The crop water requirement is a key factor in the agricultural process. It is usually estimated throughout actual evapotranspiration (ETa). This parameter is the key to develop irrigation strategies, to improve water use efficiency and to understand hydrological, climatic, and ecosystem processes. Currently, it is calculated with classical methods, which are difficult to extrapolate, or with land surface energy balance models (LSEB), such as METRIC and SEBAL, which are based on remote sensing data. This paper describes water, an open implementation of LSEB. The package provides several functions to estimate the parameters of the LSEB equation from satellite data …

Changes In R, R Core Team

The R Journal

CHANGES IN R 3.3.1 patched

CHANGES IN R 3.3.1

CHANGES IN R 3.3.0

Changes On Cran, Kurt Hornik, Achim Zeileis

The R Journal

In the past 8 months,1322 new packages were added to the CRAN package repository. 43 packages were unarchived,48 archived,1 package had to be removed.The following shows the growth of the number of active packages in the CRAN package repository:

Nonparametric Tests For The Interaction In Two-Way Factorial Designs Using R, Jos Feys

The R Journal

An increasing number of R packages include nonparametric tests for the interaction in two-way factorial designs. This paper briefly describes the different methods of testing and reports the resulting p-values of such tests on datasets for four types of designs: between, within, mixed, and pretest-posttest designs. Potential users are advised only to apply tests they are quite familiar with and not be guided by p-values for selecting packages and tests.

Using Decipher V2.0 To Analyze Big Biological Sequence Data In R, Erik S. Wright

The R Journal

In recent years, the cost of DNA sequencing has decreased at a rate that has outpaced improvements in memory capacity. It is now common to collect or have access to many gigabytes of biological sequences. This has created an urgent need for approaches that analyze sequences in subsets without requiring all of the sequences to be loaded into memory at one time. It has also opened opportunities to improve the organization and accessibility of information acquired in sequencing projects. The DECIPHER package offers solutions to these problems by assisting in the curation of large sets of biological sequences stored in …

Scmamp: Statistical Comparison Of Multiple Algorithms In Multiple Problems, Borja Calvo, Guzmán Santafé

The R Journal

Comparing the results obtained by two or more algorithms in a set of problems is a central task in areas such as machine learning or optimization. Drawing conclusions from these comparisons may require the use of statistical tools such as hypothesis testing. There are some interesting papers that cover this topic. In this manuscript we present scmamp, an R package aimed at being a tool that simplifies the whole process of analyzing the results obtained when comparing algorithms, from loading the data to the production of plots and tables.

Comparing the performance of different algorithms is an essential step …

Swmpr: An R Package For Retrieving, Organizing, And Analyzing Environmental Data For Estuaries, Marcus W. Beck

The R Journal

The System-Wide Monitoring Program (SWMP)was implemented in 1995 by the US National Estuarine Research Reserve System. This program has provided two decades of continuous monitoring data at over 140 fixed stations in 28 estuaries. However, the increasing quantity of data provided by the monitoring network has complicated broad-scale comparisons between systems and, in some cases, prevented simple trend analysis of water quality parameters at individual sites. This article describes the SWMPr package that provides several functions that facilitate data retrieval, organization, and analysis of time series data in the reserve estuaries. Previously unavailable functions for estuaries are also provided to …

Spatio-Temporal Interpolation Using Gstat, Benedikt Gräler, Edzer Pebesma, Gerard Heuvelink

The R Journal

We present new spatio-temporal geostatistical modelling and interpolation capabilities of the R package gstat. Various spatio-temporal covariance models have been implemented, such as the separable, product-sum, metric and sum-metric models. Inareal-world application we comparespatio temporal interpolations using these models with a purely spatial kriging approach. The target variable of the application is the daily mean PM10 concentration measured at rural air quality monitoring stations across Germany in 2005. R code for variogram fitting and interpolation is presented in this paper to illustrate the workflow of spatio-temporal interpolation using gstat. We conclude that the system works properly and that the …

Model Builder For Item Factor Analysis With Openmx, Joshua N. Pritikin, Karen M. Schmidt

The R Journal

We introduce a shiny web application to facilitate the construction of Item Factor Analysis (a.k.a. Item Response Theory) models using the OpenMx package. The web application assists with importing data, outcome recoding, and model specification. However, the app does not conduct any analysis but, rather, generates an analysis script. Generated Rmarkdown output serves dual purposes: to analyze a data set and demonstrate good programming practices. The app can be used as a teaching tool or as a starting point for custom analysis scripts.

Quickpsy: An R Package To Fit Psychometric Functions For Multiple Groups, Daniel Linares, Joan López-Moliner

The R Journal

quickpsy is a package to parametrically fit psychometric functions. In comparison with previous R packages, quickpsy was built to easily fit and plot data for multiple groups. Here, we describe the standard parametric model used to fit psychometric functions and the standard estimation of its parameters using maximum likelihood. We also provide examples of usage of quickpsy, including how allowing the lapse rate to vary can sometimes eliminate the bias in parameter estimation, but not in general. Finally, we describe some implementation details, such as how to avoid the problems associated to round-off errors in the maximisation of the …

Variable Clustering In High-Dimensional Linear Regression: The R Package Clere, Loïc Yengo, Julien Jacques, Christophe Biernacki, Mickael Canouil

The R Journal

Dimension reduction is one of the biggest challenges in high-dimensional regression models. We recently introduced a new methodology based on variable clustering as a means to reduce dimensionality. We present here the R package clere that implements some refinements of this methodology. An overview of the package functionalities as well as examples to run an analysis are described. Numerical experiments on real data were performed to illustrate the good predictive performance of our parsimonious method compared to standard dimension reduction approaches.

Maps, Coordinate Reference Systems And Visualising Geographic Data With Mapmisc, Patrick E. Brown

The R Journal

The mapmisc package provides functions for visualising geospatial data, including fetching background map layers, producing colour scales and legends, and adding scale bars and orientation arrows to plots. Background maps are returned in the coordinate reference system of the dataset supplied, and inset maps and direction arrows reflect the map projection being plotted. This is a “light weight” package having an emphasis on simplicity and ease of use