- #probas = rowSums(p * probas)
- X = rbind(X, newXblock)
- #Y = c( Y, vapply(probas, function(p) (ifelse(p >= .5, 1, 0)), 1) )
- Y = c( Y, vapply(probas, function(p) (rbinom(1,1,p)), 1) )
+ X <- rbind(X, newXblock)
+ Y <- c( Y, vapply(probas, function(p) (rbinom(1,1,p)), 1) )