library(dyn)
data1 = read.table("D:\\FSU-school\\statapp2\\project\\4-17\\data4.txt", header=TRUE)
attributes(data1)

attach(data1)
#TCMNOMM3 EDM3

#data1_ts1 = ts(data1$TCMNOMM3, start=c(2004, 4), frequency=365)
data1_ts2= ts(data1$EDM3, start=c(2004, 4), frequency=365)
data1_ts2

data1_ts3 = ts(spread, start=c(2006, 4), frequency=365)

data1_ts4 = ts(change, start=c(2006, 4), frequency=365)


plot.ts(data1_ts3)
plot.ts(data1_ts4)

cor(data1)
pairs(data1)

fit1 = dyn$lm(data1_ts4~data1_ts3)
fit1
summary(fit1)




plot.ts(fit2$fitted.values)


anova(fit1,fit2)





summary(fit1)
par(mfrow=c(2,2))
plot(fit1)

par(mfrow=c(2,1))
plot.ts(fit1$fitted.values)

plot.ts(data1_ts)