CompareMRF_T2: Comparison of two MRFs using Hotellings T2

In Rnets: Resistance Relationship Networks using Graphical LASSO

Description

Hotelling's T2 is the multivariate extension of Welch's t-test used to compare two sets of means. The test is used here to test for differences in the penalized partial correlations in two diffrent MRFs.

Usage

CompareMRF_T2(X, Y, L1, vertices = NULL, subset = NULL, B = 500,
  n_b = 0.5, replace = T, seed = NULL)

Arguments

X	A dataset to be used to estimate an MRF/Rnet
Y	A dataset to be used to estimate an MRF/Rnet
L1	The regularization pentalty to use for both networks
vertices	Vertices to include in MRF. Must be a subset of names(x)
subset	An expression for the subsetting X and Y based on columns not used in the MRF. Appled to both data sets
B	The number of subsamples to draw. The same number of subsamples is drawn from both data sets.
n_b	The size of the subsamples. May be provided as an integer less than the number of rows in x, or as a proportion.
replace	Logical. Is subsampling done with (T) or without (F) replacement.
seed	Random seed value for reproducibility. The same seed is applied to sampling from both data sets.

Value

An S3 object of class 'mrf_t2' containing crude test results, adjusted test results, and a table of pair-wise comparisons of unique edges

Rnets documentation built on July 23, 2019, 9:04 a.m.