給定 B1 的值，如何使用 Predict() 線性 model

Question

library(wooldridge)
library(tidyverse)
library(stargazer)
setwd("C:/Users/Charlie/Desktop/R Homework")

data(wage1)

reg_wage1 <- lm(lwage ~ female + educ + (female * educ), data = wage1)

stargazer(reg_wage1, type = "text")

female_at_zero <- data.frame(female=0)
pred <- predict(reg_wage1, female_at_zero)

stargazer(pred, type = "text")

我的問題是，當我嘗試運行此代碼時，它一直要求我輸入 educ 的值，但我不想更改教育，我只想在女性等於 0 時在 stargazer 中查看模型的結果。

Answer 1

您遇到的問題是您已經使用educ定義了 model reg_wage1 ，因此它不知道如何在沒有educ值的情況下進行預測。 lm已根據educ變量的值定義了female變量（更不用說您添加到模型中的交互變量）。 如果你想看到female = 0對lwage的影響，你將不得不重新定義你的 model。

另一種選擇是簡單地調用female_at_zero <- data.frame(female = 0, educ = 0) 。 通過將educ設置為 0，您實際上只是在查看female的效果。 但是這個解決方案是多余的，因為pred的值將只是stargazer(reg_wage1, type = "text")的 output 中Constant的值。 這是因為線性 model ( lm ) 將其截距 ( Constant ) 定義為所有自變量都等於 0 的點。