R语言学习(瑞士军刀)

158 阅读10分钟

(windows安装安装路径不要出现中文,不然RStudio找不到R软件)

1、R中的数据结构

1.1向量,标量

1.2矩阵

1.3数组

1.4列表

1.5数据框

1.6因子

1.7时间序列

2、R语言基础

2.1运算符

  R 语言赋值使用的是左箭头 <- 符号,不过一些新版本也支持等号=,也可以用->符号(最新版本)。

2.2变量的定义

  R 语言的有效的变量名称由字母,数字以及点号 . 或下划线 _ 组成。

变量名称以字母或点开头。

  注意:首字母为“.”点,或者英文字母,“.”号后面不能直接跟数字

demo

> var.1 = c(0,1,2,3) 
> print(var.1)
[1] 0 1 2 3
> var.2 = c("Hello","R")
> print(var.2)
[1] "Hello" "R"    
> var.3 <= c("Hello","Java")
错误: 找不到对象'var.3'
#符号错误 <=   应该是 <-
> var.3 <- c("Hello","Java")
> print(var.3)
[1] "Hello" "Java" 
> c(true,1)->var.4
错误: 找不到对象'true'
#布尔值大小写敏感
> c(TRUE,1) -> var.4
> print(var.4)
[1] 1 1
>

  列出已定义的变量ls()函数

> var.1 = c(0,1,2,3) 
> var.2 = c("Hello","R")
> var.3 <- c("Hello","Java")
> c(TRUE,1) -> var.4
> print(ls())
[1] "var.1" "var.2" "var.3" "var.4"
> 

  删除变量可以使用rm()函数

> print(ls())
[1] "var.1" "var.2" "var.3" "var.4"
> rm(var.1)
> print(ls())
[1] "var.2" "var.3" "var.4"
> 

  linux中进入和退出R软件

[root@VM-0-3-centos ~]# R

R version 3.6.0 (2019-04-26) -- "Planting of a Tree"
Copyright (C) 2019 The R Foundation for Statistical Computing
Platform: x86_64-redhat-linux-gnu (64-bit)

R是自由软件,不带任何担保。
在某些条件下你可以将其自由散布。
用'license()''licence()'来看散布的详细条件。

R是个合作计划,有许多人为之做出了贡献.
用'contributors()'来看合作者的详细情况
用'citation()'会告诉你如何在出版物中正确地引用R或R程序包。

用'demo()'来看一些示范程序,用'help()'来阅读在线帮助文件,或
用'help.start()'通过HTML浏览器来看帮助文件。
用'q()'退出R.

> q()
Save workspace image? [y/n/c]: n

  运行R脚本,在命令窗口使用Rscript来执行该脚本

 

#linux 创建编辑脚本文件 并运行脚本文件
[root@VM-0-3-centos zxgtest]# touch test.R
[root@VM-0-3-centos zxgtest]# vim test.R
[root@VM-0-3-centos zxgtest]# Rscript test.R
[1] 1 2 3 4 5

2.3、常见的函数

  输入输出print()函数

> print("Hello R!")
[1] "Hello R!"
> 

  拼接cat()函数

> cat("张三","今年",18,"岁",'\n')
张三 今年 18 岁 
> 

  cat()函数输出到文件中,append=TRUE输出到文件中不覆盖,追加到后面

> cat("输出到文件中",file="/zxgtest/test.txt")
> cat("追加数据不覆盖之前文件数据",file="/zxgtest/test.txt",append=TRUE)

  sink() 函数可以把控制台输出的文字直接输出到文件中去:

> sink("/zxgtest/sinkfile.txt")
#这样会将控制台输出的数据全部存储到指定的文件中,控制台不输出数据,显然这不是我们需要的效果
#split=TRUE 可以实现既将控制台输出的文件保存到文件中,又在控制台显示
> sink("/zxgtest/sinkfile.txt",split=TRUE)
#如果想取消输出到文件,可以调用无参数的 sink :
>sink()

  readLines() 从文件中读入数据

 

  获取当前工作目录,设置当前工作目录,查看当前工作目录

# 当前工作目录
print(getwd())

# 设置当前工作目录
setwd("/zxgtwst/test2")

# 查看当前工作目录
print(getwd())

3、R语言两种注释方式

3.1、单行注释 #

