Computer Sciences Commons^™

Ngspatial: A Package For Fitting The Centered Autologistic And Sparse Spatial Generalized Linear Mixed Models For Areal Data, John Hughes

The R Journal

Two important recent advances in areal modeling are the centered autologistic model and the sparse spatial generalized linear mixed model (SGLMM), both of which are reparameterizations of traditional models. The reparameterizations improve regression inference by alleviating spatial confounding, and the sparse SGLMM also greatly speeds computing by reducing the dimension of the spatial random effects. Package ngspatial (’ng’ = non-Gaussian) provides routines for fitting these new models. The package supports composite likelihood and Bayesian inference for the centered autologistic model, and Bayesian inference for the sparse SGLMM.

Coordinate-Based Meta-Analysis Of Fmri Studies With R, Andrea Stocco

The R Journal

This paper outlines how to conduct a simple meta-analysis of neuroimaging foci of activation in R. In particular, the first part of this paper reviews the nature of fMRI data, and presents a brief overview of the existing packages that can be used to analyze fMRI data in R. The second part illustrates how to handle fMRI data by showing how to visualize the results of different neuroimaging studies in a so-called orthographic view, where the spatial distribution of the foci of activation from different fMRI studies can be inspected visually.

Functional MRI (fMRI) is one of the most important …

Bshazard: A Flexible Tool For Nonparametric Smoothing Of The Hazard Function, Paola Rebora, Agus Salim, Marie Reilly

The R Journal

The hazard function is a key component in the inferential process in survival analysis and relevant for describing the pattern of failures. However, it is rarely shown in research papers due to the difficulties in nonparametric estimation. We developed the bshazard package to facilitate the computation of a nonparametric estimate of the hazard function, with data-driven smoothing. The method accounts for left truncation, right censoring and possible covariates. B-splines are used to estimate the shape of the hazard within the generalized linear mixed models frame work. Smoothness is controlled by imposing an autoregressive structure on the baseline hazard coefficients. This …

Prinsimp, Jonathan Zhang, Nancy Heckman, Davor Cubranic, Joel G. Kingsolver, Travis Gaydos, J.S. Marron

The R Journal

Principal Components Analysis (PCA) is a common way to study the sources of variation in a high-dimensional data set. Typically, the leading principal components are used to understand the variation in the data or to reduce the dimension of the data for subsequent analysis. The remaining principal components are ignored since they explain little of the variation in the data. However, the space spanned by the low variation principal components may contain interesting structure, structure that PCA cannot find. Prinsimp is an R package that looks for interesting structure of low variability. “Interesting” is defined in terms of a simplicity …

Phaser: An R Package For Phase Plane Analysis Of Autonomous Ode Systems, Michael J. Grayling

The R Journal

When modelling physical systems, analysts will frequently be confronted by differential equations which cannot be solved analytically. In this instance, numerical integration will usually be the only way forward. However, for autonomous systems of ordinary differential equations (ODEs) in one or two dimensions, it is possible to employ an instructive qualitative analysis foregoing this requirement, using so-called phase plane methods. Moreover, this qualitative analysis can even prove to be highly useful for systems that can be solved analytically, or will be solved numerically anyway. The package phaseR allows the user to perform such phase plane analyses: determining the stability of …

Taming Pitchf/X Data With Xml2r And Pitchrx, Carson Sievert

The R Journal

XML2R is a framework that reduces the effort required to transform XML content into tables in a way that preserves parent to child relationships. pitchRx applies XML2R’s grammar for XML manipulation to Major League Baseball Advanced Media (MLBAM)’s Gameday data. With pitchRx, one can easily obtain and store Gameday data in a remote database. The Gameday website hosts a wealth of XML data, but perhaps most interesting is PITCHf/x. Among other things, PITCHf/x data can be used to recreate a baseball’s flight path from a pitcher’s hand to home plate. With pitchRx, one can easily create animations …

Applying Spartan To Understand Parameter Uncertainty In Simulations, Kieran Alden, Mark Read, Paul S. Andrews, Jon Timmis, Mark Coles

The R Journal

In attempts to further understand the dynamics of complex systems, the application of computer simulation is becoming increasingly prevalent. Whereas a great deal of focus has been placed in the development of software tools that aid researchers develop simulations, similar focus has not been applied in the creation of tools that perform a rigorous statistical analysis of results generated through simulation: vital in understanding how these results offer an insight into the captured system. This encouraged us to develop spartan, a package of statistical techniques designed to assist researchers in understanding the relationship between their simulation and the real system. …

