Day1-R 01-get_start

2.1.向量生成

(1)用 c()逐一放到一起,允许重复值

> c(2,5,6,2,9) 
[1] 2 5 6 2 9
> c("a","f","md","b")
[1] "a"  "f"  "md" "b" 

(2)连续的数字用冒号“:”

> 1:15
 [1]  1  2  3  4  5  6  7  8  9 10 11 12 13 14 15

(3)有重复的用rep(),有规律的序列用seq(),等差数列随机数用rnorm

> rep("gene",times=15)  
 [1] "gene" "gene" "gene" "gene" "gene" "gene" "gene" "gene" "gene"
[10] "gene" "gene" "gene" "gene" "gene" "gene"
> seq(from = 3,to = 21,by = 3)#等差数列
[1]  3  6  9 12 15 18 21
> rnorm(n = 5, mean = 3, sd = 5)#正态曲线:平均值+-标准差
[1]  7.564325 -2.151259  5.991248  2.028028  7.973115
 #(4)通过组合,产生更为复杂的向量。
> paste0(rep("gene",times=15),1:15)#两个向量一对一连接
 [1] "gene1"  "gene2"  "gene3"  "gene4"  "gene5"  "gene6"  "gene7" 
 [8] "gene8"  "gene9"  "gene10" "gene11" "gene12" "gene13" "gene14"
[15] "gene15"

练习2: 向量生成

1.将两种不同类型的数据用c()组合在一起,看输出结果


> c(1,"a")
[1] "1" "a"
> c(T,F,0)
[1] 1 0 0
> c(T,"a")
[1] "TRUE" "a"   

2.生成1到30之间所有4的倍数,答案是


  #4,8,12,16,20,24,28
> seq(from = 4,to = 30,by = 4)
[1]  4  8 12 16 20 24 28

3.生成sample4,sample8,sample12…sample28

>  # 提示:paste0
> paste0(rep("sample",times = 7),
+        seq(from = 4,to = 30,by = 4))
[1] "sample4"  "sample8"  "sample12" "sample16" "sample20" "sample24"
[7] "sample28"

2.2对单个向量进行的操作

(1)赋值给一个变量名

> x = c(1,3,5,6,2) #随意的写法
> x
[1] 1 3 5 6 2
> x <- c(1,3,5,3) #规范的赋值符号Alt➕减号
> x#四个元素
[1] 1 3 5 3
> #赋值+输出同步运行 
> x <- c(1,3,5);x#两行命令用分号连接
[1] 1 3 5
> (x <- c(1,3,5))
[1] 1 3 5

(2)简单数学计算

> x+1
[1] 2 4 6
> log(x)
[1] 0.000000 1.098612 1.609438
> sqrt(x)
[1] 1.000000 1.732051 2.236068
> x %% 10 #取整
[1] 1 3 5
> x %/% 3 #取余
[1] 0 1 1

(3)根据某条件进行判断

> x>3
[1] FALSE FALSE  TRUE
> x==3
[1] FALSE  TRUE FALSE

(4)初级统计

> max(x) #最大值
[1] 5
> min(x) #最小值
[1] 1
> range(x) #max和min极差
[1] 1 5
> which.max(x) #最大值下标#返回的是位置(即下标,索引)
[1] 3
> which.min(x) #最小值下标
[1] 1
> sort(x,decreasing = T)
[1] 5 3 1
> mean(x) #均值
[1] 3
> median(x) #中位数
[1] 3
> var(x) #方差
[1] 4
> sd(x) #标准差
[1] 2
> sum(x) #总和
[1] 9
> length(x) #长度?
[1] 3
> unique(x) #去重复
[1] 1 3 5
> table(x) #重复值统计
x
1 3 5 
1 1 1 

2.3.对两个向量进行的操作#####

> x = c(1,3,5)
> y = c(3,2,5)

(1)逻辑比较,生成等长的逻辑向量

> x == y 
[1] FALSE FALSE  TRUE
> x %in% y #x中的每一个元素判断是否在y中,输出逻辑值
[1] FALSE  TRUE  TRUE

(2)数学计算

> x + y#一一相加
[1]  4  5 10

(3)“连接“

> paste(x,y,sep=":")#注意与paste0的区别,是否可以用sep
[1] "1:3" "3:2" "5:5"

(4)交集、并集、差集

> intersect(x,y)#交集
[1] 3 5
> union(x,y)#并集,,自动去重
[1] 1 3 5 2
> setdiff(x,y)#差集,返回前面那个有的,你有我没有,以第一个为准
[1] 1
> #当两个向量长度不一致
> x == y
[1] FALSE FALSE  TRUE

