R语言计算机转录组

R package:Mfuzz进行基因表达的时间趋势分析并划分聚

2022-12-31  本文已影响0人  佳名

安装

BiocManager::install("Mfuzz")
library("Mfuzz")
data(yeast)
##缺失值的处理
yeast.r <- filter.NA(yeast, thres=0.25)
yeast.f <- fill.NA(yeast.r,mode="mean")#还可以是knn/wknn
##标准化数据
yeast.s <- standardise(yeast.f)
## m值评估
m1 <- mestimate(yeast.s)
## 评估C聚类簇数
tmp  <- Dmin(yeast.s,m=m1,crange=seq(4,40,4),repeats=3,visu=TRUE)

##表达水平低或者波动小的数据处理
tmp <- filter.std(yeast.f,min.std=0)
聚类.png

参考资料
R语言Mfuzz包:分析时间序列的组学数据 - 墨天轮 (modb.pro)
https://mp.weixin.qq.com/s/J5hNJlNj2-fSIOkFLlfhMg

上一篇下一篇

猜你喜欢

热点阅读