# Basic Source Code for 
# "Determining Growth Determinants: Default Priors and Predictive Performance in Bayesian Model Averaging" 
# Theo Eicher, Chris Papageorgiou and Adrian Raftery
# Journal of Applied Econometrics, 26, 1, (February 2011): 30-55
#
# Programming by Adrian Raftery, Drew Creal, Amanda Cox
# LPS program by Jennifer Hoeting
# CRPS program by Tilmann Gneiting
#
# Date 4/29/2007
#
# ****BEGIN SAMPLE SCRIPT**********

require(BMA)
require(foreign)
require(rrcov)
library(MASS)

source("c:/e/source_bmacompare.txt")
source("c:/e/source_summarycompare.txt")

mydata <- read.table(file="c:/e/FLS.txt", header=T)
y <- mydata[,1]
x <- mydata[,2:42]

resultprior2 <- bma.compare(x, y, OR = 1e+12, nbest = 3, maxCol = 31, compare = TRUE, performance = TRUE,
	performance.compare=TRUE, performance.only = FALSE, define.V = FALSE, g.prior = 2,
		 delta = 0.15, gam = 0.65)

summary.compare(resultprior2, y, x,print.HPDI=FALSE,digits=4)

# ****END SAMPLE SCRIPT**********