¿Cómo dibujar gráficos B-spline mediante el procedimiento PROC GLIMMIX en SAS?

Question

Recientemente he descubierto que el PROC GLIMMIX en SAS añadió una declaración effect que puede manejar B-spline en los modelos. Lo he probado, y el programa parece:

proc glimmix data=temp;
    effect t = spline(time);
    model y = x1 x2 t / solution;
    random intercept / subject=id type=un;
run;

Está bien para obtener salidas, pero lo que realmente me interesa es el gráfico de la línea B. He revisado todo el manual, pero no he encontrado ninguna información. Su ODS GRAPH no puede generar automáticamente gráficos B-spline. Además, su ODS OUTPUT no proporciona ninguna salida para dibujar gráficos B-spline fuera. Me pregunto si alguien solía utilizar la función PROC GLIMMIX ¿procedimiento para dibujar un gráfico B-spline? Se agradece cualquier información.

Answer 1

Hago los trazados utilizando la declaración de estimación. Es decir, hacer algunos contrastes a diferentes valores de t y guardar la salida como una tabla. Entonces es sólo para trazar el valor estimado frente a t. Pero tal vez esto no es lo que usted está buscando? Si quieres una opción que te dé el gráfico como salida, no puedo ayudarte..

Answer 2

La forma más sencilla de hacerlo es añadir algunos registros adicionales al conjunto de datos temp para los valores deseados de x1, x2 y tiempo, pero con un valor perdido para la variable de resultado y. Se obtiene un conjunto de datos de salida para las predicciones utilizando una sentencia de salida como

output out=predicted pred=pred LCL=pLCL UCL=pUCL / alpha=0.05;

A continuación, puede trazar a partir del conjunto de datos creado. Este truco funciona con muchos procedimientos de SAS y, en general, es preferible a trabajar a partir de las estimaciones, porque la matriz de covarianza de los diferentes efectos del modelo se tiene en cuenta al calcular los errores estándar y los intervalos de predicción.