Help required with the stl function

["Gladwin, Philip [FI]" <philip.gladwin@citigroup.com>]

Good afternoon,

I am using S plus (windows 2000) to investigate seasonality.
I am using the command

myData.stl <- stl(myData$ts, "periodic")
ts.plot(myData$ts, myData.stl$sea, myData.stl$rem)

to decompose the data into a seasonal and remainder components and I am very 
pleased
with the plot.  

However I would like to add error bars / confidence intervals to the plot.
For the life of me I can't work out how to do it.


Can anyone out there help me?

kind regards,
Philip Gladwin