The R Project SVN R

Rev

Rev 61160 | Rev 67599 | Go to most recent revision | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed

Rev 61160 Rev 61168
Line 124... Line 124...
124 
plot(stmR, set.pars = NULL)
 124 
plot(stmR, set.pars = NULL)
125 
# mark the 'outliers' :
 125 
# mark the 'outliers' :
126 
(iO <- which(stmR $ weights  < 1e-8)) # 10 were considered outliers
 126 
(iO <- which(stmR $ weights  < 1e-8)) # 10 were considered outliers
127 
sts <- stmR$time.series
 127 
sts <- stmR$time.series
128 
points(time(sts)[iO], 0.8* sts[,"remainder"][iO], pch = 4, col = "red")
 128 
points(time(sts)[iO], 0.8* sts[,"remainder"][iO], pch = 4, col = "red")
129 
par(op)# reset
 129 
par(op)   # reset
130 
}
 130 
}
131 
\keyword{ts}
 131 
\keyword{ts}