在gnmr包中使用实例（）

我试图用R包gnm来拟合一个广义的Lee-Carter模型。 我建立了一个数据框架

EXPO <- read.table("dati/Exposures.txt",header=TRUE,skip=0)
DEATH <- read.table("dati/Deaths.txt",header=TRUE,skip=2) 
base=data.frame(
D=DEATH$Total,
E=EXPO$Total,
X=as.factor(EXPO$Age),
T=as.factor(EXPO$Year))

---

我已经试着用公式中的“base$D/base$E”替换“D/E”，但没有用…

我是

gnm

的维护者。这似乎是

实例

中的一个bug，我将做一个注释来修复它。同时，您可以指定长格式的实例：

LC2 <- gnm(D/E ~ as.factor(X) + Mult(as.factor(X), as.factor(T), inst = 1) + 
           Mult(as.factor(X), as.factor(T), inst = 2),
           weights = E, family = binomial(link = "logit"), data = base)

---

LC2如果您事先计算了
DE=D/E
，这行吗
gnm
的
实例
可能无法正确处理复杂的响应。谢谢。我预先计算了一个新变量中的比例，结果很好！
LC2 <- gnm(D/E~ as.factor(X)+instances(Mult(as.factor(X),as.factor(T)),2),
                 weights=E,family=binomial(link="logit"),data=base).

Error in eval(expr, envir, enclos) : object 'D/E' not found

LC2 <- gnm(D/E ~ as.factor(X) + Mult(as.factor(X), as.factor(T), inst = 1) + 
           Mult(as.factor(X), as.factor(T), inst = 2),
           weights = E, family = binomial(link = "logit"), data = base)