stablespec (version 0.1.3)

plotStability: Plot of edge and causal path stability.

Description

Plot each of the stability of causal path and edge including the threshold of stability and model complexity.

Usage

plotStability(listOfFronts = NULL, threshold = NULL, stableCausal = NULL, stableCausal_l1 = NULL, stableEdge = NULL, longitudinal = NULL)

Arguments

listOfFronts

list of models including their fitness and subset index, which is one of the outputs of stableSpec.

threshold

threshold of stability selection. The default is 0.6.

stableCausal

list of the stability of causal path for the whole range of model complexities, which is one of the outputs of stableSpec.

stableCausal_l1

list of the stability of causal path of length 1 for the whole range of model complexities, which is one of the outputs of stableSpec.

stableEdge

list of the stability of edge for the whole range of model complexities, which is one of the outputs of stableSpec.

longitudinal

TRUE for longitudinal data, and FALSE cross-sectional data.

Value

Plot of causal path and edge stability for every pair of variables, and a plot of altogether combined including the corresponding stability and complexity thresholds.

Examples

the_data <- adhd
numSubset <- 1
num_iteration <- 5
num_pop <- 10
mut_rate <- 0.075
cross_rate <- 0.85
longi <- FALSE
num_time <- 1
the_co <- "covariance"
# assummed that nothing causing variable Gender
cons_matrix <- matrix(c(2, 1, 3, 1, 4, 1, 5, 1, 6, 1), 5, 2, byrow=TRUE)
th <- 0.1
to_plot <- FALSE

result_adhd <- stableSpec(theData=the_data, nSubset=numSubset,
iteration=num_iteration,
nPop=num_pop, mutRate=mut_rate, crossRate=cross_rate,
longitudinal=longi, numTime=num_time,
co=the_co, consMatrix=cons_matrix, threshold=th, toPlot=to_plot)

plotStability(listOfFronts=result_adhd$listOfFronts, threshold=th,
stableCausal=result_adhd$causalStab,
stableCausal_l1=result_adhd$causalStab_l1,
stableEdge=result_adhd$edgeStab,
longitudinal=longi)