R

R语言常用语法

基本操作

数据的读入、显示、操作

setwd("C:/Users/UniversityGrass/Desktop/Rlanguage/Data")//设置路径
read.csv('SZCI.csv')->szci//读取文件，->赋值箭头指向的是被赋值变量
View(szci);names(szci)//一行里有两句话时才需要加分号。
//View可现实数据。names可输出数据的列名
szci$dailyr[1:200]<-diff(szci$Close)/szci$Close[1:200]//创建一列新数据。不输出
szci$dailyr[200]//不带分号可直接输出。$用来标识列名

一列的最小值、最大值、中位数、平均数、方差、偏度、峰度

smin<-min(szci$dailylogr);smin//最后写smin是为了输出
smax<-max(szci$dailylogr);smax
smed<-median(szci$dailylogr);smed
smean<-mean(szci$dailylogr);smean
svar<-var(szci$dailylogr);svar
sskew<-skewness(szci$dailylogr);sskew
skurt<-kurtosis(szci$dailylogr);skurt

t检验

t.test(szci$dailylogr)
t.test(szci$dailylogr,alternative = "greater")

正态性检验

tseries::jarque.bera.test(szci$dailylogr)

时间序列RW1假设检验：εt~IID(0，σ^2)

tseries::runs.test(as.factor(szci$dailylogr>0))

计算ACF

acf(diff(szci$dailylogr),6)

画图

plot(X,Y)
//以下是用来画时间序列的图。xaxt="n"是把x坐标轴设为空值
plot(pgun$PG, type="l", xaxt="n", main="P&G returns")
axis(1, at=1:264, labels = pgun$Date)
library(lubridate)
d<-c("2012-1-10","2013-5-9","2014-6-25")
year(d)
month(d)
day(d)

求和

sum(alli$res^2)

抽取时间中的年月日

library(lubridate)
d