矢島町（年次） - waseda2bun_r @ ウィキ

矢島町の年次データ

• データ出所気象庁
• 期間：1980-2009年
• ファイル名　"YSMYR8008.csv"
• 使用ソフトR

分析

• （１）データの読み込み、確認

---

YR8808 <- read.table("YSMYR8008.csv",sep=",",skip=1)

#名前を付ける

RI　<- YR8808[,1]
SU　<- YR8808[,2]
SN　<- YR8808[,3]
AT　<- YR8808[,4]
HT　<- YR8808[,5]
LT　<- YR8808[,6]

#時系列化

RI2 <- ts(RI, start=c(1980),frequency=1)
SU2 <- ts(SU, start=c(1980),frequency=1)
SN2 <- ts(SN, start=c(1980),frequency=1)
AT2 <- ts(AT, start=c(1980),frequency=1)
HT2 <- ts(HT, start=c(1980),frequency=1)
LT2 <- ts(LT, start=c(1980),frequency=1)

#時系列プロット

par(mfrow=c(3,2))
ts.plot(RI2,xlab='1980-2008',ylab='mm',main="年間降水量@yashima")
ts.plot(SU2,xlab='1980-2008',ylab='hours',main="年間日照時間@yashima")
ts.plot(SN2,xlab='1980-2008',ylab='cm',main="年間降雪量@yashima")
ts.plot(AT2,xlab='1980-2008',ylab='deg C',main="年平均気温最高気温")
ts.plot(HT2,xlab='1980-2008',ylab='deg C',main="年間最高気温@yashima")
ts.plot(LT2,xlab='1980-2008',ylab='deg C',main="年間最低気温@yashima")

---

• 結果はコチラ

• （２）データを調べる

---

#相関関係

par(mfrow=c(3,2))
acf(SN2)
ccf(AT2,SN2)
ccf(SU2,SN2)
ccf(HT2,SN2)
ccf(LT2,SN2)
ccf(AT2,SN2)

---

• 結果はコチラ