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

algo tipo isso? x = c(1,0,3,4,0,2,5,0)

x = ifelse (x>0,1,x) x [1] 1 0 1 1 0 1 1 0

Em 08/11/2013, à(s) 12:42, ASANTOS <alexandresantosbr@yahoo.com.br> 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,

-- ====================================================================== Alexandre dos Santos Proteção Florestal IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso Campus Cáceres Caixa Postal 244 Avenida dos Ramires, s/n Bairro: Distrito Industrial Cáceres - MT CEP: 78.200-000 Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos@cas.ifmt.edu.br Lattes: http://lattes.cnpq.br/1360403201088680 ======================================================================

_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br 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.