setwd("H:/myR")
DT <- read.csv("beijing.csv",header=T)
DT

library(tseries)
y=DT$y
t=2000:2014
plot(t,y,xlab="year",ylab="Average selling price")

pacf(y)
acf(y)

library(tseries)
pp.test(y)

Box.test(y,lag=6)
Box.test(y,lag=12)

t=c(1:15)
y<-ts(y)
plot(y,type="o")
abline(lm(y~t),col=2)

dt=1:14
dy=diff(y)
plot(dt,dy,type="o")
pp.test(dy)

ddt=1:13
ddy=diff(dy)
plot(ddt,ddy,type="o")
pp.test(ddy)

acf(ddy)
pacf(ddy)

Box.test(ddy,lag=6)
Box.test(ddy,lag=12)

y.fit<arima(ddy,order=c(0,2,1))
y.fit

y.fore<-forecast(y,h=5)
y.fore
plot(y.fore)


为什么在Rstudio可以运行,但是转化成PDF就不可以了?文件是Rnw格式