# 变量赋值
a <- 9

3.2、多行注释"

if(TRUE) {
    "
    这是一个多行注释的实例
    注释内容放在单引号或双引号之间
    "
}

4、基础运算

优先级符号含义
1()括号
2乘方运算
3%%整除求余
 %/%整除
4*乘法
 /除法
5+加法
 -减法
> 1 + 2*3
[1] 7
> (1+2)*3
[1] 9
> 3 / 4
[1] 0.75
> 3.4-1.2
[1] 2.2
> 1 -4 * 0.5 ^ 3
[1] 0.5
> 8 / 3 %% 2
[1] 8
> 8 / 4 %% 2
[1] Inf
> 3 %% 2^2
[1] 3
> 10 / 3 %/% 2
[1] 10

  关系运算符(>  <  ==  !=  >=  <=)

> v <- c(2,4,6,9)
> t <- c(1,4,7,9)
> print(v>t)
[1]  TRUE FALSE FALSE FALSE
> print(v < t)
[1] FALSE FALSE  TRUE FALSE
> print(v == t)
[1] FALSE  TRUE FALSE  TRUE
> print(v!=t)
[1]  TRUE FALSE  TRUE FALSE
> print(v>=t)
[1]  TRUE  TRUE FALSE  TRUE
> print(v<=t)
[1] FALSE  TRUE  TRUE  TRUE
> 

  逻辑运算

运算符描述
&元素逻辑与运算符,将第一个向量的每个元素与第二个向量的相对应元素进行组合,如果两个元素都为TRUE则记过为TRUE,否则为FALSE
|元素逻辑或运算符,将第一个向量的每个元素与第二个向量的相对应元素进行组合,如果两个元素中有一个为TRUE,则结果为TRUE,如果都为FALSE,则返回FALSE
逻辑非运算符,返回向量每个元素相反的逻辑值,如果元素为TRUE则返回FALSE,如果元素为FALSE则返回TRUE
&&逻辑与运算符,只对两个向量对一个元素进行判断,如果两个元素都为TRUE,结果为TRUE,否则为FALSE
||逻辑或运算符,只对两个向量,对第一个元素进行判断,如果两个元素中一个为TRUE,则结果为TRUE,如果都为FALSE则返回FALSE
> v = c(3,1,TRUE,2+3i)
> t = c(4,1,FALSE,2+3i)
> print(v&t)
[1]  TRUE  TRUE FALSE  TRUE
> print(v|t)
[1] TRUE TRUE TRUE TRUE
> print(!t)
[1] FALSE FALSE  TRUE FALSE
> p = c(3,0,TRUE,2+2i)
> k = c(1,3,TRUE,2+3i)
> print(p&&k)
[1] TRUE
> print(p||k)
[1] TRUE
#&&和||只对第一个元素进行比较

  赋值运算符

运算符描述
<-、=、<<-向左赋值
->、->>向右赋值

  其他运算符

运算符描述
冒号运算符,用于创建一系列数字的向量
%in%用于判断元素是否再向量里,返回布尔值,有的话返回TRUE,没有的话返回FALSE
%*%用于矩阵与它转置的矩阵相乘
> v <- 1:10  #1到10的向量
> print(v)
 [1]  1  2  3  4  5  6  7  8  9 10
> s = c(1,2,3)
> print(s)
[1] 1 2 3
> v1 = 3
> v2 = 20
> print(v1 %in% v) #判断数字是否再向量v中
[1] TRUE
> print(v2 %in% v)
[1] FALSE
> M = matrix(c(2,6,5,1,10,4),nrow = 2,ncol = 3,byrow = TRUE)#矩阵与它置换的矩阵相乘
> t = M %*% t(M)
> print(t)
     [,1] [,2]
[1,]   65   82
[2,]   82  117
> 

  常用数学函数

函数说明
sqrt(n)n的平方根
exp(n)自然常数e的n次方
log(m,n)m的对数函数,返回n的几次方等于m
log10(m)相当于log(m,10)
> sqrt(4)
[1] 2
> sqrt(16)
[1] 4
> exp(1)
[1] 2.718282
> exp(2)
[1] 7.389056
> log(2,4)
[1] 0.5
> log10(10000)
[1] 4
> log(10,100)
[1] 0.5
> 

  取整函数

