当前位置: 首页>>代码示例 >>用法及示例精选 >>正文


R sunspot.year 年度太阳黑子数据,1700–1988


R语言 sunspot.year 位于 datasets 包(package)。

说明

1700 至 1988 年间每年的太阳黑子数量(四舍五入到一位数)。

请注意,每月数字以 sunspot.month 形式提供,但开始时间稍晚。

用法

sunspot.year

格式

单变量时间序列 sunspot.year 包含 289 个观测值,属于 "ts" 类。

例子

utils::str(sm <- sunspots)# the monthly version we keep unchanged
utils::str(sy <- sunspot.year)
## The common time interval
(t1 <- c(max(start(sm), start(sy)),     1)) # Jan 1749
(t2 <- c(min(  end(sm)[1],end(sy)[1]), 12)) # Dec 1983
s.m <- window(sm, start=t1, end=t2)
s.y <- window(sy, start=t1, end=t2[1]) # {irrelevant warning}
stopifnot(length(s.y) * 12 == length(s.m),
          ## The yearly series *is* close to the averages of the monthly one:
          all.equal(s.y, aggregate(s.m, FUN = mean), tolerance = 0.0020))
## NOTE: Strangely, correctly weighting the number of days per month
##       (using 28.25 for February) is *not* closer than the simple mean:
ndays <- c(31, 28.25, rep(c(31,30, 31,30, 31), 2))
all.equal(s.y, aggregate(s.m, FUN = mean))                     # 0.0013
all.equal(s.y, aggregate(s.m, FUN = weighted.mean, w = ndays)) # 0.0017

来源

H. Tong (1996) 非线性时间序列。克拉伦登出版社,牛津,p。 471.

也可以看看

有关每月太阳黑子数量,请参阅 sunspot.monthsunspots

每年定期更新的太阳黑子数量可从 WDC-SILSO、比利时皇家天文台 http://www.sidc.be/silso/datafiles 获取

相关用法


注:本文由纯净天空筛选整理自R-devel大神的英文原创作品 Yearly Sunspot Data, 1700–1988。非经特殊声明,原始代码版权归原作者所有,本译文未经允许或授权,请勿转载或复制。