20231017思考题_Rmd笔记练习_TATA
原创
这是一道来自生信技能树的课堂练习
- 首先用
knitr::opts_chunk$set设置全局环境 echo = T表示显示结果时，一同显示代码
knitr::opts_chunk$set(echo = T,message = F,warning=F)思考题：挑出一个表达矩阵里方差最大的1000基因名
（R_04文件夹下test2.Rdata里有表达矩阵）
策略一
1、计算每个基因的方差
2、每个基因的方差降序排序
3、取前1000个数字
4、提取对应的基因名
嵌套式代码
rm(list = ls())
load("test2.Rdata")
t=names(head(sort(apply(test, 1, var),decreasing = T),1000))
t管道符传递
rm(list = ls())
load("test2.Rdata")
t=test %>%
apply(1, var) %>%
sort(decreasing = T) %>%
head(1000) %>%
names()
t策略二
1、计算每个基因的方差
2、每个基因的方差升序排序
3、取后1000个数字
4、提取对应的基因名
嵌套式代码
rm(list = ls())
load("test2.Rdata")
t=names(tail(sort(apply(test, 1, var)),1000))
t管道符传递
rm(list = ls())
load("test2.Rdata")
t=test %>%
apply(1, var) %>%
sort() %>%
tail(1000) %>%
names()
t参考：
原创声明：本文系作者授权腾讯云开发者社区发表，未经许可，不得转载。
如有侵权，请联系 cloudcommunity@tencent.com 删除。
原创声明：本文系作者授权腾讯云开发者社区发表，未经许可，不得转载。
如有侵权，请联系 cloudcommunity@tencent.com 删除。
评论
登录后参与评论
推荐阅读
目录
腾讯云开发者
