[R-br] [R_STAT] Seleção de Quadrados Médios

[Greice Laureano]

Walmes,

Muito obrigado! Funcionou perfeitamente =D

Att,
Greice.

2011/4/23 Walmes Zeviani <walmeszeviani em gmail.com>

> Greice,
>
> Segue o exemplo reproduzível de como acessar de diversas formas os valores
> do quadro de análise de variância para um modelo com Error() na fórmula.
>
> da <- expand.grid(bloco=gl(4,1,la="b"),
>                   parce=gl(3,1,la="p"),
>                   subpa=gl(5,1,la="s"))
> da$y <- rnorm(nrow(da))
>
> m0 <- aov(y~bloco+parce*subpa+Error(bloco:parce), data=da)
> summary(m0)
>
> str(m0) # apresenta a estrutura do objeto, você aprende como acessá-lo
> summ <- summary(m0)
> str(summ) # o que é? o que tem? como acesso?
> names(summ)
>
> summ[1] # porção parcela da anova
> summ[2] # porção subparcela da anova
>
> summ$"Error: bloco:parce" # acesso pelo nome do slot
> summ$"Error: Within"      # acesso pelo nome do slot
>
> summ[[1]][[1]] # acesso pelo número do slot
> summ[[2]][[1]] # acesso pelo número do slot
>
> summ$"Error: bloco:parce"[[1]][,"Df"]        # acesso por nome
> summ$"Error: bloco:parce"[[1]]["bloco","Df"] # acesso por nome
>
> summ[[1]][[1]][,"Sum Sq"]             # acesso por nome
> summ[[2]][[1]]["Residuals",]          # acesso por nome
> summ[[2]][[1]]["Residuals","Mean Sq"] # acesso por nome
>
> summ$"Error: bloco:parce"[[1]][,1]  # acesso por nome
> summ$"Error: bloco:parce"[[1]][1,1] # acesso por nome
>
> À disposição.
> Walmes.
>
> ==========================================================================
> Walmes Marques Zeviani
> LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W)
> Departamento de Estatística - Universidade Federal do Paraná
> fone: (+55) 41 3361 3573
> VoIP: (3361 3600) 1053 1173
> e-mail: walmes em ufpr.br
> twitter: @walmeszeviani
> homepage: http://www.leg.ufpr.br/~walmes
> linux user number: 531218
> ==========================================================================
>
>
> _______________________________________________
> R-br mailing list
> R-br em listas.c3sl.ufpr.br
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
>
>
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110423/c72a2a0d/attachment.html>