Generating a raster form bernoulli distiributions

Question

I have a raster file of occurrence probabilities (P) and I want to generate a binary raster from this raster (Y). So at each pixel i I generate y_i ~ Binomia