这个问题
我试图弄清楚如何去做$ \ ce{二氧化碳}$集中趋势数据中计算地球系统研究实验室提供的莫纳罗亚山天文台。数据包含每月测量,季节性修正趋势,并插入数据填补缺口。趋势和插值数据应该从原始的可计算的平均水平。我试图复制计算趋势和插值的结果,但我保持稍微偏离值。
我已经试了两个公式计算的趋势,我的电话民用和formula2下面也再现了数据准确,虽然formula2似乎符合很好。
的数据
的数据(esrl.noaa.gov)可以作为一个文本文件结构如下:
#(-99.99丢失数据;1 #每日意味着月没有数据)# # #十进制平均插值趋势#天(季节corr) 1958 1958 1958.208 315.71 315.71 314.62 1 1958.292 317.45 317.45 315.29 1 1958 1958 1958.375 317.50 317.50 314.71 1 1958.458 -99.99 317.10 314.85 1。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。在这里平均平均每月测定的$ \ ce{二氧化碳}$浓度在ppm的单位,#天是每月的有效天数的测量是否可用。
其他两个列,内插和趋势是稍微复杂一些。文件描述的方法计算值如下。
“插入”列包含从前面列平均值和插值数据在哪里失踪。插值计算在两个步骤。第一个,我们为每个月平均计算季节性周期在每个月度值7窗口。以这种方式允许季节性周期变化缓慢。然后我们确定每个月的“趋势”的价值通过消除季节性周期;这个结果显示在“趋势”列中。趋势值是线性插值失踪几个月。插值月度的意思是然后平均季节周期值之和的趋势值丢失。
公式
为了避免处理边界影响我比较结果过去第一个3.5年的记录,以便集中windows 7能正确定义。从1961年9月。
让x_t美元,美元t = 1 . . N月平均的值。
我需要首先计算”7年的季节性周期平均窗口周围每个月度值”。然后减去从数据获取的趋势。
据我所知一般季节性周期的区别是季节性过滤器和一些季节性的平均水平。
计算7季节性的过滤器t \组3 * 12 + 6美元我计算$ s_t = \压裂{1}{7}\ Sigma_ {j = 3} ^ {j = 3}间{t + 12 j} $。
现在季节平均我试着两个选项。
- 平均7年季节性过滤器$ f_t = \压裂{1}{12}\ Sigma_ j = {0} ^ {j = 11} s_ {t + j} $
- 7年平均水平$ a_t = \压裂{1}{84}\ Sigma_ {j = -42} ^ {j = 41}间{t + j} $
从而民用是x_t——(s_t-f_t)美元和formula2是x_t——(s_t-a_t)美元。
结论
在我看来,我误会是什么意思”7年的季节性周期平均窗口周围每个月度值”。既不(s_t-f_t)美元,也不(s_t-a_t)美元上面定义的健康。正确的公式是什么?
我已经把数据放进一个只读的(谷歌表快速检查)。