名称参数模型含义
round(n)对n四舍五入取整
round(n,m)对n保留m位小数四舍五入
ceiling(n)对n向上取整
floor(n)对n向下取整
> round(1.5)   #round函数谨慎使用,有时5也会舍去
[1] 2
> round(2.5)
[1] 2
> round(3.5)
[1] 4
> round(4.5)
[1] 4

  数学运算   三角函数、反三角函数

> sin(pi/6)
[1] 0.5
> cos(pi/4)
[1] 0.7071068
> tan(pi/3)
[1] 1.732051
> asin(0.5)
[1] 0.5235988
> acos(0.7071068)
[1] 0.7853981
> atan(1.732051)
[1] 1.047198
#下面几个都是正态分布的函数,名字以norm结尾,代表“正态分布”
#前缀   d   概率密度函数
#     p   概率密度积分函数(从无穷小到x的积分)
#     q   分位数函数
#     r   随机数函数(常用于概率仿真)

> dnorm(0)
[1] 0.3989423
> pnorm(0)
[1] 0.5
> qnorm(0.95)
[1] 1.644854
> rnorm(3,5,2) #产生3个平均值为5,标准差为2的正态随机数
[1] 1.351745 5.329550 2.958539

  R 语言除了含有正态分布函数以外还有泊松分布 (pois, Poisson) 等常见分布函数  --------“概率论与数理统计”

 注意:个位数字为4的话round函数

5、数据类型(数字、逻辑、文本)

****数字常量:一般型、科学计数法

  逻辑类型:常量值只有 TRUE和FALSE   区分大小写。

  对象类型:(向量、列表、矩阵、数组、因子、数据框)

下面详细介绍详细介绍对象类型

  5.1、向量(Vector)

 

  最常见的是二维向量,在平面坐标系中必然会用到。(可以看作java中的数组)。

  c()是创造向量的函数。\

  向量中的每个元素可以通过下标单独取出。(下标代表的第几个,起始值为1)(java等其他语言中的下标基本上从0开始)

  注意:二维向量和其他维度的向量运算失去数学意义,会出现警告。

 

> a=c(1,2,3,6,11,25,64,126)
> a[4]
[1] 6
#通过下标获取单独的值
> a[1:4]
[1] 1 2 3 6
#获取1到4的值
> a[c(1,3,6)]
[1]  1  3 25
#获取1,3,6位置的值
> a[c(-1,-3)]
[1]   2   6  11  25 shuchu  64 126
> 
#去掉1,3位置的值,并输出向量

  向量支持标量运算

> c(1.1,1.2,1.3) - 0.5
[1] 0.6 0.7 0.8
> a = c(1,2,4)
> a ^ 2
[1]  1  4 16

  常用的数学运算函数也适用向量运算 sqrt、exp、sort、rev、order

> a = c(1,3,5,2,4,6)
> sort(a)
[1] 1 2 3 4 5 6
> rev(a)
[1] 6 4 2 5 3 1
> order(a)
[1] 1 4 2 5 3 6
> a[order(a)]
[1] 1 2 3 4 5 6
#order(a) 函数为排序后a向量的下标

  R中有完善的统计函数

  sum、mean(求平均值)、var(方差)、sd(标准差)、min(最小值)、max(最大值)、range(取值范围)(二维向量最大值和最小值)

> sum(1:10)
[1] 55
> sd(1:10)
[1] 3.02765
> range(1:10)
[1]  1 10

  向量生成,生成等差数列,可以用seq函数

> seq(1,9,2)
[1] 1 3 5 7 9
> seq(1,8,2)
[1] 1 3 5 7
#seq(a,b,c)  从ab,步长为c
> seq(1,2,length.out=3) 
[1] 1.0 1.5 2.0
# 从m到n数列的长度为3
> rep(0, 5)
[1] 0 0 0 0 0
#重复数字0 5次

  向量中常会用到NA和NULL,介绍与区别

  NA代表的缺失,NULL代表的不存在

  NA缺失就像占位符,代表这里没有一个值,但位置存在

  NULL代表的就是数据不存在

> length(c(NA, NA, NULL))
[1] 2
> c(NA, NA, NULL, NA)
[1] NA NA NA
#NULL在R向量中没有任何意义

  which比较罕见的向量逻辑型处理函数

