Displaying 1 result from an estimated 1 matches for "stlpa".
Did you mean:
sglpa
2011 Dec 01
3
Assign name to object for each iteration in a loop.
...tart=2000
TS is a timeseries:
X2000049 X2000065 X2000081 X2000097 X2000113
1 0.2080 0.2165 0.2149 0.2314 0.2028
2 0.1578 0.1671 0.1577 0.1593 0.1672
3 0.1897 0.1948 0.2290 0.2292 0.2067
Here's the function:
STLpA<-function(TS, area, Ystart, period=23, nSG="5,5", DSG=0)
{
require (RTisean)
for(i in 1:unique(area)){
vi.metric=TS[area==i]
filt.vi<-sav_gol(vi.metric,n=nSG,D=DSG)
vi.sg<-ts(filt.vi[,1], start=Ystart,frequency=period)
stld.tmp<-stl(vi.sg, s.window="periodic", robust...