<html><body><div style="color:#000; background-color:#fff; font-family:garamond, new york, times, serif;font-size:14pt">Edson,<br>A quantidade de dados não é pequena????<br>Abs.<br><div><span><br></span></div><div><br></div> <div style="font-family: garamond, new york, times, serif; font-size: 14pt;"> <div style="font-family: times new roman, new york, times, serif; font-size: 12pt;"> <div dir="ltr"> <font face="Arial" size="2"> <hr size="1"> <b><span style="font-weight:bold;">De:</span></b> Edson Lira <edinhoestat@yahoo.com.br><br> <b><span style="font-weight: bold;">Para:</span></b> R-br Lista <r-br@listas.c3sl.ufpr.br> <br> <b><span style="font-weight: bold;">Enviadas:</span></b> Terça-feira, 23 de Outubro de 2012 12:54<br> <b><span style="font-weight: bold;">Assunto:</span></b> [R-br] Regressão Logística<br> </font> </div> <br><div id="yiv1967361227"><div><div style="color:#000;background-color:#fff;font-family:lucida console,
sans-serif;font-size:10pt;"><div><span>Estou usando as rotinas abaixo,</span></div><div style="color:rgb(0, 0, 0);font-size:13.3333px;font-family:lucida console, sans-serif;background-color:transparent;font-style:normal;"><br></div><div style="color:rgb(0, 0, 0);font-size:13.3333px;font-family:lucida console, sans-serif;background-color:transparent;font-style:normal;">PARA LER OS DADOS<br><span></span></div><div style="color:rgb(0, 0, 0);font-size:13.3333px;font-family:lucida console, sans-serif;background-color:transparent;
font-style:normal;">h<-read.csv2("C:\\Consulta\\CIESA\\bd.csv")<br>h[h==""]<-NA<br>names(h)<br>table(h$quem_elab)<br>str(h)<br>h$tempo<-2012-h$ano_ini<br>h$melh_orig<-h$niv_melhora[drop=T]<br>h$niv_melhora<-h$niv_melhora[drop=T]<br>h$bsc<-h$bsc[drop=T]<br><br><br>levels(h$niv_melhora)<br>levels(h$niv_melhora)[levels(h$niv_melhora)%in% c("elevado","muito elevado")] <- "Melhora"<br>levels(h$niv_melhora)[levels(h$niv_melhora)%in% c("nenhum","pouco",<br>"intermediário")] <- "Sem melhora"<br>#levels(h$niv_melhora)[levels(h$niv_melhora)%in% c("intermediário")] <-""<br><br>levels(h$quem_elab)<br>h$part<-h$quem_elab<br>levels(h$part)[levels(h$part)%in% c("alta direção")]<- "Direção"<br>levels(h$part)[levels(h$part)%in% c("alta direção\ngerentes",<br>"alta direção\ngerentes\nfuncionários","alta direção\ngerentes\noutros","gerentes",<br>"funcionários")] <-
"Gerentes"<br>table(h$part)<br><br><br><br><span></span></div><div style="color:rgb(0, 0, 0);font-size:13.3333px;font-family:lucida console, sans-serif;background-color:transparent;font-style:normal;"><span>#PARA SELECIONAR SOMENTE AS VARIÁVEIS QUE PENSO USAR NO MODELO<br></span></div><div style="color:rgb(0, 0, 0);font-size:13.3333px;font-family:lucida console, sans-serif;background-color:transparent;font-style:normal;"><span>nb<-data.frame(melhora=h$melh_orig,bsc=h$bsc,champion=h$champion,ap_dir=h$apoio_bsc,<br>niv_tre=h$niv_treina,cons_ext=h$cext_bsc,reu_f_est=h$reu_func_estr,<br>com_est=h$comu_estr,ajuste=h$ajustes,ind_apr=h$ind_aprop,ind_desemp=h$ind_desd_dep,<br>mapa_est=h$mapa_estr,tempo=h$tempo)<br></span></div><div style="color:rgb(0, 0, 0);font-size:13.3333px;font-family:lucida console, sans-serif;background-color:transparent;font-style:normal;"><span><br></span></div><div style="color:rgb(0, 0, 0);font-size:13.3333px;font-family:lucida
console, sans-serif;
background-color:transparent;font-style:normal;"><span><br></span></div><div style="color:rgb(0, 0, 0);font-size:13.3333px;font-family:lucida console, sans-serif;background-color:transparent;font-style:normal;"><span>levels(nb$melhora)<br>levels(nb$melhora)[levels(nb$melhora)%in% c("elevado","muito elevado")] <- 1<br>levels(nb$melhora)[levels(nb$melhora)%in% c("nenhum","pouco",<br>"intermediário")] <- 0<br></span></div><div style="color:rgb(0, 0, 0);font-size:13.3333px;font-family:lucida console, sans-serif;background-color:transparent;font-style:normal;"><br></div><div style="color:rgb(0, 0, 0);font-size:13.3333px;font-family:lucida console, sans-serif;background-color:transparent;font-style:normal;">nb$desemp<-as.numeric(as.character(nb$melhora))<br><br><span></span></div><div style="color:rgb(0, 0, 0);font-size:13.3333px;font-family:lucida console, sans-serif;background-color:transparent;
font-style:normal;"><span><br></span></div><div style="color:rgb(0, 0, 0);font-size:13.3333px;font-family:lucida console, sans-serif;background-color:transparent;font-style:normal;"><span><br></span></div><div style="color:rgb(0, 0, 0);font-size:13.3333px;font-family:lucida console, sans-serif;background-color:transparent;font-style:normal;"><span>#PARA ELIMINAR OS NA'S</span></div><div style="color:rgb(0, 0, 0);font-size:13.3333px;font-family:lucida console, sans-serif;background-color:transparent;font-style:normal;"><br><span></span></div><div style="color:rgb(0, 0, 0);font-size:13.3333px;font-family:lucida console, sans-serif;background-color:transparent;font-style:normal;"><span><br>cd <- subset(nb, subset=!is.na(desemp) <br>& !is.na(bsc)& !is.na(champion)& !is.na(ap_dir)& !is.na(niv_tre)& !is.na(cons_ext) <br>& !is.na(reu_f_est)& !is.na(com_est)& !is.na(ajuste)& !is.na(ind_apr)&
!is.na(ind_desemp)<br>& !is.na(mapa_est)& !is.na(tempo)& !is.na(desemp)) <br> dim(cd)<br></span></div><div style="color:rgb(0, 0, 0);font-size:13.3333px;font-family:lucida console, sans-serif;background-color:transparent;font-style:normal;"><br><span></span></div><div style="color:rgb(0, 0, 0);font-size:13.3333px;font-family:lucida console, sans-serif;background-color:transparent;font-style:normal;">#MODELO DE REGRESSÃO<br>glm3<-glm(desemp~bsc+champion+ap_dir+niv_tre+cons_ext+reu_f_est+com_est+<br>ajuste+ind_apr+ind_desemp+mapa_est+tempo+desemp,family=binomial,data=cd) <br><br>modelstep <- step(glm3, direction = "both")<br>logistic.display(modelstep) <br>summary(modelstep)<br><br><span></span></div><div style="color:rgb(0, 0, 0);font-size:13.3333px;font-family:lucida console, sans-serif;background-color:transparent;font-style:normal;"><span><br></span></div><div style="color:rgb(0, 0,
0);font-size:13.3333px;font-family:lucida console, sans-serif;background-color:transparent;font-style:normal;"><span>Estou tentando aplicar uma regressão logística, quero estimar um modelo aonde a variável resposta é o desempenho de determinada empresa, aonde tenho os níveis melhora=1, sem melhora=0 com outras variáveis. No link abaixo tem uma amostra do banco de dados.</span></div><div style="color:rgb(0, 0, 0);font-size:13.3333px;font-family:lucida console, sans-serif;background-color:transparent;font-style:normal;"><span><br></span></div><div><span>http://www.datafilehost.com/download-996309b6.html</span></div><div style="color:rgb(0, 0, 0);font-size:13.3333px;font-family:lucida console, sans-serif;background-color:transparent;font-style:normal;"><br></div><div style="color:rgb(0, 0, 0);font-size:13.3333px;font-family:lucida console, sans-serif;background-color:transparent;font-style:normal;">Alguém poderia
dar uma avaliada no que estou fazendo? </div><br>[ ]'s.<br><span></span><div> </div><div>Edson Lira<br>Estatístico<br>Manaus-Amazonas</div></div></div></div><br>_______________________________________________<br>R-br mailing list<br><a ymailto="mailto:R-br@listas.c3sl.ufpr.br" href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br><a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br><br> </div> </div> </div></body></html>