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
_______________________________________________
R-br mailing list
[email protected]
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código
mínimo reproduzível.