plot.abe {abe}

Plot Function

Description

Plot function for the bootstrapped version of ABE.

Usage

## S3 method for class 'abe'
plot(x, type.plot = c("coefficients", "models", "variables"),
  alpha = NULL, tau = NULL, variable = NULL, ...)

Arguments

x

an object of class "abe", an object returned by a call to abe.boot

type.plot

string which specifies the type of the plot. See details.

alpha

values of alpha for which the plot is to be made (can be a vector of length >1)

tau

values of tau for which the plot is to be made (can be a vector of length >1)

variable

variables for which the plot is to be made (can be a vector of length >1)

...

Arguments to be passed to methods, such as graphical parameters (see barplot, hist).

Details

when using type.plot="coefficients" the function plots a histogram of the estimated regression coefficients for the specified variables, alpha(s) and tau(s) obtained from different re-sampled datasets. When the variable is not included in the final model, its regression coefficient is set to zero. When using type.plot="variables" the function plots a barplot of the relative inclusion frequencies of the specified variables, for the specified values of alpha and tau. When using type.plot="models" the function plots a barplot of the relative frequencies of the final models for specified alpha(s) and tau(s).

Author(s)

Rok Blagus, [e-mail]

Sladana Babic

See Also

abe.boot, summary.abe

Examples

set.seed(1)
n=100
x1<-runif(n)
x2<-runif(n)
x3<-runif(n)
y<--5+5*x1+5*x2+ rnorm(n,sd=5)
dd<-data.frame(y=y,x1=x1,x2=x2,x3=x3)
fit<-lm(y~x1+x2+x3,x=TRUE,y=TRUE,data=dd)

fit.boot<-abe.boot(fit,data=dd,include="x1",active="x2",
tau=c(0.05,0.1),exp.beta=FALSE,exact=TRUE,
criterion="alpha",alpha=c(0.2,0.05),type.test="Chisq",
num.boot=50,type.boot="bootstrap")

plot(fit.boot,type.plot="coefficients",
alpha=0.2,tau=0.1,variable=c("x1","x3"),
col="light blue")

plot(fit.boot,type.plot="variables",
alpha=0.2,tau=0.1,variable=c("x1","x2","x3"),
col="light blue",horiz=TRUE,las=1)

par(mar=c(4,6,4,2))
plot(fit.boot,type.plot="models",
alpha=0.2,tau=0.1,col="light blue",horiz=TRUE,las=1)