Discrete dynamic Bayesian network analysis of fMRI data