> c(1,2,3) > 2
[1] FALSE FALSE  TRUE
> a = c(1,2,3)>2
> b = c(1,2,3)
> print(a)
[1] FALSE FALSE  TRUE
> print(b)
[1] 1 2 3
> which(a)
[1] 3
> b = c(1,2,3,4,5,6)
> which(a)
[1] 3
> a = b >2
> which(a)
[1] 3 4 5 6

> vector = c(10,40,78,64,53,62,69,70)
> print(vector[which(vector >= 60 & vector < 70)])
[1] 64 62 69
>
#从一个线性表中筛选大于等于 60 且小于 70 的数据  

类似的函数还有any()  all() 等

> all(c(TRUE,TRUE,TRUE))
[1] TRUE
> all(c(TRUE,FALSE,TRUE))
[1] FALSE
> any(c(TRUE,FALSE,FALSE))
[1] TRUE
> any(c(FALSE,FALSE,FALSE))
[1] FALSE

  5.2、字符串(常用的对字符串操作的函数)

> toupper("Hello R")  #转换为大写
[1] "HELLO R"
> tolower("Hello R")  #转换为小写
[1] "hello r"
> nchar("汉字",type="bytes")   #统计字节长度
[1] 6
> nchar("中文",type="byte")  
[1] 6
> nchar("汉字",type="char")  #统计字符数
[1] 2
> substr("unbelievable",1,5)  #截取字符串,从1到5
[1] "unbel"
> substring("unbelievable",5)  #截取字符串,从5到结束
[1] "lievable"
> as.numeric("24")  #将字符串转化为数字
[1] 24
> as.character(12.34)  #将数字转化为字符串
[1] "12.34"
> strsplit("hello_java_learn_R","_")  #分隔符拆分字符串
[[1]]
[1] "hello" "java"  "learn" "R"    

> gsub("/","-","2021/08/21")  #替换字符串
[1] "2021-08-21"
> 

  R支持perl语言格式的正则表达式

  5.2、矩阵

> vector = c(1,2,3,4,5,6)
> matrix(vector,2,3)
     [,1] [,2] [,3]
[1,]    1    3    5
[2,]    2    4    6
#初始化内容是由一个向量来传递的,其次要输入矩阵有多少行多少列
> vector = c(1,2,3,4,5)
> matrix(vector,2,3)
     [,1] [,2] [,3]
[1,]    1    3    5
[2,]    2    4    1
Warning message:
In matrix(vector, 2, 3) : 数据长度[5]不是矩阵行数[2]的整倍
#当初始化向量数量不足时会给出警告
> vector = c(1,2,3,4,5,6)
> matrix(vector,2,3,byrow=TRUE)
     [,1] [,2] [,3]
[1,]    1    2    3
[2,]    4    5    6
#按行来初始化矩阵
> matrix(vector,2,3,byrow=FALSE)
     [,1] [,2] [,3]
[1,]    1    3    5
[2,]    2    4    6
#按列初始化矩阵
> m1 = matrix(vector,2,3,byrow=FALSE)
> m1[1,1]
[1] 1
> m1[1,2]
[1] 3
直接访问矩阵的具体元素
> colnames(m1) = c("x","y","z")
> rownames(m1) = c("a","b")
> m1
  x y z
a 1 3 5
b 2 4 6
> #设置行列的名称,这个过程通过字符串向量批量完成

  > m1["b",]
x y z
2 4 6
> m1["b","a"]
Error in m1["b", "a"] : 下标出界
> m1["a","x"]
[1] 1
> m1["b","y"]
[1] 4

#通过设定的名称访问具体的元素

  矩阵的运算

> m1 = matrix(c(1,2,3,4),2,2)
> m2 = matrix(c(5,6,7,8),2,2)
> m1%*%m2
     [,1] [,2]
[1,]   23   31
[2,]   34   46
#矩阵的乘法运算
> A = matrix(c(1,2,3,4),2,2)
> solve(A)
     [,1] [,2]
[1,]   -2  1.5
[2,]    1 -0.5
#逆矩阵
> A
     [,1] [,2]
[1,]    1    3
[2,]    2    4
> apply(A,1,sum)
[1] 4 6
8
> apply(A,2,sum)
[1] 3 7