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


R sunspot.month 每月太阳黑子数据,从 1749 年至今


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

说明

来自世界数据中心(又名 SIDC)的每月太阳黑子数量。这是当新计数可用时偶尔会更新的数据版本。

用法

sunspot.month

格式

单变量时间序列 sunspot.yearsunspot.month 分别包含 289 个和 2988 个观测值。这些对象属于 "ts" 类。

例子

require(stats); require(graphics)
## Compare the monthly series
plot (sunspot.month,
      main="sunspot.month & sunspots [package'datasets']", col=2)
lines(sunspots) # -> faint differences where they overlap

## Now look at the difference :
all(tsp(sunspots)     [c(1,3)] ==
    tsp(sunspot.month)[c(1,3)]) ## Start & Periodicity are the same
n1 <- length(sunspots)
table(eq <- sunspots == sunspot.month[1:n1]) #>  132  are different !
i <- which(!eq)
rug(time(eq)[i])
s1 <- sunspots[i] ; s2 <- sunspot.month[i]
cbind(i = i, time = time(sunspots)[i], sunspots = s1, ss.month = s2,
      perc.diff = round(100*2*abs(s1-s2)/(s1+s2), 1))

## How to recreate the "old" sunspot.month (R <= 3.0.3):
.sunspot.diff <- cbind(
    i = c(1202L, 1256L, 1258L, 1301L, 1407L, 1429L, 1452L, 1455L,
          1663L, 2151L, 2329L, 2498L, 2594L, 2694L, 2819L),
    res10 = c(1L, 1L, 1L, -1L, -1L, -1L, 1L, -1L,
          1L, 1L, 1L, 1L, 1L, 20L, 1L))
ssm0 <- sunspot.month[1:2988]
with(as.data.frame(.sunspot.diff), ssm0[i] <<- ssm0[i] - res10/10)
sunspot.month.0 <- ts(ssm0, start = 1749, frequency = 12)

作者

R

来源

WDC-SILSO,太阳影响数据分析中心 (SIDC),比利时皇家天文台,Av. Circulaire, 3, B-1180 布鲁塞尔 目前位于http://www.sidc.be/silso/datafiles

也可以看看

sunspot.monthsunspots 的更长版本;后者运行到 1983 年并保持固定(为了作为示例数据集的再现性)。

相关用法


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