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 <c <comunello.eder@gmail.com>omunello.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.