#Random Series (with trend and cyclic)

Y<- rnorm(200) 
seq<- 1:200
 
postscript(file="RNDcom.ps",height=16,width=10,horizontal=F)
par(mfrow=c(2,1))


plot(seq,Y,type="l",lty=1,xlim=c(0,200),ylim=c(-4,4),xlab="t",
   ylab="normal random series")
abline(h=c(-2,0,2),lty=2)
title(main="Random Serie",cex=0.85)

Yt<- Y + 0.1*seq+3*sin(2*pi*seq/50)
plot(seq,Yt,type="l",lty=1,xlim=c(0,200),xlab="t",
   ylab="normal random series + trend + sinus")
abline(h=c(0,5,10,15),lty=2)
title(main="Random Series + Trend + Cyclic Component",cex=0.85)

dev.off()
rm(Y,Yt,seq)