En SAS, hay una función para generar ruido blanco y gráficos estacionarios al analizar series de tiempo ( https://go.documentation.sas.com/doc/en/pgmsascdc/9.4_3.4/etsug/etsug_tfstart_sect028.htm ), que muestran (gráficamente) las probabilidades del estadístico chi cuadrado de Ljung-Box. En JMP, la estadística Ljung-Box Q es visible de forma predeterminada, pero no parece haber una forma de obtener los gráficos correspondientes que SAS puede generar. ¿Existe una solución alternativa para generar estos gráficos y / o existe una solución que se pueda programar?
Gracias.
¿Quiere decir que desea la representación del gráfico de barras de las estadísticas de Ljung-Box Q en el informe de Residuos de un modelo de serie temporal?
No veo otra forma que no sea el guion. Pero no sé cómo escribir desde cero, pero esto es lo que puedo hacer para juntar las cosas:
- Primero necesitas hacer tu modelado de series de tiempo.
- Haga clic con el botón derecho en las estadísticas de Ljung-Box y seleccione "Convertir en tabla de datos". La tabla de datos generada tiene un script, que me dice cómo generar la tabla, junto con cómo obtenemos el análisis. (Estoy sorprendido por este compre uno y llévese otro gratis). Tenemos la primera pieza.
- Ahora use Graph Builder para dibujar el gráfico de barras como desee y genere un script después de obtener el gráfico que desea. Tenemos la tercera pieza.
Por ejemplo, aquí está el script GB
Graph Builder(
Variables(
X( :"Ljung-Box Q"n ),
Y( Transform Column( "Row", Ordinal, Formula( Row() ) ) )
),
Elements( Bar( X, Y, Legend( 6 ) ) ),
SendToReport(
Dispatch( {}, "Row", ScaleBox, {Reversed Scale} ),
Dispatch( {}, "Y title", TextEditBox, {Set Text( "Lag" )} )
)
);
que genera este gráfico
Luego coso códigos juntos. Y aqui lo que consigo
dt = open("$sample_data/time series/seriesg.jmp");
Names Default To Here( 1 );
platform = Data Table( "Seriesg" ) << Time Series(
Y( :Log Passengers ),
Variogram( 1 ),
Seasonal ARIMA( 0, 1, 1, 0, 1, 1, 12, No Intercept( 1 ) )
);
Wait( 0 );
dt2 =
Report( platform )[Outline Box( "Time Series Log Passengers" )][
Outline Box( "Model: Seasonal ARIMA(0, 1, 1)(0, 1, 1)12 No Intercept " )][
Outline Box( "Residuals" )][Table Box( 2 )] << Make Into Data Table;
Report( platform ) << Close Window;
dt2 <<
Graph Builder(
Variables(
X( :"Ljung-Box Q"n ),
Y( Transform Column( "Row", Ordinal, Formula( Row() ) ) )
),
Elements( Bar( X, Y, Legend( 6 ) ) ),
SendToReport(
Dispatch( {}, "Row", ScaleBox, {Reversed Scale} ),
Dispatch( {}, "graph title", TextEditBox, {Set Text( "Ljung-Box Q" )} ),
Dispatch( {}, "Y title", TextEditBox, {Set Text( "Lag" )} )
)
);
No puedo poner los números de línea arriba, pero puedo decirles que agregué las líneas 1, 10 y 16. Todas las demás líneas se pegan directamente en. Viendo en esta captura de pantalla.
Por el momento, la plataforma Time Series tiene un soporte limitado de pruebas de raíz unitaria y no produce residuos automáticamente. Si lo desea, primero debe guardar los residuales y luego ejecutar Time Series solo en los residuales. Varias estadísticas de prueba de raíz unitaria se encuentran en la parte superior derecha del informe. Tampoco hay valores p que los acompañen.
