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 통계를 마우스 오른쪽 버튼으로 클릭하고 "데이터 테이블로 만들기"를 선택합니다. 생성된 데이터 테이블에는 분석을 얻는 방법과 함께 테이블을 생성하는 방법을 알려주는 스크립트가 있습니다. (나는 1개 1개를 사면 깜짝 놀라게 된다.) 우리는 첫 번째 조각을 가지고 있다.
- 이제 그래프 빌더를 사용하여 원하는 대로 막대 차트를 그리고 원하는 차트를 얻은 후 스크립트를 생성하십시오. 세 번째 조각이 있습니다.
예를 들어 다음은 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-값도 없습니다.