The R Journal (December 2014) 6(2): Complete Issue, The R Foundation

The R Journal

Editorial, Deepayan Sarkar

Contributed Research Articles

Coordinate-Based Meta-Analysis of fMRI Studies with R, Andrea Stocco

Automatic Conversion of Tables to LongForm Dataframes, Jimmy Oh

Prinsimp, Jonathan Zhang, Nancy Heckman, Davor Cubranic, Joel G. Kingsolver, Travis Gaydos, and J. S. Marron

phaseR: An R Package for Phase Plane Analysis of Autonomous ODE Systems, Michael J. Grayling

Flexible R Functions for Processing Accelerometer Data, with Emphasis on NHANES 2003-2006, Dane R. Van Domelen and W. Stephen Pittard

Applying spartan to Understand Parameter Uncertainty in Simulations, Kieran Alden, Mark Read, Paul S. Andrews, Jon Timmis, and Mark Coles

ngspatial: A Package for Fitting …

Boscor: Extending R From The Desktop To The Grid, Derek J. Weitzel, Jaime Frey, Marco Mambelli, Dan Fraser, Miha Ahronovitz, David Swanson

Holland Computing Center: Faculty Publications

In this paper, we describe a framework to execute R functions on remote resources from the desktop using Bosco. The R language is attractive to researchers because of its high level programming constructs which lower the barrier of entry for use. As the use of the R programming language in HPC and High Throughput Computing (HTC) has grown, so too has the need for parallel libraries in order to utilize computing resources.

Bosco is middleware that uses common protocols to manage job submissions to a variety of remote computational platforms and resources. The researcher is able to control and monitor …

Cubic Spline Interpolation By Solving A Recurrence Equation Instead Of A Tridiagonal Matrix, Peter Revesz

School of Computing: Conference and Workshop Papers

The cubic spline interpolation method is proba- bly the most widely-used polynomial interpolation method for functions of one variable. However, the cubic spline method requires solving a tridiagonal matrix-vector equation with an O(n) computational time complexity where n is the number of data measurements. Even an O(n) time complexity may be too much in some time-ciritical applications, such as continuously estimating and updating the flight paths of moving objects. This paper shows that under certain boundary conditions the tridiagonal matrix solving step of the cubic spline method could be entirely eliminated and instead the coefficients of the unknown cubic polynomials …

Designing A Bayer Filter With Smooth Hue Transition Interpolation Using The Xilinx System Generator, Zhiqiang Li, Peter Revesz

School of Computing: Conference and Workshop Papers

This paper describes the design of a Bayer filter with smooth hue transition using the System Generator for DSP. We describe and compare experimentally two different designs, one based on a MATLAB implementation and the other based on a modification of the Bayer filter using bilinear interpolation.

A Game-Theoretic Analysis Of The Nuclear Non-Proliferation Treaty, Peter Revesz

School of Computing: Conference and Workshop Papers

Although nuclear non-proliferation is an almost universal human desire, in practice, the negotiated treaties appear unable to prevent the steady growth of the number of states that have nuclear weapons. We propose a computational model for understanding the complex issues behind nuclear arms negotiations, the motivations of various states to enter a nuclear weapons program and the ways to diffuse crisis situations.

Estimating The Flight Path Of Moving Objects Based On Acceleration Data, Peter Revesz

School of Computing: Conference and Workshop Papers

Inertial navigation is the problem of estimating the flight path of a moving object based on only acceleration measurements. This paper describes and compares two approaches for inertial navigation. Both approaches estimate the flight path of the moving object using cubic spline interpolation, but they find the coefficients of the cubic spline pieces by different methods. The first approach uses a tridiagonal matrix, while the second approach uses recurrence equations. They also require different boundary conditions. While both approaches work in O(n) time where n is the number of given acceleration measurements, the recurrence equation-based method can be easier updated …

Taming The Software Chaos: True To Its Promise, Iftdss Eases The Burden Of Fuels Treatment Planning—And Does A Lot More Besides, Gail Wells

Joint Fire Science Program Digests

A key problem reported by the fuels treatment planning community is the difficulty and inefficiency of evaluating and then applying many planning tools and applications. Fuels specialists have struggled to find, load, and learn all the different fuels and fire planning models, not to mention the interface of running, adjusting, and inputting data specific to each model without the ability to easily share inputs/outputs between models.

