\HeaderA{sunspot.month}{Monthly Sunspot Data, 1749--1997}{sunspot.month}
\keyword{datasets}{sunspot.month}
\begin{Description}\relax
Monthly numbers of sunspots.
\end{Description}
\begin{Usage}
\begin{verbatim}
sunspot.month
\end{verbatim}
\end{Usage}
\begin{Format}\relax
The univariate time series \code{sunspot.year} and
\code{sunspot.month} contain 289 and 2988 observations, respectively.
The objects are of class \code{"ts"}.
\end{Format}
\begin{Source}\relax
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}
\end{Source}
\begin{SeeAlso}\relax
\code{sunspot.month} is a longer version of \code{\LinkA{sunspots}{sunspots}}
that runs until 1988 rather than 1983.
\end{SeeAlso}
\begin{Examples}
\begin{ExampleCode}
require(stats)
## 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))
\end{ExampleCode}
\end{Examples}

