读取数据
x <- read.table('...csv', sep=',', header=T) # 需要指定sep
x <- read.csv('...csv') # 方便直接读取csv文件
x <- scan('...txt', sep='') # 读取txt为列表格式
生成循环,指定步长
c(1:5)
sep(-2, 2, 0.5) # 步长为0.5
计算
输出字符串插入变量+数字:
方法一:sprintf('P值=%f', pchisq(Q_lb, m))
方法二:paste('P值', pchisq(Q_lb,m))
在for
循环中,需要使用如下语句才能输出print(paste...)
生成指定长度的0向量 x<-vector(mode="numeric",length=m)
生成顺序数据,相当于python的range
c(1:5)
查看变量x的所有属性 str(x)
画图
plot(x, y, type, main, sub, xlab, ylab, xlim, ylim, pch, lty, lwd, col...)
- x,y: 横纵坐标数据
- type: 指定绘图的类型:
p
为点图,
l
为线图;
b
为点连线图;
o
为线穿过点;
h
为悬垂线
s
为阶梯线- main: 指定主标题
- sub : 指定副标题
- xlab: x轴标签
- ylab: y轴标签
- xlim: x轴范围
- ylim: y轴范围
- pch: 指定观察点的符号,可取从1-25的整数
- lty:指定连点类型,可取1-6的整数
- lwd:线条宽度
- col:颜色,可取整数
- family:字体
设置新罗马字体
TNM = windowsFonts(TNM = windowsFont("Times New Roman"))
plot(x, type='o', col='4', ylab='Sales Numbers', main='Timing chart of number of units sold by the company')
但是这对于acf()
函数不管用,请用以下代码。
acf(x, family="TNM", main='')
title(main='Sample autocorrelation graph of the number of units sold by the company',family="TNM")
添加标签图示
legend("topleft", inset=.05, title="Drug Type", c("A","B"),
lty=c(1, 2), pch=c(15, 17), col=c("red", "blue"))
保存图片
png(filename, width, height) ... dev.off()
通过不断调整width
、height
和res
,可以得到高清,且标题字体不小的图片。
R语言绘图、高清保存、以及颜色系统
png("type.png",width=600*3,height=3*600,res=72*3) # 这么设置,保存的图片高清
plot(GDP, type="o", xlab="年份", ylab='国内生产总值(GDP)', col=1)
dev.off() # 设置保存图片的代码的收尾
指定画布大小
- mfrow: 画子图
- pin:图形的宽和高
- mai/mar:边界的大小
转为时间序列,指定日期范围
x <- ts(data, start(2020,2), frequency=12)
# start:从2020年2月份开始算
# frequency:时间间隔
生成随机数
rpois()
rnorm(n=, mean=, std=)
生成服从正态分布的随即数。
纯随机性检验
Box.test(x, type=, lag=)
,参数type
为检验统计量类型:
type='Box-Pierce
–
Q
B
P
Q_{BP}
QBP统计量
type='Ljung-Box
–
Q
L
B
Q_{LB}
QLB统计量
参数lag=n
表示输出滞后n阶的白噪声统计量。默认lag=1
数据索引
同时提取某几列
同时提取x5、x6…、lnx13列
bo_log[c('x5','x6','lnx8','x9','x10','x11','x12','lnx13')]