<div dir="ltr">Prezados colegas,<div><br></div><div>Possuo dados de um questionário e quero estimar o efeito individual de cada variável explicativa (todas são qualitativas) sobre a variável dependente (contínua). Contudo o padrão é considerar os coeficientes estimados como a diferença entre o o grupo de referência (valor estimado do intercepto) e o grupo do coeficiente. Cada variável qualitativa diz respeito a uma característica de uma determinada atividade e a variável contínua é o tempo tomado para executar tal atividade. Destaco que o nível de referência de cada fator diz respeito ao caso em que nenhuma das opções para aquele fator foi marcado. Os dados estão em "<span style="font-family:Calibri,sans-serif;font-size:11pt;line-height:115%"><a href="http://www.datafilehost.com/d/c7f0d342" target="_blank">http://www.datafilehost.com/d/c7f0d342</a>"</span>. Não os coloquei no CMR, pois ainda não consegui (peço minhas sinceras desculpas), mas espero que isto não seja empecilho.Também não coloquei apenas uma amostra dos dados, pois estava havendo problemas de matriz singular.</div><div><br></div><div><font color="#222222">1º (e principal) Problema - Para obter o efeito individual dos níveis de cada fator, considerei que o grupo de referência possui efeito zero e fiz os seguintes passos:</font></div><div><font color="#222222"><br></font></div><div><font color="#222222">#Como o arquivo não está sendo carregado no CMR, supõe-se aqui que ele foi baixado da internet e já se encontra carregado no R.</font></div><div><font color="#222222">#Utilizarei uma regressão quantílica, por isso segue o pacote abaixo. </font></div><div><font color="#222222">install.packages(quantreg)</font></div><div><font color="#222222"><font color="#222222">library(quantreg)</font></font></div><div><font color="#222222"><font color="#222222">#Transformando as variáveis categóricas em objetos únicos</font></font></div><div><div>p_1 = table(1:length(tabela1.1$p1),as.ordered(tabela1.1$p1))</div><div>p_21 = table(1:length(tabela1.1$p21),as.factor(tabela1.1$p21))</div><div>p_22 = table(1:length(tabela1.1$p22),as.factor(tabela1.1$p22))</div><div>p_23 = table(1:length(tabela1.1$p23),as.factor(tabela1.1$p23))</div><div>p_24 = table(1:length(tabela1.1$p24),as.factor(tabela1.1$p24))</div><div>p_25 = table(1:length(tabela1.1$p25),as.factor(tabela1.1$p25))</div><div>p_34 = table(1:length(tabela1.1$p34),as.ordered(tabela1.1$p34))</div><div>p_5 = table(1:length(tabela1.1$p5),as.ordered(tabela1.1$p5))</div><div>p_6 = table(1:length(tabela1.1$p6),as.ordered(tabela1.1$p6))</div><div>p_7 = table(1:length(tabela1.1$p7),as.ordered(tabela1.1$p7))</div><div>p_8 = table(1:length(tabela1.1$p8),as.ordered(tabela1.1$p8))</div><div>p_9 = table(1:length(tabela1.1$p9),as.ordered(tabela1.1$p9))</div></div><div><font color="#222222"><font color="#222222">#Regressão do modelo sem intercepto, mas também considerando que o grupo de referência tem valor 0, haja vista que o grupo de referência significa que nenhum dos fatores foi marcado (se nenhum item foi marcado, considero que o tempo tomado para executar a atividade é praticamente zero).</font></font></div><div><font color="#222222"><font color="#222222">qrModel=rq(data=tabela1.1, pontoefetivo ~ 0 + p_1[,-1] + p_21[,-1] + p_22[,-1] + p_23[,-1] + p_24[,-1] + p_25[,-1] + p_34[,-1] + p_5[,-1] + p_6[,-1] + p_7[,-1] + p_8[,-1] + p_9[,-1], tau=0.5)<br></font></font></div><div><font color="#222222"><font color="#222222">summary(qrModel)<br></font></font></div><div><font color="#222222"><font color="#222222"><br></font></font></div><div><font color="#222222"><font color="#222222">A minha ideia foi que, dado que o grupo de referência é zero, o coeficiente estimado é justamente o efeito individual do nível escolhido da variável. A minha ideia está certa? Se não, o que faço para obter esses efeitos individuais?</font></font></div><div><font color="#222222"><font color="#222222"><br></font></font></div><div><font color="#222222"><font color="#222222">2º Problema: Supondo que fiz tudo certinho acima, os fatores ordenados não possuem efeitos crescentes [Ver summary(qrModel)]. Mas não deveriam possuir? Se deveriam, o que faço para garantir tal efeito?</font></font></div><div><font color="#222222"><font color="#222222"><br></font></font></div><div><font color="#222222"><font color="#222222">3º Problema: Mas uma vez supondo que tudo está correto, esperava que todos os efeitos sobre o tempo de execução das atividades fossem não negativos. Ou estou errado a respeito disso? Se estou certo, o que faço para garantir que todos os valores sejam não negativos.</font></font></div><div><font color="#222222"><font color="#222222"><br></font></font></div><div><font color="#222222"><font color="#222222"><font color="#222222">Desde já agradeço a todos e aguardo qualquer ajuda.<a href="http://www.datafilehost.com/get.php?file=c7f0d342" target="_blank"><br></a></font></font></font>-- <br><div><b>Rafael Carneiro da Costa<br>Doutorando em Economia<br>Centro de Pós-Graduação em Economia - CAEN<br>Universidade Federal do Ceará - UFC</b><br></div>
</div></div>