[R-br] Converter dados de contagem em binomial

[ASANTOS]

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 em 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 em yahoo.com.br
         alexandre.santos em cas.ifmt.edu.br
Lattes: http://lattes.cnpq.br/1360403201088680
======================================================================

-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20131210/21797552/attachment.html>