我有一个 R glm 对象,我想评估与新数据帧的单个输入变量相关联的系数。这是一个例子:
library(statmod)
data = tibble(X = c('a', 'b', 'c', 'd', 'e'), Z = c('a', 'c', 'c', 'c', 'a'), W = c(21:25), Y = c(0.1, 0.2, 0.3, 0.4, 0.88))
model = glm(Y ~ I(X %in% c('a', 'b')) + I(X %in% c('d', 'e')) + I(Z %in% c('a', 'c')) + I(log(W)), data=data, family=tweedie(var.power=1.44, link.power=0))
newdata = tibble(X = c('f', 'a'), Z = c('a', 'a'), W = c(5, 15))
我想提取与 X = 'f' 相关的系数(应该是基础级别,0)。我想提取与 W = 5 相关的系数(应该是 log(5) * beta)
谢谢,
约翰