这学期的《应用时间序列分析》、《R语言》、《统计案例分析》都需要使用R语言,故将课程中学到的代码以及相关补充代码分享出来,一方面促进自身学习的积极主动性,另一方面或许可以给学习R语言的同学们提供一点点帮助。
rm(list=ls())
install.packages("tseries")
install.packages("astsa")
install.packages("zoo")#插值
library(tseries)
library(astsa)
library(zoo)
##时间序列数据展示——自输
#函数ts
price = c(101,82,66,35,31,7,9,4)
price = ts(price,start = c(2005,1),frequency = 12 )
price
##时间序列数据展示——scan
price = scan()#填完数据之后要回车一下,得到Read n items
price = ts(price,start = c(2005,1),frequency = 12 )
price
##输入数据
data = read.table(".csv",sep = ",",header = T)
data = read.table(".txt",sep = "\t")
##缺失值插值
a = c(1,3,6,9,10,24,79)
a[4] = NA
a
#线性插值
y1 = na.approx(a)
y1
#样条插值
y2 = na.spline(a)
y2
##数据读出
write.table(data.frame(),file = ".csv",sep = ",",row.names = F)
#