Eu não conheço nenhuma função capaz de fazer isto o R. Em geral, o que se sugere é fazer testes de razão de verossimilhança, que são muito mais gerais e resolvem problemas mais sofisticados, até mesmo em modelos lineares generalizados.
Entretanto, nada nos impede de escrevermos nossa própria função. Afinal, temos que testar as hipóteses
H_0: \beta = \beta_0
H_1: \beta != \beta_0
em que \beta_0 é o valor de referência sob o qual desejamos realizar o teste. Por padrão, o valor de \beta_0 é igual a zero. A função TesteCoeficiente
, definida abaixo, implementa este teste de hipóteses para qualquer valor de \beta_0. Basta informar os seguintes valores para a função:
reg
: o modelo de regressão obtido com a função lm
coeficiente
: o número do coeficiente no modelo de regressão. Para y = \beta_0 + \beta_1*x
, 1 significa testar \beta_0
e 2 significa testar \beta_1
h0
: valor da hipótese nula a ser testada. Em geral, este valor é zero
Com isto definido, basta rodar o exemplo abaixo para ver como a função é executada:
x <- c(1, 2, 3, 4, 5)
y <- c(1.2, 2.4, 3.3, 4.2, 5.1)
reg <- lm(y~x)
TesteCoeficiente <- function(reg, coeficiente, h0){
estimativas <- coef(summary(reg))
estatistica <- (estimativas[coeficiente, 1]-h0)/estimativas[coeficiente, 2]
2 * pt(abs(estatistica), reg$df.residual, lower.tail = FALSE)
}
TesteCoeficiente(reg, coeficiente=2, h0=0)
[1] 2.771281e+01 1.031328e-04
summary(reg)
Call:
lm(formula = y ~ x)
Residuals:
1 2 3 4 5
-1.200e-01 1.200e-01 6.000e-02 1.457e-16 -6.000e-02
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 0.36000 0.11489 3.133 0.051929 .
x 0.96000 0.03464 27.713 0.000103 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Residual standard error: 0.1095 on 3 degrees of freedom
Multiple R-squared: 0.9961, Adjusted R-squared: 0.9948
F-statistic: 768 on 1 and 3 DF, p-value: 0.0001031
Note que o output da função TesteCoeficiente
é idêntico ao da função summary
testando h0 = 0
. Portanto, a função funciona para este caso particular. É bastante razoável supor que ela funciona para qualquer valor de \beta_0. Agora basta trocar o valor do parâmetro h0
da função para realizar o teste de hipóteses desejado.