添加系数以获得交互作用-SE怎么办?


13

我有一个多元回归,其中包括相互作用。例如,要估算最差的五分位数的治疗效果,我需要将来自治疗回归指标的系数与来自交互变量(与治疗和五分位数1相互作用)的系数相加。通过回归将两个系数相加时,如何获得标准误差?是否可以将两个系数的标准误差相加?那t统计呢?是否可以添加这些?我猜不是,但是我找不到任何指导。

提前非常感谢您的帮助!


这真的很有帮助!我希望在R中做类似的事情,但各组之间的样本大小略有不同。我是否仍可以使用相同的方程式将这两个误差结合起来以得出新的标准品。错误?预先感谢您的帮助Crystal
Crystal

1
嘿@Crystal-欢迎来到网站!这是一个很好的问题,但是您应该将其作为一个新问题(右上角的“问问题”按钮)。现在,您已经提交了它作为对此旧问题的“答案”。如果您仅将此问题的URL复制并粘贴到新问题中,我们都会理解您在说什么。
Matt Parker

Answers:


10

我认为这是的表达式:SEbnew

SE12+SE22+2Cov(b1,b2)

您可以使用这个新的标准错误来找到用于测试H o的新测试统计量β = 0Ho:β=0


您好莎拉,如果您认为已回答,则应关闭此问题。
suncoolsu

嗨-再次感谢您的回答。我忘了提到我正在使用Stata。当我将两个系数加在一起(使用Stata的输出)时,是否也可以将标准误差相加?如果是这样,那么我应该能够通过将系数的总和除以标准误差的总和来获得标准误差。你同意吗?再次感谢。
莎拉(Sarah)2010年

Sarah,在Stata中,使用“ lincom”功能。假设您有变量var1和var2,并想将var1上的系数乘以3倍,并将var2上的系数乘以2倍。输入'lincom 3 * var1 + 2 * var2'。这给出了此估计的标准误差和置信区间。
查理

5

我假设您的意思是“多变量”回归,而不是“多变量”。“多变量”是指具有多个因变量。

采取连续的预测变量并将其分成若干间隔不被认为是可接受的统计实践。这将导致残留的混淆,并使交互产生误导性的意义,因为某些交互只能反映出某些主要影响的不恰当(此处为欠拟合)。在外部五分位数中有很多无法解释的变化。另外,实际上不可能精确地解释“五分之一效应”。

对于感兴趣的比较,最容易将它们设想为预测值的差异。这是使用R rms包的示例。

require(rms)
f <- ols(y ~ x1 + rcs(x2,3)*treat)  # or lrm, cph, psm, Rq, Gls, Glm, ...
# This model allows nonlinearity in x2 and interaction between x2 and treat.
# x2 is modeled as two separate restricted cubic spline functions with 3
# knots or join points in common (one function for the reference treatment
# and one function for the difference in curves between the 2 treatments)
contrast(f, list(treat='B', x2=c(.2, .4)),
            list(treat='A', x2=c(.2, .4)))
# Provides a comparison of treatments at 2 values of x2
anova(f) # provides 2 d.f. interaction test and test of whether treatment
# is effective at ANY value of x2 (combined treat main effect + treat x x2
# interaction - this has 3 d.f. here)

2

RRβRV^RV^Rχr2r 是矩阵中的行数(假设这些行是线性独立的)。


谢谢。由于我不是统计专家,因此我将提出另一个问题,并且不确定我的问题是否明确。
萨拉(Sarah)2010年
By using our site, you acknowledge that you have read and understand our Cookie Policy and Privacy Policy.
Licensed under cc by-sa 3.0 with attribution required.