huge (version 1.0.1)

huge.NPN: NonparaNormal(NPN) transformation

Description

Implements the NonparaNormal transformation (Gausianization) to help relax the assumption of normality.

Usage

huge.NPN(x, NPN.func = "shrinkage", NPN.thresh, verbose = TRUE)

Arguments

x

The n by d data matrix representing n observations in d dimensions

NPN.func

The transformation function used in the NPN transformation. If NPN.func = "truncation", the truncated ECDF is applied. If NPN.func = "shrinkage", the shrunken ECDF is applied. The default is "shrinkage".

NPN.thresh

The truncation threshold used in NPN transformation, ONLY applicable when NPN.func = "truncation". The default value is 1/(4*(n^0.25)* sqrt(pi*log(n))).

verbose

If verbose = FALSE, tracing information printing is disabled. The default value is TRUE.

Value

• An object with S3 class "NPN" is returned:
• data
  The n by d data matrix representing n observations in d transformed dimensions
• ntdata
  The original data matrix before the NPN transformation
• NPN.func
  The transformation function used in the NPN transformation

Details

The NPN is a very helpful tool to relax the normality assumption. It assumes the data come from a Gaussian cornopula and computationally, fitting a high dimensional NPN is no more difficult than estimating a multivariate Gaussian and one can also apply to other existing estimation algorithms. The output data are already standardized as sample mean zero and unit variance.

References

1.Tuo Zhao and Han Liu. HUGE: A Package for High-dimensional Undirected Graph Estimation. Technical Report, Carnegie Mellon University, 2010 2.Han Liu, John Lafferty and Larry Wasserman. The Nonparanormal: Semiparametric Estimation of High Dimensional Undirected Graphs. Journal of Machine Learning Research (JMLR), 2009

See Also

huge and huge-package.

Examples

# generate data
L = huge.generator(graph = "cluster", g = 5)

# transform the non-Gaussian data using the shrunken ECDF
Q = huge.NPN(L$data^5)
summary(Q)
plot(Q)

# transform the non-Gaussian data using the truncated ECDF
Q = huge.NPN(5^(L$data), NPN.func = "truncation")
plot(Q)