[R-br] sugestão análise experimento no tempo

[Maurício Lordêlo]

Wagner, obrigado vou tentar verificar  o que sugeriu.
Walmes, obrigado pela atenção. Peço desculpas, mas fiz o procedimento da
mesma maneira que faço sempre aqui na lista. Não sei porque desta vez não
funcionou.
Descrevi a forma como o experimento foi montado no primeiro e-mail.
Acrescento apenas que as plantas foram cultivadas em vasos localizados em
um ambiente controlado.
Como não consegui fazer a leitura da forma como sugeriu, segue uma maneira
não muito elegante de apresentar os dados.
Agradeço caso tenha sugestões.

hormonio =
factor(c(rep("AAS",15),rep("GA3",15),rep("AJ",15),rep("H2O",15),rep("STM",15),rep("BAP",15)))
tempo = factor(rep((1:3),each=5,, times = 6))
altura = c(43.5, 44.3, 47.9, 44.8, 46.5, 44.0, 52.8, 52.0, 42.0, 45.5,
50.2, 36.4, 50.9, 50.5, 49.9,
           46.0, 47.8, 45.1, 44.6, 42.5, 50.5, 43.4, 52.2, 47.5, 45.4,
48.6, 61.5, 58.0, 56.8, 53.5,
           47.2, 40.6, 46.6, 44.7, 42.3, 43.5, 41.5, 57.0, 47.2, 51.5,
58.5, 39.6, 52.5, 59.5, 47.4,
           45.9, 42.7, 42.1, 44.5, 47.2, 54.8, 58.8, 58.0, 49.2, 50.2,
56.7, 51.2, 57.1, 63.9, 45.9,
           45.1, 42.6, 50.8, 47.6, 41.7, 44.8, 54.5, 51.7, 49.8, 57.6,
63.4, 56.6, 52.6, 62.1, 54.8,
           47.6, 42.4, 45.8, 40.7, 41.2, 44.5, 41.0, 46.0, 59.0, 48.0,
53.7, 50.6, 45.5, 46.5, 53.7)
tapply(altura,hormonio,summary)
plot(altura~hormonio)
plot(altura~tempo)
require(latticeExtra)
xyplot(altura~hormonio|tempo)
xyplot(altura~tempo|hormonio)
modelo = aov(altura ~ hormonio*tempo)
summary(modelo)
par(mfrow=c(2,2))
plot(modelo)
shapiro.test(modelo$residuals)



Em 26 de outubro de 2015 18:46, Walmes Zeviani <walmeszeviani em gmail.com>
escreveu:

> Veja o CMR.
>
> ## Isso não deu certo. Apesar da URL, seu dado não é automaticamente
> ## trazido para a sessão. Eu estava curioso então me empenhei um pouco
> ## mas mesmo usado read.table() com a URL não consegui ler. Então parei.
> da <- read.table(
>     "https://www.dropbox.com/s/lkhzb87gh1vtgas/dados.txt",
>     header=TRUE, sep="\t")
> 'data.frame':    465 obs. of  1 variable:
>  $ X..DOCTYPE.html..html.lang.en.xmlns.fb.http...ogp.me.ns.fb..xml.lang.en.class.media.desktop.xmlns.http...www.w3.org.1999.xhtml..head..script.nonce.nlNb7FY.tMIOEhehezDS.:
> Factor w/ 230 levels "    ","        ",..: 222 216 194 32 165 224 223 29 9
> 11 ...
>
> Dá para ver os dados pela URL que você indicou. Você ficou a um passo de
> deixar o seu código reproduzível. Tente fazer com que da URL seja feita a
> leitura.
>
> A estrutura é fatorial completamente cruzado de hormônios (5) com tempo
> (3) com 5 registros por cela. Vai depender de como ela montou para saber se
> é parcela subdividida ou não. A suspeita é que seja parcela subdivididas.
>
> W.
>
> _______________________________________________
> 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/20151026/992d7322/attachment.html>