Description Usage Arguments Value Examples
Produce a plot of the stability path obtained by stability selection.
1 2 3 4 
x 
output of a 
y 
used for S4 compatibility. 
xvar 
variable to plot on the Xaxis: either 
annot 
logical; should annotation be made on the graph
regarding controlled PFER (only relevant when 
main 
main title. If none given, a somewhat appropriate title is automatically generated. 
log.scale 
logical; indicates if a logscale should be used
when 
labels 
an optional vector of labels for each variable in the path (e.g., 'relevant'/'irrelevant'). See examples. 
plot 
logical; indicates if the graph should be
plotted. Default is 
sel.mode 
a character string, either 
cutoff 
value of the cutoff probability (only relevant when

PFER 
value of the perfamily error rate to control (only
relevant when 
nvar 
number of variables selected (only relevant when

... 
used for S4 compatibility. 
a list with a ggplot2 object which can be plotted
via the print
method, and a vector of selected variables
corresponding to method of choice ('rank'
or
'PFER'
)
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25  ## Not run:
## Simulating multivariate Gaussian with blockwise correlation
## and piecewise constant vector of parameters
beta < rep(c(0,1,0,1,0), c(25,10,25,10,25))
Soo < matrix(0.75,25,25) ## bloc correlation between zero variables
Sww < matrix(0.75,10,10) ## bloc correlation between active variables
Sigma < bdiag(Soo,Sww,Soo,Sww,Soo) + 0.2
diag(Sigma) < 1
n < 100
x < as.matrix(matrix(rnorm(95*n),n,95) %*% chol(Sigma))
y < 10 + x %*% beta + rnorm(n,0,10)
## Build a vector of label for true nonzeros
labels < rep("irrelevant", length(beta))
labels[beta != 0] < c("relevant")
labels < factor(labels, ordered=TRUE, levels=c("relevant","irrelevant"))
## Call to stability selection function, 200 subsampling
stab < stability(x,y, subsamples=200, lambda2=1, min.ratio=1e2)
## Build the plot an recover the selected variable
plot(stab, labels=labels)
plot(stab, xvar="fraction", labels=labels, sel.mode="PFER", cutoff=0.75, PFER=2)
## End(Not run)

Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.