SASには、時系列を分析するときにホワイトノイズと定常プロットを生成する機能があります（ https://go.documentation.sas.com/doc/en/pgmsascdc/9.4_3.4/etsug/etsug_tfstart_sect028.htm ）、これはLjung-Boxカイ2乗統計の（グラフで）確率を示します。 JMPでは、Ljung-Box Q統計がデフォルトで表示されますが、SASが生成できる対応するグラフを取得する方法はないようです。 これらのプロットを生成するための回避策はありますか、および/またはスクリプト可能な解決策はありますか？
ありがとう。
時系列モデルの残差レポートにLjung-BoxQ統計の棒グラフ表現が必要であることを意味しますか？
とにかくスクリプト以外は表示されません。 しかし、最初から書く方法はわかりませんが、物事をつなぎ合わせるためにできることは次のとおりです。
- まず、時系列モデリングを完了する必要があります。
- Ljung-Box統計を右クリックし、[データテーブルに作成]を選択します。 生成されたデータテーブルには、テーブルの生成方法と分析の取得方法を示すスクリプトがあります。 （この1つ1つ無料で購入できることに驚いています。）最初のピースがあります。
- 次に、グラフビルダーを使用して必要に応じて棒グラフを描画し、必要なグラフを取得した後にスクリプトを作成します。 3つ目があります。
たとえば、これが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値もありません。
