[R-br] Converter dados de contagem em binomial

[ASANTOS]

Obrigado Rodrigo,

     Na verdade a concatenação (c()) foi um deslize, mas eu queria fazer:

dados.bin <- dados[rep(1:nrow(dados), times=dados$y), -3]

censor<-rep(1,length(dados.bin[,1]))

dados.final<- cbind(dados.bin,censor)

Só que tenho n.indv<-100 - Número de indivíduos por tratamentos - e 
agora preciso adicionar uma repetição de zeros
no objeto censor equivalente a n.indv que restaram vivos (n.indv - soma 
de mortos por tratamento) , mas de maneira
mais automatizada que a que venho fazendo, teria alguma sugestão?

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






Em 10/12/2013 09:43, Rodrigo Coster escreveu:
> Pq tu ta usando c() e não rbind()... c() transforma tudo em vetor.
>
> Pelo que eu entendi, o que tu quer é isso:
>
> dados.bin <- dados[rep(1:nrow(dados), times=dados$y), -3]
>
>
>
> 2013/12/10 ASANTOS <alexandresantosbr em yahoo.com.br 
> <mailto:alexandresantosbr em yahoo.com.br>>
>
>     Caros Membros,
>
>        Estou tentando converter dados de contagem de insetos mortos em
>     notação binomial e estou fazendo algo errado, pois se tenho:
>
>     trat<-gl(5,10) ## 5 tratamentos com dez repetições
>
>     tempo<- gl(10,50)  ## 10 avaliações no tempo
>
>     dados<-expand.grid(trat, tempo)
>
>     colnames(dados)<-c("trat","tempo")
>
>     y<-rpois(length(dados[,1]),10) ## Simula a variável resposta
>
>     dados<-cbind(dados,y)
>
>     ## Converter os dados de contagem em binomial, por ex. se tenho 5,
>     devo converter em 1 1 1 1 1
>
>     dados.bin <- NULL
>     for (i in 1:nrow(dados))
>     dados.bin  <- c(dados.bin,
>     rep(c(1,dados[i,3])))
>
>     Por que para cada valor de contagem na coluna dados[,3] não tenho
>     a transformação que espero?
>
>     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 <tel:%28%2B55%29%2065%208132-8112> (TIM)
>     (+55) 65 9686-6970 <tel:%28%2B55%29%2065%209686-6970> (VIVO)
>     e-mails:alexandresantosbr em yahoo.com.br
>     <mailto:e-mails%3Aalexandresantosbr em yahoo.com.br>
>     alexandre.santos em cas.ifmt.edu.br
>     <mailto:alexandre.santos em cas.ifmt.edu.br>
>     Lattes: http://lattes.cnpq.br/1360403201088680
>     ======================================================================
>
>     _______________________________________________
>     R-br mailing list
>     R-br em listas.c3sl.ufpr.br <mailto: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.
>
>
>
>
> _______________________________________________
> 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.

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