Diogo, a Lista tem regra de não anexar arquivos no email para evitar vírus, dentre outros problemas, poste seus arquivos por exemplo, no Datafile Host ou outro e indique o link.
[. ]'s. Edson Lira Estatístico Ma-Am Em 17/05/2013, às 10:54, "Diogo B. Provete" <[email protected]> escreveu: > Caros, > > Tenho um problema com a análise de um experimento sobre o efeito da > diversidade no funcionamento de ecossitemas aquáticos > > O experimento manipulou 3 níveis de riqueza de espécies e analisou o seu > efeito em 3 propriedades do funcionamento do ecossistema (productivity, > decomposition, respiration). > Para evitar o confundimento entre o efeito da riqueza e composição de > espécies (efeito de amostragem), a autora aninhou 7 diferentes composições de > espécie em cada tratamento. Cada uma dessas composições de espécies foi > replicada 2 x (as linhas de poças na fig), tomando aleatoriamente spp de um > pool total de 21 espécies. > Meu objetivo é usar a diversidade filogenética (PD) como um preditor de cada > uma das propriedades do ecossistema. Como a PD é uma medida contínua que diz > respeito à composição de espécies e suas relações evolutivas, cada uma dessas > 7 composições ≠ terá um PD diferente. O problema é que pra cada tratamento eu > terei 7 PD's ≠ aninhados. Além disso, ela tomou 6 medidas ao longo do tempo > de cada uma das propriedades. > > Eu tentei rodar um modelo misto especificando o desenho aninhado para o > efeito aleatório (PD e composição de spp) no lme4, mas tá dando erro (veja o > m1 no script). > > > dados <- read.csv2("dados1.csv", header = T) > attach(dados) > m1 <- lmer(prod ~ factor(div)*factor(time) + factor(div)+ factor(time) + > (1|comp) + (factor(PD)/factor(tank)/factor(time)), data = dados) > > Então, o gde problema é conseguir fazer com que o efeito aleatório seja > aninhado e ainda adicionar as medidas repetidas, ou seja, indicar que foram > feitas várias medidas numa mesma unidade amostral (tank, na tabela abaixo) e > que o PD é um efeito aleatório, junto com a composição, e está aninhado no > efeito fixo (=tratamento) que são os 3 níveis de riqueza. > > Meus dados são assim: > tank treat comp PD div time prod > 1 18 1 1 2141.15 1 1 0.259 > 2 60 1 1 2141.15 1 1 0.269 > 3 20 4 1 2141.15 1 1 0.299 > 4 58 4 1 2141.15 1 1 0.294 > 5 18 1 1 2141.15 1 2 0.244 > 6 60 1 1 2141.15 1 2 0.295 > > > Envio em anexo os dados completos e tb o script do R, caso possam me ajudar. > O artigo do qual eu extrai esses dados é este. > > Obrigado desde já > > -- > Atenciosamente, > Diogo Borges Provete > > ============================== > Biólogo > Mestre em Biologia Animal (UNESP) > Doutorando PPG Ecologia e Evolução > Laboratório de Ecologia de Insetos (sl. 222) > Departamento de Ecologia > Instituto de Ciências Biológicas - ICB 1 > Universidade Federal de Goiás, campus II > Goiânia-GO > 74001-970 > Brazil > > Tel. Lab. +55 62 3521-1732 > Cel. +55 62 8231-5775 > > : diogoprovete > > Personal web page > > Editor Herpetology Notes > NorthWestern Journal of Zoology > ============================== > > > > > <dados.R> > <dados1.csv> > _______________________________________________ > R-br mailing list > [email protected] > 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 [email protected] 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.
