Numerical Analysis and Scientific Computing Commons^™

Editorial, Hadley Wickham

The R Journal

Welcome to volume 5, issue 2 of The R Journal. I’m very pleased to include 21 articles about R for your enjoyment.

The end of the year also brings changes to the editorial board. Martyn Plummer is leaving the board after four years. Martyn was responsible for writing up the standard operating procedures for the journal, an act which has made my life as a new editor considerably easier! We welcome Michael Lawrence, who will join the editorial board in 2014. I am stepping down as Editor-in-Chief and will be leaving this task in the capable hands of Deepayan Sarkar.

Rnetcdf: A Package For Reading And Writing Netcdf Datasets, Pavel Michna, Milton Woods

The R Journal

This paper describes the RNetCDF package (version 1.6), an interface for reading and writing files in Unidata NetCDF format, and gives an introduction to the NetCDF file format. NetCDF is a machine independent binary file format which allows storage of different types of array based data, along with short metadata descriptions. The package presented here allows access to the most important functions of the NetCDF C-interface for reading, writing, and modifying NetCDF datasets. In this paper, we present a short overview on the NetCDF file format and show usage examples of the package.

Betategarch: Simulation, Estimation And Forecasting Of Beta-Skew-T-Egarch Models, Genaro Sucarrat

The R Journal

This paper illustrates the usage of the betategarch package, a package for the simulation, estimation and forecasting of Beta-Skew-t-EGARCH models. The Beta-Skew-t-EGARCH model is a dynamic model of the scale or volatility of financial returns. The model is characterised by its robustness to jumps or outliers, and by its exponential specification of volatility. The latter enables richer dynamics, since parameters need not be restricted to be positive to ensure positivity of volatility. In addition, the model also allows for heavy tails and skewness in the conditional return (i.e. scaled return), and for leverage and a time-varying long-term component in the …

Lfe: Linear Group Fixed Effects, Simen Gaure

The R Journal

Linear models with fixed effects and many dummy variables are common in some fields. Such models are straightforward to estimate unless the factors have too many levels. The R package lfe solves this problem by implementing a generalization of the within transformation to multiple factors, tailored for large problems.

Changes To Grid For R 3.0.0, Paul Murrell

The R Journal

From R 3.0.0, there is a new recommended way to develop new grob classes in grid. In a nutshell, two new “hook” functions, makeContext() and makeContent() have been added to grid to provide an alternative to the existing hook functions preDrawDetails(), drawDetails(), and postDrawDetails(). There is also a new function called grid.force(). This article discusses why these changes have been made, provides a simple demonstration of the use of the new functions, and discusses some of the implications for packages that build on grid.

Exactcidiff: An R Package For Computing Exact Confidence Intervals For The Difference Of Two Proportions, Guogen Shan, Weizhen Wang

The R Journal

Comparing two proportions through the difference is a basic problem in statistics and has applications in many fields. More than twenty confidence intervals (Newcombe, 1998a,b) have been proposed. Most of them are approximate intervals with an asymptotic infimum coverage probability much less than the nominal level. In addition, large sample may be costly in practice. So exact optimal confidence intervals become critical for drawing valid statistical inference with accuracy and precision. Recently, Wang (2010, 2012) derived the exact smallest (optimal) one-sided 1 confidence intervals for the difference of two paired or independent proportions. His intervals, however, are computer-intensive by nature. …

Reconstructing Point Clouds Of Mid-Size Objects, Spencer Woodworth

Computer Science and Software Engineering

This project explores the use of an inexpensive 3D camera for the acquisition and reconstruction of mid-size objects. The disparity of objects between stereo image pairs are used to calculate depth and generate a depth map. The depth map is used to generate a point cloud representation of the object from a single view. Finally, point clouds are generated from several views of an object and then aligned and merged into a seamless 360-degree point cloud.

Decision Trees To Model The Impact Of Disruption And Recovery In Supply Chain Networks, Loganathan Ponnanbalam, L. Wenbin, Xiuju Fu, Xiaofeng Yin, Zhaoxia Wang, Rick S. M. Goh

Research Collection School Of Computing and Information Systems

