plot(data$q1a)

acf(data$q1a,lag.max=100)

pacf(data$q1a,lag.max=10)

q1adiff<-data$q1a[2:500]-data$q1a[1:499]
acf(q1adiff,lag.max=10)

pacf(q1adiff,lag.max=10)

ans1<- arima(data$q1a,order=c(2,1,0))
plot(data$q1a[4:500],ans1[['residuals']][1:497],xlab = 'predicted value',ylab='residual')

plot(data$q1b)

acf(data$q1b,lag.max=100)

pacf(data$q1b,lag.max=100)

ans2<-arima(data$q1b,order=c(2,0,0))
plot(data$q1b[3:500],ans2[['residuals']][1:498],xlab = 'predicted value',ylab='residual')

plot(data$q1c)

acf(data$q1c,lag.max=10)

pacf(data$q1c,lag.max=10)

q1cdiff<-data$q1c[2:500]-data$q1c[1:499]
plot(q1cdiff)

acf(q1cdiff,lag.max=10)

pacf(q1cdiff,lag.max=10)

ans3<-arima(data$q1c,order=c(0,1,2))
plot(data$q1c[2:500],ans3[['residuals']][1:499],xlab = 'predicted value',ylab='residual')

plot(data$q1d)

q1dlog<-log(data$q1d)
plot(q1dlog)

acf(q1dlog,lag.max=10)

pacf(q1dlog,lag.max=10)

ans4<-arima(q1dlog,order=c(1,0,0))
plot(data$q1d[2:500],ans4[['residuals']][1:499],xlab = 'predicted value',ylab='residual')