Obrigado pelas respostas, vocês são muitos solícitos.<div>Já tinha usado a função glht() para dados com distribuição normal, usarei então junto com a função glm().</div><div><br></div><div>Um abraço à todos,</div><div><br>
</div><div>Danilo Scorzoni Ré<br><br><div class="gmail_quote">2011/10/26 <span dir="ltr"><<a href="mailto:r-br-request@listas.c3sl.ufpr.br">r-br-request@listas.c3sl.ufpr.br</a>></span><br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
Enviar submissões para a lista de discussão R-br para<br>
<a href="mailto:r-br@listas.c3sl.ufpr.br">r-br@listas.c3sl.ufpr.br</a><br>
<br>
Para se cadastrar ou descadastrar via WWW, visite o endereço<br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
ou, via email, envie uma mensagem com a palavra 'help' no assunto ou<br>
corpo da mensagem para<br>
<a href="mailto:r-br-request@listas.c3sl.ufpr.br">r-br-request@listas.c3sl.ufpr.br</a><br>
<br>
Você poderá entrar em contato com a pessoa que gerencia a lista pelo<br>
endereço<br>
<a href="mailto:r-br-owner@listas.c3sl.ufpr.br">r-br-owner@listas.c3sl.ufpr.br</a><br>
<br>
Quando responder, por favor edite sua linha Assunto assim ela será<br>
mais específica que "Re: Contents of R-br digest..."<br>
<br>
<br>
Tópicos de Hoje:<br>
<br>
1. Re: Gráficos no R ou Gnu Plot (Fernando Colugnati)<br>
2. Comparações Múltiplas em GLM (Danilo Scorzoni Ré)<br>
3. Re: Comparações Múltiplas em GLM (Edson Lira)<br>
4. Re: Comparações Múltiplas em GLM (Benilton Carvalho)<br>
5. Re: Comparações Múltiplas em GLM (Walmes Zeviani)<br>
<br>
<br>
----------------------------------------------------------------------<br>
<br>
Message: 1<br>
Date: Tue, 25 Oct 2011 23:09:25 -0200<br>
From: Fernando Colugnati <<a href="mailto:fernando@ipti.org.br">fernando@ipti.org.br</a>><br>
To: <a href="mailto:r-br@listas.c3sl.ufpr.br">r-br@listas.c3sl.ufpr.br</a><br>
Subject: Re: [R-br] Gráficos no R ou Gnu Plot<br>
Message-ID:<br>
<<a href="mailto:CAOLzYt28ZB%2BK4we38999v_62mXJkK08inXdptTYX1SFzere%2BaA@mail.gmail.com">CAOLzYt28ZB+K4we38999v_62mXJkK08inXdptTYX1SFzere+aA@mail.gmail.com</a>><br>
Content-Type: text/plain; charset="iso-8859-1"<br>
<br>
Não devemos perpetuar práticas mesmo, e na estatística temos muitas práticas<br>
ruins sendo perpetuadas que não dizem respeito a como editar ou fazer um<br>
gráfico para que mais pessoas utilizem. Temos a horrível prática do uso do<br>
p-valor como oráculo, de testes de hipótese que testam hipóteses erradas,<br>
modelos enviesados, etc, etc...com isso perco o sono, fico angustiado<br>
etc...(e revistas de áreas usuárias da estatística pedem essas coisas, assim<br>
como pedem arquivos em xls...isso sim é o fim da picada). Com os gráficos em<br>
excel eu já realmente relaxei e passei a ter mais tempo pra estas outras<br>
questões! Há pouco tempo um cara morreu de cancer tentando mudar o<br>
mercado...e não conseguiu um grande sucesso neste quesito (mas em outros,<br>
temos de concordar)!<br>
Abs<br>
<br>
<br>
Em 25 de outubro de 2011 21:29, Eder David Borges da Silva <<a href="mailto:eder@leg.ufpr.br">eder@leg.ufpr.br</a><br>
> escreveu:<br>
<br>
> Pessoal,<br>
> Este tópico é interessante é assim como o Pedro sei como é difícil<br>
> utilizar um software que só você usar (R) dentro de uma empresa que só<br>
> utilizar pacote office (E na maioria das vezes não desfrutando das<br>
> poucas facilidades que eles te ajudam), tenho seguido a seguinte<br>
> filosofia, não gostou do gráfico faça você!! com seu pacote preferido.<br>
> Em minha análise eu creio que as pessoas de modo geral perder 50% do<br>
> seu tempo "enfrufruzando" tabelas, gráficos e mapas de um relatório,<br>
> acho isso um desperdício de vida, mas fazer oque né....<br>
> ate mais<br>
><br>
><br>
><br>
> Em 25 de outubro de 2011 19:38, Paulo Justiniano<br>
> <<a href="mailto:paulojus@leg.ufpr.br">paulojus@leg.ufpr.br</a>> escreveu:<br>
> > perpetuar a pratica nao ajuda muito....<br>
> ><br>
> ><br>
> > On Tue, 25 Oct 2011, Pedro Rafael wrote:<br>
> ><br>
> >> Não Benilton é usar gráficos mais atraentes para ser colocado em um<br>
> >> Power Point, não gráficos cegos e apagados para não dar sono no<br>
> >> pessoal rsrsrsrs. Eu particularmente acho frescura. Sofro quando passo<br>
> >> os gráficos feito no R porque tem algum fulano que quer editar tais<br>
> >> gráficos. Criei várias funções que leem arquivos DBF que é formato que<br>
> >> eu tenho aqui para calcular indicadores para estados, municípios,<br>
> >> regionais de saúde mas os gráficos o pessoal reclama porque não sabem<br>
> >> de R. Aí eu fico no dilema do que usar. Mas o Lema não é esse que você<br>
> >> citou! Fazer um gráfico com um efeito em 3D a depender do gráfico não<br>
> >> complica a cabeça de ninguém.<br>
> >><br>
> >> 2011/10/25 beniltoncarvalho [via R-br]<br>
> >> <<a href="mailto:ml-node%2Bs2285057n3936773h40@n4.nabble.com">ml-node+s2285057n3936773h40@n4.nabble.com</a>>:<br>
> >>><br>
> >>> Bem interessante a atitude da instituição: "para falar com leigos,<br>
> >>> precisamos usar gráficos que complicam o entendimento da nossa<br>
> mensagem".<br>
> >>><br>
> >>><br>
> >>> --<br>
> >>> Successful people ask better questions, and as a result, they get<br>
> better<br>
> >>> answers. (Tony Robbins)<br>
> >>><br>
> >>> _______________________________________________<br>
> >>> R-br mailing list<br>
> >>> [hidden email]<br>
> >>> <a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
> >>> Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça<br>
> >>> código<br>
> >>> mínimo reproduzível.<br>
> >>><br>
> >>> ________________________________<br>
> >>> If you reply to this email, your message will be added to the<br>
> discussion<br>
> >>> below:<br>
> >>><br>
> >>><br>
> <a href="http://r-br.2285057.n4.nabble.com/R-br-Graficos-no-R-ou-Gnu-Plot-tp3936555p3936773.html" target="_blank">http://r-br.2285057.n4.nabble.com/R-br-Graficos-no-R-ou-Gnu-Plot-tp3936555p3936773.html</a><br>
> >>> To unsubscribe from R-br, click here.<br>
> >><br>
> >><br>
> >><br>
> >> --<br>
> >> Saudações,<br>
> >> Pedro Rafael Diniz Marinho.<br>
> >> Estatístico - Secretaria de Estado da Saúde - PB.<br>
> >> _______________________________________________<br>
> >> R-br mailing list<br>
> >> <a href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br>
> >> <a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
> >> Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça<br>
> >> código mínimo reproduzível.<br>
> ><br>
> > _______________________________________________<br>
> > R-br mailing list<br>
> > <a href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br>
> > <a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
> > Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça<br>
> código<br>
> > mínimo reproduzível.<br>
> ><br>
> _______________________________________________<br>
> R-br mailing list<br>
> <a href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br>
> <a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
> Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça<br>
> código mínimo reproduzível.<br>
><br>
<br>
<br>
<br>
--<br>
Fernando A.B. Colugnati<br>
Pesquisador Associado<br>
<br>
Instituto de Pesquisas em Tecnologia e Inovação - IPTI<br>
Tel. 55 11 8704-9812<br>
<a href="http://www.ipti.org.br" target="_blank">www.ipti.org.br</a><br>
<a href="mailto:fernando@ipti.org.br">fernando@ipti.org.br</a><br>
-------------- Próxima Parte ----------<br>
Um anexo em HTML foi limpo...<br>
URL: <<a href="http://listas.inf.ufpr.br/pipermail/r-br/attachments/20111025/5a531e18/attachment-0001.html" target="_blank">http://listas.inf.ufpr.br/pipermail/r-br/attachments/20111025/5a531e18/attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Message: 2<br>
Date: Wed, 26 Oct 2011 10:11:06 -0200<br>
From: Danilo Scorzoni Ré <<a href="mailto:danilo.scorzoni@gmail.com">danilo.scorzoni@gmail.com</a>><br>
To: <a href="mailto:r-br@listas.c3sl.ufpr.br">r-br@listas.c3sl.ufpr.br</a><br>
Subject: [R-br] Comparações Múltiplas em GLM<br>
Message-ID:<br>
<CA+bJVreC=<a href="mailto:uKJ_XJs4495NcK_tUuJ8FyQ%2BsWG5ruKhD77H%2BMo3g@mail.gmail.com">uKJ_XJs4495NcK_tUuJ8FyQ+sWG5ruKhD77H+Mo3g@mail.gmail.com</a>><br>
Content-Type: text/plain; charset="iso-8859-1"<br>
<br>
Bom dia a todos,<br>
<br>
É a primeira vez que envio para essa lista e, ainda por cima, é um assunto<br>
que eu já vi comentarem por aqui, mas não tenho mais o e-mail (e não<br>
encontrei nenhum local na web que dê pra acessar as mensagens já<br>
enviadas...), mas meu problema é o seguinte.<br>
<br>
Estou analisando um conjunto de dados de um delineamento casualizado em<br>
blocos em duas áreas distintas. A primeira questão é que esse experimento<br>
está instalado com o mesmo delineamento em duas áreas distintas e gostaria<br>
de avaliar o efeito da área. Segundo o livro do Pimentel-Gomes, eu preciso<br>
checar a homogeneidade de variâncias antes de proceder a análise usual, usei<br>
então um teste de Bartlett e um de Levenne. Em ambos, a hipótese nula não<br>
foi rejeitada, ou seja, não tenho evidência de que as variâncias sejam<br>
heterogêneas. Porém, dada a natureza dos dados, não acho que seja isso<br>
(proporções de variâncias entre a menor e a maior estão acima de 7). Dessa<br>
forma, procedi a análise através do GLM com resposta com distribuição Gama.<br>
Consegui um bom ajuste com boa análise de resíduos. A minha dúvida agora é:<br>
como fazer comparações múltiplas usando o GLM? Sei que o Tukey pressupõe<br>
modelos paramétricos, mas não me recordo se ele precisa pressupor<br>
distribuição normal.<br>
<br>
Desde já agradeço a atenção de todos e desculpe pelo texto extenso.<br>
<br>
Um abraço,<br>
<br>
--<br>
Danilo Scorzoni Ré<br>
Engenheiro Florestal<br>
Mestrando em Ciência Florestal<br>
FCA / UNESP - Botucatu<br>
(14) 8180-2494<br>
-------------- Próxima Parte ----------<br>
Um anexo em HTML foi limpo...<br>
URL: <<a href="http://listas.inf.ufpr.br/pipermail/r-br/attachments/20111026/54381e6a/attachment-0001.html" target="_blank">http://listas.inf.ufpr.br/pipermail/r-br/attachments/20111026/54381e6a/attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Message: 3<br>
Date: Wed, 26 Oct 2011 05:33:17 -0700 (PDT)<br>
From: Edson Lira <<a href="mailto:edinhoestat@yahoo.com.br">edinhoestat@yahoo.com.br</a>><br>
To: "<a href="mailto:r-br@listas.c3sl.ufpr.br">r-br@listas.c3sl.ufpr.br</a>" <<a href="mailto:r-br@listas.c3sl.ufpr.br">r-br@listas.c3sl.ufpr.br</a>><br>
Subject: Re: [R-br] Comparações Múltiplas em GLM<br>
Message-ID:<br>
<<a href="mailto:1319632397.14122.YahooMailNeo@web160707.mail.bf1.yahoo.com">1319632397.14122.YahooMailNeo@web160707.mail.bf1.yahoo.com</a>><br>
Content-Type: text/plain; charset="iso-8859-1"<br>
<br>
Danilo, caso o seus dados não atenda as pressuposições,tem um teste não paramétrico que faz teste de médias. Veja o exemplo abaixo do pacote agricolae.<br>
<br>
Não esqueça de instalar o pacote agricolae<br>
<br>
library(agricolae)<br>
<br>
?kruskal<br>
<br>
data(corn)<br>
attach(corn)<br>
str(corn)<br>
comparison<-kruskal(observation,method,group=TRUE, main="corn")<br>
detach(corn)<br>
<br>
No caso da homogeneidade da variância no pacote MASS tem a transformação de BOXCOX que é muito usada.<br>
<br>
#Examples<br>
?boxcox<br>
require(MASS)<br>
boxcox(Volume ~ log(Height) + log(Girth), data = trees,<br>
lambda = seq(-0.25, 0.25, length = 10))<br>
<br>
boxcox(Days+1 ~ Eth*Sex*Age*Lrn, data = quine,<br>
lambda = seq(-0.05, 0.45, len = 20))<br>
<br>
<br>
[]'s.<br>
Edson Lira<br>
Estatístico<br>
Manaus-Amazonas<br>
<br>
<br>
________________________________<br>
De: Danilo Scorzoni Ré <<a href="mailto:danilo.scorzoni@gmail.com">danilo.scorzoni@gmail.com</a>><br>
Para: <a href="mailto:r-br@listas.c3sl.ufpr.br">r-br@listas.c3sl.ufpr.br</a><br>
Enviadas: Quarta-feira, 26 de Outubro de 2011 8:11<br>
Assunto: [R-br] Comparações Múltiplas em GLM<br>
<br>
<br>
Bom dia a todos,<br>
<br>
É a primeira vez que envio para essa lista e, ainda por cima, é um assunto que eu já vi comentarem por aqui, mas não tenho mais o e-mail (e não encontrei nenhum local na web que dê pra acessar as mensagens já enviadas...), mas meu problema é o seguinte.<br>
<br>
Estou analisando um conjunto de dados de um delineamento casualizado em blocos em duas áreas distintas. A primeira questão é que esse experimento está instalado com o mesmo delineamento em duas áreas distintas e gostaria de avaliar o efeito da área. Segundo o livro do Pimentel-Gomes, eu preciso checar a homogeneidade de variâncias antes de proceder a análise usual, usei então um teste de Bartlett e um de Levenne. Em ambos, a hipótese nula não foi rejeitada, ou seja, não tenho evidência de que as variâncias sejam heterogêneas. Porém, dada a natureza dos dados, não acho que seja isso (proporções de variâncias entre a menor e a maior estão acima de 7). Dessa forma, procedi a análise através do GLM com resposta com distribuição Gama. Consegui um bom ajuste com boa análise de resíduos. A minha dúvida agora é: como fazer comparações múltiplas usando o GLM? Sei que o Tukey pressupõe modelos paramétricos, mas não me recordo se<br>
ele precisa pressupor distribuição normal.<br>
<br>
Desde já agradeço a atenção de todos e desculpe pelo texto extenso.<br>
<br>
Um abraço,<br>
<br>
--<br>
Danilo Scorzoni Ré<br>
Engenheiro Florestal<br>
Mestrando em Ciência Florestal<br>
FCA / UNESP - Botucatu<br>
(14) 8180-2494<br>
<br>
<br>
_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br>
-------------- Próxima Parte ----------<br>
Um anexo em HTML foi limpo...<br>
URL: <<a href="http://listas.inf.ufpr.br/pipermail/r-br/attachments/20111026/dac3b546/attachment-0001.html" target="_blank">http://listas.inf.ufpr.br/pipermail/r-br/attachments/20111026/dac3b546/attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Message: 4<br>
Date: Wed, 26 Oct 2011 13:38:20 +0100<br>
From: Benilton Carvalho <<a href="mailto:beniltoncarvalho@gmail.com">beniltoncarvalho@gmail.com</a>><br>
To: <a href="mailto:r-br@listas.c3sl.ufpr.br">r-br@listas.c3sl.ufpr.br</a>, Edson Lira <<a href="mailto:edinhoestat@yahoo.com.br">edinhoestat@yahoo.com.br</a>><br>
Subject: Re: [R-br] Comparações Múltiplas em GLM<br>
Message-ID:<br>
<<a href="mailto:CAO-arWNG_7VnipqftsexdORqGrG0ZWjtt_FnKcqEgiwXZA5O5A@mail.gmail.com">CAO-arWNG_7VnipqftsexdORqGrG0ZWjtt_FnKcqEgiwXZA5O5A@mail.gmail.com</a>><br>
Content-Type: text/plain; charset=ISO-8859-1<br>
<br>
Hipoteses de distribuicao para casos assim referem-se `a distribuicao<br>
dos **estimadores**, nao dos dados... E minha (fraca) memoria nao se<br>
recorda de nenhuma hipotese de distribuicao no teste de Tukey, mas<br>
posso estar enganado....<br>
<br>
<br>
------------------------------<br>
<br>
Message: 5<br>
Date: Wed, 26 Oct 2011 10:45:07 -0200<br>
From: Walmes Zeviani <<a href="mailto:walmeszeviani@gmail.com">walmeszeviani@gmail.com</a>><br>
To: <a href="mailto:r-br@listas.c3sl.ufpr.br">r-br@listas.c3sl.ufpr.br</a><br>
Subject: Re: [R-br] Comparações Múltiplas em GLM<br>
Message-ID:<br>
<CAFU=<a href="mailto:EkYkW5EsmUR8aKG5KTsQHa%2BeVLsW4wkNmc_9LkyScO2taw@mail.gmail.com">EkYkW5EsmUR8aKG5KTsQHa+eVLsW4wkNmc_9LkyScO2taw@mail.gmail.com</a>><br>
Content-Type: text/plain; charset="iso-8859-1"<br>
<br>
Tukey supõe distribuição normal para os estimadores (como Benilton<br>
assinalou) pois é baseado da distribuição da amplitude total studentizada.<br>
Assintóticamente os estimadores nos GLMs terão distribuição normal. Então é<br>
possível usar o procedimento. Detalhe é que como nos GLMs a variância não é<br>
constante e é em geral dependente da média, as estimativas terão diferentes<br>
erros padrões, e por isso assim também serão os contrastes entre elas. A<br>
multcomp::glht() oferece recursos para isso, mas não aplica teste de Tukey,<br>
faz os contrastes de Tukey (2 a 2) e corrige p-valor por diversos métodos<br>
(opções da função).<br>
<br>
À disposição.<br>
Walmes.<br>
<br>
==========================================================================<br>
Walmes Marques Zeviani<br>
LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W)<br>
Departamento de Estatística - Universidade Federal do Paraná<br>
fone: <a href="tel:%28%2B55%29%2041%203361%203573" value="+554133613573">(+55) 41 3361 3573</a><br>
VoIP: (3361 3600) 1053 1173<br>
e-mail: <a href="mailto:walmes@ufpr.br">walmes@ufpr.br</a><br>
twitter: @walmeszeviani<br>
homepage: <a href="http://www.leg.ufpr.br/~walmes" target="_blank">http://www.leg.ufpr.br/~walmes</a><br>
linux user number: 531218<br>
==========================================================================<br>
-------------- Próxima Parte ----------<br>
Um anexo em HTML foi limpo...<br>
URL: <<a href="http://listas.inf.ufpr.br/pipermail/r-br/attachments/20111026/13fcd070/attachment-0001.html" target="_blank">http://listas.inf.ufpr.br/pipermail/r-br/attachments/20111026/13fcd070/attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
<br>
<br>
Fim da Digest R-br, volume 8, assunto 33<br>
****************************************<br>
</blockquote></div><br><br clear="all"><div><br></div>-- <br>Danilo Scorzoni Ré<br>Engenheiro Florestal<br>Mestrando em Ciência Florestal<br>FCA / UNESP - Botucatu<br>(14) 8180-2494<br><br>
</div>