Saludos Estoy utilizando el PROC GLIMMIX para ajustar un modelo conjunto data WORK.IMPORT_V;
length dist $7;
set WORK.IMPORT;
response=(VAR16=1);
dist="Binary";
output;
response = Rep;
dist="Poisson";
output;
keep APR APT Diagnostico Cursos IDs Escuela response dist ;
run;
proc glimmix data= work.import_v (where= (dist="Binary"));
class IDs Escuela Cursos;
model response(event='1') = APR APT Diagnostico / s dist=binary;
run;
proc glimmix data=work.import_v(where=(dist="Poisson"));
class IDs Escuela Cursos;
model response = APR APT Diagnostico / s dist=Poisson;
run;
proc glimmix data=work.import_v;
class dist;
class IDs Escuela Cursos;
model response(event='1') = dist dist*APR dist*APT dist*Diagnostico /
noint s dist=byobs(dist);
run;
proc glimmix data=work.import_v;
class IDs dist;
class Escuela Cursos;
model response(event='1') = dist dist*APR dist*APT dist*Diagnostico /
noint s dist=byobs(dist);
random int / subject=IDs;
run;
El modelo conjunto consta de dos submodelos, un logístico y otra poisson. La salida de dicho modelo me dice que la matriz G estima que no es claramente positiva. Adjunto imagen con la salida. La verdad no se como interpretar esta salida. Si tienes la oportunidad, también me gustaría contar con apoyo literario ya que no encuentro suficiente informacion del tema para aclarar mis dudas. Mis datos se encuentran organizados de la siguiente manera: Identificaciones Cursos Notas Reps Sexo Ano Escuela Distrito IGS APTO ABR Nivel1 Nivel2 Nacimiento Diagnostico VAR16 00BA1F4E0C MATE_3171 W 1 F 2013 Privada Mayaguez 323 544 542 0 0 1995 0 0 00D72CEB8A MATE_3171 F 2 F 2013 Publica Mayaguez 305 531 528 0 0 1995 0 0 00A362D565 MATE_3171 W 1 F 2013 Publica Mayaguez 305 531 528 0 0 1995 0 0 00AA00E246 MATE_3172 hacer 1 F 2013 Publica Mayaguez 305 531 528 0 0 1995 0 1 00CD4D51FA MATE_3171 hacer 1 METRO 2013 Publica Mayaguez 327 589 584 0 2 1994 0 1 00CD4D51FA MATE_3172 segundo 2 METRO 2013 Publica Mayaguez 327 589 584 0 2 1994 0 1 00CD4D51FA MATE_3031 W 1 METRO 2013 Publica Mayaguez 327 589 584 0 2 1994 0 0 Agradezco su amable atención y colaboración.
... View more
