Dear All,

Would anyone happen to have tips on how to do a pseudo-additive seasonal decomposition in R?

I am working on a ca. 20 year monthly time series on species abundance data, with annual peaks of varying magnitude and zero abundances between the seasonal occurrences.

I have tried to use the package "x12", which utilizes x12arima, but without luck so far. More specifically, I am trying to run the following code, using R version 2.15.0:

Cx12 <- x12(Cyanea.ts, period=12,
                x12path="C:\\x12a.exe",  use="x12",
                forecast_years=0,
                transform="none",
                x11regress=F,
                samode="pseudoadd",
                seasonalma="s3x3")

Everything seems to run ok and a number of files get saved in the source directory and folder 'gra'. However, when I try to call the resulting object Cx12, e.g. with 'summary(Cx12)', I receive an "Error: object 'Cx12' not found".

If anyone has a solution to the above problem, or alternative solutions to getting the decomposition done, I would much appreciate the help!

Regards,
Aino

Aino Hosia
Postdoc
Havforskningsinstituttet/Institute of Marine Research
PO Box 1870 Nordnes, N-5817 Bergen, Norway
(Nordnesgaten 50)
Tel: +47 55 23 53 49
E-mail: aino.hosia@imr.no<mailto:aino.hosia@imr.no>
www.imr.no<http://www.imr.no/>




	[[alternative HTML version deleted]]

