Boa tarde, Pessoal, sou novo n alista, mas já tenho minhas dúvidas. Sou usuário SAS e estou migrando para o R e consegui fazer contraste das médias, porém o parâmetro retornado é o estimate e gostaria que tivesse os mesmos dados da saída do contrats do SAS que traz a soma de quadrados, alguém sabe como saber ou algum package que faça isso? Att Daniel Lemes -- @lemes_daniel
Será o aov() ? Nunca usei o SAS. [1]Leonardo Ferreira Fontenelle Em Qui 4 set. 2014, às 12:51, Daniel Batista Lemes escreveu: Boa tarde, Pessoal, sou novo n alista, mas já tenho minhas dúvidas. Sou usuário SAS e estou migrando para o R e consegui fazer contraste das médias, porém o parâmetro retornado é o estimate e gostaria que tivesse os mesmos dados da saída do contrats do SAS que traz a soma de quadrados, alguém sabe como saber ou algum package que faça isso? Att Daniel Lemes -- @lemes_daniel _______________________________________________ R-br mailing list [2]R-br@listas.c3sl.ufpr.br [3]https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem ([4]http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível. References 1. http://lattes.cnpq.br/9234772336296638 2. mailto:R-br@listas.c3sl.ufpr.br 3. https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br 4. http://www.leg.ufpr.br/r-br-guia
Tem um livro ‘R for SAS and SPSS users’. Sugiro a leitura do mesmo De: R-br [mailto:r-br-bounces@listas.c3sl.ufpr.br] Em nome de Leonardo Ferreira Fontenelle Enviada em: quinta-feira, 4 de setembro de 2014 13:50 Para: r-br@listas.c3sl.ufpr.br Assunto: Re: [R-br] Contraste Será o aov() ? Nunca usei o SAS. Leonardo Ferreira Fontenelle <http://lattes.cnpq.br/9234772336296638> Em Qui 4 set. 2014, às 12:51, Daniel Batista Lemes escreveu: Boa tarde, Pessoal, sou novo n alista, mas já tenho minhas dúvidas. Sou usuário SAS e estou migrando para o R e consegui fazer contraste das médias, porém o parâmetro retornado é o estimate e gostaria que tivesse os mesmos dados da saída do contrats do SAS que traz a soma de quadrados, alguém sabe como saber ou algum package que faça isso? Att Daniel Lemes -- @lemes_daniel _______________________________________________ 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.
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.
Olá Walmes, Na realidade a soma de quadrados é só porque hoje é mostrado e é para fins didáticos, meu professor de delineamento que me pediu para converter os exercícios para R e trabalharmos eles na aula. no caso o SAS fornece essa saída Contrast DF Contrast SS Mean Square F Value Pr > F ab vs cd 1 1093.500000 1093.500000 10.84 0.0036 a vs b 1 243.000000 243.000000 2.41 0.1364 c vs d 1 300.000000 300.000000 2.97 0.1001 Mas se não tiver jeito deixamos assim, porque se quisermos a mesma saída do sas e do r podemos usar o estimate, obrigado. Em 4 de setembro de 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.
-- @lemes_daniel
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.
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
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.
Se forem k-1 (k=número de níveis de um fator) contrastes ortogonais, você desdobrar dentro do quadro de anova conforme tem-se nessa página http://www.leg.ufpr.br/~walmes/cursoR/ciaeear/ciaeear_08fat.html Para outros casos, tem-se outras abordagens. Bem, na minha opinião, mesmo que seja para fins didáticos, eu não posso deixar de expressar minha preocupação com o fato de as pessoas insistirem tentar fazer o R reproduzir o que outros aplicativos fazem tomando esses aplicativos como padrão ouro/referência sem antes parar para pensar se aquilo faz sentido. Não estou dizendo que o seu caso não faça sentido, ok? À disposição. Walmes.
No nosso caso usamos o SAS e damos aulas com o SAS, porém os alunos só tem acesso a eles no laboratório, então uma maneira seria ensinar os dois ou apenas o R, a ideia é fazer uma transição suave. Obrigado Em 4 de setembro de 2014 20:30, walmes . <walmeszeviani@gmail.com> escreveu:
Se forem k-1 (k=número de níveis de um fator) contrastes ortogonais, você desdobrar dentro do quadro de anova conforme tem-se nessa página
http://www.leg.ufpr.br/~walmes/cursoR/ciaeear/ciaeear_08fat.html
Para outros casos, tem-se outras abordagens. Bem, na minha opinião, mesmo que seja para fins didáticos, eu não posso deixar de expressar minha preocupação com o fato de as pessoas insistirem tentar fazer o R reproduzir o que outros aplicativos fazem tomando esses aplicativos como padrão ouro/referência sem antes parar para pensar se aquilo faz sentido. Não estou dizendo que o seu caso não faça sentido, ok?
À 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.
-- @lemes_daniel
participantes (6)
-
Daniel Batista Lemes -
Fernando Souza -
Leonard Assis -
Leonard de Assis
-
Leonardo Ferreira Fontenelle -
walmes .