[R-br] Sugestões sobre a melhor abordagens estatística

Alisson Lucrecio alissonluc em gmail.com
Segunda Junho 9 15:12:16 BRT 2014


Walmes,

Boa tarde.

Tenho o efeito fixo para tratamento com sete níveis, os quais são: com n
itrogénio; sem nitrogénio; (estirpe de rizóbioum) CIAT 899 ; UFLA 02-100;
UFLA 02-127; UFLA 02-68 e UFLA 04-195.

Tenho o efeito de Local com 10 níveis, os quais são: Lavras; Lavras2; Patos
de Minas; Patos de Minas2; PO; PO2; Bambui; GM; Luminarias e Pitangui.

Patos de Minas e Bambui foram irrigados e cultivados na estação Seca. Patos
de Minas2 e Pitangui foram irrigados e cultivados na estação Inverno.
Lavras, GM, Luminarias e PO2 foram cultivados sem irrigação e na estação
Águas. PO e Lavras2 foram cultivados sem irrigação e na estação Seca.

Existe também um efeito de nivel tecnológico onde Patos de Minas, Bambui,
Patos de Minas2 e Pitangui foi 4. Para Lavras, Lavras2, PO, Luminarias foi
2. Para GM foi 4. Para PO2 foi 3.

Além disso o Local GM foi cultivado com Plantio Direto.

Criei o modelo abaixo, mas testou com duvida em como inserir o efeito de
Irrigação, Nível Tecnológico, Safra e Preparo do Solo no Modelo.

lmer1 <- lmer(RG ~ Tratamento + ??? + (1 + Tratamento | Local ), data =
estirpe, control=lmerControl(optCtrl=list(maxfun=40000)))

Você poderia me ajudar?

Obrigado.

'data.frame':	280 obs. of  15 variables:
 $ Tratamento: Factor w/ 7 levels "T_CN","T_SN",..: 4 4 4 4 5 5 5 5 6 6 ...
 $ Local     : Factor w/ 10 levels "Lavras","Lavras2",..: 3 3 3 3 3 3
3 3 3 3 ...
 $ Bloco     : Factor w/ 4 levels "1","2","3","4": 1 2 3 4 1 2 3 4 1 2 ...
 $ Irrigacao : Factor w/ 2 levels "Sem","Com": 2 2 2 2 2 2 2 2 2 2 ...
 $ NT        : Factor w/ 3 levels "2","3","4": 3 3 3 3 3 3 3 3 3 3 ...
 $ Safra     : Factor w/ 3 levels "Seca","Aguas",..: 1 1 1 1 1 1 1 1 1 1 ...
 $ Prep_Solo : Factor w/ 2 levels "CO","PD": 1 1 1 1 1 1 1 1 1 1 ...

 $ MSPA      : num  21.5 19.2 27.1 21.4 11.1 ...
 $ NN        : int  255 136 136 186 129 79 84 130 242 106 ...
 $ MSN       : num  0.63 1.05 0.53 0.48 0.31 0.13 0.25 0.33 0.56 0.32 ...
 $ ANPA      : num  386 520 758 450 367 ...
 $ V.P       : num  13.8 9.9 6.1 9.8 3.7 5.1 9 7 7.2 7.4 ...
 $ S.V       : num  5.17 5.09 4.97 4.84 4.41 4.57 4.34 4.9 4 4.45 ...
 $ RG        : num  1362 1550 1925 1612 525 ...
 $ ANG       : num  24.5 41.9 53.9 33.9 17.3 ...




2014-05-26 14:22 GMT-03:00 walmes . <walmeszeviani em gmail.com>:

> Já antecipado que a escolha do rizóbium é dentro de cada local, poderia-se
> fazer uma análise por local com o efeito de estirpe, nitrogênio e
> tratamento (e as interações) no preditor (modelo de fatorial triplo). No
> entanto, uma análise conjunta também pode ser feita, o que considero
> melhor, pois permite "empréstimo de informação" entre os locais. Acaba-se
> caindo em modelo de efeitos mistos ao considerar efeito de parcela dentro
> de local, local e/ou bloco dentro de local como aleatórios. Para as
> respostas que são contagens (suporte nos inteiros não negativos) vale
> considerar a Poisson como candidata (candidata pois nem tudo que é contagem
> é Poisson), no entanto, análise sob uma transformação da resposta também
> pode ser conduzida. A vantagem da análise conjunta é verificar a existência
> de interação (se estirpe depende de local) e ser mais poderosa dado aumento
> do número de observações (que informam sobre os efeitos).
>
> 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.
>



-- 
Alisson Lucrecio da Costa
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20140609/8779f8dc/attachment.html>


Mais detalhes sobre a lista de discussão R-br