ภาพแบบไบนารี

# Create the vector of probabilities
p_cluster_1 <- c(0.8, 0.8, 0.2, 0.9)
# Create the sample for each pixel
pixel_1 <- sample(c(0, ___), ___, replace = TRUE, 
                  prob = c(1-p_cluster_1[___], p_cluster_1[___]))
pixel_2 <- sample(c(0, 1), 100, replace = TRUE, 
                  prob = c(1-p_cluster_1[___], p_cluster_1[___]))
pixel_3 <- sample(c(0, 1), 100, replace = TRUE, 
                  prob = c(1-p_cluster_1[___], p_cluster_1[___]))
pixel_4 <- sample(c(0, 1), 100, replace = TRUE, 
                  prob = c(1-p_cluster_1[___], p_cluster_1[___]))
# Combine the samples
sample_cluster_1 <- cbind(pixel_1, pixel_2, pixel_3, pixel_4)
# Have a look to the sample
head(___)