<div dir="ltr"><div class="gmail_default" style="font-family:monospace,monospace;font-size:large;color:#073763">Saudações amigos de R, </div><div class="gmail_default" style="font-family:monospace,monospace;font-size:large;color:#073763"><br></div><div class="gmail_default" style="font-family:monospace,monospace;font-size:large;color:#073763">Depois de muito tempo sem interagir na lista venho pedir uma luz. Estou tentando montar um modelo de previsão penalizado. Eu já tinha usado no passado remoto o pacote glmpath que regulariza/penaliza modelos glm. </div><div class="gmail_default" style="font-family:monospace,monospace;font-size:large;color:#073763">Recentemente me deparei com a função glmnet - glmnetUtils - que faz o mesmo trabalho de forma mais automatizada, inclusive com validação cruzada já implementada. Porém, ao olhar a saída do modelo fiquei com dúvidas. Fui olhar a bibliografia original do pacote e não adiantou muito. Acho que Ã© uma expectativa equivocada de minha parte que está me impedindo de entender corretamente. </div><div class="gmail_default" style="font-family:monospace,monospace;font-size:large;color:#073763"><br></div><div class="gmail_default" style="font-family:monospace,monospace;font-size:large;color:#073763">Então ajustei o seguinte modelo:</div><div class="gmail_default" style="font-family:monospace,monospace;font-size:large;color:#073763"><br></div><div class="gmail_default" style="font-family:monospace,monospace;font-size:large;color:#073763">cva.glmnet(hav ~ ano_nasc_trunc60 + idade + indicacao_outras + sexo, data = bi[setD,],alpha = seq(0, 1, len = 11), family="binomial", nfolds = 1000, grouped = FALSE, <a href="http://trace.it">trace.it</a> = TRUE, outerParallel = makeCluster(4))<br></div><div class="gmail_default" style="font-family:monospace,monospace;font-size:large;color:#073763"><br></div><div class="gmail_default" style="font-family:monospace,monospace;color:rgb(7,55,99)"><font size="1">(Inntercept) Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â 113.55782297<br>ano_nasc_trunc60 Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â -0.05781642<br>idade Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â 0.04359743<br>indicacao_outrasViajante Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â -0.85226169<br>indicacao_outrasOTHER Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  -0.53550446<br>indicacao_outrasHIV/AIDS Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  0.49319132<br>indicacao_outrasAsplenia anatômica ou funcional e doenças relacionadas (anemia falciforme, talassemia maior) Â  . Â  Â  Â  Â  <br>indicacao_outrasPneumopatias crônicas (DPOC, Bronquiectasia, Asma etc) Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  . Â  Â  Â  Â  <br>indicacao_outrasCandidatos a transplante de Ã³rgão sólido Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  0.13338936<br>indicacao_outrasCardiopatias crônicas Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  -0.05086300<br>indicacao_outrasHepatopatias crônicas, Cirrose e portadores de HCV Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  . Â  Â  Â  Â  <br>indicacao_outrasComunicantes domiciliares de imunodeprimidos Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â -0.24283444<br>indicacao_outrasImunodeficiência/Imunodepressão outras (adquirida ou congênita) Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â . Â  Â  Â  Â  <br>indicacao_outrasImunodepressão terapêutica ou devido ao câncer Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  0.15478292<br>indicacao_outrasNão informado Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  -0.18649749<br>indicacao_outrasNefropatias crônicas/hemodialisados/sindrome nefrótica Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  . Â  Â  Â  Â  <br>indicacao_outrasProfissional de saúde Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  -0.51745839<br>indicacao_outrasTransplante de Ã³rgão sólido Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â . Â  Â  Â  Â  <br>sexoFeminino Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  0.03089234<br>sexoMasculino Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  -0.01749610<br> <br></font></div><div class="gmail_default" style="font-family:monospace,monospace;color:rgb(7,55,99)"><font size="4">Reparem que há dois coeficientes para o sexo, um pra masculino e outro para Feminino. Eu esperaria que houvesse somente 1, sendo que a referencia (no caso o masculino) não apareceria ou deveria ser zero. A mesma dúvida no preditor indicação, onde a categoria Viajantes deveria ser a referência e eu não esperaria que esta aparecesse ou fosse zero. Ã‰ claro que as categorias são mutuamente exclusivas, então quando um está presente, ou outro necessariamente é zero. Mas no glm, até a minha memória vai, o efeito de uma categoria deveria ser o efeito adicional ao efeito da referência quando a categoria está presente, mas na saída como está não consigo fazer esta interpretação. Alguem pode me dar uma luz de como interpretar corretamente ester aspecto?</font></div><div class="gmail_default" style="font-family:monospace,monospace;color:rgb(7,55,99)"><font size="4"><br></font></div><div class="gmail_default" style="font-family:monospace,monospace;color:rgb(7,55,99)"><font size="4">Em análises alternativas eu explorei interações e faz sentido que a idade tivesse efeitos diferentes para diferentes anos de nascimento, sendo que efeitos não lineares não foram detectados. Quando faço um ajuste com interação a saída fica ainda mais esquisita</font></div><div class="gmail_default" style="font-family:monospace,monospace;color:rgb(7,55,99)"><font size="4"><br></font></div><div class="gmail_default" style="font-family:monospace,monospace;color:rgb(7,55,99)"><font size="4">cva.glmnet(hav ~ ano_nasc_trunc60+idade+indicacao_outras+sexo+ano_nasc_trunc60*idade, data = bi[setD,],alpha = seq(0, 1, len = 11), family="binomial", nfolds = 1000, grouped = FALSE, <a href="http://trace.it">trace.it</a> = TRUE, outerParallel = makeCluster(4))<br></font></div><div class="gmail_default" style="font-family:monospace,monospace;color:rgb(7,55,99)"><font size="4"><br></font></div><div class="gmail_default" style="font-family:monospace,monospace;color:rgb(7,55,99)">(Intercept) Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  1.056153e+02<br>ano_nasc_trunc60 Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  -2.736619e-02<br>idade Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  1.656760e-02<br>indicacao_outrasViajante Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  -8.345313e-01<br>indicacao_outrasOTHER Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â -4.632827e-01<br>indicacao_outrasHIV/AIDS Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â 4.927235e-01<br>indicacao_outrasAsplenia anatômica ou funcional e doenças relacionadas (anemia falciforme, talassemia maior) Â . Â  Â  Â  Â  Â  <br>indicacao_outrasPneumopatias crônicas (DPOC, Bronquiectasia, Asma etc) Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â . Â  Â  Â  Â  Â  <br>indicacao_outrasCandidatos a transplante de Ã³rgão sólido Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â 1.358595e-01<br>indicacao_outrasCardiopatias crônicas Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  . Â  Â  Â  Â  Â  <br>indicacao_outrasHepatopatias crônicas, Cirrose e portadores de HCV Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â . Â  Â  Â  Â  Â  <br>indicacao_outrasComunicantes domiciliares de imunodeprimidos Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  -2.076846e-01<br>indicacao_outrasImunodeficiência/Imunodepressão outras (adquirida ou congênita) Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  . Â  Â  Â  Â  Â  <br>indicacao_outrasImunodepressão terapêutica ou devido ao câncer Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â 1.432945e-01<br>indicacao_outrasNão informado Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â -1.616079e-01<br>indicacao_outrasNefropatias crônicas/hemodialisados/sindrome nefrótica Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â . Â  Â  Â  Â  Â  <br>indicacao_outrasProfissional de saúde Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â -4.942054e-01<br>indicacao_outrasTransplante de Ã³rgão sólido Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  . Â  Â  Â  Â  Â  <br>sexoFeminino Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â 3.024890e-02<br>sexoMasculino Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â -1.611223e-02<br>ano_nasc_trunc60 Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  -2.651105e-02<br>idade Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  1.603003e-02<br>ano_nasc_trunc60:idade Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â 7.512202e-06<br><br></div><div class="gmail_default" style="font-family:monospace,monospace;color:rgb(7,55,99)"><font size="4">Reparem que agora, além do sexo, o ano de nascimento e a idade têm cada uma, além do coeficiente da interação, dois coeficientes diferentes de zero, quando eu esperava que houvesse somente o coeficiente da interação. Talvez entendendo os porques da primeira questão eu consiga achar o porque da segunda. Se puderem indicar alguma leitura, eu agradeço. </font></div><div class="gmail_default" style="font-family:monospace,monospace;color:rgb(7,55,99)"><font size="4"><br></font></div><div class="gmail_default" style="font-family:monospace,monospace;color:rgb(7,55,99)"><font size="4">Olhando o desempenho do ajuste com e sem a interação, a discriminação Ã© bem pareceida, mas a calibração Ã© um tanto melhor no ajuste com a interação.</font></div><div class="gmail_default" style="font-family:monospace,monospace;color:rgb(7,55,99)"><font size="4"><br></font></div><div class="gmail_default" style="font-family:monospace,monospace;color:rgb(7,55,99)"><font size="4">Abraços fraternos a todos. </font><span style="font-size:large"> </span><br></div><div class="gmail_default" style="font-family:monospace,monospace;color:rgb(7,55,99)"><font size="4"><br></font></div><div class="gmail_default" style="font-family:monospace,monospace;color:rgb(7,55,99)"><span style="color:rgb(0,0,102);font-family:"courier new",monospace"><font size="4">Pedro Brasil</font></span><br></div></div>