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@gmail.com To: r-br@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@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.