如下图所示
除了中位数标准化之外,我们还可以使用z-score的方法来对表达谱数据进行标准化:
z-score=(表达量-均值)/标准差
那么下面小编就给大家演示一下如何使用前面讲到的☞R中的sweep...函数,使用z-score的方法来对表达谱矩阵进行标准化
#为了保证随机数保持一致,这里设置一下种子序列
set.seed(123)
#随机生成100个数,构造一个10X10的矩阵
data=matrix...(runif(100,1,10),nrow=10)
#设置行名是gene1到gene10
rownames(data)=paste0("gene",1:10)
#设置列明是sample1到sample10...colnames(data)=paste0("sample",1:10)
#计算每一行的均值
rowmean=apply(data,1,mean)
#计算每一行的标准差
rowsd=apply(data...,1,sd)
#每一行基因表达值减去这一行的均值
data1=sweep(data,1,rowmean)
#每一行基因表达值除以这一行的标准差
data2=sweep(data1,1,rowsd,'/'