421429_zh-CN

liliconcom

生成白噪声和平稳图

在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 统计量的条形图表示吗？

除了脚本之外，我看不到任何其他内容。 但我不知道如何从头开始编写，但我可以通过以下方式将事情拼凑起来：

1. 首先，您需要完成时间序列建模。
2. 右键单击 Ljung-Box 统计数据，然后选择“Make into Data Table”。 生成的数据表有一个脚本，它告诉我如何生成表，以及我们如何获得分析。 （我对买一送一感到惊讶。）我们有第一件。
3. 现在使用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 值。