calculer si le coefficient de régression est statistiquement significatif dans R

Question

J'ai des résultats d'une analyse de régression effectuée avec un autre programme et je voudrais tester avec r si elles sont importantes.Je sais que ls.diag () calcule les erreurs standard et les tests t pour les résultats de la régression, mais cela nécessite un format d'entrée très spécifique (c'est-à-dire le résultat de LSFIT ()), donc je ne pense pas pouvoir l'utiliser.Existe-t-il une fonction dans laquelle calculer les erreurs standard et les tests t pour l'analyse de régression qui me permet de donner simplement aux coefficients pertinents à la fonction à la main?

Answer 1

Je ne suis pas si sûr que c'est ce que vous cherchez, mais voici une ligne directrice

# this is a model obtained from ?lm 
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) this is our target

Supposons que nous n'avions que les coefficients de régression, ses erreurs standard et la taille de l'échantillon

beta <- coef(lm.D9)
errorBeta <- summary(lm.D9)$coefficients[,2]
n <- length(weight) # the sample size
k <- length(beta) # number of regression parameters

Je pense que c'est votre cas, si vous n'avez pas les erreurs standard de coefficient, vous devez les estimer, c'est assez facile.

Une fois que vous avez les coefficients de régression et ses erreurs standard, on peut estimer le T-Stat:

t_stats <- beta/errorBeta

La règle echolique stipule que si | t_stats |>= 2 Ensuite, le coefficient est statistiquement significatif à 5% de niveau.Mais si vous voulez connaître la valeur p, utilisez:

pt(abs(t_stats), n-k, lower.tail=FALSE)*2

Si les valeurs P> 0.05, les coefficients associés ne sont pas statistiques significatifs à ce niveau.

Tout ce dont vous avez besoin est de connaître les coefficients, ses erreurs standard et la taille de l'échantillon.Sinon, vous ne le ferez pas.