Re: Legent to the Periodogram

classic Classic list List threaded Threaded
2 messages Options
Reply | Threaded
Open this post in threaded view
|

Re: Legent to the Periodogram

Peter Maclean
How Can I add a legent (showing x1, x2, x3, x4) to the last plot?
require(TSA)
require(graphics)
require(stats)
t<-1986:2011
x1<-cos(t*1990/2011)
x2<-cos(t*2000/20011)
x3<-sin(t*1990/2011)
x4<-sin(t*2000/2011)
y<-cbind(t,x1,x2, x3,x4)
y.time = ts(y.time, start=1986, frequency=1)
y.spc<-spec.pgram(y.time, spans = c(3,3), detrend=FALSE,log="no",plot = TRUE, kernel("modified.daniell", c(5,7)))
plot(y.spc, plot.type = "marginal", main="Smoothed Periodogram")


Peter Maclean
Department of Economics
UDSM

______________________________________________
[hidden email] mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
Reply | Threaded
Open this post in threaded view
|

Re: Legent to the Periodogram

Mikhail Titov-2
plot.spec uses matplot. see ?matplot for default col & lty and use
legend as usual.

P.S. You can add plot=FALSE to spec.pgram to prevent it from plotting

On 08/27/2011 05:39 PM, Peter Maclean wrote:

> How Can I add a legent (showing x1, x2, x3, x4) to the last plot?
> require(TSA)
> require(graphics)
> require(stats)
> t<-1986:2011
> x1<-cos(t*1990/2011)
> x2<-cos(t*2000/20011)
> x3<-sin(t*1990/2011)
> x4<-sin(t*2000/2011)
> y<-cbind(t,x1,x2, x3,x4)
> y.time = ts(y.time, start=1986, frequency=1)
> y.spc<-spec.pgram(y.time, spans = c(3,3), detrend=FALSE,log="no",plot = TRUE, kernel("modified.daniell", c(5,7)))
> plot(y.spc, plot.type = "marginal", main="Smoothed Periodogram")
>
>
> Peter Maclean
> Department of Economics
> UDSM
>
> ______________________________________________
> [hidden email] mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.

______________________________________________
[hidden email] mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.