Additional changes to fitBMA()

I did the following
- Added libraries necessary to run parallel
- Changed the default parallel argument in the function from FALSE to
TRUE (because in most cases people using this function will need to run
it in parallel)
- Added additional argument that specifies the number of cores to run
parallel in (I’ve set the default as 10)