sunspot.month            package:datasets            R Documentation

Monthly Sunspot Data, 1749-1997

Description:

     Monthly numbers of sunspots.

Usage:

     sunspot.month

Format:

     The univariate time series 'sunspot.year' and 'sunspot.month'
     contain 289 and 2988 observations, respectively. The objects are
     of class '"ts"'.

Source:

     World Data Center-C1 For Sunspot Index Royal Observatory of
     Belgium, Av. Circulaire, 3, B-1180 BRUSSELS <URL:
     http://www.oma.be/KSB-ORB/SIDC/sidc_txt.html>

See Also:

     'sunspot.month' is a longer version of 'sunspots' that runs until
     1988 rather than 1983.

Examples:

     require(stats); require(graphics)
     ## Compare the monthly series 
     plot (sunspot.month, main = "sunspot.month [stats]", col = 2)
     lines(sunspots) # "very barely" see something

     ## 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, sunspots = s1, ss.month = s2,
           perc.diff = round(100*2*abs(s1-s2)/(s1+s2), 1))

