start writing script to perform exact model selection for hierarchical models