The Interagency Fuels Treatment Decision Support System (IFTDSS) was conceived as a way for users to learn one interface, access a variety of data and models all in one place, and pass data (inputs …

Across-Speaker Articulatory Normalization For Speaker-Independent Silent Speech Recognition, Jun Wang, Ashok Samal, Jordan Green

School of Computing: Conference and Workshop Papers

Silent speech interfaces (SSIs), which recognize speech from articulatory information (i.e., without using audio information), have the potential to enable persons with laryngectomy or a neurological disease to produce synthesized speech with a natural sounding voice using their tongue and lips. Current approaches to SSIs have largely relied on speaker-dependent recognition models to minimize the negative effects of talker variation on recognition accuracy. Speaker-independent approaches are needed to reduce the large amount of training data required from each user; only limited articulatory samples are often available for persons with moderate to severe speech impairments, due to the logistic difficulty of …

A Comparative Study Of Underwater Robot Path Planning Algorithms For Adaptive Sampling In A Network Of Sensors, Sreeja Banerjee

School of Computing: Dissertations, Theses, and Student Research

Monitoring lakes, rivers, and oceans is critical to improving our understanding of complex large-scale ecosystems. We introduce a method of underwater monitoring using semi-mobile underwater sensor networks and mobile underwater robots in this thesis. The underwater robots can move freely in all dimension while the sensor nodes are anchored to the bottom of the water column and can move only up and down along the depth of the water column. We develop three different algorithms to optimize the path of the underwater robot and the positions of the sensors to improve the overall quality of sensing of an area of …

Improving Preference Recommendation And Customization In Real World Highly Configurable Software Systems, Dongpu Jin

School of Computing: Dissertations, Theses, and Student Research

Highly configurable software systems, such as web browsers or office applications, may have a large number of preferences that the user can customize. When faced with the task of trying to identify which configuration option should be modified to change a particular system behavior, the user, tester or debugger may have to search through hundreds or thou- sands of options, and documentation may be scarce. Simple pattern matching utilities exist, but these searches are sensitive to using the right keyword. Static analysis may help, but will require access to source code. Alternatively a user may ask questions on help forums, …

Ultimate Codes: Near-Optimal Mds Array Codes For Raid-6, Zhijie Huang, Hong Jiang, Chong Wang, Ke Zhou, Yuhong Zhao

School of Computing: Technical Reports

As modern storage systems have grown in size and complexity, RAID-6 is poised to replace RAID-5 as the dominant form of RAID architectures due to its ability to protect against double disk failures. Many excellent erasure codes specially designed for RAID-6 have emerged in recent years. However, all of them have limitations. In this paper, we present a class of near perfect erasure codes for RAID-6, called the Ultimate codes. These codes encode, update and decode either optimally or nearly optimally, regardless of what the code length is. This implies that utilizing these codes we can build highly efficient and …

S-Code: Lowest Density Mds Array Codes For Raid-6, Zhijie Huang, Hong Jiang, Ke Zhou, Yuhong Zhao, Chong Wang

School of Computing: Technical Reports

RAID, a storage architecture designed to exploit I/O parallelism and provide data reliability, has been deployed widely in computing systems as a storage building block. In large scale storage systems, in particular, RAID-6 is gradually replacing RAID-5 as the dominant form of disk arrays due to its capability of tolerating concurrent failures of any two disks. MDS (maximum distance separable) array codes are the most popular erasure codes that can be used for implementing RAID-6, since they enable optimal storage efficiency and efficient encoding and decoding algorithms. In this paper, we propose a new class of MDS array codes called …

Preliminary Test Of A Real-Time, Interactive Silent Speech Interface Based On Electromagnetic Articulograph, Jun Wang, Ashok Samal, Jordan R. Green

School of Computing: Conference and Workshop Papers

A silent speech interface (SSI) maps articulatory movement data to speech output. Although still in experimental stages, silent speech interfaces hold significant potential for facilitating oral communication in persons after laryngectomy or with other severe voice impairments. Despite the recent efforts on silent speech recognition algorithm development using offline data analysis, online test of SSIs have rarely been conducted. In this paper, we present a preliminary, online test of a real-time, interactive SSI based on electromagnetic motion tracking. The SSI played back synthesized speech sounds in response to the user’s tongue and lip movements. Three English talkers participated in this …

Simexplorer: A Testing Framework To Detect Elusive Software Faults, Tingting Yu

School of Computing: Dissertations, Theses, and Student Research

Modern computer software systems are prone to various classes of runtime faults due to their reliance on features such as concurrency and peripheral devices such as sensors. Testing remains a common method for uncovering faults in these systems. However, commonly used testing techniques that execute the program with test inputs and inspect program outputs to detect failures are often ineffective on these systems. In addition, existing testing techniques focus primarily on single applications, neglecting elusive faults that occur at the whole-system level due to complex system event interactions.

This research provides a framework, SIMEXPLORER, that allows engineers to effectively test …

Power Management In The Cluster System, Leping Wang

School of Computing: Dissertations, Theses, and Student Research

With growing cost of electricity, the power management of server clusters has become an important problem. However, most previous researchers have only addressed the challenge in traditional homogeneous environments. Considering the increasing popularity of heterogeneous and virtualized systems, this thesis develops a series of efficient algorithms respectively for power management of heterogeneous soft real-time clusters and a virtualized cluster system. It is built on simple but effective mathematical models. When deployed to a new platform, the software incurs low configuration cost because no extensive performance measurements and profiling are required. Built upon optimization, queuing theory and control theory techniques, our …

Changes On Cran, Kurt Hornik, Achim Zeileis

The R Journal

New packages in CRAN task views

New contributed packages

Other changes

The R Journal (June 2014) 6(1): Complete Issue, The R Foundation

The R Journal

Editorial, Deepayan Sarkar

Contributed Research Articles

Taming PITCHf/x Data with XML2R and pitchRx, Carson Sievert

A Multiscale Test of Spatial Stationarity for Textured Images in R, Matthew A. Nunes, Sarah L. Taylor, and Idris A. Eckley

Stratified Weibull Regression Model for Interval-Censored Data, Xiangdong Gu, David Shapiro, Michael D. Hughes, and Raji Balasubramanian

brainR: Interactive 3 and 4D Images of High Resolution Neuroimage Data, John Muschelli, Elizabeth Sweeney, and Ciprian Crainiceanu

The RWiener Package: an R Package Providing Distribution Functions for the Wiener Diffusion Model, Dominik Wabersich, and Joachim Vandekerckhove

PivotalR: A Package for Machine Learning on Big Data, Hai …

R Foundation News, Kurt Hornik

The R Journal

Donations and new members

• Donations
• New benefactors
• New supporting institutions
• New supporting members

Web Technologies Task View, Patrick Mair, Scott Chamberlain

The R Journal

This article presents the CRAN Task View on Web Technologies. We describe the most important aspects of Web Technologies and Web Scraping and list some of the packages that are currently available on CRAN. Finally, we plot the network of Web Technology related package dependencies.

Oligomask: A Framework For Assessing And Removing The Effect Of Genetic Variants On Microarray Probes, Daniel Bottomly, Beth Wilmot, Shannon K. Mcweeney

The R Journal

As expression microarrays are typically designed relative to a reference genome, any individual genetic variant that overlaps a probe’s genomic position can possibly cause a reduction in hybridization due to the probe no longer being a perfect match to a given sample’s mRNA at that locus. If the samples or groups used in a microarray study differ in terms of genetic variants, the results of the microarray experiment can be negatively impacted. The oligoMask package is an R/SQLite framework which can utilize publicly available genetic variants and works in conjunction with the oligo package to read in the expression data …

Archiving Reproducible Research With R And Dataverse, Thomas J. Leeper

The R Journal

Reproducible research and data archiving are increasingly important issues in research involving statistical analyses of quantitative data. This article introduces the dvn package, which allows R users to publicly archive datasets, analysis files, codebooks, and associated metadata in Dataverse Network online repositories, an open-source data archiving project sponsored by Harvard University. In this article I review the importance of data archiving in the context of reproducible research, introduce the Dataverse Network, explain the implementation of the dvn package, and provide example code for archiving and releasing data using the package

Investr: An R Package For Inverse Estimation, Brandon M. Greenwell, Christine M. Schubert Kabban

The R Journal

Inverse estimation is a classical and well-known problem in regression. In simple terms, it involves the use of an observed value of the response to make inference on the corresponding unknown value of the explanatory variable. To our knowledge, however, statistical software is somewhat lacking the capabilities for analyzing these types of problems. In this paper, we introduce investr (which stands for inverse estimation in R), a package for solving inverse estimation problems in both linear and nonlinear regression models.1

Rotations: An R Package For So(3) Data, Bryan Stanfill, Heike Hofmann, Ulrike Genschel

The R Journal

In this article we introduce the rotations package which provides users with the ability to simulate, analyze and visualize three-dimensional rotation data. More specifically it includes four commonly used distributions from which to simulate data, four estimators of the central orientation, six confidence region estimation procedures and two approaches to visualizing rotation data. All of these features are available for two different parameterizations of rotations: three-by-three matrices and quaternions. In addition, two datasets are included that illustrate the use of rotation data in practice