[R-br] Teste de Hipótese de Coeficientes de uma regressão

Quarta Junho 11 11:05:17 BRT 2014

Não é o teste apresentado no final do summary?

> ctl <- c(4.17,5.58,5.18,6.11,4.50,4.61,5.17,4.53,5.33,5.14)
> trt <- c(4.81,4.17,4.41,3.59,5.87,3.83,6.03,4.89,4.32,4.69)
> group <- gl(2, 10, 20, labels = c("Ctl","Trt"))
> weight <- c(ctl, trt)
> lm.D9 <- lm(weight ~ group)
> summary(lm.D9)

Call:
lm(formula = weight ~ group)

Residuals:
    Min      1Q  Median      3Q     Max
-1.0710 -0.4938  0.0685  0.2462  1.3690

Coefficients:
            Estimate Std. Error t value Pr(>|t|)
(Intercept)   5.0320     0.2202  22.850 9.55e-15 ***
groupTrt     -0.3710     0.3114  -1.191    0.249
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Residual standard error: 0.6964 on 18 degrees of freedom
Multiple R-squared:  0.07308,   Adjusted R-squared:  0.02158
*F-statistic: 1.419 on 1 and 18 DF,  p-value: 0.249*

2014-06-11 11:01 GMT-03:00 Adriano Borges Costa <adrianobfc em gmail.com>:

> Benilton, obrigado pela resposta, mas o teste para cada parâmetro isolado
> está ok. O que preciso é testar os parâmetros juntos. Vamos ver se consigo
> me explicar melhor.
>
> Em uma regressão linear como esta,
>
> modelo<-lm(y ~ α + β1X1 + β2X2 + β3X3 + ui)
>
> digamos que por meio do summary(modelo) vejo que os três parâmetros são
> estatisticamente significantes, mas com valores próximos de zero.
>
> Então gostaria de fazer um teste F ou uma anova para verificar, por
> exemplo, se β1=β2. Ou se β1+β2=0. Ou se β1=β2=β3=0. Ou se 2*β1*β2=β3. Ou
> testar qualquer outra relação entre os parâmetros. Eu encontrei como se faz
> este teste no Stata, conforme este tutorial
> http://www.ats.ucla.edu/stat/stata/library/anova_comp.htm, mas não estou
> achando no R. Este teste está disponível no summary(modelo)?
>
> Ficou mais claro agora? Me desculpe se estou sendo meio confuso.
>
> Obrigado
>
> Adriano
>
>
>
> Em 10 de junho de 2014 08:31, Benilton Carvalho <
> beniltoncarvalho em gmail.com> escreveu:
>
> vc aparentemente deixou escapar a secao "Coefficients" do resultado de
>> summary(Modelo)... esta contem sim os testes de hipoteses para cada
>> coeficiente de modo isolado.
>>
>> Uma vez dado este passo inicial, vc podera' procurar como ajustar os
>> contrastes de interesse... (para  a comparacao entre os coeficientes -
>> assumindo que sejam padronizados, claro).
>>
>> Mas, ao inves de ficarmos 'chutando' sugestoes, um exemplo reprodutivel
>> e' boa pratica....
>>
>>
>> Em 9 de junho de 2014 11:39, Adriano Borges Costa <adrianobfc em gmail.com>
>> escreveu:
>>
>> Obrigado pela resposta Paulo, mas no summary não me apresenta este teste.
>>>
>>> Eu preciso fazer algo como este teste que pode ser feito no stata:
>>> http://www.ats.ucla.edu/stat/stata/library/anova_comp.htm
>>>
>>> Você pode me ajudar? Desde já obrigado
>>>
>>> Adriano
>>>
>>>
>>> Em 6 de junho de 2014 14:12, Paulo Justiniano <paulojus em leg.ufpr.br>
>>> escreveu:
>>>
>>>> para comecar veja o resultado de
>>>>
>>>> summary(Modelo)
>>>>
>>>>
>>>> On Fri, 6 Jun 2014, Adriano Borges Costa wrote:
>>>>
>>>>  Boa Tarde,
>>>>> Sou iniciante no uso do R e estou com uma dívida bastante simples.
>>>>> Estou fazendo a análise de uma política pública e quero medir a relação
>>>>> entre o
>>>>> grau de execução da política e características das gestões municipais.
>>>>> Para tanto, estou usando como variáveis de controle informações
>>>>> socioeconômicas das cidades (população, renda, IDH, etc).
>>>>>
>>>>> Para tanto, estou rodando uma regressão linerar simples mais ou menos
>>>>> desta forma:
>>>>>
>>>>> Modelo<-lm(ExecuçãodoPrograma~VariáveisdeGestãoPública+Controles)
>>>>>
>>>>> O Modelo deu certo e estou conseguindo fazer a análise dos parâmetros.
>>>>> No entanto, preciso fazer um teste de hipótese para verificar se alguns
>>>>> parâmetros são iguais a zero e se uns são iguais aos outros.
>>>>>
>>>>> Como faço este teste no R? Eu devo usar as funções test.t() e
>>>>> var.test()?
>>>>>
>>>>> Obrigado
>>>>>
>>>>> --
>>>>> Adriano Borges Costa
>>>>>
>>>>>
>>>> _______________________________________________
>>>> R-br mailing list
>>>> R-br em listas.c3sl.ufpr.br
>>>> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
>>>> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
>>>> código mínimo reproduzível.
>>>>
>>>
>>>
>>>
>>> --
>>> Adriano Borges Costa
>>>
>>> _______________________________________________
>>> R-br mailing list
>>> R-br em listas.c3sl.ufpr.br
>>> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
>>> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
>>> código mínimo reproduzível.
>>>
>>
>>
>>
>> --
>>
>> [image: Benilton Carvalho on about.me]
>>
>> Benilton Carvalho
>> about.me/benilton
>>   <http://about.me/benilton>
>>
>> _______________________________________________
>> R-br mailing list
>> R-br em listas.c3sl.ufpr.br
>> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
>> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
>> código mínimo reproduzível.
>>
>
>
>
> --
> Adriano Borges Costa
>
> _______________________________________________
> R-br mailing list
> R-br em listas.c3sl.ufpr.br
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
> código mínimo reproduzível.
>
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20140611/8607f013/attachment.html>