[R-br] RES: Contraste

Fernando Souza nandodesouza em gmail.com
Quinta Setembro 4 20:03:23 BRT 2014


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 em gmail.com 
> <mailto:assis.leonard em 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 em gmail.com
>     <mailto:walmeszeviani em 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 em listas.c3sl.ufpr.br <mailto:R-br em 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 em listas.c3sl.ufpr.br <mailto:R-br em 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 em 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.

-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20140904/88d84d62/attachment.html>


Mais detalhes sobre a lista de discussão R-br