p-ts(data$水位，start=c(1860 )，frequency=1) #Draw a sequence diagramplot(p ) p， type='o ' ) #drawautocorrelation lag.max=25 ) # checkwhetherthesequenceisapurerandomsequencebox.test (p，lag=12，p )
执行结果
Box-Ljung test的糟糕唇膏pdata: p
X-squared=329.2，df=12，p-value 2.2e-16