Version imprimable
bonjour,
je pars d'une matrice 4x2 composée de couples de variables issues d'une regression;
A cette matrice je voudrais ajouter une troisième colonne contenant le resultat d'un test d'égalité des coefficients genre linear.hypothesis ou wald.test pour chaque couple de variable (chaque ligne sur les 2 colonnes)
je pense utiliser un mapply avec comme syntaxe qqch du genre
mat_col3<-(mat,function(col1,col2) linear.hypothesis(model,"col1-col2=0"))
pensez vous que cela est concevable?
merci
Bonsoir,
Peux tu nous mettre un exemple reproductible pour y voir en plus clair
Merci
Manoir
une regression nommée "Reg" type lm par exemple peut sortir les résultats suivants;
Estimate Std. error
(Intercept) 3.95811722 0.000719656
exerciceavant 0.178572325 0.000285748
anc_souAsou0-2 0.131144251 0.000420326
anc_souBsou3-5 0.005896919 0.000388712
anc_souDsou11-20 -0.076310332 0.000395076
je voudrais mener 2 à 2 un test d'égalité des coefficients via la fonction wald.test par exemple du pack aod.
exerciceavant est le 2ème terme, (intercept) le premier, anc_souD... le 5ème.
Ainsi deux à deux nous avons les couples de termes suivant c(1,2),c(1,3),c(1,4),c(2,3),c(2,4),c(3,4),c(4,4) qu'on peut mettre sous forme de data.frame "couple": col1=1er élément du couple, col2=2ème élt du couple.
sur ce data.frame, je veux ajouter une troisième colonne qui ligne par ligne (couple par couple) contient le résultat (la pvalue) de la fonction suivante :
wald.test(Sigma=vcov(Reg),b=coef(Reg),Terms=col1,H0=coef(Reg)[col2]).
Je pensais donc appliquer à mon data.frame "couple" la fonction ci dessus via un apply ...