Teste de Hipótese de Coeficientes de uma regressão
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
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
Ola pessoal tudo bem. Existe algum pacote no Rque faça a árvore de regressão? Abraços Enviado do Yahoo Mail no Android
rpart e tree Leonardo Ferreira Fontenelle [1]http://lattes.cnpq.br/9234772336296638 Em Seg 9 jun. 2014, às 10:41, [2]geovanecb@yahoo.com.br escreveu: Ola pessoal tudo bem. Existe algum pacote no R que faça a árvore de regressão? Abraços [3]Enviado do Yahoo Mail no Android __________________________________________________________________ From: Paulo Justiniano <paulojus@leg.ufpr.br>; To: <r-br@listas.c3sl.ufpr.br>; Subject: Re: [R-br] Teste de Hipótese de Coeficientes de uma regressão Sent: Fri, Jun 6, 2014 5:12:41 PM 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 [4]R-br@listas.c3sl.ufpr.br [5]https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem ([6]http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível. References 1. http://lattes.cnpq.br/9234772336296638 2. mailto:geovanecb@yahoo.com.br 3. https://br.overview.mail.yahoo.com/mobile/?.src=Android 4. mailto:R-br@listas.c3sl.ufpr.br 5. https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br 6. http://www.leg.ufpr.br/r-br-guia
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@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@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
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@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@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@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@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>
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@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@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@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@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@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@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
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@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@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@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@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@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@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@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@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.
Para testar hipóteses sobre os parâmetros tem-se várias situações. Seja o seu modelo um linear do tipo b_0+b_1 x_1+...+b_k x_k, portanto com um vetor de parâmetros com k+1 elementos. Temos os seguintes grupos de hipótese: 1) Se você quer testar a hipótese nula **conjunta** H0: [b_1, b_2, ..., b_k] = [0,0,...,0], sou seja, todos os b (exceto o b_0) serem **simultaneamente** 0, observe o valor do teste F que saí como última linha do output do summary() do modelo. Uma forma equivalente é fazer a razão de verossimilhanças entre o modelo com as cováriaveis e outro com apenas o intercepto. 2) Se você quer testar a hipótese nula **individual e marginal** H0: b_i = 0, onde i=0,1,...,k, esse teste saí na tabela de estimativas do summary() do modelo. Esse teste é para um b_i marginal aos demais, ou seja, não é feita restrição algumas sobre os b restantes. Sendo assim, se você aceitar que b_1 e b_2 são, cada um, 0 pelo teste t não significa que eles são simultaneamente 0. Reforçando, aceitar duas hipóteses marginais não implica em aceitar uma hipótese conjunta. 3) Se você quer testar uma hipótese que é uma **função linear** dos b, H0: Lb = 0, como um contraste envolvendo dois deles (comum em experimentos com fatores categóricos) então basta você especificar o vetor que representa a hipótese e passar, por exemplo, para a multcomp::glht(). Instale peça a documentação para ver exemplos de como usar. 4) Se você quer testar uma hipótese que é uma **função não linear** entre os b, por exemplo H0: b1/b2 = 10, então você tem que usar o método delta, obter o erro padrão, em seguida aplicar o teste. Encontre funções para método delta buscando com o comando RSiteSearch("delta method") em uma sessão do R. À disposição. Walmes.
Walmes, Era exatamente isso que eu precisava. No caso, meu interesse especificamente era saber como fazer o teste do ítem 3. Muito obrigado Adriano Em 11 de junho de 2014 13:31, walmes . <walmeszeviani@gmail.com> escreveu:
Para testar hipóteses sobre os parâmetros tem-se várias situações. Seja o seu modelo um linear do tipo b_0+b_1 x_1+...+b_k x_k, portanto com um vetor de parâmetros com k+1 elementos. Temos os seguintes grupos de hipótese: 1) Se você quer testar a hipótese nula **conjunta** H0: [b_1, b_2, ..., b_k] = [0,0,...,0], sou seja, todos os b (exceto o b_0) serem **simultaneamente** 0, observe o valor do teste F que saí como última linha do output do summary() do modelo. Uma forma equivalente é fazer a razão de verossimilhanças entre o modelo com as cováriaveis e outro com apenas o intercepto. 2) Se você quer testar a hipótese nula **individual e marginal** H0: b_i = 0, onde i=0,1,...,k, esse teste saí na tabela de estimativas do summary() do modelo. Esse teste é para um b_i marginal aos demais, ou seja, não é feita restrição algumas sobre os b restantes. Sendo assim, se você aceitar que b_1 e b_2 são, cada um, 0 pelo teste t não significa que eles são simultaneamente 0. Reforçando, aceitar duas hipóteses marginais não implica em aceitar uma hipótese conjunta. 3) Se você quer testar uma hipótese que é uma **função linear** dos b, H0: Lb = 0, como um contraste envolvendo dois deles (comum em experimentos com fatores categóricos) então basta você especificar o vetor que representa a hipótese e passar, por exemplo, para a multcomp::glht(). Instale peça a documentação para ver exemplos de como usar. 4) Se você quer testar uma hipótese que é uma **função não linear** entre os b, por exemplo H0: b1/b2 = 10, então você tem que usar o método delta, obter o erro padrão, em seguida aplicar o teste. Encontre funções para método delta buscando com o comando RSiteSearch("delta method") em uma sessão do R.
À disposição. Walmes.
_______________________________________________ R-br mailing list R-br@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
Deixa eu ver se entendi sua necessidade: Para uma variável nominal, portanto com vários parâmetros: modelo1 <- lm(y ~ a + b + c) modelo2 <- update(modelo1, . ~ . - c) anova(modelo1, modelo2) Para várias variáveis ao mesmo tempo: modelo1 <- lm(y ~ a + b + c) modelo2 <- update(modelo1, . ~ . - b - c) anova(modelo1, modelo2) É isso? Leonardo Ferreira Fontenelle [1]http://lattes.cnpq.br/9234772336296638 Em Sex 6 jun. 2014, às 12:41, Adriano Borges Costa escreveu: 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 [2]R-br@listas.c3sl.ufpr.br [3]https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem ([4]http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível. References 1. http://lattes.cnpq.br/9234772336296638 2. mailto:R-br@listas.c3sl.ufpr.br 3. https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br 4. http://www.leg.ufpr.br/r-br-guia
participantes (7)
-
Adriano Borges Costa -
Benilton Carvalho -
geovanecb＠yahoo.com.br -
Leonardo Ferreira Fontenelle -
Paulo Justiniano -
Rodrigo Coster -
walmes .