<div dir="ltr">Obrigado pelas respostas Leonardo e Daniel. <div><br></div><div>Concordo que a supervisão de campo tem que garantir os protocolos de aleatoriedade. No entanto, no exercício estou mostrando que de fato o tratamento e controle eram iguais no momento inicial em termos de renda e que após o treinamento eles se diferenciam. Mas quero mostrar também que o aumento da renda não está relacionada a nenhuma das características observáveis do beneficiários, apenas ao fato de terem participado do treinamento.</div><div><br></div><div>Daniel, você entendeu o que eu quero fazer. Queria apenas mostrar que a inclusão de variáveis no segundo modelo não altera significativamente o efeito estimado do programa, justamente por estarmos tratando de um experimento aleatório simples.</div><div><br></div><div>Você sabe se posso fazer um teste F para verificar se os dois parâmetros estimados são estatisticamente diferentes usando a função <span style="font-family:arial,sans-serif;font-size:13px">linearHypothesis?</span></div><div><br></div><div>Obrigado</div><div><br></div><div>Adriano</div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">Em 22 de outubro de 2014 00:14, Leonardo Ferreira Fontenelle <span dir="ltr"><<a href="mailto:leonardof@leonardof.med.br" target="_blank">leonardof@leonardof.med.br</a>></span> escreveu:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><u></u>




