生信 R语言

11.芯片表达矩阵下游分析

R 复制代码
​rm(list = ls())#清除所有变量
options(stringsAsFactors = F)
#BiocManager::install("CLL")
suppressPackageStartupMessages(library(CLL))
data("sCLLex")
sCLLex
R 复制代码
## ExpressionSet (storageMode: lockedEnvironment)
## assayData: 12625 features, 22 samples 
##   element names: exprs 
## protocolData: none
## phenoData
##   sampleNames: CLL11.CEL CLL12.CEL ... CLL9.CEL (22 total)
##   varLabels: SampleID Disease
##   varMetadata: labelDescription
## featureData: none
## experimentData: use 'experimentData(object)'
## Annotation: hgu95av2
R 复制代码
exprSet=exprs(sCLLex)
samples=sampleNames(sCLLex)
pData=pData(sCLLex)
#没有代码,暂未做

12.RNA-seq 表达矩阵差异分析

R 复制代码
rm(list = ls())
options(stringsAsFactors = F)
library(airway)
data("airway") #加载数据
exprSet=assay(airway)
colnames(exprSet)
R 复制代码
group_list=colData(airway)[,3]
exprSet=exprSet[apply(exprSet, 1,function(x)sum(x>1)>5),]
table(group_list)
R 复制代码
group_list
R 复制代码
exprSet[1:4,1:4]
R 复制代码
boxplot(exprSet)
R 复制代码
boxplot(log(exprSet+1))
dev.off()
相关推荐
是阿建吖!6 分钟前
【Linux | 网络】网络基础
linux·网络
mmoyula6 分钟前
【RK3568 驱动开发:实现一个最基础的网络设备】
android·linux·驱动开发
阿杆27 分钟前
😡同事查日志太慢,我现场教他一套 grep 组合拳!
linux·后端
墨小傲39 分钟前
基于Linux下的vscode c/c++开发环境搭建详细教程
linux·c语言·vscode
cui_win1 小时前
【网络】Linux 内核优化实战 - net.ipv4.tcp_timestamps
linux·网络·tcp/ip
czhc11400756631 小时前
Linux 77 FTP
linux·运维·服务器
魂尾ac1 小时前
Centos Docker 安装(100%成功)
linux·docker·centos
天若有情6731 小时前
【Linux】02_CentOS 7 开机运行级别详解:从基础概念到实战配置
linux·运维·centos
ajassi20001 小时前
开源 C# .net mvc 开发(八)IIS Express轻量化Web服务器的配置和使用
linux·开源·c#·mvc·.net
張萠飛1 小时前
Linux中程序的limits中的Max open files的配置由哪些参数决定
linux·运维·服务器