楼主: superhugo
9419 10

[问答] 如果时间序列不存在明显趋势Trend,是否还能用STL分解?其分解效果如何? [推广有奖]

  • 0关注
  • 0粉丝

硕士生

50%

还不是VIP/贵宾

-

威望
0
论坛币
7145 个
通用积分
3.3519
学术水平
8 点
热心指数
8 点
信用等级
7 点
经验
2592 点
帖子
131
精华
0
在线时间
144 小时
注册时间
2011-4-8
最后登录
2023-2-27

相似文件 换一批

+2 论坛币
k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群
赵安豆老师微信:zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

感谢您参与论坛问题回答

经管之家送您两个论坛币!

+2 论坛币
R中有STL分解,其函数描述为:
Decompose a time series into seasonal, trend and irregular components using loess, acronym STL.

感谢之前ryusukekenji 的热心回帖,解决了时间序列的定义问题。再把进一步的问题列示如下:

定义一个时间序列y,这个时间序列由若干年若干天组成,中间不存在缺失值。其周期为年,包含365天。所以定义了n(p)=365,对其进行stl分解。得到以下结果:

截图00.jpg

从结果来看,感觉STL分解效果很差。参数的设置修改了好几次,效果都不理想。想问一下,对于不存在趋势项的情况,只是想去除序列在各年同一天的季节项,这个STL分解还适用吗?
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

关键词:trend 时间序列 不存在 STL End 如何

沙发
ryusukekenji 发表于 2012-2-25 20:24:02 |只看作者 |坛友微信交流群
时间变量有误,楼主需要仔细查看...

使用道具

藤椅
superhugo 发表于 2012-2-27 11:07:36 |只看作者 |坛友微信交流群
ryusukekenji 发表于 2012-2-25 20:24
时间变量有误,楼主需要仔细查看...
stl(x, s.window, s.degree = 0,
    t.window = NULL, t.degree = 1,
    l.window = nextodd(period), l.degree = t.degree,
    s.jump = ceiling(s.window/10),
    t.jump = ceiling(t.window/10),
    l.jump = ceiling(l.window/10),
    robust = FALSE,
    inner = if(robust)  1 else 2,
    outer = if(robust) 15 else 0,
    na.action = na.fail)

对变量x的说明如下:
x  univariate time series to be decomposed. This should be an object of class "ts" with a frequency greater than one.

能解释一下后面那句话的意思吗:“This should be an object of class "ts" with a frequency greater than one.”

使用道具

板凳
superhugo 发表于 2012-2-27 11:08:19 |只看作者 |坛友微信交流群
请用过 stl 函数的朋友说一下这个函数中的x应该怎么定义?

使用道具

报纸
ryusukekenji 发表于 2012-2-27 13:54:26 |只看作者 |坛友微信交流群
superhugo 发表于 2012-2-27 11:08
请用过 stl 函数的朋友说一下这个函数中的x应该怎么定义?
x变量必须是ts时间序列属性,并且超过一次频率、亦重复时间(月份、周日等变量)
This should be an object of class "ts" with a frequency greater than one.

例如:
            Sun Mon Tue Wed Thu Fri Sat
week 1
week 2
...


         Jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov Dec
1920
1921
1922
...
  1. >nottem
  2. >attributes(nottem)
  3. >class(nottem)
  4. >stl(nottem, s.window=2)
复制代码
已有 1 人评分学术水平 热心指数 信用等级 收起 理由
superhugo + 1 + 1 + 1 热心帮助其他会员

总评分: 学术水平 + 1  热心指数 + 1  信用等级 + 1   查看全部评分

使用道具

地板
superhugo 发表于 2012-2-27 21:56:59 |只看作者 |坛友微信交流群
感觉STL分解s.window的设置比较不确定,论文中是说的该值≥7,且为奇数。如果不考虑趋势项 还可以用STL吗?

使用道具

7
superhugo 发表于 2012-2-29 10:21:26 |只看作者 |坛友微信交流群
用R软件实现了co2浓度序列的STL分解,见图。效果应该还是不错的,见图1。考虑robust与不考虑robust情况下得到的结果基本一致。
但STL分解对其它趋势trend不显著或很小、季节性seasonal很强的数据,似乎效果不是很好。

截图00.png (20.81 KB)

截图00.png

使用道具

8
zppn123 发表于 2014-6-18 10:42:42 |只看作者 |坛友微信交流群
superhugo 发表于 2012-2-29 10:21
用R软件实现了co2浓度序列的STL分解,见图。效果应该还是不错的,见图1。考虑robust与不考虑robust情况下得 ...
楼主还在不?我现在也在用这个,我分解的是季节性很强的数据,有什么指标可以判断分解的好不好吗?

使用道具

9
lsc989818 在职认证  发表于 2015-3-16 15:38:14 |只看作者 |坛友微信交流群
同样也是这个问题,如何判断STL分解出来结果的好坏了?

使用道具

10
lsc989818 在职认证  发表于 2015-3-20 14:41:59 |只看作者 |坛友微信交流群
还有,STL方法应该不能用来分解年度数据吧

使用道具

您需要登录后才可以回帖 登录 | 我要注册

本版微信群
加好友,备注cda
拉您进交流群

京ICP备16021002-2号 京B2-20170662号 京公网安备 11010802022788号 论坛法律顾问:王进律师 知识产权保护声明   免责及隐私声明

GMT+8, 2024-4-28 20:19