https://www.zybuluo.com/zhoujj2013/note/847578
https://bioconductor.org/packages/release/bioc/html/Mfuzz.html
https://www.ncbi.nlm.nih.gov/pmc/articles/PMC2139991/
library(Mfuzz)
data(yeast)
yeast.r <- filter.NA(yeast, thres=0.25)
yeast.f <- fill.NA(yeast.r,mode="mean")
tmp <- filter.std(yeast.f,min.std=0)
yeast.s <- standardise(yeast.f)
cl <- mfuzz(yeast.s,c=16,m=1.25)
pdf("aa.pdf")
mfuzz.plot(yeast.s,cl=cl,mfrow=c(4,4),time.labels=seq(0,160,10), new.window=FALSE)
dev.off()
aa = data.frame(cl)
write.table(aa,"aa.txt",sep="\t", row.names=T, col.names=F, quote=F)
网友评论