異なる回帰の係数を比較する

Question

私は同じ変数で2つの線形回帰の係数を比較しようとしていますが、異なるサブグループに対して実行しています。モデル1の係数がモデル2の係数と等しいかどうかを確認したいと思います。それぞれの係数について知る必要があります。

マイ再現性のあるデータ：

Data <- data.frame(
gender = sample (c("men", "women"), 2000, replace = TRUE), 
var1 = sample (c("value1", "value2"), 2000, replace = TRUE), 
var2 = sample (c("valueA", "valueB"), 2000, replace = TRUE), 
y = sample(0:10, 2000, replace = TRUE) 
) 

私は2つの回帰を実行します。基本的には

men <- subset(Data, gender =="men") 
women <- subset(Data, gender =="women") 

lm.men <- lm(y~var1+var2, data = men) 
summary(lm.men) 
lm.women <- lm(y~var1+var2, data = women) 
summary(lm.women) 

を、私がテストしたい場合：

• 係数var1lm.men =係数var1中で
• 係数lm.menでvar2 =係数、私の2つの標本が異なっているので、私は、anova()機能を使用することはできませんlm.women

でvar2。私はFテストを適用すべきだと思うが、私はこのテストのための機能を見つけることができない。

誰も私の問題を解決する方法を知っていますか？

Answer 1

@Barkerがコメントで指摘しているように、この質問の統計的部分はalready answered on CrossValidatedです。いくつかのRコーディングの詳細をここに追加します。変数ごとの男女相互作用モデルに適合し、（var1とvar2の効果は、男性と女性の間で大きく異なるのですか？」）これらの質問に答えると相互作用項を試験するために

。

Data <- data.frame(
    gender = sample (c("men", "women"), 2000, replace = TRUE), 
    var1 = sample (c("value1", "value2"), 2000, replace = TRUE), 
    var2 = sample (c("valueA", "valueB"), 2000, replace = TRUE), 
    y = sample(0:10, 2000, replace = TRUE) 
) 
mm <- lm(y~(var1+var2)*gender,Data) 

ここ

は、相互作用項である：

interax <- c("var1value2:genderwomen","var2valueB:genderwomen") 
printCoefmat(coef(summary(mm))[interax,]) 
##      Estimate Std. Error t value Pr(>|t|) 
## var1value2:genderwomen 0.20144 0.28241 0.7133 0.4758 
## var2valueB:genderwomen -0.15423 0.28266 -0.5456 0.5854