【问题标题】:R lm Capture interaction terms, but not categorical variableR lm 捕获交互项,但不是分类变量
【发布时间】:2020-08-07 11:14:19
【问题描述】:

我想估计以下回归模型: y = b0 + b1 * X + b2 * x * 虚拟

其中 y 和 x 是连续的,dummy 是分类(虚拟变量)。

换句话说,我希望我的估计模型估计三个系数:bo、b1 和 b2。

我已经尝试了以下...

lm(y ~ x + x * dummy, data)

但它在模型中添加了变量dummy,并估计了dummy的系数。

以下内容接近我想要做的,但它将交互项转换为二进制变量(真/假)。

lm(y ~ x + I(!x * dummy), data)

对于复制,请考虑以下示例:

data <- tibble(y=rnorm(10), x=runif(10), dummy=ifelse(x>.5,1,0))
lm(y ~ x + x * dummy, data)
lm(y ~ x + I(!x * dummy), data)

谢谢

【问题讨论】:

    标签: r lm


    【解决方案1】:

    这里:

    > summary(lm(y ~ x+ x : dummy, data))
    
    Call:
    lm(formula = y ~ x + x:dummy, data = data)
    
    Residuals:
         Min       1Q   Median       3Q      Max 
    -0.61312 -0.15558 -0.00354  0.23965  0.47351 
    
    Coefficients:
                Estimate Std. Error t value Pr(>|t|)
    (Intercept)  0.06755    0.36162   0.187    0.857
    x            0.94953    1.18299   0.803    0.449
    x:dummy     -1.10220    0.88112  -1.251    0.251
    
    Residual standard error: 0.4148 on 7 degrees of freedom
    Multiple R-squared:  0.2645,    Adjusted R-squared:  0.05438 
    F-statistic: 1.259 on 2 and 7 DF,  p-value: 0.3412
    

    【讨论】:

      猜你喜欢
      • 2018-06-11
      • 2019-02-27
      • 1970-01-01
      • 2021-08-07
      • 2014-02-14
      • 2015-01-28
      • 1970-01-01
      • 1970-01-01
      • 2018-03-27
      相关资源
      最近更新 更多