
8 Nov
2013
8 Nov
'13
14:52
Em sexta-feira, 8 de novembro de 2013 11:42:22, ASANTOS escreveu:
Boa tarde pessoal,
No seguinte exemplo:
da <- expand.grid(trat=gl(2,1), tempo=1:10) X <- model.matrix(~trat+tempo, da); ncol(X) betas <- c(0.1,0.1,0.3) eta <- X%*%betas y1 <- rpois(da$trat, lambda=exp(eta)) y2 <- rbinom(y1, size=1, prob=0.7) da$y <- y1*y2 head(da)
Gostaria de criar um novo vetor ybinom, através da substituição de valores em da$y >0 por 1, e os zeros continuam zeros mesmo. Alguém poderia me dar um help?
Obrigado,
Você pode fazer da seguinte forma: da$ybinom<-ifelse(da$y>0,1,0) da