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)
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 alexandresanto...@yahoo.com.br
Caros Membros,
Estou tentando converter dados de contagem de insetos
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.
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
Eu não sei exatamente porque você precisa dos dados dessa maneira, mas se é
para correr um glm() binomial, basta que você tenha um vetor com o número
de mortos e outro com o número de vivos. Não precisar estar em binário
{0,1} não.
tot - 212
da - expand.grid(trat=gl(3,4), tempo=1:5)
da$y -
Caros membros,
Vou usar da lista para possivelmente redimir uma dúvida conceitual,
pois não vejo melhor local para fazê-lo, dado ao perfil dos usuários.
Muitos autores classificam os modelos de regressão linear assim:
1- Simples (y = Bo + B1X)
2- Polinomial (y = Bo + B1X + B2X^2 + + B3X^3...)
3-
Por sinal, pegando carona na última colocação do Walmes, acho que a
abordagem de análise de sobrevivência (que neste caso sobrevivência passa
a ser o antônimo de germinação, mas vamos pensar que ainda assim seria uma
metáfora), seria mais intuitiva para o tempo de germinação.
Meus 0,019
Walmes,
Solicitei desta maneira porque eu achava que em análises de
sobrevivência, como no meu caso que utilizo o modelo de Weibull com a
função survreg() do pacote survival, fosse necessário informar o número
de indivíduos vivos (zeros) ao final do tempo experimental para se
efetuar