Fmri data simulation

Is there a way to simulate BOLD activation data in fmri?
Preferably, I would like to provide some roi, and an effect size, and then I could simulate activity based on the design.