ENH load.besa:  add dat_file, dat_set and roi readers
ENH Factor:  .isin() and related methods
ENH Dataset: allow adding factor with slicing

>>> ds['name', :] = 'value'
ENH plot.UTSStat:  add xlim option
 WIP:  FIX stats:  Nested effects