Rの情報

10：白書さん：2010/10/15(金) 21:07:23 HOST:wcache2.waseda.ac.jp: a <- read.table("GDP1980-2010.csv",sep=",",skip=1)
x <- a[,1]
GDP <- ts(x,start=c(1980,1),frequency=4)

bunkai1 <- decompose(GDP,type="multiplicative")
GDPTC <- bunkai1$trend
GDPTC2 <- GDPTC[is.na(GDPTC)==FALSE]
GDPTC3 <- ts(GDPTC2,start=c(1980,3),frequency=4)
n <- length(GDPTC3)
x <- rep(0,13);b <- 0.7;p <- 20
w <- (1-b)*(cumprod(c(1,rep(b,p-1))));w
sum(w)

GDPINT <- rev(GDPTC3[(n-p+1):n])
YOSOKU <- filter(x,filter=w,method="recursive",init=GDPINT)
YOSOKU2 <- ts(YOSOKU,start=c(2005,1),frequency=4)
ttl <- "指数平滑法による予測"
ts.plot(GDPTC3,YOSOKU2,type="l",lty=c(1,2),main=ttl)

新着レスの表示

※書き込む際の注意事項はこちら

※画像アップローダーはこちら

（画像を表示できるのは「画像リンクのサムネイル表示」がオンの掲示板に限ります）

掲示板管理者へ連絡無料レンタル掲示板