In SAS, esiste una funzione per generare rumore bianco e grafici stazionari durante l'analisi delle serie temporali ( https://go.documentation.sas.com/doc/en/pgmsascdc/9.4_3.4/etsug/etsug_tfstart_sect028.htm ), che mostrano (graficamente) le probabilità della statistica del chi quadrato di Ljung-Box. In JMP, la statistica Ljung-Box Q è visibile come predefinita, ma non sembra esserci un modo per ottenere i grafici corrispondenti che SAS può generare. Esiste una soluzione alternativa per generare questi grafici e/o esiste una soluzione scriptabile?
Grazie.
Intendi dire che vuoi la rappresentazione del grafico a barre delle statistiche Ljung-Box Q nel rapporto sui residui di un modello di serie temporale?
Non vedo comunque altro che scripting. Ma non so come scrivere da zero, ma ecco cosa posso fare per mettere insieme le cose:
- Per prima cosa devi completare la modellazione delle serie temporali.
- Fare clic con il pulsante destro del mouse sulle statistiche di Ljung-Box e scegliere "Crea tabella dati". La tabella dei dati generata ha uno script, che mi dice come generare la tabella, insieme a come otteniamo l'analisi. (Sono sorpreso da questo compra uno prendi uno gratis.) Abbiamo il primo pezzo.
- Ora usa Graph Builder per disegnare il grafico a barre come desideri e produrre uno script dopo aver ottenuto il grafico desiderato. Abbiamo il terzo pezzo.
Ad esempio, ecco lo 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" )} )
)
);
che genera questo grafico
Poi cucio insieme i codici. Ed ecco cosa ottengo
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" )} )
)
);
Non riesco a mettere i numeri di riga sopra, ma posso dirti che ho aggiunto la riga 1, 10 e 16. Tutte le altre linee vengono incollate direttamente. Vedendo in questo screenshot.
Al momento, la piattaforma Time Series ha un supporto limitato per i test di root delle unità e non produce automaticamente i residui. Se lo desideri, devi prima salvare i residui, quindi eseguire Time Series solo sui residui. Diverse statistiche sui test della radice dell'unità si trovano in alto a destra nel report. Non ci sono nemmeno valori p che li accompagnano.
