我在R中运行OLS回归,从中得到一对系数。以下是部分代码:最小二乘回归系数非线性函数的标准误和置信区间
Attacks <- Treat.Terr.Dataset$Attacks[2:30]
Attackslag <- Treat.Terr.Dataset$Attacks[1:29]
TreatmentEffect <- Treat.Terr.Dataset$TreatmentEffect[2:30]
TreatmentEffectlag <- Treat.Terr.Dataset$TreatmentEffect[1:29]
olsreg <- lm(TreatmentEffect ~ TreatmentEffectlag + Attacks + Attackslag)
coeffs<-olsreg$coefficients
然后我需要计算:(Attacks + Attackslag)/(1 - TreatmentEffectlag)
。问题是我可以使用(coeffs[3] + coeffs[4])/(1 - coeffs[2])
在R上做这个,但结果是没有任何p值或置信区间的固定数字,就像计算器会显示我一样。
有谁知道是否有任何函数可以用来计算这个置信区间?
编者注
如果目标量是回归系数的线性函数,那么问题降低到一般的线性假设检验,其中精确推断是可能的。
'bootstrap'它。 – user20650
欢迎来到StackOverflow!请阅读关于[如何提出一个好问题](http://stackoverflow.com/help/how-to-ask)以及如何给出[可重现的示例]的信息(http://stackoverflow.com/questions/ 5963269 /如何对化妆一个伟大-R-重复性,例如/ 5963610)。这会让其他人更容易帮助你。 – Jaap
你打算如何使用这些?这将决定最佳的回应。 – Elin