非平衡面板数据回归

连续观察患者治疗后1，2，3，6，12月某些指标(A,B,C等,为连续变量)的变化。25例患者。
结果变量为结局outcome，为二分类。如果能够同时考虑结局发生时间更好，不能考虑就是简单的二分类结果变量。
除了A,B,C等指标外，有其他的患者其他变量也会影响结果，这些变量不随时间变化，比如年龄，性别等。

我们的目的是，分析A,B,C是否影响结果变量，A，B,C之间有无相关性。其次年龄，性别等是否影响A，B，C的动态变化，从而影响最后结果outcome。

想请教，这种用什么方法比较合适，
广义估计 GEE 或者 广义线性混合效应GLMMs，或者其他更好办法?


请提供分析过程和命令。
谢谢！