[R-br] Anova multivariada.

[Edson Lira]

Vale Walmes, boa noite!


[]'s

Edson Lira

Estatístico

Manaus-Amazonas

--- Em ter, 12/4/11, Walmes Zeviani <walmeszeviani em gmail.com> escreveu:

De: Walmes Zeviani <walmeszeviani em gmail.com>
Assunto: Re: [R-br] Anova multivariada.
Para: r-br em listas.c3sl.ufpr.br
Data: Terça-feira, 12 de Abril de 2011, 21:55

Edson,

Combine os resultados dos objetos "an" e "m1" abaixo em um novo data.frame. Ou você pode extrair a matriz do modelo e usa-la novamente dentro de outro lm(), pois a lm() pode receber fórmula (em que cada termo possui seus graus de liberdade reconhecidos) ou pode receber uma matriz (em que o número de colunas-1 da matriz é o grau de liberdade). Veja o CMR.


n <- 100
da <- data.frame(x=runif(n), z=3*runif(n))

da$y <- da$x+da$z+rnorm(n)
m0 <- lm(y~1, da)

m1 <- lm(y~x+z, da)
an <- anova(m0, m1)

str(an)
str(m1)


reg <- model.matrix(m1)
m2 <- lm(y~reg, da)

anova(m2)

À 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
==========================================================================


-----Anexo incorporado-----

_______________________________________________
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/20110412/dacb2481/attachment.html>