2.4.利用循环补齐简化代码

> paste0(rep("gene",15),1:15)#按tab键,然后选择下一类
 [1] "gene1"  "gene2"  "gene3"  "gene4"  "gene5"  "gene6"  "gene7" 
 [8] "gene8"  "gene9"  "gene10" "gene11" "gene12" "gene13" "gene14"
[15] "gene15"
> paste0("gene",1:15)#循环补齐原理
 [1] "gene1"  "gene2"  "gene3"  "gene4"  "gene5"  "gene6"  "gene7" 
 [8] "gene8"  "gene9"  "gene10" "gene11" "gene12" "gene13" "gene14"
[15] "gene15"
> x <- 8:12
> #根据逻辑值取子集
> x[x==10]
[1] 10
> x[x<12]
[1]  8  9 10 11
> x[x %in% c(9,13)]
[1] 9
> #根据位置取子集
> x[4]
[1] 11
> x[2:4]
[1]  9 10 11
> x[c(1,5)]
[1]  8 12
> x[-4]
[1]  8  9 10 12
> x[-(2:4)]
[1]  8 12

2.5.修改向量中的某个/某些元素:取子集+赋值

> x[4] <- 40
> x
[1]  8  9 10 40 12
> x[x>3] <- 3
> x
[1] 3 3 3 3 3

练习3

1.将基因名"ACTR3B","ANLN","BAG1","BCL2","BIRC5","RAB","ABCT","ANF","BAD","BCF","BARC7","BALV"组成一个向量,赋值给x

> x = c("ACTR3B","ANLN","BAG1","BCL2","BIRC5","RAB","ABCT","ANF","BAD","BCF","BARC7","BALV")

2.用函数计算向量长度

> length(x)
[1] 12

3.用向量取子集的方法,选出第1,3,5,7,9,11个基因名。

> seq(1,11,2)
[1]  1  3  5  7  9 11
> x[seq(1,11,2)]
[1] "ACTR3B" "BAG1"   "BIRC5"  "ABCT"   "BAD"    "BARC7" 
> x[c(1,3,5,7,9)]
[1] "ACTR3B" "BAG1"   "BIRC5"  "ABCT"   "BAD"   

4.用向量取子集的方法,选出除倒数第2个以外所有的基因名。

> x[-11]
 [1] "ACTR3B" "ANLN"   "BAG1"   "BCL2"   "BIRC5"  "RAB"    "ABCT"  
 [8] "ANF"    "BAD"    "BCF"    "BALV"  
> length(x)-1#倒数第二个
[1] 11
> x[-(length(x)-1)]#反向取值
 [1] "ACTR3B" "ANLN"   "BAG1"   "BCL2"   "BIRC5"  "RAB"    "ABCT"  
 [8] "ANF"    "BAD"    "BCF"    "BALV"  

5.用向量取子集的方法,选出出在c("ANLN", "BCL2","TP53")中有的基因名。

> # 提示:%in%
> y = c("ANLN", "BCL2","TP53")
> x %in% y
 [1] FALSE  TRUE FALSE  TRUE FALSE FALSE FALSE FALSE FALSE FALSE FALSE
[12] FALSE
> x[x %in% y]
[1] "ANLN" "BCL2"

6.修改第6个基因名为"a"并查看是否成功

> x[6] <- "a"

7.生成100个随机数: rnorm(n=100,mean=0,sd=18)

> z = rnorm(n=100,mean=0,sd=18)
> #将小于-2的统一改为-2,将大于2的统一改为2
> z[z< (-2)] <- -2#加括号shift+9或者空格  不然成了赋值号
> z[z>2] <- 2
> 

1.两个不同类型的向量合并在一起会怎么样?
1.两个不同类型的向量合并在一起会怎么样?
2.如何将两个向量合到一起,组成一个长向量?
1.两个不同类型的向量合并在一起会怎么样?
2.如何将两个向量合到一起,组成一个长向量?
3.如何在向量首/尾/中间某位置增加一个元素?

©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 203,772评论 6 477
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 85,458评论 2 381
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 150,610评论 0 337
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 54,640评论 1 276
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 63,657评论 5 365
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 48,590评论 1 281
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 37,962评论 3 395
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 36,631评论 0 258
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 40,870评论 1 297
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 35,611评论 2 321
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 37,704评论 1 329
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 33,386评论 4 319
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 38,969评论 3 307
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 29,944评论 0 19
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 31,179评论 1 260
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 44,742评论 2 349
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 42,440评论 2 342

推荐阅读更多精彩内容