Prezados listeiros Bom dia. Eu tenho dados de produtividade de grãos em milho como variável resposta e dados de estande (número de plantas por parcela) como outra variável. Como nós tivemos problemas com irregularidade de estande, eu quero usar o número de plantas como covariável no no meu modelo. Abaixo o modelo: exp01.anova.2<-lm(gy ~ prog+rep+bloc%in%rep, data = exp01[exp01$local==2,]) Como eu coloco essa variável como covariável? O número de plantas ideal é 25 plantas por parcela, mas tenho parcelas com 16, 18 etc. Como coloco isso no aquivo? Eu devo subtrair o número de plantas ideal do real e, aquelas parcelas que tem o número ideal ficariam com zero, aquelas com mais de 25 plantas com número positivo e aquelas com menos de 25 plantas com número negativo? Isso mesmo? Ele vai para o modelo como fator ou numero? Muito obrigado. Att, Rodrigo. -- Prof. Rodrigo Oliveira de Lima Federal University of Viçosa Department of Plant Sciences Viçosa, MG, Brazil, 36570-900, Phone: +55(31)3899-2781 Fax: +55(31)3899-2614 Email: rodrigoodelima@ufv.br CV: http://lattes.cnpq.br/2761299616351806
Olá Rodrigo, O número de plantas não deve ser covariável, mas algo como offset, Ntrials ou exposure. Uma opção ingênua seria trabalhar com produtividade por planta. Porém, produtividade por planta igual em parcelas com número diferente de plantas provavelmente têm variância diferente. Elias
Complementando a sugestão do Elias, Você pode centrar o stand em 25 o no valor médio do experimento. Usar ela como covariável na lm(). Seria então uma análise de covariância. Talvez sirva de inspiração esse documento que analisa peso de leitões tendo como covariáveis peso inicial e idade inicial. Lembre-se de um detalhe, para uma análise de covariância supõe-se que as covariáveis não dependam dos fatores estudados. Eu suponho que no seu modelo "prog" seja progênie e pode ser que o stand dependa de progênie porque umas germinam com maior taxa ou têm menor caimento por vento, etc. Se o stand é função da cultivar e você usa o stand para corrigir a produtividade, o que você está fazendo em outras palavras e comparando elas numa situação de stand padronizado. A pergunta é, será que isso deve ser feito? No exemplo dos leitões, peso e idade foram medidos antes do experimento começar, antes mesmo de casualizar os animais, e por isso antecipamos que não haverá dependência dos tratamentos aplicações. À disposição. Walmes.
O link acabei esquecendo: http://www.leg.ufpr.br/~walmes/cursoR/cpaa/cpaa11anc.html. W.
participantes (3)
-
Elias Teixeira Krainski -
Rodrigo Lima -
Walmes Zeviani