1

我使用 delta 方法从 glm 拟合估计两个系数之间的差异（下面附上代码）。现在，我想将此估计值与零进行比较（即，没有差异的零假设）。一篇文章提到使用单样本位置、双向近似 Z 检验来计算这种差异。

但是，我找不到使用增量差异在 R 中执行此操作的简单方法。我查看了两个样本的 Z 测试文档，并可能考虑在 z-stat 公式中使用差异作为替代......但我不确定这是否是最好的方法。

##GENERATE DATA SET
y <- c(1:12)
x1 <- rep(c(1000, 4000, 0), each = 4)
x2 <- rep(c(0, 1000, 4000), each = 4)
df <- data.frame(y, x1, x2)

##RUN GLM
library(lmerTest)
g1 <- glm(log(y) ~ x1 + x2, data = df)

##Use delta-method to estimate the difference between coefficients of x1 and x2 (Ritz & Streibig 2008)
library(car)
g1.delta <- deltaMethod(g1,"(-x1) - (-x2)")


              Estimate   SE         2.5 %       97.5 % 
(-x1) - (-x2) 2.3217e-04 7.3180e-05 8.8738e-05  4e-04
4

0 回答 0