Um for resolve isso:

dados <- read.table('C:/Users/rcoster/Downloads/RCBD_dad_B.RData.txt', header=T)

mat <- matrix(0, 20, 20)

for (i in 1:20) {

 for (j in i:20) {

  mat[i,j] <- mat[j,i] <- t.test(subset(dados, Treat == i)$y, y = subset(dados, Treat == j)$y, var.equal = FALSE)$p.value

 }

}

Mas não é o caso de se usar ANOVA e depois Tukey?

2014-03-11 15:37 GMT-03:00 Luiz Roberto Martins Pinto <luizroberto.uesc@gmail.com>:

Éder,

Agradeço a sugestão, mas esta função usa 'The pool.sd switch calculates a common SD for all groups'. E eu preciso de uma função que viabilize a utilização de variancias diferentes para os tratamentos.

Luiz Roberto

Luiz Roberto Martins Pinto
Prof. Pleno/DCET/UESC

Laboratório de Estatística Computacional

Universidade Estadual de Santa Cruz

Ilhéus-Bahia

luizroberto.uesc@gmail.com
skype: lrmpinto

http://lattes.cnpq.br/2732314327604831

Em 11 de março de 2014 15:31, Éder Comunello <comunello.eder@gmail.com> escreveu:

Luiz Roberto, boa tarde!

Já consultou a função pairwise.t.test() também do pacote stats?

Éder Comunello <comunello.eder@gmail.com>
Dourados, MS - [22 16.5'S, 54 49'W]

_______________________________________________
R-br mailing list
R-br@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.

_______________________________________________
R-br mailing list
R-br@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.

_______________________________________________
R-br mailing list
R-br@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.