Bootstrap code in dic.fit relies heavily on data frames and lists. Could probably be sped up at least mildly by using matrices