Re: [R-br] Testes de Comparação Múltipla! Análise Multivariada!

Walmes; Obrigado pela resposta!Mas será que eu não posso utilizar o teste de Bonferroni para a comparação múltipla duas a duas? Vc não tem algum exemplo que possa ser reproduzido? Utilizando a comparação de modelos aninhados como vc sugeriu? Novamente agradeço-lhe! Thiago de Paula Protásio Acadêmico de Engenharia Florestal Universidade Federal de Lavras Laboratório de Energia da Biomassa Florestal Laboratório de Biomateriais (035) 9183-2246

Thiago, Confesso que nunca vi alguém apresentar em um artigo uma comparação múltipla multivariada. Acredito que você pode usar um algoritmo k-means classificar os centróides em grupos homogeneos e tals, da mesma forma como SkottKnott faz. De qualquer forma, a comparação de modelos que eu me referi está representada no código abaixo da <- data.frame(tr=gl(3,8)) da$y1 <- rnorm(da$tr) da$y2 <- rnorm(da$tr) da$tr12vs3 <- as.factor(ifelse(as.numeric(da$tr)<=1,1,2)) m0 <- aov(cbind(y1,y2)~1, da) m1 <- aov(cbind(y1,y2)~tr12vs3, da) m2 <- aov(cbind(y1,y2)~tr, da) anova(m0, m1) anova(m1, 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 ==========================================================================

Para reflexão: se w é uma obs Bivariada com componentes (x1, y1), que critérios usar p determinar se ela é maior / menor / igual q z = (x2, y2)? - from a mobile device On 30 Jun 2011 23:32, "Walmes Zeviani" <walmeszeviani@gmail.com> wrote:
Thiago,
Confesso que nunca vi alguém apresentar em um artigo uma comparação múltipla multivariada. Acredito que você pode usar um algoritmo k-means classificar os centróides em grupos homogeneos e tals, da mesma forma como SkottKnott faz. De qualquer forma, a comparação de modelos que eu me referi está representada no código abaixo
da <- data.frame(tr=gl(3,8)) da$y1 <- rnorm(da$tr) da$y2 <- rnorm(da$tr) da$tr12vs3 <- as.factor(ifelse(as.numeric(da$tr)<=1,1,2))
m0 <- aov(cbind(y1,y2)~1, da) m1 <- aov(cbind(y1,y2)~tr12vs3, da) m2 <- aov(cbind(y1,y2)~tr, da)
anova(m0, m1) anova(m1, 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 ==========================================================================

Em 30/6/2011 20:51, Benilton Carvalho escreveu:
Para reflexão: se w é uma obs Bivariada com componentes (x1, y1), que critérios usar p determinar se ela é maior / menor / igual q z = (x2, y2)?
- from a mobile device Benilton,
Antes do critério, você precisa dizer qual a norma (em termos matemáticos estatísticos) ou "métrica" (em termos de linguagem de negócio) você usa para ordenar as componentes. Um dos grandes desafios está exatamente aí. -- Cesar Rabak GNU/Linux User 52247. Get counted: http://counter.li.org/

O 'truque' da norma, utilizado por todos (e eu incluso), e' apenas por nao termos uma definicao de ordenacao em R^n... e, assim, acaba-se por colapsar R^n -> R... mas, sim, a 'solucao' corrente... e ainda fica o tema para reflexao. b 2011/7/1 Cesar Rabak <cesar.rabak@gmail.com>:
Em 30/6/2011 20:51, Benilton Carvalho escreveu:
Para reflexão: se w é uma obs Bivariada com componentes (x1, y1), que critérios usar p determinar se ela é maior / menor / igual q z = (x2, y2)?
- from a mobile device
Benilton,
Antes do critério, você precisa dizer qual a norma (em termos matemáticos estatísticos) ou "métrica" (em termos de linguagem de negócio) você usa para ordenar as componentes.
Um dos grandes desafios está exatamente aí.
-- Cesar Rabak
participantes (4)
-
Benilton Carvalho
-
Cesar Rabak
-
Thiago De paula protásio
-
Walmes Zeviani