<div><div>Veja se isso faz sentido para o seu problema: talvez você possa testar a associação entre o tratamento e os demais preditores, em vez de comparar os dois modelos de regressão já apresentados.</div>
<div> </div>
<div>Mas isso se você realmente achar importante testar a associação. Por exemplo, num ensaio clínico randomizado, não faz sentido testar a hipótese de alocação aleatória. Em princípio, qualquer associação entre o tratamento alocado e as variáveis de exposição <b>é</b> aleatório, por menor que seja o valor de p. A supervisão de campo é quem vai dizer se houve violação do protocolo. E diferenças grandes entre o grupo de intervenção e o de controle podem precisar ser levadas em consideração nas análises, independentemente delas terem aparecido ao acaso ou não.</div>
<div> </div>
<div><div><a href="http://lattes.cnpq.br/9234772336296638" title="Currículo Lattes" target="_blank">Leonardo Ferreira Fontenelle</a><br></div>
</div>
<div> </div>
<div> </div>
<div>Em Ter 21 out. 2014, às 23:52, Adriano Borges Costa escreveu:<br></div>
<blockquote type="cite"><div><div class="h5"><div dir="ltr"><div style="font-family:arial,sans-serif;font-size:13.3333339691162px">Caros,<br></div>
<div style="font-family:arial,sans-serif;font-size:13.3333339691162px"> </div>
<div style="font-family:arial,sans-serif;font-size:13.3333339691162px">Estou fazendo um exercício teórico de regressão para estimar o efeito de um programa a partir da base de dados disponível no código baixo.<br></div>
<div style="font-family:arial,sans-serif;font-size:13.3333339691162px"> </div>
<div style="font-family:arial,sans-serif;font-size:13.3333339691162px">library(foreign)<br></div>
<div style="font-family:arial,sans-serif;font-size:13.3333339691162px">lalonde <- read.dta("<a href="http://dl.dropboxusercontent.com/u/7453595/BD_LAB/lalonde.dta" target="_blank">http://dl.dropboxusercontent.com/u/7453595/BD_LAB/lalonde.dta</a>") #carrega a base de dados<br></div>
<div style="font-family:arial,sans-serif;font-size:13.3333339691162px"> </div>
<div style="font-family:arial,sans-serif;font-size:13.3333339691162px">str(lalonde) # variáveis<br></div>
<div style="font-family:arial,sans-serif;font-size:13.3333339691162px"> </div>
<div style="font-family:arial,sans-serif;font-size:13.3333339691162px"># As variáveis re74, re75 e re78 contêm informações sobre o rendimento para os anos de 1974, 1975 e 1978, respectivamente. <br></div>
<div style="font-family:arial,sans-serif;font-size:13.3333339691162px"># A variável de tratamento t indica se o indivíduo participou do programa (t=1) ou se não participou do programa (t=0).<br></div>
<div style="font-family:arial,sans-serif;font-size:13.3333339691162px"># As demais variáveis contém características dos trabalhadores<br></div>
<div style="font-family:arial,sans-serif;font-size:13.3333339691162px"> </div>
<div style="font-family:arial,sans-serif;font-size:13.3333339691162px">Trata-se de um experimento aleatório e quero mostrar que há aleatoriedade no experimento. Assim, estou rodando a seguinte regressão:<br></div>
<div style="font-family:arial,sans-serif;font-size:13.3333339691162px"> </div>
<div style="font-family:arial,sans-serif;font-size:13.3333339691162px">regressao1 <- lm(re78 ~ t, data=lalonde)<br></div>
<div style="font-family:arial,sans-serif;font-size:13.3333339691162px">summary(regressao1)<br></div>
<div style="font-family:arial,sans-serif;font-size:13.3333339691162px"> </div>
<div style="font-family:arial,sans-serif;font-size:13.3333339691162px">E em seguida incluindo os controles em uma outra regressão e gostaria de mostrar que o estimador de t (tratamento) varia pouco (cerca de <span style="font-family:Calibri,sans-serif"><span style="font-size:11pt">U$ 100,00)</span></span>.<br></div>
<div style="font-family:arial,sans-serif;font-size:13.3333339691162px"> </div>
<div style="font-family:arial,sans-serif;font-size:13.3333339691162px">regressao2 <- lm(re78 ~ t+age+educ+nodegree+black+hisp+married, data=lalonde)<br></div>
<div style="font-family:arial,sans-serif;font-size:13.3333339691162px">summary(regressao2)<br></div>
<div style="font-family:arial,sans-serif;font-size:13.3333339691162px"> </div>
<div style="font-family:arial,sans-serif;font-size:13.3333339691162px">Assim, gostaria de fazer um teste para verificar se o estimador de t na regressão 1 é estatisticamente diferente da mesmo estimador na regressão 2. <br></div>
<div style="font-family:arial,sans-serif;font-size:13.3333339691162px"> </div>
<div style="font-family:arial,sans-serif;font-size:13.3333339691162px">Está certo eu fazer da seguinte forma?<br></div>
<div style="font-family:arial,sans-serif;font-size:13.3333339691162px"> </div>
<div style="font-family:arial,sans-serif;font-size:13.3333339691162px">install.packages("car")<br></div>
<div style="font-family:arial,sans-serif;font-size:13.3333339691162px">library(car)<br></div>
<div style="font-family:arial,sans-serif;font-size:13.3333339691162px">linearHypothesis(regressao2,"t=1794.3", test = "F")<br></div>
<div style="font-family:arial,sans-serif;font-size:13.3333339691162px"> </div>
<div style="font-family:arial,sans-serif;font-size:13.3333339691162px">Desde já obrigado pelas ajudas.<br></div>
<div style="font-family:arial,sans-serif;font-size:13.3333339691162px"> </div>
<div style="font-family:arial,sans-serif;font-size:13.3333339691162px">At.<br></div>
<div> </div>
<div>-- <br></div>
<div>Adriano Borges Costa<br></div>
</div>
</div></div><div><u>_______________________________________________</u><br></div>
<div>R-br mailing list<br></div>
<div><a href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">R-br@listas.c3sl.ufpr.br</a><br></div>
<div><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></div>
<div>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></div>
</blockquote><div> </div>
</div>

<br>_______________________________________________<br>
R-br mailing list<br>
<a 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></blockquote></div><br><br clear="all"><div><br></div>-- <br>Adriano Borges Costa
</div>