Increase in the frequency of disruptions in the recent times and their impact have increased the attention in supply chain disruption management research. The objective of this paper is to understand as to how a disruption might affect the supply chain network - depending upon the network structure, the node that is disrupted, the disruption in production capacity of the disrupted node and the period of the disruption - via decision trees. To this end, we first developed a 5-tier agent-based supply chain model and then simulated it for various what-if disruptive scenarios for 3 different network structures (80 trials …

The Influence Of Online Word-Of-Mouth On Long Tail Formation, Bin Gu, Qian Tang, Andrew B. Whinston

Research Collection School Of Computing and Information Systems

The long tail phenomenon has been attributed to both supply side and demand side economies. While the cause on the supply side is well-known, research on the demand side has largely focused on the awareness effect of online information that helps consumers discover new and often niche products. This study expands the demand side factors by showing that online information also influences the long tail phenomenon through the informative effect, which affects consumers' evaluation of product quality. We examine the informative effect in the context of online WOM. Two sets of theories suggest opposite directions for the implication of the …

An Efficient Two Sample Capture-Recapture Model With High Recaptures, Danjuma Jibasen, Yusuf J. Adams

CBN Journal of Applied Statistics (JAS)

This paper proposed an efficient two sample capture-recapture model (Ma) with high recaptures and compared it with the existing models like the model of no factor effect (Mo), behavioral response model (Mb) and the Petersen model (Ms), using simulated data. We found that the Petersen model provides a better estimate of the population size when the observations follow a hypergeometric distribution and the population is overestimated when recapture is high. It was also found that the proposed model provides a better estimator of the population size than the existing ones when the recapture is high. This model is particularly useful …

Causal Relationship Between Stock Market Index And Exchange Rate: Evidence From Nigeria, Abdulrasheed Zubair

CBN Journal of Applied Statistics (JAS)

This paper uses Johansen’s cointegration to test for the possibility of cointegration and Granger-causality to estimate the causal relationship between stock market index and monetary indicators (exchange rate and M2) before and during the global financial crisis for Nigeria, using monthly data for the period 2001–2011. Results suggest absence of long-run relationship before and during the crisis. The Granger-causality tests show a uni-directional causality running from M2 to ASI before the crisis while during the period of the crisis there is absence of causality between the variables. This suggests that ASI show responsiveness to M2. Thus, absence of the direct …

Investigating Chaos In The Nigerian Asset And Resource Management (Arm) Discovery Fund, Ibiyinka A. Fuwape, Samuel T. Ogunjo

CBN Journal of Applied Statistics (JAS)

This paper investigates chaos in a Nigerian mutual fund, Asset and Resource Management Company Limited (ARM) for a period of eleven years. The existence of chaotic signals in the data was identified by the reconstruction of the phase space of the daily closing price of the fund and the delay time was quantified using mutual information function and the embedding dimension by the false nearest neighbours, where the values were identified to be 15 and 20 respectively. The presence of chaotic signals in the ARM data was further confirmed by the correlation dimension method which yielded a dimension of 2.2 …

Modeling The Nigerian Inflation Rates Using Periodogram And Fourier Series Analysis, Chukwuemeka O. Omekara,, Emmanuel J. Ekpenyong, Micheal P. Ekerete

CBN Journal of Applied Statistics (JAS)

This work considers the application of Periodogram and Fourier Series Analysis to model all-items monthly inflation rates in Nigeria from 2003 to 2011. The main objectives are to identify inflation cycles, fit a suitable model to the data and make forecasts of future values. To achieve these objectives, monthly all-items inflation rates for the period were obtained from the Central Bank of Nigeria (CBN) website. Periodogram and Fourier series methods of analysis are used to analyze the data. Based on the analysis, it was found that inflation cycle within the period was fifty one (51) months, which coincides with the …

Nigerian Stock Index: A Search For Optimal Garch Model Using High Frequency Data, Olaoluwa Simon Yaya

CBN Journal of Applied Statistics (JAS)

This paper attempts to fit the best Generalized Autoregressive Conditional Heteroscedastic (GARCH) model for All Share Index (ASI) of Nigerian Stock Exchange (NSE) returns. A search is made on various GARCH variants specified on the assumptions of stationarity and asymmetry. Fractionally integrated types are also considered to capture the possibility of return series having property of long range dependency. The parameter estimations are carried out on the assumptions of normality and non-normality of GARCH innovations, with models and forecasts evaluated using information criteria and loss functions respectively. Under normality assumption, Hyperbolic GARCH (HYGARCH(1,d,1)) model is selected and Integrated GARCH (IGARCH(1,1)) …

Time Series Modeling Of Nigeria External Reserves, Iheanyichukwu S. Iwueze, Eleazar C. Nwogu, Valentine U. Nlebedim

CBN Journal of Applied Statistics (JAS)

This paper discusses the levels and trend of external reserves in Nigeria. The relevance of this lies in the fact that it could help to monitor the reserves and throw early warning signal about any economic crisis. Monthly data on Nigeria external reserves for the period January 1999 to December, 2008 derived from the 2008 CBN Statistical Bulletin was analyzed using ARIMA model. Results of the analyses show that (i) the data requires logarithmic transformation to stabilize the variance and make the distribution normal (ii) the appropriate model that best describes the pattern in the transformed data is the Autoregressive- …

A Markov Decision Process Approach To Optimal Control Of A Multi-Level Hierarchical Manpower System, Akaninyene U. Udom

CBN Journal of Applied Statistics (JAS)

A recurrent problem in manpower control is how to attain the desired structural configuration in an optimal way, since it is possible to reach a desired structural configuration using different control inputs. The major aim of this paper is to develop a Markov Decision Process for optimal control of a Multi-level Hierarchical Manpower System (MHMS) by promotion and interdepartmental transfers. This is examined under control by intervention and contraction cost Markov Decision Process.

Modeling Preferences With Availability Constraints, Bingtian Dai, Hady W. Lauw

Research Collection School Of Computing and Information Systems

User preferences are commonly learned from historical data whereby users express preferences for items, e.g., through consumption of products or services. Most work assumes that a user is not constrained in their selection of items. This assumption does not take into account the availability constraint, whereby users could only access some items, but not others. For example, in subscription-based systems, we can observe only those historical preferences on subscribed (available) items. However, the objective is to predict preferences on unsubscribed (unavailable) items, which do not appear in the historical observations due to their (lack of) availability. To model preferences in …

Modeling Temporal Adoptions Using Dynamic Matrix Factorization, Freddy Chong-Tat Chua, Richard Jayadi Oentaryo, Ee Peng Lim

Research Collection School Of Computing and Information Systems

The problem of recommending items to users is relevant to many applications and the problem has often been solved using methods developed from Collaborative Filtering (CF). Collaborative Filtering model-based methods such as Matrix Factorization have been shown to produce good results for static rating-type data, but have not been applied to time-stamped item adoption data. In this paper, we adopted a Dynamic Matrix Factorization (DMF) technique to derive different temporal factorization models that can predict missing adoptions at different time steps in the users' adoption history. This DMF technique is an extension of the Non-negative Matrix Factorization (NMF) based on …

Topicsketch: Real-Time Bursty Topic Detection From Twitter, Wei Xie, Feida Zhu, Jing Jiang, Ee Peng Lim, Ke Wang

Research Collection School Of Computing and Information Systems

Twitter has become one of the largest platforms for users around the world to share anything happening around them with friends and beyond. A bursty topic in Twitter is one that triggers a surge of relevant tweets within a short time, which often reflects important events of mass interest. How to leverage Twitter for early detection of bursty topics has therefore become an important research problem with immense practical value. Despite the wealth of research work on topic modeling and analysis in Twitter, it remains a huge challenge to detect bursty topics in real-time. As existing methods can hardly scale …

Efficient Index-Based Approaches For Skyline Queries In Location-Based Applications, Ken C. K. Lee, Baihua Zheng, Cindy Chen, Chi-Yin Chow

Research Collection School Of Computing and Information Systems

Enriching many location-based applications, various new skyline queries are proposed and formulated based on the notion of locational dominance, which extends conventional one by taking objects' nearness to query positions into account additional to objects' nonspatial attributes. To answer a representative class of skyline queries for location-based applications efficiently, this paper presents two index-based approaches, namely, augmented R-tree and dominance diagram. Augmented R-tree extends R-tree by including aggregated nonspatial attributes in index nodes to enable dominance checks during index traversal. Dominance diagram is a solution-based approach, by which each object is associated with a precomputed nondominance scope wherein query points …