2019年2月17日日曜日

plot abline eps GSPC




plot.default(diff(eps_year_xts["2007::2018"])[-1],type='h',axes=F)
par(new=T)
qtr <- seq(as.Date("2007-04-01"),as.Date("2018-12-31"),by='quarters')
plot.default(qtr,to.quarterly(GSPC["2007-04::2018"])[,4],type='l')
abline(v=as.Date("2018-01-01"),col = "gray60",lty=3)
abline(v=as.Date("2015-01-01"),col = "gray60",lty=3)
abline(v=as.Date("2017-01-01"),col = "gray60",lty=3)
abline(v=as.Date("2016-01-01"),col = "gray60",lty=3)
abline(h=2000,col="gray60",lty=3)
abline(v=as.Date("2014-01-01"),col = "gray60",lty=3)


For the longer period, use SP5 instead of GSPC

plot.default(diff(eps_year_xts["1992::2018"])[-1],type='h',axes=F,col=3)
par(new=T)
qtr <- seq(as.Date("1992-04-01"),as.Date("2018-12-31"),by='quarters')
plot.default(qtr,to.quarterly(SP5["1992-04::2018"])[,4],type='l')
abline(v=as.Date("2018-01-01"),col = "gray60",lty=3)
abline(v=as.Date("2015-01-01"),col = "gray60",lty=3)
abline(v=as.Date("2017-01-01"),col = "gray60",lty=3)
abline(v=as.Date("2016-01-01"),col = "gray60",lty=3)
abline(h=2000,col="gray60",lty=3)
abline(v=as.Date("2014-01-01"),col = "gray60",lty=3)




To see the relation between S&P500 and payroll. Change frequency from quarterly to monthly.

plot.default(diff(PA["1992::2018"])[-1],type='h',axes=F,col=3)
par(new=T)
mnt <- seq(as.Date("1992-02-01"),as.Date("2018-12-31"),by='months')
plot.default(mnt,to.monthly(SP5["1992-02::2018"])[,4],type='l')
abline(v=as.Date("2018-01-01"),col = "gray60",lty=3)
abline(v=as.Date("2015-01-01"),col = "gray60",lty=3)
abline(v=as.Date("2017-01-01"),col = "gray60",lty=3)
abline(v=as.Date("2016-01-01"),col = "gray60",lty=3)
abline(h=2000,col="gray60",lty=3)
abline(v=as.Date("2014-01-01"),col = "gray60",lty=3)
par(new=T)
plot.default(mnt,as.vector(cli_xts$oecd["1992-02::2018"]),axes=F,col=4,type='l')
par(new=T)
plot.default(mnt,na.trim(diff(cli_xts$oecd["1991-08::2018"],lag=6)),axes=F,col=6,type='l')




0 件のコメント: