09/09 fixed gamma0 model to handle cases where there are no 0 forecasts

06/09 modified estimation in gamma0 and gamma to do true EM 
      updated bias correction strategy to handle 0 forc in gamma and gamma0 

05/09 fixed crps calculation for gamma0 and gamma models
      changed call to cdfBMAgamma0 to cdfBMAgamma in cdf.ensembleBMAgamma

03/09 changed nonunique forecast error exit to warning

01/09 added subsetting (via date) for fitted BMA models
      changed ensembleDates to ensembleValidDates
      added forecastHour and initializationTime attrributes to model fits  
      changed crps to CRPS and added modeified crps; changed MAE to mae
      added verifRankHist function

12/08 modified ensembleData function for efficiency    
      added forecastHour requirement for ensembleData
      trainingRule changed to trainingDays
      added wind speed data modeling
      changed eps to tol in control functions

11/08 added pit function

10/08 fixed brierScore so that it works with one threshold
      fixed help file syntax