Faça através da função glht {multcomp} dados<-structure(list(Variedade = structure(c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L), .Label = c("RB92579 ", "RB931003", "RB931011", "RB951541", "RB98710", "RB99395"), class = "factor"), Rep = c(1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L), Coleta = structure(c(1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L), .Label = c("1", "2"), class = "factor"), MO = c(97.25, 96.2, 95.19, 96.06, 98.79, 98.66, 98.92, 98.85, 92.47, 96.33, 94.39, 96.55, 99.13, 99.02, 98.93, 99.06, 95.3, 94.07, 95, 95.34, 99.18, 99.71, 99.48, 98.63, 96.38, 97.9, 96.78, 97.46, 98.72, 99.26, 98.97, 98.8, 95.42, 93.7, 93.38, 94.99, 98.12, 99.83, 99.03, 99.1, 94.18, 94.89, 94.98, 94.19, 98.87, 99.59, 98.73, 98.76)), .Names = c("Variedade", "Rep", "Coleta", "MO"), class = "data.frame", row.names = c(NA, -48L)) analise<-lm(MO~Variedade,data=dados[,]) comparacaoVar<-glht(analise,linfct=mcp(Variedade="Tukey")) #onde esta escrito tukey você pode colocar uma matriz ou vetor de contrastes de seu interesse. Neste exemplo o argumento "Tukey" cria uma matriz de contraste de todos os níveis do fator Variedade summary(comparacaoVar) On 04-09-2014 19:33, Daniel Batista Lemes wrote:
Obrigado, estou engatinhando no R, mas pro que eu preciso não preciso de muita coisa.
Abraço
Em 4 de setembro de 2014 19:32, Leonard Assis <assis.leonard@gmail.com <mailto:assis.leonard@gmail.com>> escreveu:
Se não me engano o car tem funções que retornam resultado comparável ao SAS
Em 04/09/2014 19:05, "walmes ." <walmeszeviani@gmail.com <mailto:walmeszeviani@gmail.com>> escreveu:
Se o contrate for sobre uma hipótese (H0: A-B==0) então a estatística t ao quadrado equivale à F que vem dessa saída que mostra soma de quadrados. O p-valor será o mesmo, entende? Então você de fato não precisa conhecer a soma de quadrados para testar essa hipótese. Se o contraste corresponde à uma hipótese conjunta (H0: A-B==0 & A-C==0 & A-D==0) então as coisas diferem e aí tem-se mais de uma forma de avaliar essa hipótese, assim como no caso de uma hipótese simples. Meu workhorse para testes de hipótese linear é o pacote multcomp. O pacote car também fornce algumas funções úteis.
À disposição. Walmes.
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br <mailto: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 <mailto: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.
--
@lemes_daniel
_______________________________________________ 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.