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