Boa noite a todos. Estou replicando um exemplo de um livro que verifica a dependência entre o tempo de viagem(y) ~ número de entregas(x1) +número de milhas percorridas. Quando peço a anova no R. anova(mod.cam), o R me dá a saída abaixo. Analysis of Variance Table Response: cam$tem_via Df Sum Sq Mean Sq F value Pr(>F) cam$n_entrega 1 9.0449 9.0449 27.535 0.001190 ** cam$milhas 1 12.5556 12.5556 38.222 0.000453 *** Residuals 7 2.2994 0.3285 --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 Eu gostaria da ANOVA com o Df da Regressão e do Resíduo, ou seja a anova reduzida. Alguém poderia me indicar qual a função ou rotina para isso? Em anexo o banco de dados em csv e a rotina. []' s Edson Lira Estatístico Manaus-Amazonas
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@ufpr.br twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================
Vale Walmes, boa noite! []'s Edson Lira Estatístico Manaus-Amazonas --- Em ter, 12/4/11, Walmes Zeviani <walmeszeviani@gmail.com> escreveu: De: Walmes Zeviani <walmeszeviani@gmail.com> Assunto: Re: [R-br] Anova multivariada. Para: r-br@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@ufpr.br twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ========================================================================== -----Anexo incorporado----- _______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
participantes (2)
-
Edson Lira -
Walmes Zeviani