[R-br] Extraindo formula do preditor linear GEE

[Pedro Emmanuel Alvarenga Americano do Brasil]

Amigos de R,

Venho procurar uma luz. Tenho certeza que é alguma coisa ridícula que não
consigo ver. Eu gostaria de adicionar num relatório a fórmula do preditor
linear de um modelo GEE. No pacote Hmisc há uma função Formula que faz
isso, mas não funciona para objetos geeglm. Além disso, a formula possui um
termo spline e depois de um tempo fuçando por ai eu não consegui fazer a
formula funcionar.

fitar1_11 <- geeglm(outcome ~ v1 + v2 + v4_recoded + v5_trunc + bs(weeks,
df = 3), family = binomial(), data = M2, id = f1p0,  corstr = "ar1")
fitar1_11
geeglm(formula = outcome ~ v1 + v2 + v4_recoded + v5_trunc +
    bs(weeks, df = 3), family = binomial(), data = M2, id = f1p0,
    corstr = "ar1")
Coefficients:
                           (Intercept)
v1
                         -1.586386e+00
8.363458e-02
                                    v2 v4_recodedNET, NGT, JF, OET, GF or
TPN
                         -3.763166e-02
4.503599e+15
                              v5_trunc                     bs(weeks, df =
3)1
                          7.261799e-02
3.469451e-03
                    bs(weeks, df = 3)2                     bs(weeks, df =
3)3
                          3.756985e-01
4.542696e+00

Degrees of Freedom: 472 Total (i.e. Null);  464 Residual
> nd # newdata to predict at
  v1 v2  v4_recoded v5_trunc weeks
1  8  3 No catheter       36     2
> predict(fitar1_11, type = "response", newdata = nd) # prediction from the
model
        1
0.8345324
# Not working prediction from the hand made formula
> plogis(-1.586386 + 0.08363458*nd$v1 + -0.03763166*nd$v2 +
4503599142291968*0 + 0.07261799*nd$v5_trunc + 0.003469451*nd$weeks +
0.3756985*pmax(nd$weeks-1,0)^3 + 4.542696*pmax(nd$weeks-8,0)^3)
[1] 0.8772541

Alguem consegue me dar uma luz o que está errado para eu acertar esta
previsão? Se puder indicar uma leitura que resolva eu gostaria também.

Abraço forte e que a força esteja com voces,

Pedro Brasil
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20210618/e8ddc57c/attachment.html>