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


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.