建站如何注重内容建设,网站建设服务费应该做到什么科目,海安建设银行网站,接订单去哪个网站目录 1.纯随机序列的定义
2.性质
3.纯随机性检验 1.纯随机序列的定义
纯随机序列也称为白噪声序列#xff0c;满足如下性质#xff1a;2.性质
纯随机性#xff08;无记忆性#xff09;方差齐性举例#xff0c;随机生成1000个白噪声序列
用正态分布序列 rnorm(数量满足如下性质2.性质
纯随机性无记忆性方差齐性举例随机生成1000个白噪声序列
用正态分布序列 rnorm(数量均值方差)如下为1000个标准正态的分布图形
a-rnorm(1000)
x-ts(a)
plot(x)时序图如下 自相关图
acf(x)返回 如图可以看出自相关系数基本分布在二倍标准差之间但由于数据是随机的所以有一定的误差
3.纯随机性检验
Bartlett定理如果一个时间序列是纯随机的得到一个观察期数为n的观察序列那么该序列的延迟非零期的样本自相关系数将近似服从均值为0方差为序列观察期数倒数的正态分布 原假设 备择假设 Q统计量Box和Pierce服从自由度为m的卡方分布其对大样本检验效果较好 LB统计量Box和LjungQ统计量的修正现在使用较普遍 拒绝域 R语言白噪声检验
Box.test(x,type,lag6)
其中
typeBox-Pierce Q统计量 默认
typeLjung-Box LB统计量
举例1对上面随机生成的白噪声序列进行6阶和12阶的LB统计量
for(i in 1:2)print( Box.test(x,typeLjung-Box,lag6*i))返回 举例2对1900年到1998年全球7级以上地震法伤次数序列进行平稳性和纯随机性检验
读取数据
a-read.table(D:/桌面/E2_5.csv,sep,,headerT)
a
返回 选择变量序列
x-ts(a$number,start1900)
x返回 绘制时序图
plot(x)绘制自相关图
acf(x)6阶LB统计量
Box.test(x,typeLjung-Box,lag6)返回
