itemConfirm: Replicability of Items from EGA

Description

Based on the bootEGA results, this function computes and plots the number of times an item (variable) is estimated in the same factor/dimension as originaly estimated by EGA. The output also each item's likelihood (i.e., proprotion of bootstraps that it appeared in) for every possible dimension (defined by the maximum number of dimensions across the bootstrap samples).

Usage

itemConfirm(bootega.obj, confirm, item.rep = 0.1, plot.ic = TRUE)

Arguments

bootega.obj

A bootEGA object

confirm

A vector with community numbers or labels for each item

item.rep

A value for lowest likelihood allowed in item.likelihood output. Removes noise from table to allow for easier interpretation. Defaults to .10

plot.ic

Should the plot be produced? Defaults to TRUE

Value

Returns a list containing:

item.confirm

The proporton of times each item replicated within the defined dimension

item.likelihood

The proportion of times each item replicated within each possible dimension. Dimensions greater than the maximum number used in the confirm argument are labeled based on the largest remaining components after the dimensions used to confirm.

A matrix containing the community membership values for each bootstrapped sample. The values correspond to the values input for the confrim argument.

References

Danon, L., Diaz-Guilera, A., Duch, J., & Arenas, A. (2005). Comparing community structure identification. Journal of Statistical Mechanics: Theory and Experiment, 9, P09008. <doi:10.1088/1742-5468/2005/09/P09008>

Examples

Run this code

# NOT RUN {
ega.wmt <- EGA(data = wmt2[,7:24], model = "glasso")

boot.wmt <- bootEGA(data = wmt2[,7:24], n = 100, typicalStructure = TRUE,
plot.typicalStructure = TRUE, model = "glasso",
type = "parametric", ncores = 4)

itemConfirm(boot.wmt, confirm = ega.wmt$wc)
# }
# NOT RUN {
# }