時間序列1時序圖的繪制和純隨機性檢驗#R


一、繪制時序圖的參數設置

plot(x,tpye=,pch=,lty=,lwd=,col=,main=,sub=,xlab=,ylab=,xlim=,ylim=)

   type=來控制點線輸出結構。type=“l”為線,type=“o”為線穿過點。

參數取值

描述

type=“p”

type=“l”

type=“b”

點連線

type=“o”

線穿過點

type=“h”

懸垂線

type=“s”

階梯線

   pch=選項設置觀察點的符號。pch=1為空心圓,pch=2為空心三角形。

 

 

   lty=選項設置線的類型

 

  lwd=選項設置線的寬度

  col=選項設置點線顏色

  xlim=選項指定橫坐標的范圍,用ylim=選項指定縱坐標的范圍

  abline()函數為圖形添加參考線

     (1)v= 指定垂直參照線,如果需要多條參照線,可以使用v=c(a1,a2,)指定多條垂直參照線。例:

abline(h=c(15.5,16.6), lty=2)

    (2)h= 指定水平參照線,h=c(a1,a2,)指定多條水平參照線。

    (3)lty= 指定線條的類型。

    (4)lwd= 指定線條的寬度。

    (5)col= 指定線條的顏色。

二、繪制序列自相關圖

acf(x, lag=)
  
 

 x:變量名

 lag:延遲階數,用戶不特殊指定的話,系統會根據序列長度自動指定延遲階數。例:

acf(yield)

三、純隨機性檢驗

Box.test(x, type= ,lag= )

   x:變量名

   type:檢驗統計量類型。type=“Box-Pierce”輸出白噪聲檢驗的Q統計量。該統計量為系統默認輸出結果。type=“Ljung-Box”輸出白噪聲檢驗的LB統計量。

   lag:延遲階數。lag=n 表示輸出滯后n階的白噪聲檢驗統計量。忽略該選項時,默認輸出滯后1階的檢驗統計量結果。

   可以編寫一個循環命令。最簡單的循環函數是for函數,它的命令格式如下:

for(i in n1:n2 ) print( Box.test(yield,lag=i))

四、產生隨機序列

    正態分布隨機數生成函數:

rnorm(n=,mean=,sd= )
   -n:隨機數個數

  -mean:均值

  -sd:標准差

    如果要產生n個服從標准正態分布的隨機數,可以簡寫為rnorm(n )

 


免責聲明!

本站轉載的文章為個人學習借鑒使用,本站對版權不負任何法律責任。如果侵犯了您的隱私權益,請聯系本站郵箱yoyou2525@163.com刪除。



 
粵ICP備18138465號   © 2018-2025 CODEPRJ.COM