R中nlcom(Stata)的等价物?回归系数的非线性变换
Equivalent of nlcom (Stata) in R? Nonlinear transformations of regression coefficients
我想对回归系数进行非线性变换。例如:
,或
。
Stata 有一个方便的实现 nlcom
this,它使用 delta 方法来估计标准误差和相应的置信区间。我知道发布的简单转换可以通过直接处理模型中的感兴趣系数来简单地完成。但是,如果我们对几个线性和非线性组合的比率感兴趣,那么在这样的转换中产生置信区间的有效方法是什么?此外,当系数有一个完整的 co-variance 矩阵时,标准误差与它们一起估计。
为了回答我自己的问题,我发现 library(msm)
包的功能 deltamethod()
很好地满足了我的要求。加州大学洛杉矶分校对这种方法的描述非常好,所以我为可能有类似需求的任何人提供 link。
Using the delta method for nonlinear transformations of regression coefficients.
包 car
中的 deltaMethod()
函数也完成相同的工作,提供估计值、标准误差和 95% 置信区间作为其输出。
我想对回归系数进行非线性变换。例如:
,或
。
Stata 有一个方便的实现 nlcom
this,它使用 delta 方法来估计标准误差和相应的置信区间。我知道发布的简单转换可以通过直接处理模型中的感兴趣系数来简单地完成。但是,如果我们对几个线性和非线性组合的比率感兴趣,那么在这样的转换中产生置信区间的有效方法是什么?此外,当系数有一个完整的 co-variance 矩阵时,标准误差与它们一起估计。
为了回答我自己的问题,我发现 library(msm)
包的功能 deltamethod()
很好地满足了我的要求。加州大学洛杉矶分校对这种方法的描述非常好,所以我为可能有类似需求的任何人提供 link。
Using the delta method for nonlinear transformations of regression coefficients.
包 car
中的 deltaMethod()
函数也完成相同的工作,提供估计值、标准误差和 95% 置信区间作为其输出。