R语言学习笔记

读取数据

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()                                                # 设置保存图片的代码的收尾

指定画布大小

R语言_par()函数用法

• 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_{BP}$统计量
type='Ljung-Box – $Q_{LB}$统计量
参数lag=n表示输出滞后n阶的白噪声统计量。默认lag=1

数据索引

同时提取某几列

同时提取x5、x6…、lnx13列

bo_log[c('x5','x6','lnx8','x9','x10','x11','x12','lnx13')]