<div dir="ltr"><div class="gmail_default" style="font-family:monospace,monospace;font-size:large;color:#073763"><div class="gmail_default">Amigos de R, </div><div class="gmail_default"><br></div><div class="gmail_default">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. </div><div class="gmail_default"><br></div><div class="gmail_default"><font size="1">fitar1_11 <- geeglm(outcome ~ v1 + v2 + v4_recoded + v5_trunc + bs(weeks, df = 3), family = binomial(), data = M2, id = f1p0,  corstr = "ar1")<br></font></div><div class="gmail_default"><font size="1">fitar1_11 </font>   <font size="1"><br></font></div><div class="gmail_default"><font size="1">geeglm(formula = outcome ~ v1 + v2 + v4_recoded + v5_trunc + <br>    bs(weeks, df = 3), family = binomial(), data = M2, id = f1p0, <br>    corstr = "ar1")</font></div><div class="gmail_default"><font size="1">Coefficients:<br>                           (Intercept)                                     v1 <br>                         -1.586386e+00                           8.363458e-02 <br>                                    v2 v4_recodedNET, NGT, JF, OET, GF or TPN <br>                         -3.763166e-02                           4.503599e+15 <br>                              v5_trunc                     bs(weeks, df = 3)1 <br>                          7.261799e-02                           3.469451e-03 <br>                    bs(weeks, df = 3)2                     bs(weeks, df = 3)3 <br>                          3.756985e-01                           4.542696e+00 </font><br><br><font size="1">Degrees of Freedom: 472 Total (i.e. Null);  464 Residual<br></font></div><div class="gmail_default"><font size="1">> nd # newdata to predict at<br>  v1 v2  v4_recoded v5_trunc weeks<br>1  8  3 No catheter       36     2<br>> predict(fitar1_11, type = "response", newdata = nd) # prediction from the model<br>        1 <br>0.8345324</font></div><div class="gmail_default"><font size="1"># Not working prediction from the hand made formula<br>> 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)<br>[1] 0.8772541</font><br></div><div class="gmail_default"><br></div><div class="gmail_default">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. </div><div class="gmail_default"><br></div><div class="gmail_default">Abraço forte e que a força esteja com voces,  </div><div class="gmail_default"><br></div><div><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><font face="'courier new', monospace" color="#000066" size="4">Pedro Brasil</font></div></div></div></div></div></div></div></div></div></div>