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> 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> 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
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.

--

@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.