Tengo la siguiente situación: General practitionar (gp), paciente (pat) y de consulta (contras).
Cada gp tiene varias paciente y de cada paciente puede tener 1 de consultas y más con las características específicas de cons_x. El resultado es y=0/1 en el nivel de consulta.
Un posible modelo podría ser (en lmer/R notación)
y ~ gp_sex + pat_sex + pat_age + cons_x + (1|gp) + (1|gp:pat)
gp es una variable aleatoria y pat está anidado en el gp. Dado que el resultado es 0/1 (derivación a especialista sí/no) el modelo adecuado es un multinivel de regresión logística.
(1) El problema: el 50% de los pacientes tienen sólo una consulta, es decir, no hay repetición, otros tienen 2, 3 hasta 6 consultas.Me gustaría saber cómo lidiar con este tipo de situación. Con otras palabras: he grupos (pacientes) con una única consulta. Así, dentro del grupo es perfecta correlación debido a que el único miembro en el grupo. ¿Los modelos mixtos marco de lidiar con esto?
(2) Otra pregunta que surge es: El número de consultas pesos el efecto en el paciente, ya que el mismo paciente tiene tantas entradas como se ha consultas. El más consulta un paciente tiene más características de los pacientes son considerados. Este no sería el caso, si todos los pacientes tienen el mismo número de consultas.
Cualquier ayuda y sugerencia se agradece.
Actualización:
Para simplificar el análisis podría agregado en el nivel del paciente: Si el paciente tiene varias consultas y al menos una consulta tiene 1 entonces el paciente ha resultado 1 caso contrario 0. Entonces tengo un modelo más simple:
y ~ gp_sex + pat_sex + pat_age + (1|gp)
En este caso el problema (1) será obsoleto pero que todavía tengo problema (2), ya que cada gp tiene diferente número de pacientes.