Plot histogram of AIC values of top MAR models

Plots a histogram of the AIC values of the best-fit models contained in the $top.bestfit component of an object of class MAR

Usage

# S3 method for MARtop
hist(x, ...)

Arguments

x: The $top.bestfit component of an object of class MAR
...: Further arguments passed to or from other methods

Details

The AIC values of the top best-fit models are plotted as a frequency histogram. A blue asterisk indicates the AIC value of the selected best-fit model.

Value

A graphics device containing the histogram of top best-fit model AIC values

Author

LP Scheef

Examples

data(L4.mar)
myvar <- c(0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 2)
myres <- matrix(0.5,
  nrow = length(which(myvar == 1)),
  ncol = length(which(myvar != 0))
) 
small <- L4.mar[1:20,]
run1 <- run.mar(small, variables = myvar, restrictions = myres, search = "random", boot=100)
#>  
#>  
#> searching for best-fit model...
#>        ...BEST-FIT MODEL SELECTED
#> 	           ( search time:  0 minutes 5 seconds )
#>  
#>  
#>  
#> identifying 10 lowest AIC models...
#>        ...TOP MODELS RETAINED 
#>  
#>  
#>  
#> bootstrapping best-fit model...
#>        ...BOOTSTRAPPING COMPLETE
#>  
#>  
#> ════════════════════════════════════════════════════════════════════════════════ 
#>  
#>  

hist(run1$top.bestfit)