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