<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40"><head><meta http-equiv=Content-Type content="text/html; charset=utf-8"><meta name=Generator content="Microsoft Word 15 (filtered medium)"><!--[if !mso]><style>v\:* {behavior:url(#default#VML);}
o\:* {behavior:url(#default#VML);}
w\:* {behavior:url(#default#VML);}
.shape {behavior:url(#default#VML);}
</style><![endif]--><style><!--
/* Font Definitions */
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
@font-face
{font-family:"dejavu sans mono";
panose-1:0 0 0 0 0 0 0 0 0 0;}
@font-face
{font-family:Consolas;
panose-1:2 11 6 9 2 2 4 3 2 4;}
@font-face
{font-family:Verdana;
panose-1:2 11 6 4 3 5 4 4 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
margin-bottom:.0001pt;
font-size:12.0pt;
font-family:"Times New Roman",serif;}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:purple;
text-decoration:underline;}
pre
{mso-style-priority:99;
mso-style-link:"Pré-formatação HTML Char";
margin:0cm;
margin-bottom:.0001pt;
font-size:10.0pt;
font-family:"Courier New";}
p.msonormal0, li.msonormal0, div.msonormal0
{mso-style-name:msonormal;
mso-margin-top-alt:auto;
margin-right:0cm;
mso-margin-bottom-alt:auto;
margin-left:0cm;
font-size:12.0pt;
font-family:"Times New Roman",serif;}
span.Pr-formataoHTMLChar
{mso-style-name:"Pré-formatação HTML Char";
mso-style-priority:99;
mso-style-link:"Pré-formatação HTML";
font-family:Consolas;
mso-fareast-language:PT-BR;}
span.m8358487097371014076gmail-gghfmyibcpb
{mso-style-name:m_8358487097371014076gmail-gghfmyibcpb;}
span.m8358487097371014076gmail-gghfmyibcob
{mso-style-name:m_8358487097371014076gmail-gghfmyibcob;}
span.m8358487097371014076gmail-gghfmyibgob
{mso-style-name:m_8358487097371014076gmail-gghfmyibgob;}
span.m8358487097371014076m6541542930473827578m6389802118341864428m799196135501887637m-6160994631685705737gmail-gghfmyibcpb
{mso-style-name:m_8358487097371014076m_6541542930473827578m6389802118341864428m799196135501887637m-6160994631685705737gmail-gghfmyibcpb;}
span.m8358487097371014076m6541542930473827578m6389802118341864428m799196135501887637m-6160994631685705737gmail-gghfmyibcob
{mso-style-name:m_8358487097371014076m_6541542930473827578m6389802118341864428m799196135501887637m-6160994631685705737gmail-gghfmyibcob;}
span.EstiloDeEmail26
{mso-style-type:personal-reply;
font-family:"Calibri",sans-serif;
color:windowtext;}
.MsoChpDefault
{mso-style-type:export-only;
font-family:"Calibri",sans-serif;
mso-fareast-language:EN-US;}
@page WordSection1
{size:612.0pt 792.0pt;
margin:70.85pt 3.0cm 70.85pt 3.0cm;}
div.WordSection1
{page:WordSection1;}
/* List Definitions */
@list l0
{mso-list-id:113015061;
mso-list-template-ids:-1147741792;}
@list l1
{mso-list-id:1129477554;
mso-list-template-ids:-538422238;}
ol
{margin-bottom:0cm;}
ul
{margin-bottom:0cm;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]--></head><body lang=PT-BR link=blue vlink=purple><div class=WordSection1><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri",sans-serif;mso-fareast-language:EN-US'>Marcos,<o:p></o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri",sans-serif;mso-fareast-language:EN-US'><o:p> </o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri",sans-serif;mso-fareast-language:EN-US'>Não vou discutir escolhas de Factor e etc pois não é o assunto. Sugiro que revise o assunto.<o:p></o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri",sans-serif;mso-fareast-language:EN-US'><o:p> </o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri",sans-serif;mso-fareast-language:EN-US'><o:p> </o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri",sans-serif;mso-fareast-language:EN-US'>O erro no primeiro ajuste(binomial(link=log) é que o modelo não está conseguindo convergir.<o:p></o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri",sans-serif;mso-fareast-language:EN-US'>O modelo Poisson, (Poisson, quase Poisson, binomial negativo, etc) não faz sentido em seus dados, vide o que expliquei nos e-mails anteriores.<o:p></o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri",sans-serif;mso-fareast-language:EN-US'><o:p> </o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri",sans-serif;mso-fareast-language:EN-US'>O que você pode fazer?<o:p></o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri",sans-serif;mso-fareast-language:EN-US'><o:p> </o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri",sans-serif;mso-fareast-language:EN-US'>Você tem dúzias de artigos comentando que o melhor seria binomial(link=log), mas você tem problemas de convergência neste caso. Como você já tem uma estimativa dos parâmetros com link logit, utilize estes coeficientes como chute inicial do modelo.<o:p></o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri",sans-serif;mso-fareast-language:EN-US'><o:p> </o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri",sans-serif;mso-fareast-language:EN-US'>Não gosto destes testes gráficos de ajuste, prefiro olhar cross-validation, roc curve ou algo assim.<o:p></o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri",sans-serif;mso-fareast-language:EN-US'><o:p> </o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri",sans-serif;mso-fareast-language:EN-US'>Leonard<o:p></o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri",sans-serif;mso-fareast-language:EN-US'><o:p> </o:p></span></p><p class=MsoNormal><b><span style='font-size:11.0pt;font-family:"Calibri",sans-serif'>De:</span></b><span style='font-size:11.0pt;font-family:"Calibri",sans-serif'> R-br [mailto:r-br-bounces@listas.c3sl.ufpr.br] <b>Em nome de </b>Marcos Bissoli via R-br<br><b>Enviada em:</b> quinta-feira, 9 de fevereiro de 2017 10:35<br><b>Para:</b> a lista Brasileira oficial de discussão do programa R. <r-br@listas.c3sl.ufpr.br><br><b>Assunto:</b> [R-br] Fwd: Distribuição para regressão de resposta binária<o:p></o:p></span></p><p class=MsoNormal><o:p> </o:p></p><div><p class=MsoNormal><o:p> </o:p></p><div><p class=MsoNormal style='margin-bottom:12.0pt'>---------- Mensagem encaminhada ----------<br>De: <b>Marcos Bissoli</b> <<a href="mailto:mbissoli@gmail.com">mbissoli@gmail.com</a>><br>Data: 9 de fevereiro de 2017 09:38<br>Assunto: Re: [R-br] Distribuição para regressão de resposta binária<br>Para: Leonard Mendonça de Assis <<a href="mailto:assis.leonard@gmail.com">assis.leonard@gmail.com</a>><br><br><o:p></o:p></p><div><p class=MsoNormal>Bom dia,<o:p></o:p></p><div><p class=MsoNormal><o:p> </o:p></p></div><div><p class=MsoNormal>Mais uma vez, muito obrigado pela oportunidade de diálogo.<o:p></o:p></p></div><div><p class=MsoNormal><o:p> </o:p></p></div><div><p class=MsoNormal>Seguem códigos aplicados ao meu banco de dados. Modelei com diversas técnicas. Perceba que os gráficos de diagnóstico não diferem muito entre si, nestas que estou testando. E a dúvida, talvez, central, seja essa: esses gráficos são, de fato, tão relevantes? Pois, como eu disse, alguns testes de ajuste que fiz deram resultados satisfatórios. Até que ponto esse diagnóstico visual de ajuste é um problema? Percebam, também, que não há grandes diferenças inclusive com o gráfico da regressão logística. Meu problema com a logística, repito, é que ela me retorna odds ratios, e essa medida não tem sido bem aceita entre epidemiologistas. Por isso eu prefiro famílias com função de ligação "log", pois assim eu consigo facilmente as razões de prevalência.<o:p></o:p></p></div><div><p class=MsoNormal><o:p> </o:p></p></div><div><p class=MsoNormal>Percebam, ainda, que a log binomial retorna erro e não gera um objeto. Há algum ajuste a mais que eu possa fazer no código para ela funcionar? Na literatura que consultei, quando isso ocorre, a sugestão é que se use a Poisson com variância robusta. Percebam, também, que os coeficientes significantes não mudam de modelo para modelo, à exceção da Quasi-Poisson que é mais "sensível" (estou usando um termo epidemiológico aqui) para detectar fatores potencialmente causais.<o:p></o:p></p></div><div><p class=MsoNormal><o:p> </o:p></p></div><div><p class=MsoNormal>Tabagismo é a variável binária de resposta: 1 para fumantes, 0 para não fumantes. As demais marcadas com "factor" são também binárias. As sem "factor" são contínuas ou discretas.<o:p></o:p></p></div><div><p class=MsoNormal><o:p> </o:p></p></div><pre style='line-height:12.0pt;word-break:break-all;outline:none;white-space:pre-wrap' id="m_8358487097371014076gmail-rstudio_console_output"><span class=m8358487097371014076gmail-gghfmyibcpb><span style='font-family:"dejavu sans mono",serif;color:blue'>> </span></span><span class=m8358487097371014076gmail-gghfmyibcob><span style='font-family:"dejavu sans mono",serif;color:blue'>TabModelagem <- data.frame(Tabagismo,factor(SexoDic),factor(Branca),factor(Negra),factor(Parda),<o:p></o:p></span></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span class=m8358487097371014076gmail-gghfmyibcpb><span style='font-family:"dejavu sans mono",serif;color:blue'>+ </span></span><span class=m8358487097371014076gmail-gghfmyibcob><span style='font-family:"dejavu sans mono",serif;color:blue'> factor(Amarela),factor(SemReligiao),factor(Catolica),<o:p></o:p></span></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span class=m8358487097371014076gmail-gghfmyibcpb><span style='font-family:"dejavu sans mono",serif;color:blue'>+ </span></span><span class=m8358487097371014076gmail-gghfmyibcob><span style='font-family:"dejavu sans mono",serif;color:blue'> factor(Espirita),factor(Evangelica),factor(AfroBrasileira),<o:p></o:p></span></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span class=m8358487097371014076gmail-gghfmyibcpb><span style='font-family:"dejavu sans mono",serif;color:blue'>+ </span></span><span class=m8358487097371014076gmail-gghfmyibcob><span style='font-family:"dejavu sans mono",serif;color:blue'> factor(Turno),factor(Aposentado),factor(OcupaEstDiApenasDesemp),<o:p></o:p></span></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span class=m8358487097371014076gmail-gghfmyibcpb><span style='font-family:"dejavu sans mono",serif;color:blue'>+ </span></span><span class=m8358487097371014076gmail-gghfmyibcob><span style='font-family:"dejavu sans mono",serif;color:blue'> factor(ComFamilia),factor(ComOutParentes),factor(Republica),<o:p></o:p></span></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span class=m8358487097371014076gmail-gghfmyibcpb><span style='font-family:"dejavu sans mono",serif;color:blue'>+ </span></span><span class=m8358487097371014076gmail-gghfmyibcob><span style='font-family:"dejavu sans mono",serif;color:blue'> factor(Sozinho),factor(Pensao),factor(OutroMoradia),factor(RU),<o:p></o:p></span></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span class=m8358487097371014076gmail-gghfmyibcpb><span style='font-family:"dejavu sans mono",serif;color:blue'>+ </span></span><span class=m8358487097371014076gmail-gghfmyibcob><span style='font-family:"dejavu sans mono",serif;color:blue'> factor(praec4),IdadeA,escola,RendaPC,Dist,PraecSoma)<o:p></o:p></span></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span class=m8358487097371014076gmail-gghfmyibcpb><span style='font-family:"dejavu sans mono",serif;color:blue'>> </span></span><span class=m8358487097371014076gmail-gghfmyibcob><span style='font-family:"dejavu sans mono",serif;color:blue'>ModeloLogBinomial <- glm(Tabagismo~.,data = TabModelagem,family = binomial(link = log))<o:p></o:p></span></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span class=m8358487097371014076gmail-gghfmyibgob><span style='font-family:"dejavu sans mono",serif;color:#C5060B'>Erro: nenhum jogo válido de coeficientes tem sido encontrado: por favor, fornece valores iniciais<o:p></o:p></span></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span class=m8358487097371014076gmail-gghfmyibcpb><span style='font-family:"dejavu sans mono",serif;color:blue'>> </span></span><span class=m8358487097371014076gmail-gghfmyibcob><span style='font-family:"dejavu sans mono",serif;color:blue'>ModeloLogistico <- glm(Tabagismo~.,data = TabModelagem,family = binomial(link = logit))<o:p></o:p></span></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span class=m8358487097371014076gmail-gghfmyibcpb><span style='font-family:"dejavu sans mono",serif;color:blue'>> </span></span><span class=m8358487097371014076gmail-gghfmyibcob><span style='font-family:"dejavu sans mono",serif;color:blue'>summary(ModeloLogistico)<o:p></o:p></span></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'><o:p> </o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>Call:<o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>glm(formula = Tabagismo ~ ., family = binomial(link = logit), <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'> data = TabModelagem)<o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'><o:p> </o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>Deviance Residuals: <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'> Min 1Q Median 3Q Max <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>-2.0034 -0.8824 -0.5459 1.0052 2.3721 <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'><o:p> </o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>Coefficients:<o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'> Estimate Std. Error z value Pr(>|z|) <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>(Intercept) 1.215e+01 8.827e+02 0.014 0.989021 <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.SexoDic.1 1.042e+00 1.430e-01 7.281 3.31e-13 ***<o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.Branca.1 -1.481e+01 8.827e+02 -0.017 0.986618 <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.Negra.1 -1.474e+01 8.827e+02 -0.017 0.986677 <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.Parda.1 -1.487e+01 8.827e+02 -0.017 0.986563 <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.Amarela.1 -1.545e+01 8.827e+02 -0.018 0.986037 <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.SemReligiao.1 -1.956e-01 4.151e-01 -0.471 0.637543 <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.Catolica.1 -1.009e+00 4.026e-01 -2.506 0.012225 * <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.Espirita.1 -3.190e-01 4.624e-01 -0.690 0.490209 <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.Evangelica.1 -1.759e+00 4.694e-01 -3.747 0.000179 ***<o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.AfroBrasileira.1 1.438e+01 4.219e+02 0.034 0.972808 <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.Turno.1 3.093e-02 1.920e-01 0.161 0.872056 <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.Aposentado.1 -9.971e-02 1.901e-01 -0.524 0.600027 <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.OcupaEstDiApenasDesemp.1 1.799e-01 2.670e-01 0.674 0.500536 <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.ComFamilia.1 -6.624e-01 4.098e-01 -1.616 0.106055 <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.ComOutParentes.1 -7.401e-01 5.805e-01 -1.275 0.202332 <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.Republica.1 6.300e-02 3.837e-01 0.164 0.869588 <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.Sozinho.1 -5.343e-01 4.299e-01 -1.243 0.213914 <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.Pensao.1 -1.324e+00 6.812e-01 -1.944 0.051938 . <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.OutroMoradia.1 -9.557e-01 5.861e-01 -1.631 0.102934 <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.RU.1 -3.592e-01 1.911e-01 -1.880 0.060152 . <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.praec4.1 -2.845e-01 4.985e-01 -0.571 0.568191 <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>IdadeA 7.917e-02 2.043e-02 3.875 0.000107 ***<o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>escola 1.200e-01 6.289e-02 1.909 0.056297 . <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>RendaPC 1.699e-04 4.093e-05 4.151 3.31e-05 ***<o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>Dist 1.634e-04 2.525e-04 0.647 0.517573 <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>PraecSoma 5.281e-02 5.452e-02 0.969 0.332692 <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>---<o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1<o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'><o:p> </o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>(Dispersion parameter for binomial family taken to be 1)<o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'><o:p> </o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'> Null deviance: 1516.6 on 1135 degrees of freedom<o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>Residual deviance: 1279.3 on 1109 degrees of freedom<o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'> (587 observations deleted due to missingness)<o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>AIC: 1333.3<o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'><o:p> </o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>Number of Fisher Scoring iterations: 13<o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'><o:p> </o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span class=m8358487097371014076gmail-gghfmyibcpb><span style='font-family:"dejavu sans mono",serif;color:blue'>> </span></span><span class=m8358487097371014076gmail-gghfmyibcob><span style='font-family:"dejavu sans mono",serif;color:blue'>ModeloPoisson <- glm(Tabagismo~.,data = TabModelagem,family = poisson(link = log))<o:p></o:p></span></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span class=m8358487097371014076gmail-gghfmyibcpb><span style='font-family:"dejavu sans mono",serif;color:blue'>> </span></span><span class=m8358487097371014076gmail-gghfmyibcob><span style='font-family:"dejavu sans mono",serif;color:blue'>summary(ModeloPoisson)<o:p></o:p></span></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'><o:p> </o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>Call:<o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>glm(formula = Tabagismo ~ ., family = poisson(link = log), data = TabModelagem)<o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'><o:p> </o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>Deviance Residuals: <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'> Min 1Q Median 3Q Max <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>-1.4867 -0.7821 -0.5889 0.5349 1.6624 <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'><o:p> </o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>Coefficients:<o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'> Estimate Std. Error z value Pr(>|z|) <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>(Intercept) -1.245e+00 1.125e+00 -1.107 0.26845 <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.SexoDic.1 5.800e-01 1.065e-01 5.447 5.11e-08 ***<o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.Branca.1 -8.332e-01 1.009e+00 -0.826 0.40870 <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.Negra.1 -8.210e-01 1.028e+00 -0.799 0.42446 <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.Parda.1 -9.009e-01 1.012e+00 -0.890 0.37337 <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.Amarela.1 -1.089e+00 1.092e+00 -0.998 0.31852 <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.SemReligiao.1 -9.670e-02 2.430e-01 -0.398 0.69062 <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.Catolica.1 -4.813e-01 2.396e-01 -2.009 0.04459 * <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.Espirita.1 -1.235e-01 2.806e-01 -0.440 0.65987 <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.Evangelica.1 -9.177e-01 3.126e-01 -2.936 0.00332 ** <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.AfroBrasileira.1 6.068e-01 5.538e-01 1.096 0.27325 <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.Turno.1 1.534e-03 1.331e-01 0.012 0.99081 <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.Aposentado.1 -4.516e-02 1.357e-01 -0.333 0.73937 <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.OcupaEstDiApenasDesemp.1 7.249e-02 1.816e-01 0.399 0.68972 <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.ComFamilia.1 -4.323e-01 2.739e-01 -1.578 0.11447 <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.ComOutParentes.1 -5.029e-01 4.527e-01 -1.111 0.26655 <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.Republica.1 8.985e-03 2.522e-01 0.036 0.97157 <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.Sozinho.1 -2.475e-01 2.878e-01 -0.860 0.38987 <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.Pensao.1 -8.439e-01 5.148e-01 -1.639 0.10117 <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.OutroMoradia.1 -5.262e-01 4.316e-01 -1.219 0.22275 <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.RU.1 -1.937e-01 1.363e-01 -1.421 0.15517 <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.praec4.1 -1.583e-01 3.432e-01 -0.461 0.64469 <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>IdadeA 3.787e-02 1.207e-02 3.136 0.00171 ** <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>escola 8.576e-02 4.429e-02 1.936 0.05281 . <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>RendaPC 4.045e-05 1.690e-05 2.393 0.01670 * <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>Dist 2.605e-05 1.668e-04 0.156 0.87586 <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>PraecSoma 2.419e-02 3.972e-02 0.609 0.54262 <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>---<o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1<o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'><o:p> </o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>(Dispersion parameter for poisson family taken to be 1)<o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'><o:p> </o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'> Null deviance: 834.67 on 1135 degrees of freedom<o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>Residual deviance: 706.16 on 1109 degrees of freedom<o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'> (587 observations deleted due to missingness)<o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>AIC: 1640.2<o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'><o:p> </o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>Number of Fisher Scoring iterations: 5<o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'><o:p> </o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span class=m8358487097371014076gmail-gghfmyibcpb><span style='font-family:"dejavu sans mono",serif;color:blue'>> </span></span><span class=m8358487097371014076gmail-gghfmyibcob><span style='font-family:"dejavu sans mono",serif;color:blue'>ModeloQuasiPoisson <- glm(Tabagismo~.,data = TabModelagem,family = quasipoisson(link = log))<o:p></o:p></span></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span class=m8358487097371014076gmail-gghfmyibcpb><span style='font-family:"dejavu sans mono",serif;color:blue'>> </span></span><span class=m8358487097371014076gmail-gghfmyibcob><span style='font-family:"dejavu sans mono",serif;color:blue'>summary(ModeloQuasiPoisson)<o:p></o:p></span></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'><o:p> </o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>Call:<o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>glm(formula = Tabagismo ~ ., family = quasipoisson(link = log), <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'> data = TabModelagem)<o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'><o:p> </o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>Deviance Residuals: <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'> Min 1Q Median 3Q Max <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>-1.4867 -0.7821 -0.5889 0.5349 1.6624 <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'><o:p> </o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>Coefficients:<o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'> Estimate Std. Error t value Pr(>|t|) <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>(Intercept) -1.245e+00 8.738e-01 -1.424 0.154644 <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.SexoDic.1 5.800e-01 8.273e-02 7.011 4.11e-12 ***<o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.Branca.1 -8.332e-01 7.836e-01 -1.063 0.287863 <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.Negra.1 -8.210e-01 7.987e-01 -1.028 0.304185 <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.Parda.1 -9.009e-01 7.863e-01 -1.146 0.252163 <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.Amarela.1 -1.089e+00 8.481e-01 -1.284 0.199466 <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.SemReligiao.1 -9.670e-02 1.888e-01 -0.512 0.608566 <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.Catolica.1 -4.813e-01 1.862e-01 -2.585 0.009863 ** <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.Espirita.1 -1.235e-01 2.181e-01 -0.566 0.571230 <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.Evangelica.1 -9.177e-01 2.429e-01 -3.779 0.000166 ***<o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.AfroBrasileira.1 6.068e-01 4.303e-01 1.410 0.158794 <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.Turno.1 1.534e-03 1.034e-01 0.015 0.988169 <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.Aposentado.1 -4.516e-02 1.055e-01 -0.428 0.668597 <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.OcupaEstDiApenasDesemp.1 7.249e-02 1.411e-01 0.514 0.607474 <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.ComFamilia.1 -4.323e-01 2.128e-01 -2.031 0.042444 * <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.ComOutParentes.1 -5.029e-01 3.517e-01 -1.430 0.153011 <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.Republica.1 8.985e-03 1.959e-01 0.046 0.963429 <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.Sozinho.1 -2.475e-01 2.236e-01 -1.107 0.268673 <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.Pensao.1 -8.439e-01 4.000e-01 -2.110 0.035106 * <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.OutroMoradia.1 -5.262e-01 3.353e-01 -1.569 0.116880 <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.RU.1 -1.937e-01 1.059e-01 -1.830 0.067589 . <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.praec4.1 -1.583e-01 2.666e-01 -0.594 0.552951 <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>IdadeA 3.787e-02 9.381e-03 4.037 5.79e-05 ***<o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>escola 8.576e-02 3.441e-02 2.492 0.012836 * <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>RendaPC 4.045e-05 1.313e-05 3.080 0.002119 ** <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>Dist 2.605e-05 1.296e-04 0.201 0.840689 <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>PraecSoma 2.419e-02 3.086e-02 0.784 0.433427 <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>---<o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1<o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'><o:p> </o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>(Dispersion parameter for quasipoisson family taken to be 0.6036898)<o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'><o:p> </o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'> Null deviance: 834.67 on 1135 degrees of freedom<o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>Residual deviance: 706.16 on 1109 degrees of freedom<o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'> (587 observations deleted due to missingness)<o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>AIC: NA<o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'><o:p> </o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>Number of Fisher Scoring iterations: 5<o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'><o:p> </o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span class=m8358487097371014076gmail-gghfmyibcpb><span style='font-family:"dejavu sans mono",serif;color:blue'>> </span></span><span class=m8358487097371014076gmail-gghfmyibcob><span style='font-family:"dejavu sans mono",serif;color:blue'>library(MASS)<o:p></o:p></span></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span class=m8358487097371014076gmail-gghfmyibcpb><span style='font-family:"dejavu sans mono",serif;color:blue'>> </span></span><span class=m8358487097371014076gmail-gghfmyibcob><span style='font-family:"dejavu sans mono",serif;color:blue'>ModeloBinomialNegativa <- glm.nb(Tabagismo~.,data = TabModelagem,link = log)<o:p></o:p></span></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span class=m8358487097371014076gmail-gghfmyibgob><span style='font-family:"dejavu sans mono",serif;color:#C5060B'>Warning messages:<o:p></o:p></span></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span class=m8358487097371014076gmail-gghfmyibgob><span style='font-family:"dejavu sans mono",serif;color:#C5060B'>1: In <a href="http://theta.ml" target="_blank">theta.ml</a>(Y, mu, sum(w), w, limit = control$maxit, trace = control$trace > :<o:p></o:p></span></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span class=m8358487097371014076gmail-gghfmyibgob><span style='font-family:"dejavu sans mono",serif;color:#C5060B'> iteration limit reached<o:p></o:p></span></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span class=m8358487097371014076gmail-gghfmyibgob><span style='font-family:"dejavu sans mono",serif;color:#C5060B'>2: In <a href="http://theta.ml" target="_blank">theta.ml</a>(Y, mu, sum(w), w, limit = control$maxit, trace = control$trace > :<o:p></o:p></span></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span class=m8358487097371014076gmail-gghfmyibgob><span style='font-family:"dejavu sans mono",serif;color:#C5060B'> iteration limit reached<o:p></o:p></span></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span class=m8358487097371014076gmail-gghfmyibcpb><span style='font-family:"dejavu sans mono",serif;color:blue'>> </span></span><span class=m8358487097371014076gmail-gghfmyibcob><span style='font-family:"dejavu sans mono",serif;color:blue'>summary(ModeloBinomialNegativa)<o:p></o:p></span></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'><o:p> </o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>Call:<o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>glm.nb(formula = Tabagismo ~ ., data = TabModelagem, link = log, <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'> init.theta = 9955.862378)<o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'><o:p> </o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>Deviance Residuals: <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'> Min 1Q Median 3Q Max <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>-1.4866 -0.7821 -0.5889 0.5349 1.6624 <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'><o:p> </o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>Coefficients:<o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'> Estimate Std. Error z value Pr(>|z|) <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>(Intercept) -1.245e+00 1.125e+00 -1.107 0.26847 <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.SexoDic.1 5.800e-01 1.065e-01 5.447 5.12e-08 ***<o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.Branca.1 -8.332e-01 1.009e+00 -0.826 0.40872 <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.Negra.1 -8.210e-01 1.028e+00 -0.799 0.42448 <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.Parda.1 -9.009e-01 1.012e+00 -0.890 0.37339 <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.Amarela.1 -1.089e+00 1.092e+00 -0.997 0.31853 <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.SemReligiao.1 -9.670e-02 2.430e-01 -0.398 0.69063 <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.Catolica.1 -4.813e-01 2.396e-01 -2.008 0.04459 * <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.Espirita.1 -1.235e-01 2.807e-01 -0.440 0.65989 <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.Evangelica.1 -9.177e-01 3.126e-01 -2.936 0.00333 ** <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.AfroBrasileira.1 6.068e-01 5.539e-01 1.096 0.27327 <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.Turno.1 1.530e-03 1.331e-01 0.011 0.99083 <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.Aposentado.1 -4.516e-02 1.358e-01 -0.333 0.73937 <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.OcupaEstDiApenasDesemp.1 7.249e-02 1.816e-01 0.399 0.68973 <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.ComFamilia.1 -4.323e-01 2.739e-01 -1.578 0.11449 <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.ComOutParentes.1 -5.029e-01 4.527e-01 -1.111 0.26658 <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.Republica.1 8.986e-03 2.522e-01 0.036 0.97157 <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.Sozinho.1 -2.475e-01 2.878e-01 -0.860 0.38988 <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.Pensao.1 -8.439e-01 5.149e-01 -1.639 0.10120 <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.OutroMoradia.1 -5.262e-01 4.316e-01 -1.219 0.22277 <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.RU.1 -1.937e-01 1.363e-01 -1.421 0.15519 <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.praec4.1 -1.583e-01 3.432e-01 -0.461 0.64469 <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>IdadeA 3.787e-02 1.207e-02 3.136 0.00171 ** <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>escola 8.576e-02 4.429e-02 1.936 0.05282 . <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>RendaPC 4.045e-05 1.690e-05 2.393 0.01670 * <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>Dist 2.605e-05 1.668e-04 0.156 0.87585 <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>PraecSoma 2.419e-02 3.972e-02 0.609 0.54264 <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>---<o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1<o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'><o:p> </o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>(Dispersion parameter for Negative Binomial(9955.862) family taken to be 1)<o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'><o:p> </o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'> Null deviance: 834.65 on 1135 degrees of freedom<o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>Residual deviance: 706.13 on 1109 degrees of freedom<o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'> (587 observations deleted due to missingness)<o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>AIC: 1642.2<o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'><o:p> </o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>Number of Fisher Scoring iterations: 1<o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'><o:p> </o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'><o:p> </o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'> Theta: 9956 <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'> Std. Err.: 36769 <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>Warning while fitting theta: iteration limit reached <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'><o:p> </o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'> 2 x log-likelihood: -1586.177 <o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span class=m8358487097371014076gmail-gghfmyibcpb><span style='font-family:"dejavu sans mono",serif;color:blue'>> </span></span><span class=m8358487097371014076gmail-gghfmyibcob><span style='font-family:"dejavu sans mono",serif;color:blue'>par(mfrow=c(2,2))<o:p></o:p></span></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span class=m8358487097371014076gmail-gghfmyibcpb><span style='font-family:"dejavu sans mono",serif;color:blue'>> </span></span><span class=m8358487097371014076gmail-gghfmyibcob><span style='font-family:"dejavu sans mono",serif;color:blue'>qqnorm(resid(ModeloLogistico, type = "deviance"),<o:p></o:p></span></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span class=m8358487097371014076gmail-gghfmyibcpb><span style='font-family:"dejavu sans mono",serif;color:blue'>+ </span></span><span class=m8358487097371014076gmail-gghfmyibcob><span style='font-family:"dejavu sans mono",serif;color:blue'> pch = 20, main = "Logística", las = 1)<o:p></o:p></span></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span class=m8358487097371014076gmail-gghfmyibcpb><span style='font-family:"dejavu sans mono",serif;color:blue'>> </span></span><span class=m8358487097371014076gmail-gghfmyibcob><span style='font-family:"dejavu sans mono",serif;color:blue'>qqline(resid(ModeloLogistico, type = "deviance"))<o:p></o:p></span></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span class=m8358487097371014076gmail-gghfmyibcpb><span style='font-family:"dejavu sans mono",serif;color:blue'>> </span></span><span class=m8358487097371014076gmail-gghfmyibcob><span style='font-family:"dejavu sans mono",serif;color:blue'>qqnorm(resid(ModeloPoisson, type = "deviance"),<o:p></o:p></span></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span class=m8358487097371014076gmail-gghfmyibcpb><span style='font-family:"dejavu sans mono",serif;color:blue'>+ </span></span><span class=m8358487097371014076gmail-gghfmyibcob><span style='font-family:"dejavu sans mono",serif;color:blue'> pch = 20, main = "Poisson", las = 1)<o:p></o:p></span></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span class=m8358487097371014076gmail-gghfmyibcpb><span style='font-family:"dejavu sans mono",serif;color:blue'>> </span></span><span class=m8358487097371014076gmail-gghfmyibcob><span style='font-family:"dejavu sans mono",serif;color:blue'>qqline(resid(ModeloPoisson, type = "deviance"))<o:p></o:p></span></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span class=m8358487097371014076gmail-gghfmyibcpb><span style='font-family:"dejavu sans mono",serif;color:blue'>> </span></span><span class=m8358487097371014076gmail-gghfmyibcob><span style='font-family:"dejavu sans mono",serif;color:blue'>qqnorm(resid(ModeloQuasiPoisson, type = "deviance"),<o:p></o:p></span></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span class=m8358487097371014076gmail-gghfmyibcpb><span style='font-family:"dejavu sans mono",serif;color:blue'>+ </span></span><span class=m8358487097371014076gmail-gghfmyibcob><span style='font-family:"dejavu sans mono",serif;color:blue'> pch = 20, main = "Quasi-Poisson", las = 1)<o:p></o:p></span></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span class=m8358487097371014076gmail-gghfmyibcpb><span style='font-family:"dejavu sans mono",serif;color:blue'>> </span></span><span class=m8358487097371014076gmail-gghfmyibcob><span style='font-family:"dejavu sans mono",serif;color:blue'>qqline(resid(ModeloQuasiPoisson, type = "deviance"))<o:p></o:p></span></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span class=m8358487097371014076gmail-gghfmyibcpb><span style='font-family:"dejavu sans mono",serif;color:blue'>> </span></span><span class=m8358487097371014076gmail-gghfmyibcob><span style='font-family:"dejavu sans mono",serif;color:blue'>qqnorm(resid(ModeloBinomialNegativa, type = "deviance"),<o:p></o:p></span></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span class=m8358487097371014076gmail-gghfmyibcpb><span style='font-family:"dejavu sans mono",serif;color:blue'>+ </span></span><span class=m8358487097371014076gmail-gghfmyibcob><span style='font-family:"dejavu sans mono",serif;color:blue'> pch = 20, main = "Binomial Negativa", las = 1)<o:p></o:p></span></span></pre><pre style='line-height:12.0pt;word-break:break-all'><span class=m8358487097371014076gmail-gghfmyibcpb><span style='font-family:"dejavu sans mono",serif;color:blue'>> </span></span><span class=m8358487097371014076gmail-gghfmyibcob><span style='font-family:"dejavu sans mono",serif;color:blue'>qqline(resid(ModeloBinomialNegativa, type = "deviance"))</span></span><span style='font-family:"dejavu sans mono",serif;color:black'><o:p></o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all;outline:none;white-space:pre-wrap' id="m_8358487097371014076gmail-rstudio_console_output"><span style='font-family:"dejavu sans mono",serif;color:black'><o:p> </o:p></span></pre><pre style='line-height:12.0pt;word-break:break-all;outline:none;white-space:pre-wrap' id="m_8358487097371014076gmail-rstudio_console_output"><span style='font-family:"dejavu sans mono",serif;color:blue'><img border=0 width=562 height=391 style='width:5.8541in;height:4.0729in' id="_x0000_i1025" src="cid:image003.png@01D282C6.0FA4F4D0" alt="Imagem inline 1"></span><span style='font-family:"dejavu sans mono",serif;color:black'><o:p></o:p></span></pre></div><div><div><div><p class=MsoNormal><o:p> </o:p></p><div><p class=MsoNormal>Em 8 de fevereiro de 2017 20:35, Leonard Mendonça de Assis <<a href="mailto:assis.leonard@gmail.com" target="_blank">assis.leonard@gmail.com</a>> escreveu:<o:p></o:p></p><blockquote style='border:none;border-left:solid #CCCCCC 1.0pt;padding:0cm 0cm 0cm 6.0pt;margin-left:4.8pt;margin-right:0cm'><div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'><span style='font-size:11.0pt;font-family:"Calibri",sans-serif'>Exato,</span><o:p></o:p></p><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'><span style='font-size:11.0pt;font-family:"Calibri",sans-serif'> </span><o:p></o:p></p><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'><span style='font-size:11.0pt;font-family:"Calibri",sans-serif'>Existem várias funções de ligação pre definidas para cada família. Não tenho certeza se ‘log’ é um link válido para binomial, nunca o usei com binomial.</span><o:p></o:p></p><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'><span style='font-size:11.0pt;font-family:"Calibri",sans-serif'> </span><o:p></o:p></p><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'><span style='font-size:11.0pt;font-family:"Calibri",sans-serif'>Caso este link não exista, você terá que cria-lo antes, veja um exemplo neste local: <a href="https://stat.ethz.ch/pipermail/r-help/2013-November/362787.html" target="_blank">https://stat.ethz.ch/pipermail/r-help/2013-November/362787.html</a></span><o:p></o:p></p><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'><span style='font-size:11.0pt;font-family:"Calibri",sans-serif'> </span><o:p></o:p></p><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'><span style='font-size:11.0pt;font-family:"Calibri",sans-serif'>Neste exemplo, é criado o link ‘clog’, muito próximo ao que você precisa.</span><o:p></o:p></p><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'><span style='font-size:11.0pt;font-family:"Calibri",sans-serif'> </span><o:p></o:p></p><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'><span style='font-size:11.0pt;font-family:"Calibri",sans-serif'>Mas o que você chama de ‘não consigo criar o modelo’? Eu já ajustei modelo logístico com mais de 200 variáveis (Contínuas, discretas, fatores ...) sem problemas, a não ser problemas de algumas variáveis serem linearmente dependentes e, com isto, me gerou problemas de estimação. Outro problema que tive (não o tenho a mais de ano) foi estourar memória do computador.</span><o:p></o:p></p><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'><span style='font-size:11.0pt;font-family:"Calibri",sans-serif'> </span><o:p></o:p></p><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'><span style='font-size:11.0pt;font-family:"Calibri",sans-serif'>Mande pra mim o erro especifico que está dando ao executar o comando. Só assim posso ser mais preciso em te ajudar</span><o:p></o:p></p><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'><span style='font-size:11.0pt;font-family:"Calibri",sans-serif'> </span><o:p></o:p></p><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'><b><span style='font-size:11.0pt;font-family:"Calibri",sans-serif'>De:</span></b><span style='font-size:11.0pt;font-family:"Calibri",sans-serif'> Marcos Bissoli [mailto:<a href="mailto:mbissoli@gmail.com" target="_blank">mbissoli@gmail.com</a>] <br><b>Enviada em:</b> quarta-feira, 8 de fevereiro de 2017 19:35<br><b>Para:</b> Leonard Mendonça de Assis <<a href="mailto:assis.leonard@gmail.com" target="_blank">assis.leonard@gmail.com</a>></span><o:p></o:p></p><div><div><p class=MsoNormal><br><b>Assunto:</b> Re: [R-br] Distribuição para regressão de resposta binária<o:p></o:p></p></div></div><div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'> <o:p></o:p></p><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>Prezado Leonard e amigos,<o:p></o:p></p><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'> <o:p></o:p></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>Agradeço muito pelo debate. Venho tentando cada vez mais dialogar com estatísticos, pois respeito muito o trabalho de vocês, embora admita ainda ser "um menino" na arte.<o:p></o:p></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'> <o:p></o:p></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>Em relação à função de ligação para logística, talvez eu não tenha sido claro. Eu tentei uma função "ln(y)" que, ao menos no material que venho consultando, seria uma regressão log-binomial. Este seria um modelo ideal, e foi minha primeira tentativa. Ou seja, usei um código semelhante a:<o:p></o:p></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'> <o:p></o:p></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>> Modelo -> glm(y~., data = Dados, family = binomial(link = "log"))<o:p></o:p></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'> <o:p></o:p></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>O problema é que o R não consegue criar tal modelo. Tenho muitas variáveis explicativas (isso é bastante comum em estudos epidemiológicos descritivos), incluindo cinco contínuas, se é que esse seja o motivo. O fato é que em uma das referências que citei em e-mail anterior, os autores tratam deste problema. Veja os resultados apresentados no resumo de Coutinho et al:<o:p></o:p></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'> <o:p></o:p></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>"<b><span style='font-family:"Verdana",sans-serif;color:black'>RESULTADOS:</span></b><span style='font-family:"Verdana",sans-serif;color:black'> As estimativas por ponto e por intervalo [das razões de prevalência] obtidas pelas regressões de Cox e Poisson foram semelhantes à obtida pela estratificação de Mantel-Haenszel [considerada 'prova-ouro' para a Epidemiologia], independentemente da prevalência do desfecho [variável resposta] e das covariáveis [variáveis explicativas, pode-se dizer; talvez, para que tem uma formação mais voltada para análises experimentais, poderíamos dizer que covariáveis referem-se mais a 'blocos'] do modelo. <b>O modelo log-binomial apresentou dificuldade de convergência</b> quando o desfecho tinha prevalência alta e havia covariável contínua no modelo. A regressão logística [valendo-se de logito como função de ligação] produziu estimativas por ponto e por intervalo maiores do que as obtidas pelos outros métodos, principalmente para os desfechos com maiores prevalências iniciais. <b>Se interpretados como estimativas de RP, os OR superestimariam as associações</b> para os desfechos com prevalência inicial baixa, intermediária e alta em 13%, quase 100% e quatro vezes mais, respectivamente."</span><o:p></o:p></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'> <o:p></o:p></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>[notas minhas] [grifos meus]<o:p></o:p></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'> <o:p></o:p></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>Portanto, os autores indicam (e vi isso em outros artigos) regressão de Poisson, mesmo admitindo a variável resposta como sendo binária, variando de 0 a 1, que representa a probabilidade de ocorrência do desfecho (doença). Quase todos são unânimes em recomendar apenas que se use ajuste de variância robusta para sanar problemas nos intervalos de confiança dos coeficientes. Em Epidemiologia, mais importante que os valores p são esses intervalos de confiança, pois há muitos desdobramentos inferenciais que são feitos a partir deles. Portanto, creio que a justificativa para adoção de Poisson seja esta: a não convergência da log-binomial.<o:p></o:p></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'> <o:p></o:p></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>Talvez uma outra função de ligação em família binomial possa ser uma solução, então? Como disse, tentei "log" e o próprio "logit". "Log" não deu convergência e o modelo nem foi gerado. O "logit" eu também tentei, e o gráfico de resíduos do modelo foi praticamente idêntico a este de quasi-poisson que postei na primeira mensagem. Seria grato caso pudesse me indicar algum referencial sobre outras funções de ligação, preferencialmente com aplicações. Mas acho que isso se tornará um problema, pois os coeficientes gerados certamente me resultarão indicadores não reconhecidos na área da Epidemiologia. Não sei até que ponto eu posso "converter" coeficientes livremente aplicando pura e simplesmente aritméticas a estes coeficientes. Como acho que já expliquei acima: "logit" me devolve razão de chances (odds ratio, OR) e "log" me retorna a tão desejada razão de prevalência (razão de riscos, RR). Se eu não tiver como converter meus coeficientes em uma dessas razões (e preferencialmente a RR) eu vou apanhar tanto da banca que vou sair dali roxo e sem título. :D<o:p></o:p></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'> <o:p></o:p></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>(Um parênteses: essas questões, como a brincadeira acima, tem me motivado muito a ingressar em um doutorado em Estatística assim que eu concluir esse. Creio que há muito há se propor de novo para a Epidemiologia, a partir de um conhecimento mais profundo em Estatística. A Epidemiologia é, de fato, bastante limitada naquilo que ela "aceita" como técnica válida para suas análises inferenciais)<o:p></o:p></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'> <o:p></o:p></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>Quando o senhor diz:<o:p></o:p></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'> <o:p></o:p></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>"<span style='font-size:11.0pt;font-family:"Calibri",sans-serif'>Uma segunda forma de analisar, seria termos uma resposta composta de número de ocorrências do evento em um total possível",</span><o:p></o:p></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'> <o:p></o:p></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>eu posso interpretar que isso também seria adaptável a uma probabilidade? Valeria a pena eu tentar ajuste com Gamma ou Beta? A binomial negativa eu até tentei, mais por curiosidade, pois tenho seguido o material do curso que fiz com o Walmes em Varginha e ele preconizou que ela seria recomendada para casos de superdispersão, e minha variável resposta aparenta subdispersão. Mas a binomial negativa também gerou o mesmo gráfico de resíduos. Eu posso usar uma variável binária como resposta num modelo Gamma ou Beta? Ou teria que dar algum "tratamento" na variável antes de aplicar estes modelos. Confesso que estes ainda não tentei com essa resposta com a qual estou enfrentando o problema. Em outras variáveis eu experimentei modelos Gamma, mas eles não apresentaram melhor ajuste que o gaussiano, a ponto de justificar eu enfrentar tamanha novidade com a banca de epidemiologistas. :D<o:p></o:p></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'> <o:p></o:p></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>Agradeço, já, e muito, o diálogo estabelecido.<o:p></o:p></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'> <o:p></o:p></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>Abraços fraternos,<o:p></o:p></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'> <o:p></o:p></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>Marcos<o:p></o:p></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'> <o:p></o:p></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'> <o:p></o:p></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'> <o:p></o:p></p></div></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'> <o:p></o:p></p><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>Em 8 de fevereiro de 2017 17:28, Leonard Mendonça de Assis <<a href="mailto:assis.leonard@gmail.com" target="_blank">assis.leonard@gmail.com</a>> escreveu:<o:p></o:p></p><blockquote style='border:none;border-left:solid #CCCCCC 1.0pt;padding:0cm 0cm 0cm 6.0pt;margin-left:4.8pt;margin-top:5.0pt;margin-right:0cm;margin-bottom:5.0pt'><div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'><span style='font-size:11.0pt;font-family:"Calibri",sans-serif'>Marcos,</span><o:p></o:p></p><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'><span style='font-size:11.0pt;font-family:"Calibri",sans-serif'>Não tenho acesso aos artigos, mas ... baseado em 27 anos de experiência ajustando regressão logística, vamos aos meus pitacos:</span><o:p></o:p></p><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'><span style='font-size:11.0pt;font-family:"Calibri",sans-serif'> </span><o:p></o:p></p><ol start=1 type=1><li class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:0cm;mso-list:l0 level1 lfo1'><span style='font-size:11.0pt;font-family:"Calibri",sans-serif'>Quando eu tenho os Y em forma binária (Presente/ausente), isto, estatisticamente falando, é uma distribuição de Bernoulli</span><o:p></o:p></li><li class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:0cm;mso-list:l0 level1 lfo1'><span style='font-size:11.0pt;font-family:"Calibri",sans-serif'>Esta distribuição de bernoulli tem como parâmetro, a proporção. Esta proporção varia de 0 a 1.</span><o:p></o:p></li><li class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:0cm;mso-list:l0 level1 lfo1'><span style='font-size:11.0pt;font-family:"Calibri",sans-serif'>Uma forma de ajustar este tipo de dados é a regressão logística, esta pode assumir vários tipos de ligação. Aqui, de cabeça, eu lembro uns 5, mas deve ter muito mais.</span><o:p></o:p></li><li class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:0cm;mso-list:l0 level1 lfo1'><span style='font-size:11.0pt;font-family:"Calibri",sans-serif'>Estas funções de ligação se ajustam melhor a determinados tipos de dados e algumas áreas de conhecimento às vezes preferem um em detrimento de outros.</span><o:p></o:p></li></ol><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'><span style='font-size:11.0pt;font-family:"Calibri",sans-serif'> </span><o:p></o:p></p><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'><span style='font-size:11.0pt;font-family:"Calibri",sans-serif'>Baseado neste cenário acima, acho “estranho” utilizar outro tipo de modelo (com dados na característica acima), sem uma justificativa bastante forte.</span><o:p></o:p></p><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'><span style='font-size:11.0pt;font-family:"Calibri",sans-serif'> </span><o:p></o:p></p><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'><span style='font-size:11.0pt;font-family:"Calibri",sans-serif'>Uma segunda forma de analisar, seria termos uma resposta composta de número de ocorrências do evento em um total possível. Neste caso, teríamos uma variedade de distribuições que provavelmente se encaixariam nos documentos que você apresentou. Neste caso, o modelo seria ou binomial negativa, ou Gama, ou Beta, ou qualquer outra similar.</span><o:p></o:p></p><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'><span style='font-size:11.0pt;font-family:"Calibri",sans-serif'> </span><o:p></o:p></p><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'><span style='font-size:11.0pt;font-family:"Calibri",sans-serif'>Ajustados os conceitos, vamos agora à minha opinião sobre seu problema.</span><o:p></o:p></p><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'><span style='font-size:11.0pt;font-family:"Calibri",sans-serif'> </span><o:p></o:p></p><ol start=1 type=1><li class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:0cm;mso-list:l1 level1 lfo2'><span style='font-size:11.0pt;font-family:"Calibri",sans-serif'>Pelo que consegui entender de seu texto inicial e seu código fornecido, você está ajustando algo que é 0 ou 1 como sendo uma quase-poisson (esta é uma das opções de ajuste para o que expliquei acima, onde existe uma determinada quantidade.</span><o:p></o:p></li><li class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:0cm;mso-list:l1 level1 lfo2'><span style='font-size:11.0pt;font-family:"Calibri",sans-serif'>Como seus dados são 0/1 (suposição esta que faço baseado em sua explicação), o ajuste estar deficiente é algo bem esperado</span><o:p></o:p></li></ol><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'><span style='font-size:11.0pt;font-family:"Calibri",sans-serif'>Seria bom se você informasse do que se trata a variável tabagismo, se ela é 0/1 ou se é uma quantidade. Se for 0/1, certamente o problema dos dados é esperado, por serem oriundos de uma distribuição diversa da que você está ajustando.</span><o:p></o:p></p><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'><span style='font-size:11.0pt;font-family:"Calibri",sans-serif'> </span><o:p></o:p></p><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'><span style='font-size:11.0pt;font-family:"Calibri",sans-serif'>No aguardo</span><o:p></o:p></p><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'><span style='font-size:11.0pt;font-family:"Calibri",sans-serif'> </span><o:p></o:p></p><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'><span style='font-size:11.0pt;font-family:"Calibri",sans-serif'>Leonard</span><o:p></o:p></p><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'><span style='font-size:11.0pt;font-family:"Calibri",sans-serif'> </span><o:p></o:p></p><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'><b><span style='font-size:11.0pt;font-family:"Calibri",sans-serif'>De:</span></b><span style='font-size:11.0pt;font-family:"Calibri",sans-serif'> Marcos Bissoli [mailto:<a href="mailto:mbissoli@gmail.com" target="_blank">mbissoli@gmail.com</a>] <br><b>Enviada em:</b> quarta-feira, 8 de fevereiro de 2017 14:57<br><b>Para:</b> Leonard Assis <<a href="mailto:assis.leonard@gmail.com" target="_blank">assis.leonard@gmail.com</a>>; a lista Brasileira oficial de discussão do programa R. <<a href="mailto:r-br@listas.c3sl.ufpr.br" target="_blank">r-br@listas.c3sl.ufpr.br</a>><br><b>Assunto:</b> Re: [R-br] Distribuição para regressão de resposta binária</span><o:p></o:p></p><div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'> <o:p></o:p></p><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>Olá Leonard,<o:p></o:p></p><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'> <o:p></o:p></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>Muito obrigado pelo interesse no debate.<o:p></o:p></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'> <o:p></o:p></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>Minha afirmação baseada em epidemiologistas é, na verdade, baseada numa série de artigos que venho estudando recentemente. Seguem algumas referências, dentre outras várias, sobre as quais venho me fundamentando.<o:p></o:p></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'> <o:p></o:p></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>ZOU, G. A Modified Poisson Regression Approach to Prospective Studies with Binary Data. <b>American Journal of Epidemiology</b>, v. 159, n. 7, p. 702–706, 1 abr. 2004. <o:p></o:p></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'> <o:p></o:p></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>COUTINHO, L. M. S.; SCAZUFCA, M.; MENEZES, P. R. Métodos para estimar razão de prevalência em estudos de corte transversal. <b>Revista de Saúde Pública</b>, v. 42, n. 6, p. 992–998, dez. 2008. <o:p></o:p></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'> <o:p></o:p></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>BARROS, A. J.; HIRAKATA, V. N. Alternatives for logistic regression in cross-sectional studies: an empirical comparison of models that directly estimate the prevalence ratio. <b>BMC Medical Research Methodology</b>, v. 3, n. 1, p. 21, 20 dez. 2003. <o:p></o:p></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'> <o:p></o:p></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>FRANCISCO, P. M. S. B. et al. Medidas de associação em estudo transversal com delineamento complexo: razão de chances e razão de prevalência. <b>Revista Brasileira de Epidemiologia</b>, v. 11, n. 3, p. 347–355, set. 2008. <o:p></o:p></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'> <o:p></o:p></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>WILLIAMSON, T.; ELIASZIW, M.; FICK, G. H. Log-binomial models: exploring failed convergence. <b>Emerging themes in epidemiology</b>, v. 10, n. 1, p. 14, 13 dez. 2013. <o:p></o:p></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'> <o:p></o:p></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>Abraços fraternos,<o:p></o:p></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'> <o:p></o:p></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>Marcos<o:p></o:p></p></div></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'> <o:p></o:p></p><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>Em 7 de fevereiro de 2017 22:22, Leonard Assis via R-br <<a href="mailto:r-br@listas.c3sl.ufpr.br" target="_blank">r-br@listas.c3sl.ufpr.br</a>> escreveu:<o:p></o:p></p><blockquote style='border:none;border-left:solid #CCCCCC 1.0pt;padding:0cm 0cm 0cm 6.0pt;margin-left:4.8pt;margin-top:5.0pt;margin-right:0cm;margin-bottom:5.0pt'><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>Tem ruído aí nesta explicação. Na verdade, o que o "epidemiologista" alegou, não me convenceu.<o:p></o:p></p><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'> <o:p></o:p></p></div></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'> <o:p></o:p></p><div><div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>Em 7 de fev de 2017 9:14 PM, "Marcos Bissoli via R-br" <<a href="mailto:r-br@listas.c3sl.ufpr.br" target="_blank">r-br@listas.c3sl.ufpr.br</a>> escreveu:<o:p></o:p></p></div></div><blockquote style='border:none;border-left:solid #CCCCCC 1.0pt;padding:0cm 0cm 0cm 6.0pt;margin-left:4.8pt;margin-top:5.0pt;margin-right:0cm;margin-bottom:5.0pt'><div><div><div><div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>Prezados,<o:p></o:p></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'> <o:p></o:p></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>De antemão peço desculpas se desvio o tópico da lista. Mas creio que o tema da mensagem é minimamente transversal aos aqui tratados.<o:p></o:p></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'> <o:p></o:p></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>Tenho uma variável resposta binária. Como a frequência da resposta é alta (38,11%), teóricos da Estatística aplicada à Epidemiologia sugerem que não seja usada uma regressão logística. Neste caso (de alta prevalência do desfecho), a primeira opção deveria ser uma log-binomial. Mas (e isso não é raro de ocorrer), minha log-binomial não apresentou convergência.<o:p></o:p></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'> <o:p></o:p></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>Quando não há convergência, os teóricos sugerem uma regressão de Poisson com variância robusta. Entretanto, como meus dados sugerem subdispersão, optei por um modelo de quasi-poisson. Isso já deu certo em outras análises que fiz para terceiros. Inclusive, tenho conseguido adaptar a variância robusta ao modelo de quasi-poisson. Mas justamente agora, com os dados de minha tese...<o:p></o:p></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'> <o:p></o:p></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>O diagnóstico visual está, ao meu ver, péssimo, para ajuste. A imagem anexa é do modelo de quasi-poisson. Mas experimentei todos os acima citados (logística e Poisson) e o gráfico não diferiu muito.<o:p></o:p></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'> <o:p></o:p></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'><img border=0 width=544 height=379 style='width:5.6666in;height:3.9479in' id="m_8358487097371014076m_6541542930473827578m_6389802118341864428_x005f_x0000_i1025" src="cid:image004.png@01D282C6.0FA4F4D0" alt="Imagem inline 1"><o:p></o:p></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'> <o:p></o:p></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>A dúvida é... Há alguma outra alternativa de técnica de regressão que eu poderia tentar? Minhas variáveis explicativas são diversas, em quantidade e tipo (há contínuas, ordinais e binárias). Ou será (embora eu ache pouco provável) que este gráfico não significa um grande incômodo?<o:p></o:p></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'> <o:p></o:p></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>Fiz o teste de qui-quadrado da deviance residual e estranhamente o valor p está resultando em 1, tanto para Poisson quanto para quasi-Poisson. Um outro fato estranho é o pseudo R² de Nagelkerke ter acusado 20%: todas as outras minhas variáveis resposta não passaram de 12%. Não sei se é correto (consultei bibliografia que sugeria isso para a regressão logística), mas apliquei um teste de Hosmer e Lemeshow e ele acusou um bom ajuste do modelo, também (p = 0,2718). Até uma curva de ROC eu fiz e a área está grande no gráfico (mais uma técnica que não sei se deve ser aplicada além da regressão logística,).<o:p></o:p></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'> <o:p></o:p></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>Seguem alguns resultados, caso possa ajudar em algo.<o:p></o:p></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'> <o:p></o:p></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>Desde já agradeço qualquer comentário. E reforço minhas desculpas caso eu tenha desviado do tópico além do esperado, e desde já acato qualquer negativa em prosseguir o debate. Nesse caso, se possível, aceitaria sugestões de boas listas para debates nesse nível onde eu pudesse me inscrever.<o:p></o:p></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'> <o:p></o:p></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>Há braços,<o:p></o:p></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'> <o:p></o:p></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>Marcos Bissoli<o:p></o:p></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>Faculdade de Nutrição<o:p></o:p></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>Unifal-MG<o:p></o:p></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'> <o:p></o:p></p></div><pre style='line-height:12.0pt;word-break:break-all;outline:none;white-space:pre-wrap' id="m_8358487097371014076m_6541542930473827578m_6389802118341864428m_799196135501887637m_-6160994631685705737gmail-rstudio_console_output"><span class="m8358487097371014076m6541542930473827578m6389802118341864428m799196135501887637m-6160994631685705737gmail-gghfmyibcpb"><span style='font-family:"dejavu sans mono",serif;color:blue'>> </span></span><span class="m8358487097371014076m6541542930473827578m6389802118341864428m799196135501887637m-6160994631685705737gmail-gghfmyibcob"><span style='font-family:"dejavu sans mono",serif;color:blue'>Mod1 <- glm(Tabagismo~.,data = TabModelagem,family = quasipoisson)</span></span><o:p></o:p></pre><pre style='line-height:12.0pt;word-break:break-all'><span class="m8358487097371014076m6541542930473827578m6389802118341864428m799196135501887637m-6160994631685705737gmail-gghfmyibcpb"><span style='font-family:"dejavu sans mono",serif;color:blue'>> </span></span><span class="m8358487097371014076m6541542930473827578m6389802118341864428m799196135501887637m-6160994631685705737gmail-gghfmyibcob"><span style='font-family:"dejavu sans mono",serif;color:blue'>summary(Mod1)</span></span><o:p></o:p></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'> </span><o:p></o:p></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>Call:</span><o:p></o:p></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>glm(formula = Tabagismo ~ ., family = quasipoisson, data = TabModelagem)</span><o:p></o:p></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'> </span><o:p></o:p></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>Deviance Residuals: </span><o:p></o:p></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'> Min 1Q Median 3Q Max </span><o:p></o:p></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>-1.4867 -0.7821 -0.5889 0.5349 1.6624 </span><o:p></o:p></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'> </span><o:p></o:p></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>Coefficients:</span><o:p></o:p></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'> Estimate Std. Error t value Pr(>|t|) </span><o:p></o:p></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>(Intercept) -1.245e+00 8.738e-01 -1.424 0.154644 </span><o:p></o:p></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.SexoDic.1 5.800e-01 8.273e-02 7.011 4.11e-12 ***</span><o:p></o:p></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.Branca.1 -8.332e-01 7.836e-01 -1.063 0.287863 </span><o:p></o:p></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.Negra.1 -8.210e-01 7.987e-01 -1.028 0.304185 </span><o:p></o:p></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.Parda.1 -9.009e-01 7.863e-01 -1.146 0.252163 </span><o:p></o:p></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.Amarela.1 -1.089e+00 8.481e-01 -1.284 0.199466 </span><o:p></o:p></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.SemReligiao.1 -9.670e-02 1.888e-01 -0.512 0.608566 </span><o:p></o:p></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.Catolica.1 -4.813e-01 1.862e-01 -2.585 0.009863 ** </span><o:p></o:p></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.Espirita.1 -1.235e-01 2.181e-01 -0.566 0.571230 </span><o:p></o:p></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.Evangelica.1 -9.177e-01 2.429e-01 -3.779 0.000166 ***</span><o:p></o:p></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.AfroBrasileira.1 6.068e-01 4.303e-01 1.410 0.158794 </span><o:p></o:p></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.Turno.1 1.534e-03 1.034e-01 0.015 0.988169 </span><o:p></o:p></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.Aposentado.1 -4.516e-02 1.055e-01 -0.428 0.668597 </span><o:p></o:p></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.OcupaEstDiApenasDesemp.1 7.249e-02 1.411e-01 0.514 0.607474 </span><o:p></o:p></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.ComFamilia.1 -4.323e-01 2.128e-01 -2.031 0.042444 * </span><o:p></o:p></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.ComOutParentes.1 -5.029e-01 3.517e-01 -1.430 0.153011 </span><o:p></o:p></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.Republica.1 8.985e-03 1.959e-01 0.046 0.963429 </span><o:p></o:p></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.Sozinho.1 -2.475e-01 2.236e-01 -1.107 0.268673 </span><o:p></o:p></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.Pensao.1 -8.439e-01 4.000e-01 -2.110 0.035106 * </span><o:p></o:p></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.OutroMoradia.1 -5.262e-01 3.353e-01 -1.569 0.116880 </span><o:p></o:p></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.RU.1 -1.937e-01 1.059e-01 -1.830 0.067589 . </span><o:p></o:p></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>factor.praec4.1 -1.583e-01 2.666e-01 -0.594 0.552951 </span><o:p></o:p></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>IdadeA 3.787e-02 9.381e-03 4.037 5.79e-05 ***</span><o:p></o:p></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>escola 8.576e-02 3.441e-02 2.492 0.012836 * </span><o:p></o:p></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>RendaPC 4.045e-05 1.313e-05 3.080 0.002119 ** </span><o:p></o:p></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>Dist 2.605e-05 1.296e-04 0.201 0.840689 </span><o:p></o:p></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>PraecSoma 2.419e-02 3.086e-02 0.784 0.433427 </span><o:p></o:p></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>---</span><o:p></o:p></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1</span><o:p></o:p></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'> </span><o:p></o:p></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>(Dispersion parameter for quasipoisson family taken to be 0.6036898)</span><o:p></o:p></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'> </span><o:p></o:p></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'> Null deviance: 834.67 on 1135 degrees of freedom</span><o:p></o:p></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>Residual deviance: 706.16 on 1109 degrees of freedom</span><o:p></o:p></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>AIC: NA</span><o:p></o:p></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'> </span><o:p></o:p></pre><pre style='line-height:12.0pt;word-break:break-all'><span style='font-family:"dejavu sans mono",serif;color:black'>Number of Fisher Scoring iterations: 5</span><o:p></o:p></pre><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'> <o:p></o:p></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'> <o:p></o:p></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'> <o:p></o:p></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>-- <o:p></o:p></p></div><div><div><div><div><p class=MsoNormal style='mso-margin-top-alt:auto;margin-bottom:12.0pt'>MARCOS BISSOLI<o:p></o:p></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>Faculdade de Nutrição<br>Universidade Federal de Alfenas<br><br>Blog: <span style='font-size:9.5pt'><a href="http://bocademiamaldita.blogspot.com/" target="_blank">bocademiamaldita.blogspot.com/</a></span><o:p></o:p></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>E-mail: <a href="mailto:mbissoli@gmail.com" target="_blank">mbissoli@gmail.com</a><br>Twitter: #mbissoli<o:p></o:p></p><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'><br>Alfenas, Minas Gerais, Brasil<br><br><br>*****Pense na Natureza antes de Imprimir*****<br>Divulgue ON-LINE<br><br>Eu apoio a ENEN "na luta por um Brasil sem fome"<br><br>"por ĉiu popolo ties propran lingvon, por ĉiuj popoloj la esperantan"<br>(para cada povo sua própria língua, para todos os povos o Esperanto)<o:p></o:p></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'> <o:p></o:p></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>E nunca votarei no PSDB/DEM!<o:p></o:p></p></div></div></div></div></div></div></div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'> <o:p></o:p></p></div></div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>_______________________________________________<br>R-br mailing list<br><a href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">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.<o:p></o:p></p></blockquote></div></div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'><br>_______________________________________________<br>R-br mailing list<br><a href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">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.<o:p></o:p></p></blockquote></div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'><br><br clear=all><o:p></o:p></p><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'> <o:p></o:p></p></div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>-- <o:p></o:p></p><div><div><div><div><p class=MsoNormal style='mso-margin-top-alt:auto;margin-bottom:12.0pt'>MARCOS BISSOLI<o:p></o:p></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>Faculdade de Nutrição<br>Universidade Federal de Alfenas<br><br>Blog: <span style='font-size:9.5pt'><a href="http://bocademiamaldita.blogspot.com/" target="_blank">bocademiamaldita.blogspot.com/</a></span><o:p></o:p></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>E-mail: <a href="mailto:mbissoli@gmail.com" target="_blank">mbissoli@gmail.com</a><br>Twitter: #mbissoli<o:p></o:p></p><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'><br>Alfenas, Minas Gerais, Brasil<br><br><br>*****Pense na Natureza antes de Imprimir*****<br>Divulgue ON-LINE<br><br>Eu apoio a ENEN "na luta por um Brasil sem fome"<br><br>"por ĉiu popolo ties propran lingvon, por ĉiuj popoloj la esperantan"<br>(para cada povo sua própria língua, para todos os povos o Esperanto)<o:p></o:p></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'> <o:p></o:p></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>E nunca votarei no PSDB/DEM!<o:p></o:p></p></div></div></div></div></div></div></div></div></div></div></blockquote></div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'><br><br clear=all><o:p></o:p></p><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'> <o:p></o:p></p></div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>-- <o:p></o:p></p><div><div><div><div><p class=MsoNormal style='mso-margin-top-alt:auto;margin-bottom:12.0pt'>MARCOS BISSOLI<o:p></o:p></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>Faculdade de Nutrição<br>Universidade Federal de Alfenas<br><br>Blog: <span style='font-size:9.5pt'><a href="http://bocademiamaldita.blogspot.com/" target="_blank">bocademiamaldita.blogspot.com/</a></span><o:p></o:p></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>E-mail: <a href="mailto:mbissoli@gmail.com" target="_blank">mbissoli@gmail.com</a><br>Twitter: #mbissoli<o:p></o:p></p><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'><br>Alfenas, Minas Gerais, Brasil<br><br><br>*****Pense na Natureza antes de Imprimir*****<br>Divulgue ON-LINE<br><br>Eu apoio a ENEN "na luta por um Brasil sem fome"<br><br>"por ĉiu popolo ties propran lingvon, por ĉiuj popoloj la esperantan"<br>(para cada povo sua própria língua, para todos os povos o Esperanto)<o:p></o:p></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'> <o:p></o:p></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>E nunca votarei no PSDB/DEM!<o:p></o:p></p></div></div></div></div></div></div></div></div></div></div></blockquote></div><p class=MsoNormal><br><br clear=all><o:p></o:p></p><div><p class=MsoNormal><o:p> </o:p></p></div><p class=MsoNormal>-- <o:p></o:p></p><div><div><div><div><p class=MsoNormal style='margin-bottom:12.0pt'>MARCOS BISSOLI<o:p></o:p></p></div><div><p class=MsoNormal>Faculdade de Nutrição<br>Universidade Federal de Alfenas<br><br>Blog: <span style='font-size:9.5pt'><a href="http://bocademiamaldita.blogspot.com/" target="_blank">bocademiamaldita.blogspot.com/</a></span><o:p></o:p></p></div><div><p class=MsoNormal>E-mail: <a href="mailto:mbissoli@gmail.com" target="_blank">mbissoli@gmail.com</a><br>Twitter: #mbissoli<o:p></o:p></p><div><p class=MsoNormal><br>Alfenas, Minas Gerais, Brasil<br><br><br>*****Pense na Natureza antes de Imprimir*****<br>Divulgue ON-LINE<br><br>Eu apoio a ENEN "na luta por um Brasil sem fome"<br><br>"por ĉiu popolo ties propran lingvon, por ĉiuj popoloj la esperantan"<br>(para cada povo sua própria língua, para todos os povos o Esperanto)<o:p></o:p></p></div><div><p class=MsoNormal><o:p> </o:p></p></div><div><p class=MsoNormal>E nunca votarei no PSDB/DEM!<o:p></o:p></p></div></div></div></div></div></div></div></div></div><p class=MsoNormal><br><br clear=all><o:p></o:p></p><div><p class=MsoNormal><o:p> </o:p></p></div><p class=MsoNormal>-- <o:p></o:p></p><div><div><div><div><p class=MsoNormal style='margin-bottom:12.0pt'>MARCOS BISSOLI<o:p></o:p></p></div><div><p class=MsoNormal>Faculdade de Nutrição<br>Universidade Federal de Alfenas<br><br>Blog: <span style='font-size:9.5pt'><a href="http://bocademiamaldita.blogspot.com/" target="_blank">bocademiamaldita.blogspot.com/</a></span><o:p></o:p></p></div><div><p class=MsoNormal>E-mail: <a href="mailto:mbissoli@gmail.com" target="_blank">mbissoli@gmail.com</a><br>Twitter: #mbissoli<o:p></o:p></p><div><p class=MsoNormal><br>Alfenas, Minas Gerais, Brasil<br><br><br>*****Pense na Natureza antes de Imprimir*****<br>Divulgue ON-LINE<br><br>Eu apoio a ENEN "na luta por um Brasil sem fome"<br><br>"por ĉiu popolo ties propran lingvon, por ĉiuj popoloj la esperantan"<br>(para cada povo sua própria língua, para todos os povos o Esperanto)<o:p></o:p></p></div><div><p class=MsoNormal><o:p> </o:p></p></div><div><p class=MsoNormal>E nunca votarei no PSDB/DEM!<o:p></o:p></p></div></div></div></div></div></div></div></body></html>