2

一般化された Lee-Carter モデルを R パッケージ gnm に当てはめようとしています。データフレームを構築しました

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))

さて、これが機能している間

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

これに移動すると、エラーが発生します

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

数式で「D/E」を「base$D/base$E」に置き換えようとしましたが、うまくいきませんでした...

4

1 に答える 1

1

私は のメンテナーですgnm。instancesこれは、私が修正するメモを作成するバグのようです。それまでの間、インスタンスを長い形式で指定できます。

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)

または、Josh が提案するようにプロポーションを事前に計算することもできます。

于 2015-08-04T08:17:58.607 に答える