Walmes, Perfeito, mas agora preciso achar um jeito para registrar os indivíduos que permanecem vivos no último tempo, através da declaração do número de indivíduos (tot) por tratamento, sendo: dados <- expand.grid(trat=gl(3,4), tempo=1:5) dados$y <- rpois(nrow(dados), lambda=10) dados$z <- 1 dados2 <- dados[rep(1:nrow(dados), dados$y),] xtabs(y~trat+tempo, data=dados) xtabs(z~trat+tempo, data=dados2) tot.dead<-with(dados2, tapply(z, list(trat), sum, na.rm = T)) tot<-212 ### Número de insetos por tratamento dif<-tot-tot.dead dif # 1 2 3 10 36 9 Essa diferença (dif) precisa ser representada por zeros no max(tempo) e fundido a variável z no objeto dados2, com isso consigo criar a variável reposta para realizar análises de sobrevivência, como eu poderia fazer isto? Tenho tentado criar um rep(0, dif[,i]) no max(dados2$tempo), mas sem sucesso, Obrigado, Em 10/12/2013 13:33, walmes . escreveu:
Acredito que o jeito mais simples seja assim
dados <- expand.grid(trat=gl(3,4), tempo=1:5) dados$y <- rpois(nrow(dados), lambda=10) dados$z <- 1
dados2 <- dados[rep(1:nrow(dados), dados$y),]
xtabs(y~trat+tempo, data=dados) xtabs(z~trat+tempo, data=dados2)
À disposição. Walmes.
_______________________________________________ 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.
-- ====================================================================== 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 ======================================================================