# read data (may need to include full path)
y1<-scan("wine.dat")
plot(y1,type='l',xlab='Month',ylab='Wine Sales')

# various Box-Cox transforms
y2<-y1^0.5
plot(y2,type='l',xlab='Month',ylab='Transformed Wine Sales')
y3<-y1^0.25
plot(y3,type='l',xlab='Month',ylab='Transformed Wine Sales')
y4<-log(y1)
plot(y4,type='l',xlab='Month',ylab='Transformed Wine Sales')

# Plot differenced data at lag 12
y3dif<-y3[13:length(y3)]-y3[1:(length(y3)-12)]
plot(y3dif,type='l',xlab='Month',ylab='Transformed Wine Sales')

# Show acf and pacf of differenced transformed data
acf(y3dif)
pacf(y3dif)