我使用包delthamethod中的msm来推导转换变量的标准错误。
示例代码：
require(msm)
x1 <- 1:10
x2 <- c(0,0,0,0,0,0,0,0,0,0)
y <- c(1,3,3,4,5,7,7,8,9,10)
m1 <- lm(y~x1+x2)
summary(m1)
deltamethod(~ (1-x1), coef(m1), vcov(m1))
我得到的错误是＃34;协方差应该是3x3矩阵＆＃34;。原因是1变量没有任何变化（x2总是为零）并且有＃34; NA＆＃34;在回归输出中。
这有一个简单的解决方法吗？我知道我可以将变量排除在外，但是我运行了超过1.000的回归，每个参数估计大约有15个，并且NA变量（没有变化）每次都是不同的变量。
答案 0 :(得分：1)
怎么样：
deltamethod(~(1-x1), na.omit(coef(m1)), vcov(m1))
# [1] 0.2949063