在SAS中，有一個功能可以在分析時間序列時生成白噪聲和平穩圖（ https://go.documentation.sas.com/doc/en/pgmsascdc/9.4_3.4/etsug/etsug_tfstart_sect028.htm )，其中顯示了 Ljung-Box 卡方統計量的（以圖形方式）概率。 在 JMP 中，Ljung-Box Q 統計是默認可見的，但似乎沒有辦法獲得 SAS 可以生成的相應圖形。 是否有生成這些圖的解決方法和/或是否有可編寫腳本的解決方案？
謝謝。
您的意思是您想要時間序列模型的殘差報告中 Ljung-Box Q 統計量的條形圖表示嗎？
除了腳本之外，我看不到任何其他內容。 但我不知道如何從頭開始編寫，但我可以通過以下方式將事情拼湊起來：
- 首先，您需要完成時間序列建模。
- 右鍵單擊 Ljung-Box 統計數據，然後選擇“Make into Data Table”。 生成的數據表有一個腳本，它告訴我如何生成表，以及我們如何獲得分析。 （我對買一送一感到驚訝。）我們有第一件。
- 現在使用Graph Builder根據需要繪製條形圖，並在獲得所需圖表後生成腳本。 我們有第三件。
例如，這裡是 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" )} )
)
);
生成此圖表
然後我將代碼拼接在一起。 在這裡我得到了什麼
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" )} )
)
);
我不能把行號放在上面，但我可以告訴你我添加了第 1、10 和 16 行。 所有其他行直接粘貼到。 在這個截圖中看到。
目前，時間序列平台對單位根檢驗的支持有限，並且不會自動生成殘差。 如果需要，您需要先保存殘差，然後僅對殘差運行時間序列。 報告的右上角有幾個單位根檢驗統計量。 它們也沒有 p 值。
