added simulation study R scripts