<html><body><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></body></html>