fit <-glm(y ~ log( x_i +1) , family = binomial(link="logit"),
          data = gdata, maxit=50, trace=TRUE, na.action =  "na.omit")

for (k in 1:NumberOfDays){
   runis     = runif(10,0,1)
   ytest[,k] = ifelse(runis < proba[k],1,0)
   }