[R-br] Parcela subdividida: problema na declaração do erro associado a parcela

Jose Claudio Faria joseclaudio.faria em gmail.com
Quarta Dezembro 4 18:46:47 BRST 2013


Verdade!

Como deu erro eu não tentei rodar o summary.

Muito obrigado Tiago!

2013/12/4 Tiago Souza Marçal <tiagosouzamarcal em hotmail.com>:
> Mesmo com o aviso use o comando summary(m0) e o quadro será formado.
>
> Att.
>
> Tiago.
>
>
> #################################################################
>
> Tiago de Souza Marçal - Graduando em Agronomia pelo CCA-UFES
>
> Bolsista de Iniciação Científica da área de Genética e Melhoramento de
> Plantas
>
> #################################################################
>
>
>> Date: Wed, 4 Dec 2013 18:20:19 -0200
>> From: joseclaudio.faria em gmail.com
>> To: r-br em listas.c3sl.ufpr.br
>> Subject: [R-br] Parcela subdividida: problema na declaração do erro
>> associado a parcela
>
>>
>> Pessoal,
>>
>> O CMR abaixo é da apostila do Walmes.
>>
>>
>> #------------------------------------------------------------------------------------------
>> # dados
>> ps <- expand.grid(BL=c("I", "II", "III", "IV"),
>> ES=c("e1", "e2"),
>> AD=c("a1", "a2", "a3"))
>> ps$alt <-
>> c(58,77,38,52,44,59,30,34,85,90,73,77,59,68,45,55,66,93,67,64,54,75,53,48)
>> str(ps)
>> #
>>
>> #------------------------------------------------------------------------------------------
>> # análise de variância (erro A = BL x AD)
>> m0 <- aov(alt ~ BL + AD + Error(BL:AD) + ES + AD:ES, data=ps) #
>> termo Error para declarar erro A
>> m0 <- aov(alt ~ BL + AD * ES + Error(BL:AD), data=ps) # o mesmo com
>> fórmula comprimida
>>
>> summary(m0)
>>
>> #------------------------------------------------------------------------------------------
>>
>> Encontrei muitos outros exemplos (na web e em livros) com essa
>> notação: Error(BL:AD)
>>
>> Contudo o R (3.0.2patched) não está aceitando essa notação!
>> Fornece sempre a mensagem de erro:
>> In aov(alt ~ BL + AD + Error(BL:AD) + ES + AD:ES, data = ps) :
>> modelo Error() é singular
>>
>> Contudo o R está aceitando esse modelo:
>> m0 <- aov(alt ~ BL + AD + Error(BL/AD) + ES + AD:ES, data=ps)
>>
>> Mas nesse caso o resíduo A (associado a bloco e parcela - AD) não fica
>> correto.
>>
>> Alguém que anda trabalhando com esses modelos poderia ajudar?
>>
>> Ab,
>> --
>> ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\
>> Jose Claudio Faria
>> Estatistica
>> UESC/DCET/Brasil
>> joseclaudio.faria at gmail.com
>> Telefones:
>> 55(73)3680.5545 - UESC
>> 55(73)9100.7351 - TIM
>> 55(73)8817.6159 - OI
>> ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\
>> _______________________________________________
>> 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.
>
> _______________________________________________
> 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.



-- 
///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\
Jose Claudio Faria
Estatistica
UESC/DCET/Brasil
joseclaudio.faria at gmail.com
Telefones:
55(73)3680.5545 - UESC
55(73)9100.7351 - TIM
55(73)8817.6159 - OI
///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\


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