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

[FHRB Toledo]

Cláudio... isso não seria um warnings!

É apenas para avisar que o a solução do sistema está vindo por uma inversa
generalizada pois sua matriz de incidência não tem posto completo!


2013/12/4 Jose Claudio Faria <joseclaudio.faria em gmail.com>

> 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
> ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\
> _______________________________________________
> 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.
>
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20131204/4a7ac6e3/attachment.html>